mmg3d
mmg3d.h
Go to the documentation of this file.
1 /* =============================================================================
2 ** This file is part of the mmg software package for the tetrahedral
3 ** mesh modification.
4 ** Copyright (c) Bx INP/CNRS/Inria/UBordeaux/UPMC, 2004-
5 **
6 ** mmg is free software: you can redistribute it and/or modify it
7 ** under the terms of the GNU Lesser General Public License as published
8 ** by the Free Software Foundation, either version 3 of the License, or
9 ** (at your option) any later version.
10 **
11 ** mmg is distributed in the hope that it will be useful, but WITHOUT
12 ** ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
13 ** FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
14 ** License for more details.
15 **
16 ** You should have received a copy of the GNU Lesser General Public
17 ** License and of the GNU General Public License along with mmg (in
18 ** files COPYING.LESSER and COPYING). If not, see
19 ** <http://www.gnu.org/licenses/>. Please read their terms carefully and
20 ** use this copy of the mmg distribution only if you accept them.
21 ** =============================================================================
22 */
23 
24 #ifndef MMG3D_H
25 #define MMG3D_H
26 
27 #include "libmmg3d.h"
28 #include "libmmgcommon.h"
29 
30 #ifdef __cplusplus
31 extern "C" {
32 #endif
33 
35 #define MMG5_RETURN_AND_FREE(mesh,met,ls,disp,val)do \
36  { \
37  if ( !MMG3D_Free_all(MMG5_ARG_start, \
38  MMG5_ARG_ppMesh,&mesh,MMG5_ARG_ppMet,&met, \
39  MMG5_ARG_ppLs,&ls,MMG5_ARG_ppDisp,&disp, \
40  MMG5_ARG_end) ) { \
41  return MMG5_LOWFAILURE; \
42  } \
43  return val; \
44  }while(0)
45 
48 #define MMG3D_POINT_REALLOC(mesh,sol,ip,wantedGap,law,o,tag,src ) do \
49  { \
50  int klink,oldnpmax; \
51  assert ( mesh && mesh->point ); \
52  \
53  oldnpmax = mesh->npmax; \
54  MMG5_TAB_RECALLOC(mesh,mesh->point,mesh->npmax,wantedGap,MMG5_Point, \
55  "larger point table",law); \
56  \
57  mesh->npnil = mesh->np+1; \
58  for (klink=mesh->npnil; klink<mesh->npmax-1; klink++) \
59  mesh->point[klink].tmp = klink+1; \
60  \
61  /* solution */ \
62  if ( sol ) { \
63  if ( sol->m ) { \
64  MMG5_ADD_MEM(mesh,(sol->size*(mesh->npmax-sol->npmax))*sizeof(double), \
65  "larger solution", \
66  MMG5_SAFE_RECALLOC(mesh->point,mesh->npmax+1,oldnpmax+1,MMG5_Point,,); \
67  mesh->memCur -= (mesh->npmax - oldnpmax)*sizeof(MMG5_Point); \
68  mesh->npmax = oldnpmax; \
69  mesh->np = mesh->npmax-1; \
70  mesh->npnil = 0; \
71  law); \
72  MMG5_SAFE_REALLOC(sol->m,sol->size*(sol->npmax+1), \
73  sol->size*(mesh->npmax+1), \
74  double,"larger solution", \
75  MMG5_SAFE_RECALLOC(mesh->point,mesh->npmax+1,oldnpmax+1,MMG5_Point,,); \
76  mesh->memCur -= (mesh->npmax - oldnpmax)*sizeof(MMG5_Point); \
77  mesh->npmax = oldnpmax; \
78  mesh->np = mesh->npmax-1; \
79  mesh->npnil = 0; \
80  law); \
81  } \
82  sol->npmax = mesh->npmax; \
83  } \
84  \
85  /* We try again to add the point */ \
86  ip = MMG3D_newPt(mesh,o,tag,src); \
87  if ( !ip ) { law; } \
88  }while(0)
89 
90 
93 #define MMG3D_TETRA_REALLOC(mesh,jel,wantedGap,law ) do \
94  { \
95  int klink,oldSiz; \
96  \
97  oldSiz = mesh->nemax; \
98  MMG5_TAB_RECALLOC(mesh,mesh->tetra,mesh->nemax,wantedGap,MMG5_Tetra, \
99  "larger tetra table",law); \
100  \
101  mesh->nenil = mesh->ne+1; \
102  for (klink=mesh->nenil; klink<mesh->nemax-1; klink++) \
103  mesh->tetra[klink].v[3] = klink+1; \
104  \
105  if ( mesh->adja ) { \
106  /* adja table */ \
107  MMG5_ADD_MEM(mesh,4*(mesh->nemax-oldSiz)*sizeof(int), \
108  "larger adja table",law); \
109  MMG5_SAFE_RECALLOC(mesh->adja,4*oldSiz+5,4*mesh->nemax+5,int \
110  ,"larger adja table",law); \
111  } \
112  \
113  /* We try again to add the point */ \
114  jel = MMG3D_newElt(mesh); \
115  if ( !jel ) {law;} \
116  }while(0)
117 
118 /* numerical accuracy */
119 #define MMG3D_ALPHAD 20.7846096908265 /* 12*sqrt(3) */
120 #define MMG3D_LLONG 2.5
121 #define MMG3D_LSHRT 0.3
122 #define MMG3D_LOPTL 1.3
123 #define MMG3D_LOPTS 0.6
124 
125 #define MMG3D_SWAP06 0.0288675 /* 0.6/MMG3D_ALPHAD */
126 #define MMG3D_SSWAPIMPROVE 1.053
127 #define MMG3D_LSWAPIMPROVE 1.1
128 #define MMG3D_DET2VOL 0.1666666666666667 /* 1/6 */
129 
130 #define MMG3D_BADKAL 0.2
131 #define MMG3D_MAXKAL 1.
132 
133 
134 #define MMG3D_NPMAX 1000000 //200000
135 #define MMG3D_NAMAX 200000 //40000
136 #define MMG3D_NTMAX 2000000 //400000
137 #define MMG3D_NEMAX 6000000 //1200000
138 
139 #define MMG3D_SHORTMAX 0x7fff
140 
141 #define MMG3D_VOLFRAC 1.e-5
142 
144 static const uint8_t MMG5_inxt3[7] = { 1,2,3,0,1,2,3 };
146 static const uint8_t MMG5_iprv3[7] = { 3,0,1,2,3,0,1 };
148 static const uint8_t MMG5_idir[4][3] = { {1,2,3}, {0,3,2}, {0,1,3}, {0,2,1} };
149 /* \var idirinv[i][j]: num of the jth point in the ith face */
150 static const int8_t MMG5_idirinv[4][4] = {{-1,0,1,2},{0,-1,2,1},{0,1,-1,2},{0,2,1,-1}};
152 static const int8_t MMG5_iarf[4][3] = { {5,4,3}, {5,1,2}, {4,2,0}, {3,0,1} };
154 static const int8_t MMG5_iarfinv[4][6] = { {-1,-1,-1,2,1,0}, {-1,1,2,-1,-1,0},{2,-1,1,-1,0,-1},{1,2,-1,0,-1,-1}};
156 static const uint8_t MMG5_iare[6][2] = { {0,1}, {0,2}, {0,3}, {1,2}, {1,3}, {2,3} };
158 static const uint8_t MMG5_ifar[6][2] = { {2,3}, {1,3}, {1,2}, {0,3}, {0,2}, {0,1} };
160 static const uint8_t MMG5_isar[6][2] = { {2,3}, {3,1}, {1,2}, {0,3}, {2,0}, {0,1} };
162 static const uint8_t MMG5_arpt[4][3] = { {0,1,2}, {0,4,3}, {1,3,5}, {2,5,4} };
163 
165 static const uint8_t MMG5_idir_pr[5][4] = { {0,1,2,0},{3,5,4,3},{1,4,5,2},{0,2,5,3},{0,3,4,1} };
167 static const uint8_t MMG5_iarf_pr[5][5] = { {0,1,3,0}, {6,8,7,6}, {3,5,8,4}, {5,1,2,7},{0,4,6,2} };
168 
185 static const uint8_t MMG5_permedge[12][6] = {
186  {0,1,2,3,4,5}, {1,2,0,5,3,4}, {2,0,1,4,5,3}, {0,4,3,2,1,5},
187  {3,0,4,1,5,2}, {4,3,0,5,2,1}, {1,3,5,0,2,4}, {3,5,1,4,0,2},
188  {5,1,3,2,4,0}, {2,5,4,1,0,3}, {4,2,5,0,3,1}, {5,4,2,3,1,0} };
189 
195 typedef struct MMG3D_PROctree_s
196 {
198  int* v;
199  int nbVer;
200  int depth;
202 
207 typedef struct
208 {
209  int nv;
210  int nc;
214 
215 
216 /* PROctree */
221 int MMG3D_isCellIncluded(double* cellCenter, double l, double* zoneCenter, double l0);
222 void MMG3D_placeInListDouble(double*, double, int, int);
224 int MMG3D_seekIndex (double* distList, double dist, int indexMin, int indexMax);
225 int MMG3D_intersectRect(double *rectin, double *rectinout);
226 int MMG3D_getListSquareRec(MMG3D_PROctree_s*,double*,double*,
227  MMG3D_PROctree_s***,double*,double*,double, int, int, int*);
229 int MMG3D_addPROctreeRec(MMG5_pMesh,MMG3D_PROctree_s*,double*, const int, int);
232 int MMG3D_movePROctree(MMG5_pMesh, MMG3D_pPROctree,int, double*, double*);
234 void MMG3D_mergeBranches(MMG5_pMesh mesh,MMG3D_PROctree_s* q, int dim, int nv);
235 int MMG3D_delPROctreeRec(MMG5_pMesh,MMG3D_PROctree_s*,double*, const int,const int);
236 int MMG3D_delPROctree(MMG5_pMesh mesh, MMG3D_pPROctree q, const int no);
237 void MMG3D_printArbreDepth(MMG3D_PROctree_s* q, int depth, int nv, int dim);
239 void MMG3D_sizeArbreRec(MMG3D_PROctree_s* q, int nv, int dim, int*,int*);
240 int* MMG3D_sizeArbre(MMG3D_PROctree* q, int dim);
243 int64_t MMG3D_getPROctreeCoordinate(MMG3D_pPROctree q, double* ver, int dim);
244 
245 /* prototypes */
246 int MMG3D_tetraQual(MMG5_pMesh mesh, MMG5_pSol met,int8_t metRidTyp);
247 extern int MMG5_directsurfball(MMG5_pMesh mesh, int ip, int *list, int ilist, double n[3]);
248 
249 int MMG3D_Init_mesh_var( va_list argptr );
250 int MMG3D_Free_all_var( va_list argptr );
251 int MMG3D_Free_structures_var( va_list argptr );
252 int MMG3D_Free_names_var( va_list argptr );
254 int MMG3D_newPt(MMG5_pMesh mesh,double c[3],int16_t tag,int src);
256 int MMG3D_delElt(MMG5_pMesh mesh,int iel);
257 void MMG3D_delPt(MMG5_pMesh mesh,int ip);
264 extern int MMG5_BezierTgt(double c1[3],double c2[3],double n1[3],double n2[3],double t1[3],double t2[3]);
265 extern double MMG5_BezierGeod(double c1[3], double c2[3], double t1[3], double t2[3]);
266 int MMG3D_bezierInt(MMG5_pBezier pb,double uv[2],double o[3],double no[3],double to[3]);
267 extern int MMG5_BezierReg(MMG5_pMesh mesh,int ip0, int ip1, double s, double v[3], double *o, double *no);
268 extern int MMG5_BezierRef(MMG5_pMesh mesh,int ip0, int ip1, double s, double *o, double *no, double *to);
269 extern int MMG5_BezierEdge(MMG5_pMesh mesh,int ip0, int ip1, double b0[3], double b1[3],int8_t isrid, double v[3]);
270 extern int MMG5_BezierRidge(MMG5_pMesh mesh,int ip0, int ip1, double s, double *o, double *no1, double *no2, double *to);
271 extern int MMG5_BezierNom(MMG5_pMesh mesh,int ip0,int ip1,double s,double *o,double *no,double *to);
272 int MMG5_norface(MMG5_pMesh mesh ,int k, int iface, double v[3]);
273 int MMG3D_findEdge(MMG5_pMesh,MMG5_pTetra,int,int,int,int,int8_t*,int8_t* );
274 int MMG5_boulernm (MMG5_pMesh mesh,MMG5_Hash *hash, int start, int ip, int *ng, int *nr);
275 int MMG5_boulenm(MMG5_pMesh mesh, int start, int ip, int iface, double n[3],double t[3]);
276 int MMG5_boulenmInt(MMG5_pMesh mesh, int start, int ip, double t[3]);
277 int MMG5_boulevolp(MMG5_pMesh mesh, int start, int ip, int * list);
278 int MMG5_boulesurfvolpNom(MMG5_pMesh mesh,int start,int ip,int iface,int *listv,
279  int *ilistv,int *lists,int*ilists,int*refmin,int*refplus,int isnm);
280 int MMG5_boulesurfvolp(MMG5_pMesh mesh,int start,int ip,int iface,int *listv,
281  int *ilistv,int *lists,int*ilists, int isnm);
282 int MMG5_bouletrid(MMG5_pMesh,int,int,int,int *,int *,int *,int *,int *,int *);
283 int MMG5_startedgsurfball(MMG5_pMesh mesh,int nump,int numq,int *list,int ilist);
284 int MMG5_srcbdy(MMG5_pMesh mesh,int start,int ia);
285 int MMG5_coquil(MMG5_pMesh mesh, int start, int ia, int * list);
286 int MMG5_coquilface(MMG5_pMesh mesh, int start,int8_t iface,int,int*,int*,int*,int);
287 int MMG3D_coquilFaceFirstLoop(MMG5_pMesh mesh,int start,int na,int nb,int8_t iface,
288  int8_t ia,int *list,int *ilist,int *it1,int *it2,
289  int *piv,int *adj,int8_t *hasadja,int *nbdy,int silent);
290 void MMG3D_coquilFaceSecondLoopInit(MMG5_pMesh mesh,int piv,int8_t *iface,int8_t *i,
291  int *list,int *ilist,int *it1,int *pradj,
292  int *adj);
293 void MMG5_coquilFaceErrorMessage(MMG5_pMesh mesh, int k1, int k2);
294 int16_t MMG5_coquilTravel(MMG5_pMesh,int,int,int*,int*,int8_t*,int8_t*);
295 int16_t MMG5_openCoquilTravel(MMG5_pMesh,int,int,int*,int*,int8_t*,int8_t*);
296 int MMG5_settag(MMG5_pMesh,int,int,int16_t,int);
297 int MMG5_deltag(MMG5_pMesh,int,int,int16_t);
299 int MMG5_chkcol_int(MMG5_pMesh,MMG5_pSol,int,int8_t,int8_t,int*,int,int8_t);
300 int MMG5_chkcol_bdy(MMG5_pMesh,MMG5_pSol,int,int8_t,int8_t,int*,int,int*,int,int,int,int8_t,int);
301 int MMG5_chkcol_nomint(MMG5_pMesh,MMG5_pSol,int,int8_t,int8_t,int*,int,int8_t);
302 int MMG5_chkmanicoll(MMG5_pMesh,int,int,int,int,int,int,int,int8_t,int8_t);
304 int MMG5_colver(MMG5_pMesh,MMG5_pSol,int *,int,int8_t,int8_t);
314 int MMG5_hashFace(MMG5_pMesh,MMG5_Hash*,int,int,int,int);
315 int MMG5_hashGetFace(MMG5_Hash*,int,int,int);
318 int MMG5_hashPop(MMG5_Hash *hash,int a,int b);
319 int MMG5_hPop(MMG5_HGeom *hash,int a,int b,int *ref,int16_t *tag);
320 int MMG5_hTag(MMG5_HGeom *hash,int a,int b,int ref,int16_t tag);
321 int MMG5_hGet(MMG5_HGeom *hash,int a,int b,int *ref,int16_t *tag);
322 int MMG5_hEdge(MMG5_pMesh mesh,MMG5_HGeom *hash,int a,int b,int ref,int16_t tag);
323 int MMG5_hNew(MMG5_pMesh mesh,MMG5_HGeom *hash,int hsiz,int hmax);
330 int MMG5_cntbdypt(MMG5_pMesh mesh, int nump);
331 size_t MMG5_memSize(void);
339 int MMG5_mmg3dChkmsh(MMG5_pMesh,int,int);
340 int MMG3D_setMeshSize_initData(MMG5_pMesh,int,int,int,int,int,int);
342 int MMG3D_split1_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
343 int MMG5_split1(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t metRidTyp);
344 int MMG5_split1b(MMG5_pMesh,MMG5_pSol,int*,int,int,int,int8_t,int8_t);
345 int MMG5_splitedg(MMG5_pMesh mesh, MMG5_pSol met,int iel, int iar, double crit);
346 int MMG3D_split2sf_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
347 int MMG5_split2sf(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t);
348 int MMG3D_split2_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
349 int MMG5_split2(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t);
350 int MMG3D_split3_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
351 int MMG5_split3(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t);
352 int MMG3D_split3cone_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
353 int MMG5_split3cone(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t);
354 int MMG3D_split3op_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
355 int MMG5_split3op(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6],int8_t);
356 int MMG3D_split4sf_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
357 int MMG5_split4sf(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t);
358 int MMG3D_split4op_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
359 int MMG5_split4op(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t);
360 int MMG3D_split5_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
361 int MMG5_split5(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t);
362 int MMG3D_split6_sim(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6]);
363 int MMG5_split6(MMG5_pMesh mesh,MMG5_pSol met,int k,int vx[6],int8_t);
364 int MMG5_split4bar(MMG5_pMesh mesh,MMG5_pSol met,int k,int8_t);
365 int MMG3D_simbulgept(MMG5_pMesh mesh,MMG5_pSol met, int *list, int ilist,int);
366 void MMG5_nsort(int ,double *,int8_t *);
368 int MMG5_movintpt_iso(MMG5_pMesh ,MMG5_pSol,MMG3D_pPROctree, int *, int , int);
369 int MMG3D_movnormal_iso(MMG5_pMesh ,MMG5_pSol ,int ,int );
373  int*, int, int*, int, int ,int);
375  int*, int, int*, int, int ,int);
377  int*, int ,int);
379  int*, int ,int);
381  int*, int ,int);
384  int*, int ,int);
386  int*, int ,int);
388  int*, int ,int);
389 int MMG3D_movv_ani(MMG5_pMesh ,MMG5_pSol ,int ,int );
390 int MMG3D_movv_iso(MMG5_pMesh ,MMG5_pSol ,int ,int );
391 int MMG3D_normalAdjaTri(MMG5_pMesh,int,int8_t,int,double n[3]);
392 int MMG5_chkswpbdy(MMG5_pMesh, MMG5_pSol,int*, int, int, int,int8_t);
393 int MMG5_swpbdy(MMG5_pMesh,MMG5_pSol,int*,int,int,MMG3D_pPROctree,int8_t);
394 int MMG5_swpgen(MMG5_pMesh,MMG5_pSol,int, int, int*,MMG3D_pPROctree,int8_t);
395 int MMG5_chkswpgen(MMG5_pMesh,MMG5_pSol,int,int,int*,int*,double,int8_t);
396 int MMG3D_swap23(MMG5_pMesh,MMG5_pSol,int,int8_t,int,int,int,int);
397 int MMG5_srcface(MMG5_pMesh mesh,int n0,int n1,int n2);
398 int MMG5_chkptonbdy(MMG5_pMesh,int);
399 double MMG5_orcal_poi(double a[3],double b[3],double c[3],double d[3]);
400 int MMG5_countelt(MMG5_pMesh mesh,MMG5_pSol sol, double *weightelt, long *npcible);
401 /*function for agressive optimization*/
404 int MMG3D_splitItem(MMG5_pMesh , MMG5_pSol ,MMG3D_pPROctree ,int ,int ,double );
407 
408 int MMG5_trydisp(MMG5_pMesh,double *,short);
411 int MMG5_ppgdisp(MMG5_pMesh,double *);
413 int MMG3D_displayQualHisto(int,double,double,double,int,int,int,
414  int his[5],int,int,int);
415 int MMG3D_displayQualHisto_internal(int,double,double,double,int,int,int,
416  int his[5],int,int,int);
417 void MMG3D_computeInqua(MMG5_pMesh,MMG5_pSol,int*,double*,double*,double*,int*,int*,
418  int*,int his[5],int);
420 void MMG3D_computeOutqua(MMG5_pMesh,MMG5_pSol,int*,double*,double*,double*,int*,int*,
421  int*,int his[5],int*,int);
423 void MMG3D_computeLESqua(MMG5_pMesh,MMG5_pSol,int*,double*,double*,double*,int*,int*,
424  int*,int his[5],int);
425 int MMG3D_computePrilen(MMG5_pMesh,MMG5_pSol,double*,double*,double*,int*,int*,int*,
426  int*,int*,int*,int8_t,double**, int [9] );
427 int MMG3D_prilen(MMG5_pMesh mesh,MMG5_pSol met,int8_t);
429 int MMG5_intridmet(MMG5_pMesh,MMG5_pSol,int,int,double,double*,double*);
430 int MMG5_intregmet(MMG5_pMesh,MMG5_pSol,int,int8_t,double, double*);
431 int MMG5_intvolmet(MMG5_pMesh,MMG5_pSol,int,int8_t,double, double*);
432 int MMG3D_localParamReg(MMG5_pMesh,int,int*,int,int*,int,double*,double*,double*);
433 int MMG3D_localParamNm(MMG5_pMesh,int,int,int,double*,double*,double*);
434 int MMG3D_mark_packedPoints(MMG5_pMesh mesh,int *np,int *nc);
445 
446 /* rmc option */
447 double MMG3D_vfrac(MMG5_pMesh ,MMG5_pSol ,int ,int );
449 int MMG5_isbr(MMG5_pMesh ,int );
450 
451 
452 /* tools_3d.c */
453 void MMG3D_keep_only1Subdomain ( MMG5_pMesh mesh,int nsd );
454 
455 /* useful functions to debug */
456 int MMG3D_indElt(MMG5_pMesh mesh,int kel);
457 int MMG3D_indPt(MMG5_pMesh mesh,int kp);
458 void MMG5_printTetra(MMG5_pMesh mesh,char* fileName);
459 
460 
461 #ifdef USE_SCOTCH
462 int MMG5_mmg3dRenumbering(int,MMG5_pMesh,MMG5_pSol,MMG5_pSol,int*);
463 #endif
464 
465 int MMG5_meancur(MMG5_pMesh mesh,int np,double c[3],int ilist,int *list,double h[3]);
466 double MMG5_surftri(MMG5_pMesh,int,int);
470 
471 /* Lagrangian mode functions */
472 double MMG5_estavglen(MMG5_pMesh);
473 int MMG5_stiffelt(MMG5_pMesh,int,double*,double*);
477 
478 /* Delaunay functions*/
479 int MMG5_delone(MMG5_pMesh mesh,MMG5_pSol sol,int ip,int *list,int ilist);
480  int MMG5_cavity_iso(MMG5_pMesh mesh,MMG5_pSol sol,int iel,int ip,int *list,int lon,double volmin);
481  int MMG5_cavity_ani(MMG5_pMesh mesh,MMG5_pSol sol,int iel,int ip,int *list,int lon,double volmin);
482 int MMG5_cenrad_iso(MMG5_pMesh mesh,double *ct,double *c,double *rad);
483 int MMG5_cenrad_ani(MMG5_pMesh mesh,double *ct,double *m,double *c,double *rad);
484 
485 /* mmg3d1.c */
486 void MMG5_tet2tri(MMG5_pMesh mesh,int k,int8_t ie,MMG5_Tria *ptt);
487 int MMG3D_dichoto(MMG5_pMesh mesh,MMG5_pSol met,int k,int *vx);
488 int MMG3D_dichoto1b(MMG5_pMesh mesh,MMG5_pSol met,int *list,int ret,int);
489 int8_t MMG5_chkedg(MMG5_pMesh mesh,MMG5_Tria *pt,int8_t ori,double,double,int);
491  int8_t,int8_t,int* );
492 int MMG5_anatet(MMG5_pMesh mesh,MMG5_pSol met, int8_t typchk, int patternMode) ;
494  double clickSurf,double clickVol,int moveVol,int improveSurf,int improveVolSurf,
495  int improveVol,int maxit,int testmark);
496 int MMG5_swpmsh(MMG5_pMesh mesh,MMG5_pSol met,MMG3D_pPROctree PROctree, int);
497  int MMG5_swptet(MMG5_pMesh mesh,MMG5_pSol met,double,double,MMG3D_pPROctree, int,int);
498 
499 /* pointers */
500 /* init structures */
502 /* iso/aniso computations */
504 extern double MMG5_lenedgCoor_iso(double*, double*, double*, double*);
505 int MMG5_intmet_iso(MMG5_pMesh,MMG5_pSol,int,int8_t,int, double);
506 int MMG5_intmet_ani(MMG5_pMesh,MMG5_pSol,int,int8_t,int, double);
507 int MMG3D_intmet33_ani(MMG5_pMesh,MMG5_pSol,int,int8_t,int, double);
517 double MMG5_meansizreg_iso(MMG5_pMesh,MMG5_pSol,int,int*,int,double,double);
519 extern int MMG5_moymet(MMG5_pMesh ,MMG5_pSol ,MMG5_pTetra ,double *);
522 
523 extern double (*MMG5_lenedg)(MMG5_pMesh ,MMG5_pSol ,int, MMG5_pTetra );
524 extern double (*MMG5_lenedgspl)(MMG5_pMesh ,MMG5_pSol ,int, MMG5_pTetra );
525 extern double (*MMG5_caltet)(MMG5_pMesh mesh,MMG5_pSol met,MMG5_pTetra pt);
526 extern double (*MMG5_caltri)(MMG5_pMesh mesh,MMG5_pSol met,MMG5_pTria ptt);
527 extern int (*MMG3D_defsiz)(MMG5_pMesh ,MMG5_pSol );
528 extern int (*MMG3D_gradsiz)(MMG5_pMesh ,MMG5_pSol );
529 extern int (*MMG3D_gradsizreq)(MMG5_pMesh ,MMG5_pSol );
530 extern int (*MMG5_intmet)(MMG5_pMesh,MMG5_pSol,int,int8_t,int, double);
531 extern int (*MMG5_interp4bar)(MMG5_pMesh,MMG5_pSol,int,int,double *);
532 extern int (*MMG5_movintpt)(MMG5_pMesh ,MMG5_pSol, MMG3D_pPROctree ,int *, int , int );
533 extern int (*MMG5_movbdyregpt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree ,int*, int, int*, int, int ,int);
534 extern int (*MMG5_movbdyrefpt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree ,int*, int, int*, int ,int);
535 extern int (*MMG5_movbdynompt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree ,int*, int, int*, int ,int);
536 extern int (*MMG5_movbdyridpt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree ,int*, int, int*, int ,int);
537 extern int (*MMG5_cavity)(MMG5_pMesh ,MMG5_pSol ,int ,int ,int *,int ,double);
538 extern int (*MMG3D_PROctreein)(MMG5_pMesh ,MMG5_pSol ,MMG3D_pPROctree ,int,double );
539 
540 /* input */
541 int MMG3D_openMesh(int imprim,const char *filename,FILE **inm,int *bin,char*,char*);
542 int MMG3D_loadMesh_opened(MMG5_pMesh mesh,FILE *inm,int bin);
543 
550 static inline
552  if ( mesh->xt ) {
553  if ( mesh->xt != mesh->ne ) {
554  fprintf(stderr,"\n ## Warning: %s: %d tetra on %d reoriented.\n",
555  __func__,mesh->xt,mesh->ne);
556  fprintf(stderr," Your mesh may be non-conform.\n");
557  }
558  else {
559  fprintf(stderr,"\n ## Warning: %s: all tetra reoriented.\n",__func__);
560  }
561  }
562  mesh->xt = 0;
563 }
564 
569 static inline
577 #ifdef USE_SCOTCH
578  MMG5_renumbering = MMG5_mmg3dRenumbering;
579 #endif
580 }
581 
582 #ifdef __cplusplus
583 }
584 #endif
585 
586 #endif
MMG5_pMesh char * filename
Definition: API_functionsf_3d.c:1015
MMG5_pMesh * mesh
Definition: API_functionsf_3d.c:65
const int va_list argptr
Definition: API_functionsf_3d.c:977
int MMG5_grad2metSurfreq(MMG5_pMesh mesh, MMG5_pSol met, MMG5_pTria pt, int npmaster, int npslave)
Definition: anisosiz.c:1448
int MMG5_grad2metSurf(MMG5_pMesh mesh, MMG5_pSol met, MMG5_pTria pt, int np1, int np2)
Definition: anisosiz.c:910
API headers for the mmg3d library.
API header for the common part of the MMG libraries.
MMG5_Sol * MMG5_pSol
Definition: libmmgtypes.h:623
MMG5_Mesh * MMG5_pMesh
Definition: libmmgtypes.h:604
MMG5_Tetra * MMG5_pTetra
Definition: libmmgtypes.h:381
int MMG5_paktet(MMG5_pMesh mesh)
Definition: hash_3d.c:49
void MMG3D_sizeArbreRec(MMG3D_PROctree_s *q, int nv, int dim, int *, int *)
Definition: PRoctree_3d.c:1075
int MMG5_intregmet(MMG5_pMesh, MMG5_pSol, int, int8_t, double, double *)
int MMG5_hTag(MMG5_HGeom *hash, int a, int b, int ref, int16_t tag)
Definition: hash_3d.c:854
int MMG5_lapantilap(MMG5_pMesh, double *)
int MMG5_boulesurfvolpNom(MMG5_pMesh mesh, int start, int ip, int iface, int *listv, int *ilistv, int *lists, int *ilists, int *refmin, int *refplus, int isnm)
Definition: boulep_3d.c:755
int MMG5_hashGetFace(MMG5_Hash *, int, int, int)
Definition: hash_3d.c:143
int MMG5_dichodisp(MMG5_pMesh, double *)
int MMG3D_chk4ridVertices(MMG5_pMesh mesh, MMG5_pTetra pt)
Definition: anisosiz_3d.c:39
int MMG5_split1(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t metRidTyp)
Definition: split_3d.c:116
int MMG3D_bezierInt(MMG5_pBezier pb, double uv[2], double o[3], double no[3], double to[3])
Definition: bezier_3d.c:571
int MMG5_cavity_iso(MMG5_pMesh mesh, MMG5_pSol sol, int iel, int ip, int *list, int lon, double volmin)
Definition: delaunay_3d.c:741
int MMG5_srcface(MMG5_pMesh mesh, int n0, int n1, int n2)
void MMG3D_placeInListDouble(double *, double, int, int)
Definition: PRoctree_3d.c:296
static void MMG3D_Set_commonFunc(void)
Definition: mmg3d.h:570
double MMG5_lenedgCoor_iso(double *, double *, double *, double *)
Compute edge length from edge's coordinates.
Definition: isosiz_3d.c:60
int MMG5_moymet(MMG5_pMesh, MMG5_pSol, MMG5_pTetra, double *)
Definition: anisosiz_3d.c:70
int MMG5_cntbdypt(MMG5_pMesh mesh, int nump)
Definition: chkmsh_3d.c:363
int MMG3D_gradsizreq_iso(MMG5_pMesh, MMG5_pSol)
Definition: isosiz_3d.c:1143
int MMG5_intmet_ani(MMG5_pMesh, MMG5_pSol, int, int8_t, int, double)
Definition: intmet_3d.c:51
int(* MMG3D_defsiz)(MMG5_pMesh, MMG5_pSol)
Definition: mmg3dexterns.c:9
int MMG5_ppgdisp(MMG5_pMesh, double *)
int MMG3D_movnormal_iso(MMG5_pMesh, MMG5_pSol, int, int)
Definition: movpt_3d.c:2066
void MMG3D_set_reqBoundaries(MMG5_pMesh mesh)
Definition: analys_3d.c:45
size_t MMG5_memSize(void)
Definition: tools.c:547
int MMG3D_pack_sol(MMG5_pMesh mesh, MMG5_pSol sol)
Definition: libmmg3d.c:420
static const int8_t MMG5_idirinv[4][4]
Definition: mmg3d.h:150
int MMG3D_simbulgept(MMG5_pMesh mesh, MMG5_pSol met, int *list, int ilist, int)
Definition: split_3d.c:322
int MMG5_movbdynompt_ani(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: anisomovpt_3d.c:1028
int MMG3D_split4sf_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:3655
int MMG3D_gradsizreq_ani(MMG5_pMesh, MMG5_pSol)
Definition: anisosiz_3d.c:2025
void MMG3D_mergeBranches(MMG5_pMesh mesh, MMG3D_PROctree_s *q, int dim, int nv)
Definition: PRoctree_3d.c:862
int MMG5_split3(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:1793
int MMG3D_seekIndex(double *distList, double dist, int indexMin, int indexMax)
Definition: PRoctree_3d.c:335
int MMG5_cavity_ani(MMG5_pMesh mesh, MMG5_pSol sol, int iel, int ip, int *list, int lon, double volmin)
Definition: delaunay_3d.c:582
void MMG3D_coquilFaceSecondLoopInit(MMG5_pMesh mesh, int piv, int8_t *iface, int8_t *i, int *list, int *ilist, int *it1, int *pradj, int *adj)
Definition: boulep_3d.c:1727
int MMG5_hGeom(MMG5_pMesh mesh)
Definition: hash_3d.c:1047
int MMG3D_hashTria(MMG5_pMesh mesh, MMG5_Hash *)
Definition: hash_3d.c:779
int MMG3D_defsiz_iso(MMG5_pMesh, MMG5_pSol)
Definition: isosiz_3d.c:641
int MMG3D_swap23(MMG5_pMesh, MMG5_pSol, int, int8_t, int, int, int, int)
Definition: swap_3d.c:566
void MMG5_coquilFaceErrorMessage(MMG5_pMesh mesh, int k1, int k2)
Definition: boulep_3d.c:1559
void MMG3D_keep_only1Subdomain(MMG5_pMesh mesh, int nsd)
Definition: tools_3d.c:1314
int MMG3D_localParamReg(MMG5_pMesh, int, int *, int, int *, int, double *, double *, double *)
Definition: tools_3d.c:950
int MMG5_velextLS(MMG5_pMesh, MMG5_pSol)
int MMG5_hGet(MMG5_HGeom *hash, int a, int b, int *ref, int16_t *tag)
Definition: hash_3d.c:945
int MMG3D_set_metricAtPointsOnReqEdges(MMG5_pMesh, MMG5_pSol, int8_t)
Definition: isosiz_3d.c:544
int(* MMG5_movbdynompt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: mmg3dexterns.c:17
int(* MMG5_movbdyrefpt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: mmg3dexterns.c:16
int MMG3D_update_xtetra(MMG5_pMesh mesh)
Definition: mmg3d2.c:1296
int MMG5_movintptLES_iso(MMG5_pMesh mesh, MMG5_pSol met, MMG3D_pPROctree, int *, int, int)
Definition: movpt_3d.c:194
int MMG5_intridmet(MMG5_pMesh, MMG5_pSol, int, int, double, double *, double *)
int MMG3D_delPROctreeRec(MMG5_pMesh, MMG3D_PROctree_s *, double *, const int, const int)
Definition: PRoctree_3d.c:894
MMG3D_PROctree * MMG3D_pPROctree
Definition: mmg3d.h:213
int MMG3D_addPROctree(MMG5_pMesh mesh, MMG3D_PROctree *q, const int no)
int MMG5_interp4bar_iso(MMG5_pMesh, MMG5_pSol, int, int, double *)
int MMG5_deltag(MMG5_pMesh, int, int, int16_t)
Definition: boulep_3d.c:1266
static const uint8_t MMG5_isar[6][2]
isar[i][]: vertices of extremities of the edge opposite to the ith edge
Definition: mmg3d.h:160
int MMG5_movbdyregpt_ani(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int, int)
Definition: anisomovpt_3d.c:195
int MMG3D_rmc(MMG5_pMesh, MMG5_pSol)
Definition: mmg3d2.c:767
void MMG3D_freePROctree(MMG5_pMesh, MMG3D_PROctree **q)
int MMG5_mmg3d1_pattern(MMG5_pMesh, MMG5_pSol, int *)
Definition: mmg3d1_pattern.c:454
void MMG3D_Free_topoTables(MMG5_pMesh mesh)
Definition: libmmg3d.c:65
int MMG5_swptet(MMG5_pMesh mesh, MMG5_pSol met, double, double, MMG3D_pPROctree, int, int)
Definition: mmg3d1.c:611
int MMG3D_optbdry(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int)
Definition: optbdry_3d.c:250
int MMG3D_PROctreein_ani(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int, double)
Definition: PRoctree_3d.c:1220
int MMG3D_mark_packedPoints(MMG5_pMesh mesh, int *np, int *nc)
Definition: libmmg3d.c:197
struct MMG3D_PROctree_s MMG3D_PROctree_s
int MMG5_BezierTgt(double c1[3], double c2[3], double n1[3], double n2[3], double t1[3], double t2[3])
Definition: bezier_3d.c:53
int MMG3D_swpItem(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int, int)
Definition: opttyp_3d.c:326
double MMG5_meansizreg_iso(MMG5_pMesh, MMG5_pSol, int, int *, int, double, double)
Definition: isosiz_3d.c:462
int MMG3D_intersectRect(double *rectin, double *rectinout)
Definition: PRoctree_3d.c:383
double(* MMG5_lenedg)(MMG5_pMesh, MMG5_pSol, int, MMG5_pTetra)
Definition: mmg3dexterns.c:5
void MMG3D_placeInListPROctree(MMG3D_PROctree_s **, MMG3D_PROctree_s *, int, int)
Definition: PRoctree_3d.c:314
int MMG3D_indElt(MMG5_pMesh mesh, int kel)
Definition: tools_3d.c:868
int MMG5_interp4bar_ani(MMG5_pMesh, MMG5_pSol, int, int, double *)
int MMG3D_delPROctreeVertex(MMG5_pMesh, MMG3D_PROctree_s *q, int no)
Definition: PRoctree_3d.c:798
int MMG5_split1b(MMG5_pMesh, MMG5_pSol, int *, int, int, int, int8_t, int8_t)
Definition: split_3d.c:512
int MMG5_swpmsh(MMG5_pMesh mesh, MMG5_pSol met, MMG3D_pPROctree PROctree, int)
Definition: mmg3d1.c:542
int MMG3D_indPt(MMG5_pMesh mesh, int kp)
Definition: tools_3d.c:884
int MMG5_norface(MMG5_pMesh mesh, int k, int iface, double v[3])
Definition: tools_3d.c:69
int MMG5_hNew(MMG5_pMesh mesh, MMG5_HGeom *hash, int hsiz, int hmax)
Definition: hash_3d.c:1019
void MMG3D_initPROctree_s(MMG3D_PROctree_s *q)
Definition: PRoctree_3d.c:46
int MMG5_meancur(MMG5_pMesh mesh, int np, double c[3], int ilist, int *list, double h[3])
int MMG5_anatet(MMG5_pMesh mesh, MMG5_pSol met, int8_t typchk, int patternMode)
Definition: mmg3d1.c:2689
void MMG3D_printArbreDepth(MMG3D_PROctree_s *q, int depth, int nv, int dim)
Definition: PRoctree_3d.c:1003
int MMG3D_splsurfedge(MMG5_pMesh, MMG5_pSol, int, MMG5_pTetra, MMG5_pxTetra, int8_t, int8_t, int8_t, int *)
Definition: mmg3d1.c:1449
int MMG3D_analys(MMG5_pMesh mesh)
Definition: analys_3d.c:761
int MMG3D_hashPrism(MMG5_pMesh mesh)
Definition: hash_3d.c:297
int MMG3D_dichoto1b(MMG5_pMesh mesh, MMG5_pSol met, int *list, int ret, int)
Definition: mmg3d1.c:240
int MMG5_split6(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:4834
void MMG5_freeXPrisms(MMG5_pMesh mesh)
Definition: zaldy_3d.c:346
int MMG5_directsurfball(MMG5_pMesh mesh, int ip, int *list, int ilist, double n[3])
Definition: tools_3d.c:83
int MMG3D_pack_tetra(MMG5_pMesh mesh)
Definition: libmmg3d.c:332
int(* MMG5_cavity)(MMG5_pMesh, MMG5_pSol, int, int, int *, int, double)
Definition: mmg3dexterns.c:19
int MMG5_boulenm(MMG5_pMesh mesh, int start, int ip, int iface, double n[3], double t[3])
Definition: boulep_3d.c:155
int MMG5_chkcol_int(MMG5_pMesh, MMG5_pSol, int, int8_t, int8_t, int *, int, int8_t)
Definition: colver_3d.c:42
int MMG5_movtet(MMG5_pMesh mesh, MMG5_pSol met, MMG3D_pPROctree PROctree, double clickSurf, double clickVol, int moveVol, int improveSurf, int improveVolSurf, int improveVol, int maxit, int testmark)
Definition: mmg3d1.c:673
int MMG5_swpbdy(MMG5_pMesh, MMG5_pSol, int *, int, int, MMG3D_pPROctree, int8_t)
Definition: swap_3d.c:434
int MMG3D_split3cone_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:2005
int MMG3D_Free_all_var(va_list argptr)
Definition: variadic_3d.c:271
void MMG5_defaultValues(MMG5_pMesh)
int MMG5_hashPop(MMG5_Hash *hash, int a, int b)
Definition: hash_3d.c:791
void MMG5_freeXTets(MMG5_pMesh mesh)
Definition: zaldy_3d.c:327
int MMG5_swpgen(MMG5_pMesh, MMG5_pSol, int, int, int *, MMG3D_pPROctree, int8_t)
Definition: swapgen_3d.c:246
int MMG5_chkcol_nomint(MMG5_pMesh, MMG5_pSol, int, int8_t, int8_t, int *, int, int8_t)
Definition: colver_3d.c:728
int(* MMG3D_gradsizreq)(MMG5_pMesh, MMG5_pSol)
Definition: mmg3dexterns.c:11
int MMG3D_Free_structures_var(va_list argptr)
Definition: variadic_3d.c:459
int(* MMG5_movintpt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int)
Definition: mmg3dexterns.c:14
int MMG5_hPop(MMG5_HGeom *hash, int a, int b, int *ref, int16_t *tag)
Definition: hash_3d.c:886
static const uint8_t MMG5_permedge[12][6]
Definition: mmg3d.h:185
void MMG3D_computeOutqua(MMG5_pMesh, MMG5_pSol, int *, double *, double *, double *, int *, int *, int *, int his[5], int *, int)
Definition: quality_3d.c:661
void MMG5_Init_parameters(MMG5_pMesh mesh)
Definition: API_functions.c:51
int MMG5_bdyMCF(MMG5_pMesh)
void MMG3D_computeInqua(MMG5_pMesh, MMG5_pSol, int *, double *, double *, double *, int *, int *, int *, int his[5], int)
Definition: quality_3d.c:542
int MMG3D_gradsiz_ani(MMG5_pMesh, MMG5_pSol)
Definition: anisosiz_3d.c:1886
int MMG3D_splitItem(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int, int, double)
Definition: opttyp_3d.c:390
int MMG3D_memOption_memRepartition(MMG5_pMesh mesh)
Definition: zaldy_3d.c:172
int16_t MMG5_coquilTravel(MMG5_pMesh, int, int, int *, int *, int8_t *, int8_t *)
Definition: boulep_3d.c:1892
int MMG3D_delElt(MMG5_pMesh mesh, int iel)
Definition: zaldy_3d.c:121
int MMG3D_gradsiz_iso(MMG5_pMesh, MMG5_pSol)
Definition: isosiz_3d.c:1059
int MMG5_chkmani(MMG5_pMesh mesh)
Definition: mmg3d2.c:1508
int MMG3D_opttyp(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int)
Definition: opttyp_3d.c:467
int MMG5_setdhd(MMG5_pMesh mesh)
Definition: analys_3d.c:288
int MMG5_norver(MMG5_pMesh mesh)
Definition: analys_3d.c:499
int MMG3D_split3op_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:2720
int MMG5_interp4bar33_ani(MMG5_pMesh, MMG5_pSol, int, int, double *)
int MMG5_split3op(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:2852
int MMG3D_displayQualHisto_internal(int, double, double, double, int, int, int, int his[5], int, int, int)
Definition: quality_3d.c:477
double MMG5_volint(MMG5_pMesh)
double MMG3D_vfrac(MMG5_pMesh, MMG5_pSol, int, int)
Definition: mmg3d2.c:105
static const int8_t MMG5_iarfinv[4][6]
num of the j^th edge in the i^th face
Definition: mmg3d.h:154
int MMG5_chkswpgen(MMG5_pMesh, MMG5_pSol, int, int, int *, int *, double, int8_t)
Definition: swapgen_3d.c:55
int MMG3D_split3_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:1722
int MMG3D_mmg3d2(MMG5_pMesh, MMG5_pSol, MMG5_pSol)
Definition: mmg3d2.c:2188
void MMG3D_mark_pointsOnReqEdge_fromTetra(MMG5_pMesh mesh)
Definition: isosiz_3d.c:1024
int MMG5_movbdyrefpt_ani(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: anisomovpt_3d.c:678
double MMG5_timestepMCF(MMG5_pMesh, double)
static const uint8_t MMG5_arpt[4][3]
arpt[i]: edges passing through vertex i
Definition: mmg3d.h:162
int MMG5_chkswpbdy(MMG5_pMesh, MMG5_pSol, int *, int, int, int, int8_t)
Definition: swap_3d.c:55
void MMG3D_Free_arrays(MMG5_pMesh *, MMG5_pSol *, MMG5_pSol *, MMG5_pSol *, MMG5_pSol *)
Definition: variadic_3d.c:369
int MMG3D_optlap(MMG5_pMesh, MMG5_pSol)
Definition: optlap_3d.c:44
double(* MMG5_lenedgspl)(MMG5_pMesh, MMG5_pSol, int, MMG5_pTetra)
Definition: mmg3dexterns.c:6
int MMG5_chkfemtopo(MMG5_pMesh mesh)
Definition: chkmsh_3d.c:411
double(* MMG5_caltri)(MMG5_pMesh mesh, MMG5_pSol met, MMG5_pTria ptt)
Definition: mmg3dexterns.c:8
int(* MMG3D_gradsiz)(MMG5_pMesh, MMG5_pSol)
Definition: mmg3dexterns.c:10
int MMG5_chkBdryTria(MMG5_pMesh mesh)
Definition: hash_3d.c:1452
int MMG3D_addPROctreeRec(MMG5_pMesh, MMG3D_PROctree_s *, double *, const int, int)
Definition: PRoctree_3d.c:648
int MMG5_BezierReg(MMG5_pMesh mesh, int ip0, int ip1, double s, double v[3], double *o, double *no)
Definition: tools_3d.c:702
int MMG5_delone(MMG5_pMesh mesh, MMG5_pSol sol, int ip, int *list, int ilist)
Definition: delaunay_3d.c:140
int MMG5_BezierRidge(MMG5_pMesh mesh, int ip0, int ip1, double s, double *o, double *no1, double *no2, double *to)
Definition: tools_3d.c:160
int MMG3D_packMesh(MMG5_pMesh, MMG5_pSol, MMG5_pSol)
Definition: libmmg3d.c:900
int MMG5_boulernm(MMG5_pMesh mesh, MMG5_Hash *hash, int start, int ip, int *ng, int *nr)
Definition: boulep_3d.c:451
int MMG3D_defsiz_ani(MMG5_pMesh, MMG5_pSol)
Definition: anisosiz_3d.c:1254
int MMG5_bdryIso(MMG5_pMesh)
int MMG5_bouletrid(MMG5_pMesh, int, int, int, int *, int *, int *, int *, int *, int *)
Definition: boulep_3d.c:935
int MMG3D_split6_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:4765
void MMG5_tet2tri(MMG5_pMesh mesh, int k, int8_t ie, MMG5_Tria *ptt)
Definition: mmg3d1.c:52
int MMG5_trydisp(MMG5_pMesh, double *, short)
int MMG5_coquilface(MMG5_pMesh mesh, int start, int8_t iface, int, int *, int *, int *, int)
Definition: boulep_3d.c:1783
int MMG5_startedgsurfball(MMG5_pMesh mesh, int nump, int numq, int *list, int ilist)
Definition: tools_3d.c:107
int MMG3D_initPROctree(MMG5_pMesh, MMG3D_pPROctree *q, int nv)
Definition: PRoctree_3d.c:63
int MMG5_countelt(MMG5_pMesh mesh, MMG5_pSol sol, double *weightelt, long *npcible)
Definition: quality_3d.c:786
int MMG3D_tetraQual(MMG5_pMesh mesh, MMG5_pSol met, int8_t metRidTyp)
Definition: quality_3d.c:48
int MMG5_boulevolp(MMG5_pMesh mesh, int start, int ip, int *list)
Definition: boulep_3d.c:53
int MMG3D_memOption_memSet(MMG5_pMesh mesh)
Definition: zaldy_3d.c:157
int MMG5_BezierEdge(MMG5_pMesh mesh, int ip0, int ip1, double b0[3], double b1[3], int8_t isrid, double v[3])
Definition: bezier_3d.c:152
int MMG5_chkptonbdy(MMG5_pMesh, int)
Definition: chkmsh_3d.c:301
int MMG3D_PROctreein_iso(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int, double)
Definition: PRoctree_3d.c:1140
int(* MMG5_movbdyridpt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: mmg3dexterns.c:18
int MMG3D_displayQualHisto(int, double, double, double, int, int, int, int his[5], int, int, int)
Definition: quality_3d.c:443
int MMG3D_dichoto(MMG5_pMesh mesh, MMG5_pSol met, int k, int *vx)
Definition: mmg3d1.c:88
double MMG5_surftri(MMG5_pMesh, int, int)
int MMG3D_update_eltsVertices(MMG5_pMesh mesh)
Definition: libmmg3d.c:803
int MMG5_cenrad_iso(MMG5_pMesh mesh, double *ct, double *c, double *rad)
Definition: cenrad_3d.c:45
int MMG3D_pack_points(MMG5_pMesh mesh)
Definition: libmmg3d.c:849
int MMG5_settag(MMG5_pMesh, int, int, int16_t, int)
Definition: boulep_3d.c:1148
int MMG5_split4op(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:4116
int MMG3D_normalAdjaTri(MMG5_pMesh, int, int8_t, int, double n[3])
Definition: split_3d.c:463
int MMG5_split5(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:4527
int8_t MMG5_chkedg(MMG5_pMesh mesh, MMG5_Tria *pt, int8_t ori, double, double, int)
Definition: mmg3d1.c:309
static const int8_t MMG5_iarf[4][3]
iarf[i]: edges of face opposite to vertex i
Definition: mmg3d.h:152
double MMG5_caltet33_ani(MMG5_pMesh mesh, MMG5_pSol met, MMG5_pTetra pt)
Definition: quality_3d.c:95
int * MMG3D_sizeArbre(MMG3D_PROctree *q, int dim)
int MMG5_split4bar(MMG5_pMesh mesh, MMG5_pSol met, int k, int8_t)
Definition: split_3d.c:3302
int MMG3D_Free_names_var(va_list argptr)
Definition: variadic_3d.c:557
int MMG3D_openMesh(int imprim, const char *filename, FILE **inm, int *bin, char *, char *)
Definition: inout_3d.c:52
static const uint8_t MMG5_iare[6][2]
vertices of extremities of the edges of the tetra
Definition: mmg3d.h:156
void MMG5_nsort(int, double *, int8_t *)
Definition: tools_3d.c:51
int MMG5_movbdyridpt_iso(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: movpt_3d.c:1612
int MMG5_denoisbdy(MMG5_pMesh)
int MMG5_movbdyregpt_iso(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int, int)
Definition: movpt_3d.c:366
int MMG3D_isCellIncluded(double *cellCenter, double l, double *zoneCenter, double l0)
Definition: PRoctree_3d.c:263
int MMG3D_newPt(MMG5_pMesh mesh, double c[3], int16_t tag, int src)
Definition: zaldy_3d.c:39
int MMG5_movbdynompt_iso(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: movpt_3d.c:1135
int MMG3D_bdryBuild(MMG5_pMesh)
Definition: libmmg3d.c:97
int MMG3D_outqua(MMG5_pMesh mesh, MMG5_pSol met)
Definition: quality_3d.c:744
double MMG5_BezierGeod(double c1[3], double c2[3], double t1[3], double t2[3])
Definition: bezier_3d.c:111
int MMG5_movintpt_ani(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int)
Definition: anisomovpt_3d.c:55
int MMG3D_intmet33_ani(MMG5_pMesh, MMG5_pSol, int, int8_t, int, double)
Definition: intmet_3d.c:101
int MMG5_colver(MMG5_pMesh, MMG5_pSol, int *, int, int8_t, int8_t)
Definition: colver_3d.c:839
void MMG5_printTetra(MMG5_pMesh mesh, char *fileName)
Definition: tools_3d.c:900
int MMG3D_newElt(MMG5_pMesh mesh)
Definition: zaldy_3d.c:98
int(* MMG5_movbdyregpt)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int, int)
Definition: mmg3dexterns.c:15
static const uint8_t MMG5_ifar[6][2]
ifar[i][]: faces sharing the ith edge of the tetra
Definition: mmg3d.h:158
int MMG5_setadj(MMG5_pMesh mesh)
Definition: analys_3d.c:96
int MMG3D_localParamNm(MMG5_pMesh, int, int, int, double *, double *, double *)
Definition: tools_3d.c:1070
static const uint8_t MMG5_iprv3[7]
previous vertex of tetra: {3,0,1,2,3,0,1}
Definition: mmg3d.h:146
int MMG5_bdryUpdate(MMG5_pMesh)
Definition: hash_3d.c:2042
int MMG3D_split4op_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:3995
int MMG5_split2sf(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:1225
static const uint8_t MMG5_iarf_pr[5][5]
iarf[i]: edges of face i for a prism
Definition: mmg3d.h:167
int MMG5_boulenmInt(MMG5_pMesh mesh, int start, int ip, double t[3])
Definition: boulep_3d.c:316
int MMG5_splitedg(MMG5_pMesh mesh, MMG5_pSol met, int iel, int iar, double crit)
Definition: split_3d.c:5264
int MMG5_movbdyrefpt_iso(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: movpt_3d.c:782
int MMG3D_zaldy(MMG5_pMesh mesh)
Definition: zaldy_3d.c:314
void MMG3D_mergeBranchesRec(MMG3D_PROctree_s *, MMG3D_PROctree_s *, int, int, int *)
Definition: PRoctree_3d.c:833
int MMG5_movintpt_iso(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int)
Definition: movpt_3d.c:55
int MMG3D_getListSquareRec(MMG3D_PROctree_s *, double *, double *, MMG3D_PROctree_s ***, double *, double *, double, int, int, int *)
Definition: PRoctree_3d.c:441
int MMG5_setNmTag(MMG5_pMesh mesh, MMG5_Hash *hash)
Definition: hash_3d.c:757
int MMG5_chkVertexConnectedDomains(MMG5_pMesh mesh)
Definition: analys_3d.c:355
int(* MMG3D_PROctreein)(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int, double)
Definition: mmg3dexterns.c:20
int MMG3D_findEdge(MMG5_pMesh, MMG5_pTetra, int, int, int, int, int8_t *, int8_t *)
Definition: boulep_3d.c:110
double(* MMG5_caltet)(MMG5_pMesh mesh, MMG5_pSol met, MMG5_pTetra pt)
Definition: mmg3dexterns.c:7
static const uint8_t MMG5_idir_pr[5][4]
idir[i]: vertices of face i for a prism
Definition: mmg3d.h:165
int MMG5_saveDisp(MMG5_pMesh, MMG5_pSol)
Definition: mmg3d3.c:609
int MMG3D_split1_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:50
int MMG5_mmg3dChkmsh(MMG5_pMesh, int, int)
Definition: chkmsh_3d.c:111
int MMG5_BezierNom(MMG5_pMesh mesh, int ip0, int ip1, double s, double *o, double *no, double *to)
Definition: tools_3d.c:535
int MMG5_mmg3d1_delone(MMG5_pMesh, MMG5_pSol, int *)
Definition: mmg3d1_delone.c:1298
int MMG3D_split5_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:4436
int64_t MMG3D_getPROctreeCoordinate(MMG3D_pPROctree q, double *ver, int dim)
Definition: PRoctree_3d.c:181
int MMG5_bdryPerm(MMG5_pMesh)
Definition: hash_3d.c:2132
int MMG3D_Init_mesh_var(va_list argptr)
Definition: variadic_3d.c:177
int MMG5_boulesurfvolp(MMG5_pMesh mesh, int start, int ip, int iface, int *listv, int *ilistv, int *lists, int *ilists, int isnm)
Definition: boulep_3d.c:594
int MMG3D_movv_iso(MMG5_pMesh, MMG5_pSol, int, int)
Definition: movpt_3d.c:2176
int MMG3D_coquilFaceFirstLoop(MMG5_pMesh mesh, int start, int na, int nb, int8_t iface, int8_t ia, int *list, int *ilist, int *it1, int *it2, int *piv, int *adj, int8_t *hasadja, int *nbdy, int silent)
Definition: boulep_3d.c:1632
double MMG5_orcal_poi(double a[3], double b[3], double c[3], double d[3])
static const uint8_t MMG5_inxt3[7]
next vertex of tetra: {1,2,3,0,1,2,3}
Definition: mmg3d.h:144
int MMG5_BezierRef(MMG5_pMesh mesh, int ip0, int ip1, double s, double *o, double *no, double *to)
Definition: tools_3d.c:362
double MMG5_estavglen(MMG5_pMesh)
Definition: mmg3d3.c:44
int MMG3D_split2sf_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:1115
int MMG3D_loadMesh_opened(MMG5_pMesh mesh, FILE *inm, int bin)
Definition: inout_3d.c:112
int MMG3D_pack_tetraAndAdja(MMG5_pMesh mesh)
Definition: libmmg3d.c:272
int MMG5_movbdyridpt_ani(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int *, int, int)
Definition: anisomovpt_3d.c:1375
int MMG5_singul(MMG5_pMesh mesh)
Definition: analys_3d.c:425
int MMG5_split4sf(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:3758
int MMG3D_delPROctree(MMG5_pMesh mesh, MMG3D_pPROctree q, const int no)
Definition: PRoctree_3d.c:975
int MMG3D_getListSquare(MMG5_pMesh, double *, MMG3D_PROctree *, double *, MMG3D_PROctree_s ***)
int MMG3D_nmgeom(MMG5_pMesh mesh)
Definition: analys_3d.c:671
int MMG3D_inqua(MMG5_pMesh mesh, MMG5_pSol met)
Definition: quality_3d.c:613
int MMG5_mmg3dBezierCP(MMG5_pMesh mesh, MMG5_Tria *pt, MMG5_pBezier pb, int8_t ori)
Definition: bezier_3d.c:327
int MMG5_stiffelt(MMG5_pMesh, int, double *, double *)
int MMG3D_setMeshSize_initData(MMG5_pMesh, int, int, int, int, int, int)
Definition: API_functions_3d.c:214
int(* MMG5_interp4bar)(MMG5_pMesh, MMG5_pSol, int, int, double *)
Definition: mmg3dexterns.c:13
int MMG5_cenrad_ani(MMG5_pMesh mesh, double *ct, double *m, double *c, double *rad)
Definition: cenrad_3d.c:142
int MMG5_hEdge(MMG5_pMesh mesh, MMG5_HGeom *hash, int a, int b, int ref, int16_t tag)
Definition: hash_3d.c:977
int MMG5_hashFace(MMG5_pMesh, MMG5_Hash *, int, int, int, int)
Definition: hash_3d.c:95
void MMG3D_freePROctree_s(MMG5_pMesh, MMG3D_PROctree_s *q, int nv)
Definition: PRoctree_3d.c:111
int MMG3D_movePROctree(MMG5_pMesh, MMG3D_pPROctree, int, double *, double *)
Definition: PRoctree_3d.c:223
int16_t MMG5_openCoquilTravel(MMG5_pMesh, int, int, int *, int *, int8_t *, int8_t *)
Definition: boulep_3d.c:1944
int MMG5_movbdynomintpt_iso(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int *, int, int)
Definition: movpt_3d.c:1482
int MMG3D_split2_sim(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6])
Definition: split_3d.c:1457
void MMG3D_computeLESqua(MMG5_pMesh, MMG5_pSol, int *, double *, double *, double *, int *, int *, int *, int his[5], int)
Definition: quality_3d.c:359
int MMG5_chkmanicoll(MMG5_pMesh, int, int, int, int, int, int, int, int8_t, int8_t)
Definition: mmg3d2.c:1653
int MMG3D_memOption(MMG5_pMesh mesh)
Definition: zaldy_3d.c:239
void MMG3D_delPt(MMG5_pMesh mesh, int ip)
Definition: zaldy_3d.c:79
int MMG3D_computePrilen(MMG5_pMesh, MMG5_pSol, double *, double *, double *, int *, int *, int *, int *, int *, int *, int8_t, double **, int[9])
Definition: quality_3d.c:206
int MMG5_mmg3d3(MMG5_pMesh, MMG5_pSol, MMG5_pSol, int **)
Definition: mmg3d3.c:651
int MMG3D_setMeshSize_alloc(MMG5_pMesh)
Definition: zaldy_3d.c:256
int MMG3D_prilen(MMG5_pMesh mesh, MMG5_pSol met, int8_t)
Definition: quality_3d.c:324
int MMG3D_movetetrapoints(MMG5_pMesh, MMG5_pSol, MMG3D_pPROctree, int)
Definition: optbdry_3d.c:45
int MMG5_split2(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:1519
static const uint8_t MMG5_idir[4][3]
idir[i]: vertices of face opposite to vertex i
Definition: mmg3d.h:148
int MMG3D_movv_ani(MMG5_pMesh, MMG5_pSol, int, int)
Definition: movpt_3d.c:1951
int MMG3D_pack_prismsAndQuads(MMG5_pMesh mesh)
Definition: libmmg3d.c:376
int MMG3D_pack_pointArray(MMG5_pMesh mesh)
Definition: libmmg3d.c:468
static void MMG5_warnOrientation(MMG5_pMesh mesh)
Definition: mmg3d.h:551
int(* MMG5_intmet)(MMG5_pMesh, MMG5_pSol, int, int8_t, int, double)
Definition: mmg3dexterns.c:12
int MMG5_coquil(MMG5_pMesh mesh, int start, int ia, int *list)
Definition: boulep_3d.c:1362
int MMG5_split3cone(MMG5_pMesh mesh, MMG5_pSol met, int k, int vx[6], int8_t)
Definition: split_3d.c:2186
int MMG5_isbr(MMG5_pMesh, int)
Definition: tools_3d.c:41
int MMG5_bdrySet(MMG5_pMesh)
Definition: hash_3d.c:1757
int MMG5_intmet_iso(MMG5_pMesh, MMG5_pSol, int, int8_t, int, double)
Definition: intmet_3d.c:131
void MMG3D_unset_reqBoundaries(MMG5_pMesh mesh)
Definition: libmmg3d.c:871
void MMG3D_printArbre(MMG3D_PROctree *q)
int MMG5_intvolmet(MMG5_pMesh, MMG5_pSol, int, int8_t, double, double *)
int MMG5_srcbdy(MMG5_pMesh mesh, int start, int ia)
Definition: boulep_3d.c:1500
int MMG5_chkcol_bdy(MMG5_pMesh, MMG5_pSol, int, int8_t, int8_t, int *, int, int *, int, int, int, int8_t, int)
Definition: colver_3d.c:363
int(* MMG5_chkmsh)(MMG5_pMesh, int, int)
Definition: mmgexterns.c:27
int(* MMG5_bezierCP)(MMG5_pMesh, MMG5_Tria *, MMG5_pBezier, int8_t)
Definition: mmgexterns.c:28
int(* MMG5_indPt)(MMG5_pMesh mesh, int kp)
Definition: mmgexterns.c:31
int(* MMG5_grad2metreq_ani)(MMG5_pMesh, MMG5_pSol, MMG5_pTria, int, int)
Definition: mmgexterns.c:33
int(* MMG5_indElt)(MMG5_pMesh mesh, int kel)
Definition: mmgexterns.c:30
int(* MMG5_grad2met_ani)(MMG5_pMesh, MMG5_pSol, MMG5_pTria, int, int)
Definition: mmgexterns.c:32
Definition: mmg3d.h:196
struct MMG3D_PROctree_s * branches
Definition: mmg3d.h:197
int depth
Definition: mmg3d.h:200
int * v
Definition: mmg3d.h:198
int nbVer
Definition: mmg3d.h:199
Definition: mmg3d.h:208
int nv
Definition: mmg3d.h:209
MMG3D_PROctree_s * q0
Definition: mmg3d.h:211
int nc
Definition: mmg3d.h:210
Definition: mmgcommon.h:557
Hash table to store geometric edges.
Definition: libmmgtypes.h:521
Identic as MMG5_HGeom but use MMG5_hedge to store edges instead of MMG5_hgeom (memory economy).
Definition: libmmgtypes.h:542
MMG mesh structure.
Definition: libmmgtypes.h:552
int xt
Definition: libmmgtypes.h:568
int ne
Definition: libmmgtypes.h:559
Definition: libmmgtypes.h:610
Definition: libmmgtypes.h:370
Definition: libmmgtypes.h:301
Structure to store the surface tetrahedra of a MMG mesh.
Definition: libmmgtypes.h:387