45 #define M_MAX(a,b) (((a) > (b)) ? (a) : (b))
46 #define M_MIN(a,b) (((a) < (b)) ? (a) : (b))
51 #define MMG2D_EPSD 1.e-10
52 #define MMG2D_EPSA 1.e-12
54 #define MMG2D_PRECI 1.
55 #define MMG2D_SIZE 0.75
56 #define MMG2D_ALPHA 0.28867513459
57 #define MMG2D_ALPHAD 3.464101615137755
58 #define MMG2D_LONMAX 1024
59 #define MMG2D_BADKAL 0.2
60 #define MMG2D_NULKAL 1.e-6
61 #define MMG2D_ANGCORN -1.e-6
62 #define MMG2D_SHORTMAX 0x7fff
63 #define MMG2D_LPARMAX 200
64 #define MMG2D_VOLFRAC 1.e-5
66 #define MMG2D_LLONG 2.0
67 #define MMG2D_LSHRT 0.3
68 #define MMG2D_LOPTL 1.4
69 #define MMG2D_LOPTS 0.71
71 #define MMG2D_NPMAX 50000
72 #define MMG2D_NEDMAX 100000
73 #define MMG2D_NEMAX 100000
76 static const uint8_t
MMG2D_idir_q[4][2] = { {0,1},{0,3},{1,2},{2,3} };
79 #define MMG2D_RETURN_AND_FREE(mesh,met,ls,disp,val)do \
81 if ( !MMG2D_Free_all(MMG5_ARG_start, \
82 MMG5_ARG_ppMesh,&mesh,MMG5_ARG_ppMet,&met, \
83 MMG5_ARG_ppLs,&ls,MMG5_ARG_ppDisp,&disp, \
85 return MMG5_LOWFAILURE; \
98 fprintf(stdout,
"\n Unexpected error:"); fflush(stdout);
101 fprintf(stdout,
" Abnormal stop\n");
break;
103 fprintf(stdout,
" Floating-point exception\n");
break;
105 fprintf(stdout,
" Illegal instruction\n");
break;
107 fprintf(stdout,
" Segmentation fault\n");
break;
110 fprintf(stdout,
" Program killed\n");
break;
136 #define MMG2D_POINT_REALLOC(mesh,sol,ip,wantedGap,law,o,tag ) do \
140 assert ( mesh && mesh->point ); \
141 MMG5_TAB_RECALLOC(mesh,mesh->point,mesh->npmax,wantedGap,MMG5_Point, \
142 "larger point table",law); \
144 mesh->npnil = mesh->np+1; \
145 for (klink=mesh->npnil; klink<mesh->npmax-1; klink++) \
146 mesh->point[klink].tmp = klink+1; \
151 MMG5_ADD_MEM(mesh,(sol->size*(mesh->npmax-sol->npmax))*sizeof(double), \
152 "larger solution",law); \
153 MMG5_SAFE_REALLOC(sol->m,sol->size*(sol->npmax+1), \
154 sol->size*(mesh->npmax+1), \
155 double,"larger solution",law); \
157 sol->npmax = mesh->npmax; \
161 ip = MMG2D_newPt(mesh,o,tag); \
167 #define MMG2D_TRIA_REALLOC(mesh,jel,wantedGap,law ) do \
171 oldSiz = mesh->ntmax; \
172 MMG5_TAB_RECALLOC(mesh,mesh->tria,mesh->ntmax,wantedGap,MMG5_Tria, \
173 "larger tria table",law); \
175 mesh->nenil = mesh->nt+1; \
176 for (klink=mesh->nenil; klink<mesh->ntmax-1; klink++) \
177 mesh->tria[klink].v[2] = klink+1; \
179 if ( mesh->adja ) { \
181 MMG5_ADD_MEM(mesh,3*(mesh->ntmax-oldSiz)*sizeof(int), \
182 "larger adja table",law); \
183 MMG5_SAFE_RECALLOC(mesh->adja,3*oldSiz+5,3*mesh->ntmax+5,int \
184 ,"larger adja table",law); \
188 jel = MMG2D_newElt(mesh); \
345 double long_ani(
double *ca,
double *cb,
double *ma,
double *mb);
346 double long_iso(
double *ca,
double *cb,
double *ma,
double *mb);
MMG5_pMesh MMG5_pSol * sol
Definition: API_functionsf_2d.c:63
MMG5_pMesh * mesh
Definition: API_functionsf_2d.c:63
const int va_list argptr
Definition: API_functionsf_2d.c:728
API headers for the mmg2d library.
API header for the common part of the MMG libraries.
MMG5_Tria * MMG5_pTria
Definition: libmmgtypes.h:325
MMG5_Sol * MMG5_pSol
Definition: libmmgtypes.h:647
MMG5_Mesh * MMG5_pMesh
Definition: libmmgtypes.h:627
int MMG2D_chkmanimesh(MMG5_pMesh)
Definition: mmg2d6.c:478
int optlen_iso(MMG5_pMesh mesh, MMG5_pSol sol, double declic, int base)
Bucket * pBucket
Definition: mmg2d.h:125
int MMG2D_mmg2d1n(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d1.c:808
int MMG2D_split3_sim(MMG5_pMesh, MMG5_pSol, int, int vx[3])
Definition: split_2d.c:587
int MMG2D_setadj(MMG5_pMesh)
Definition: analys_2d.c:48
int MMG2D_cuttri_ls(MMG5_pMesh, MMG5_pSol, MMG5_pSol)
Definition: mmg2d6.c:796
int MMG2D_split2(MMG5_pMesh, MMG5_pSol, int, int vx[3])
Definition: split_2d.c:500
int MMG2D_setMeshSize_alloc(MMG5_pMesh)
Definition: zaldy_2d.c:237
int MMG2D_cenrad_iso(MMG5_pMesh, double *, double *, double *)
Definition: cenrad_2d.c:42
int MMG2D_movintpt_ani(MMG5_pMesh, MMG5_pSol, int, int *, int8_t)
Definition: anisomovpt_2d.c:37
int MMG2D_memOption(MMG5_pMesh mesh)
Definition: zaldy_2d.c:218
pQueue MMG2D_kiuini(MMG5_pMesh mesh, int nbel, double declic, int base)
void MMG2D_Init_parameters(MMG5_pMesh mesh)
Definition: API_functions_2d.c:80
int MMG2D_movedgpt(MMG5_pMesh, MMG5_pSol, int, int *, int8_t)
Definition: movpt_2d.c:52
double MMG2D_quickcal(MMG5_pMesh, MMG5_pTria)
Definition: quality_2d.c:45
size_t MMG5_memSize(void)
Definition: tools.c:547
int MMG2D_kiupop(pQueue q)
int MMG2D_delone(MMG5_pMesh, MMG5_pSol, int, int *, int)
Definition: delone_2d.c:247
int MMG2D_locateEdge(MMG5_pMesh, int, int, int *, int *)
Definition: locate_2d.c:323
int MMG2D_colver2(MMG5_pMesh, int *)
Definition: colver_2d.c:418
int MMG2D_findtrianglestate(MMG5_pMesh, int, int, int, int, int, int)
Definition: mmg2d2.c:132
int interp_iso(double *, double *, double *, double)
int MMG5_mmg2dChkmsh(MMG5_pMesh, int, int)
Definition: chkmsh_2d.c:34
static const int MMG2D_iopp[3][2]
Definition: mmg2d.h:129
int MMG2D_bezierCurv(MMG5_pMesh, int, int8_t, double, double *, double *)
Definition: bezier_2d.c:118
int MMG2D_scaleMesh(MMG5_pMesh, MMG5_pSol, MMG5_pSol)
Definition: scalem_2d.c:48
int MMG2D_intmet_ani(MMG5_pMesh, MMG5_pSol, int, int8_t, int, double)
Definition: intmet_2d.c:209
int(* MMG2D_optlen)(MMG5_pMesh, MMG5_pSol, double, int)
double MMG2D_vfrac(MMG5_pMesh, MMG5_pSol, int, int)
Definition: mmg2d6.c:69
int MMG2D_zaldy(MMG5_pMesh mesh)
Definition: zaldy_2d.c:288
int MMG2D_colelt(MMG5_pMesh, MMG5_pSol, int)
Definition: mmg2d1.c:435
int MMG2D_parsop(MMG5_pMesh, MMG5_pSol)
Definition: libmmg2d_tools.c:84
int lissmet_iso(MMG5_pMesh mesh, MMG5_pSol sol)
int MMG2D_defsiz_iso(MMG5_pMesh, MMG5_pSol)
Definition: isosiz_2d.c:130
void MMG5_delEdge(MMG5_pMesh mesh, int iel)
Definition: zaldy_2d.c:70
int MMG2D_defaultmet_2d(MMG5_pMesh, MMG5_pSol, int, int8_t)
Definition: anisosiz_2d.c:50
int MMG2D_savemesh_db(MMG5_pMesh, char *, int8_t)
Definition: inout_2d.c:1686
double(* MMG2D_lencurv)(MMG5_pMesh, MMG5_pSol, int, int)
Definition: mmg2dexterns.c:6
int MMG2D_gradsiz_ani(MMG5_pMesh, MMG5_pSol)
int interp_ani(double *, double *, double *, double)
void MMG2D_solTruncatureForOptim(MMG5_pMesh mesh, MMG5_pSol met)
Definition: libmmg2d.c:56
int MMG2D_unscaleMesh(MMG5_pMesh, MMG5_pSol, MMG5_pSol)
int MMG2D_boulep(MMG5_pMesh, int, int, int *)
Definition: boulep_2d.c:34
int MMG2D_bdryEdge(MMG5_pMesh)
Definition: hash_2d.c:433
int MMG2D_chkmsh(MMG5_pMesh)
Definition: chkmsh_2d.c:204
int lissmet_ani(MMG5_pMesh mesh, MMG5_pSol sol)
Definition: lissmet_2d.c:47
int MMG2D_settagtriangles(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d2.c:99
int MMG2D_splitbar(MMG5_pMesh, int, int)
Definition: split_2d.c:716
double MMG2D_lencurv_ani(MMG5_pMesh, MMG5_pSol, int, int)
Definition: length_2d.c:81
int MMG2D_defsiz_ani(MMG5_pMesh, MMG5_pSol)
Definition: anisosiz_2d.c:360
int MMG2D_colver(MMG5_pMesh, int, int *)
Definition: colver_2d.c:270
int MMG2D_freeLocalPar(MMG5_pMesh)
Definition: libmmg2d_tools.c:198
int MMG2D_assignEdge(MMG5_pMesh)
Definition: hash_2d.c:333
double MMG2D_caltri_ani(MMG5_pMesh mesh, MMG5_pSol sol, MMG5_pTria)
Definition: quality_2d.c:95
void MMG2D_kiufree(pQueue q)
int MMG2D_kiuput(pQueue q, int iel)
int MMG2D_hashQuad(MMG5_pMesh mesh)
Definition: hash_2d.c:151
Queue * pQueue
Definition: mmg2d.h:118
int MMG2D_gradsizreq_ani(MMG5_pMesh, MMG5_pSol)
int MMG2D_savemet_db(MMG5_pMesh, MMG5_pSol, char *, int8_t)
Definition: inout_2d.c:1779
static const unsigned int MMG2D_idir[5]
Definition: mmg2d.h:130
int MMG2D_prilen(MMG5_pMesh, MMG5_pSol)
Definition: length_2d.c:124
int MMG2D_sum_reqEdgeLengthsAtPoint(MMG5_pMesh, MMG5_pSol, MMG5_pTria, int8_t)
Definition: isosiz_2d.c:50
int MMG2D_resetRef(MMG5_pMesh)
Definition: mmg2d6.c:191
int MMG2D_cutEdge(MMG5_pMesh, MMG5_pTria, MMG5_pPoint, MMG5_pPoint)
Definition: locate_2d.c:90
int MMG2D_ismaniball(MMG5_pMesh, MMG5_pSol, int, int8_t)
Definition: mmg2d6.c:221
int MMG2D_mmg2d9(MMG5_pMesh, MMG5_pSol, MMG5_pSol, int **)
Definition: mmg2d9.c:571
int MMG2D_adpcol(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d1.c:690
int MMG2D_swpmsh(MMG5_pMesh, MMG5_pSol, int)
Definition: mmg2d1.c:509
int MMG2D_grad2met_ani(MMG5_pMesh, MMG5_pSol, MMG5_pTria, int, int)
Definition: anisosiz_2d.c:574
static void MMG2D_Set_commonFunc(void)
Definition: mmg2d.h:375
static const unsigned int MMG2D_inxt[5]
Definition: mmg2d.h:131
static void MMG2D_excfun(int sigid)
Definition: mmg2d.h:97
int MMG2D_findTria(MMG5_pMesh, int)
Definition: locate_2d.c:215
int MMG2D_movintpt(MMG5_pMesh, MMG5_pSol, int, int *, int8_t)
Definition: movpt_2d.c:212
int MMG2D_snapval(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d6.c:312
int MMG2D_hashTria(MMG5_pMesh)
Definition: hash_2d.c:35
int MMG2D_analys(MMG5_pMesh)
Definition: analys_2d.c:761
int MMG2D_split3(MMG5_pMesh, MMG5_pSol, int, int vx[3])
Definition: split_2d.c:625
int MMG2D_indPt(MMG5_pMesh mesh, int kp)
Definition: tools_2d.c:69
int MMG2D_swapdelone(MMG5_pMesh, MMG5_pSol, int, int8_t, double, int *)
Definition: swapar_2d.c:39
double MMG2D_voltri(MMG5_pMesh, int, int, int)
Definition: mmg2d6.c:45
int MMG2D_coorbary(MMG5_pMesh, MMG5_pTria, double c[2], double *, double *, double *)
Definition: locate_2d.c:37
static const int MMG2D_iare[3][2]
Definition: mmg2d.h:128
void MMG2D_delPt(MMG5_pMesh mesh, int ip)
Definition: zaldy_2d.c:57
int MMG2D_bouleendp(MMG5_pMesh, int, int8_t, int *, int *)
Definition: boulep_2d.c:282
int MMG2D_norver(MMG5_pMesh, int)
Definition: analys_2d.c:400
int MMG2D_indElt(MMG5_pMesh mesh, int kel)
Definition: tools_2d.c:45
int MMG2D_cutEdgeTriangle(MMG5_pMesh, int, int, int)
Definition: locate_2d.c:129
int MMG2D_chkspl(MMG5_pMesh, MMG5_pSol, int, int8_t)
Definition: split_2d.c:50
int MMG2D_set_metricAtPointsOnReqEdges(MMG5_pMesh, MMG5_pSol, int8_t)
Definition: isosiz_2d.c:75
int MMG2D_split1(MMG5_pMesh, MMG5_pSol, int, int vx[3])
Definition: split_2d.c:385
int(* MMG2D_gradsizreq)(MMG5_pMesh, MMG5_pSol)
Definition: mmg2dexterns.c:7
int MMG2D_adpspl(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d1.c:634
double MMG2D_caltri_iso(MMG5_pMesh mesh, MMG5_pSol sol, MMG5_pTria)
Definition: quality_2d.c:59
int MMG2D_bdryenforcement(MMG5_pMesh, MMG5_pSol)
Definition: enforcement_2d.c:34
int MMG2D_split1b(MMG5_pMesh, int, int8_t, int)
Definition: split_2d.c:239
int MMG2D_regnor(MMG5_pMesh)
Definition: analys_2d.c:505
double MMG2D_lencurv_iso(MMG5_pMesh, MMG5_pSol, int, int)
Definition: length_2d.c:62
int MMG2D_mmg2d2(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d2.c:495
int MMG2D_mmg2d1(MMG5_pMesh, MMG5_pSol)
int MMG2D_chkedg(MMG5_pMesh, int)
Definition: bezier_2d.c:28
int MMG2D_baseBdry(MMG5_pMesh mesh)
int MMG2D_grad2metreq_ani(MMG5_pMesh, MMG5_pSol, MMG5_pTria, int, int)
Definition: anisosiz_2d.c:637
int MMG2D_isInTriangle(MMG5_pMesh, int, double c[2])
Definition: locate_2d.c:69
int MMG2D_anaelt(MMG5_pMesh, MMG5_pSol, int)
Definition: mmg2d1.c:110
int MMG2D_delElt(MMG5_pMesh mesh, int iel)
Definition: zaldy_2d.c:105
int MMG2D_cavity(MMG5_pMesh, MMG5_pSol, int, int *)
Definition: delone_2d.c:147
double long_ani(double *ca, double *cb, double *ma, double *mb)
Definition: length_2d.c:46
int MMG2D_split1_sim(MMG5_pMesh, MMG5_pSol, int, int vx[3])
Definition: split_2d.c:340
int MMG2D_singul(MMG5_pMesh, int)
Definition: analys_2d.c:267
int MMG2D_chkor(MMG5_pMesh)
Definition: chkmsh_2d.c:404
int(* MMG2D_intmet)(MMG5_pMesh, MMG5_pSol, int, int8_t, int, double)
Definition: mmg2dexterns.c:5
int(* MMG2D_defsiz)(MMG5_pMesh, MMG5_pSol)
Definition: mmg2dexterns.c:4
int MMG2D_setref_ls(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d6.c:945
int MMG2D_Free_all_var(va_list argptr)
Definition: variadic_2d.c:268
int MMG5_getnElt(MMG5_pMesh mesh, int n)
Definition: zaldy_2d.c:128
int MMG2D_velextLS(MMG5_pMesh, MMG5_pSol)
int MMG2D_anatri(MMG5_pMesh, MMG5_pSol, int8_t)
Definition: mmg2d1.c:39
int MMG2D_2dMeshCheck(MMG5_pMesh mesh)
Definition: inout_2d.c:595
int MMG2D_newElt(MMG5_pMesh mesh)
Definition: zaldy_2d.c:85
int MMG2D_newPt(MMG5_pMesh mesh, double c[2], int16_t tag)
Definition: zaldy_2d.c:38
int MMG2D_Free_structures_var(va_list argptr)
Definition: variadic_2d.c:377
int MMG2D_chkswp(MMG5_pMesh, MMG5_pSol, int, int8_t, int8_t)
Definition: swapar_2d.c:128
int MMG2D_boulet(MMG5_pMesh, int, int8_t, int *)
Definition: boulep_2d.c:229
int MMG2D_kiudel(pQueue q, int iel)
int MMG2D_swapar(MMG5_pMesh, int, int8_t)
Definition: swapar_2d.c:221
int optlen_iso_bar(MMG5_pMesh mesh, MMG5_pSol sol, double declic, int base)
int MMG2D_savedisp_db(MMG5_pMesh mesh, MMG5_pSol, char *, int8_t)
Definition: inout_2d.c:1894
int MMG2D_rmc(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d6.c:552
int MMG2D_outqua(MMG5_pMesh, MMG5_pSol)
Definition: quality_2d.c:156
int MMG5_interpmet22(MMG5_pMesh, double *, double *, double, double *)
Definition: intmet_2d.c:54
int MMG2D_Free_names_var(va_list argptr)
Definition: variadic_2d.c:492
int(* MMG2D_gradsiz)(MMG5_pMesh, MMG5_pSol)
Definition: mmg2dexterns.c:9
double long_iso(double *ca, double *cb, double *ma, double *mb)
Definition: length_2d.c:26
int optlen_ani(MMG5_pMesh mesh, MMG5_pSol sol, double declic, int base)
static const uint8_t MMG2D_idir_q[4][2]
idir[i]: vertices of edge i for a quad
Definition: mmg2d.h:76
int MMG2D_boulen(MMG5_pMesh, int, int8_t, int *, int *, double *)
Definition: boulep_2d.c:112
int MMG2D_Init_mesh_var(va_list argptr)
Definition: variadic_2d.c:175
int MMG2D_chkcol(MMG5_pMesh, MMG5_pSol, int, int8_t, int *, int8_t)
Definition: colver_2d.c:41
int MMG2D_movtri(MMG5_pMesh, MMG5_pSol, int, int8_t)
Definition: mmg2d1.c:747
void MMG2D_keep_only1Subdomain(MMG5_pMesh mesh, int nsd)
Definition: tools_2d.c:91
int MMG2D_pack(MMG5_pMesh, MMG5_pSol, MMG5_pSol)
Definition: hash_2d.c:514
double(* MMG2D_caltri)(MMG5_pMesh, MMG5_pSol, MMG5_pTria)
Definition: mmg2dexterns.c:8
int MMG2D_intmet_iso(MMG5_pMesh, MMG5_pSol, int, int8_t, int, double)
Definition: intmet_2d.c:38
int MMG2D_savenor_db(MMG5_pMesh, char *, int8_t)
Definition: inout_2d.c:1839
int MMG2D_adptri(MMG5_pMesh, MMG5_pSol)
Definition: mmg2d1.c:544
int MMG2D_mmg2d6(MMG5_pMesh, MMG5_pSol, MMG5_pSol)
Definition: mmg2d6.c:1009
int MMG2D_dichoto(MMG5_pMesh, MMG5_pSol, int, int *)
Definition: mmg2d1.c:359
int MMG2D_split2_sim(MMG5_pMesh, MMG5_pSol, int, int vx[3])
Definition: split_2d.c:451
int MMG2D_colver3(MMG5_pMesh, int *)
Definition: colver_2d.c:356
int MMG2D_defmetbdy_2d(MMG5_pMesh, MMG5_pSol, int, int8_t)
Definition: anisosiz_2d.c:85
int(* MMG5_chkmsh)(MMG5_pMesh, int, int)
Definition: mmgexterns.c:27
int(* MMG5_grad2metreq_ani)(MMG5_pMesh, MMG5_pSol, MMG5_pTria, int, int)
Definition: mmgexterns.c:33
int(* MMG5_grad2met_ani)(MMG5_pMesh, MMG5_pSol, MMG5_pTria, int, int)
Definition: mmgexterns.c:32
int size
Definition: mmg2d.h:121
int * head
Definition: mmg2d.h:122
int * link
Definition: mmg2d.h:123
MMG mesh structure.
Definition: libmmgtypes.h:575
Structure to store points of a MMG mesh.
Definition: libmmgtypes.h:252
Definition: libmmgtypes.h:633
Definition: libmmgtypes.h:312
int * stack
Definition: mmg2d.h:116
int cur
Definition: mmg2d.h:116