libmmgs.h

Go to the documentation of this file.

/* =============================================================================
**  This file is part of the mmg software package for the tetrahedral
**  mesh modification.
**  Copyright (c) Bx INP/CNRS/Inria/UBordeaux/UPMC, 2004-
**
**  mmg is free software: you can redistribute it and/or modify it
**  under the terms of the GNU Lesser General Public License as published
**  by the Free Software Foundation, either version 3 of the License, or
**  (at your option) any later version.
**
**  mmg is distributed in the hope that it will be useful, but WITHOUT
**  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
**  FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
**  License for more details.
**
**  You should have received a copy of the GNU Lesser General Public
**  License and of the GNU General Public License along with mmg (in
**  files COPYING.LESSER and COPYING). If not, see
**  <http://www.gnu.org/licenses/>. Please read their terms carefully and
**  use this copy of the mmg distribution only if you accept them.
** =============================================================================
*/
 
#ifndef MMGSLIB_H
#define MMGSLIB_H
 
 
#ifdef __cplusplus
extern "C" {
#endif
 
#include "mmg/mmgs/libmmgtypes.h"
 
#define MMGS_LMAX      1024
 
enum MMGS_Param {
  MMGS_IPARAM_verbose,           
  MMGS_IPARAM_mem,               
  MMGS_IPARAM_debug,             
  MMGS_IPARAM_angle,             
  MMGS_IPARAM_iso,               
  MMGS_IPARAM_keepRef,           
  MMGS_IPARAM_optim,             
  MMGS_IPARAM_noinsert,          
  MMGS_IPARAM_noswap,            
  MMGS_IPARAM_nomove,            
  MMGS_IPARAM_nreg,              
  MMGS_IPARAM_numberOfLocalParam,
  MMGS_IPARAM_numsubdomain,      
  MMGS_IPARAM_renum,             
  MMGS_IPARAM_anisosize,        
  MMGS_IPARAM_nosizreq,          
  MMGS_DPARAM_angleDetection,    
  MMGS_DPARAM_hmin,              
  MMGS_DPARAM_hmax,              
  MMGS_DPARAM_hsiz,              
  MMGS_DPARAM_hausd,             
  MMGS_DPARAM_hgrad,             
  MMGS_DPARAM_hgradreq,          
  MMGS_DPARAM_ls,                
  MMGS_PARAM_size,               
};
 
/*----------------------------- functions header -----------------------------*/
/* Initialization functions */
/* init structures */
int MMGS_Init_mesh(const int starter,...);
 
void  MMGS_Init_fileNames(MMG5_pMesh mesh, MMG5_pSol sol);
void  MMGS_Init_parameters(MMG5_pMesh mesh);
 
/* init file names */
int  MMGS_Set_inputMeshName(MMG5_pMesh mesh, const char* meshin);
int  MMGS_Set_outputMeshName(MMG5_pMesh mesh, const char* meshout);
int  MMGS_Set_inputSolName(MMG5_pMesh mesh,MMG5_pSol sol, const char* solin);
int  MMGS_Set_outputSolName(MMG5_pMesh mesh,MMG5_pSol sol, const char* solout);
 
/* init structure sizes */
int  MMGS_Set_solSize(MMG5_pMesh mesh, MMG5_pSol sol, int typEntity, int np, int typSol);
  int MMGS_Set_solsAtVerticesSize(MMG5_pMesh mesh, MMG5_pSol *sol,int nsols,
                                  int nentities, int *typSol);
int  MMGS_Set_meshSize(MMG5_pMesh mesh, int np, int nt, int na);
 
/* init structure datas */
int  MMGS_Set_vertex(MMG5_pMesh mesh, double c0, double c1,
                     double c2, int ref,int pos);
 int  MMGS_Set_vertices(MMG5_pMesh mesh, double *vertices,int *refs);
int  MMGS_Set_triangle(MMG5_pMesh mesh, int v0, int v1,
                       int v2, int ref,int pos);
  int  MMGS_Set_triangles(MMG5_pMesh mesh, int *tria, int *refs);
int  MMGS_Set_edge(MMG5_pMesh mesh, int v0, int v1, int ref,int pos);
int  MMGS_Set_corner(MMG5_pMesh mesh, int k);
int  MMGS_Unset_corner(MMG5_pMesh mesh, int k);
int  MMGS_Set_requiredVertex(MMG5_pMesh mesh, int k);
  int  MMGS_Unset_requiredVertex(MMG5_pMesh mesh, int k);
 
int  MMGS_Set_requiredTriangle(MMG5_pMesh mesh, int k);
 
  int  MMGS_Unset_requiredTriangle(MMG5_pMesh mesh, int k);
 
int  MMGS_Set_ridge(MMG5_pMesh mesh, int k);
 
int  MMGS_Unset_ridge(MMG5_pMesh mesh, int k);
 
int  MMGS_Set_requiredEdge(MMG5_pMesh mesh, int k);
 
  int  MMGS_Unset_requiredEdge(MMG5_pMesh mesh, int k);
 
  int MMGS_Set_edges(MMG5_pMesh mesh, int *edges, int* refs);
  int MMGS_Get_edges(MMG5_pMesh mesh,int *edges,int* refs,
                     int *areRidges,int *areRequired);
 
int  MMGS_Set_normalAtVertex(MMG5_pMesh mesh, int k, double n0, double n1, double n2) ;
 
  double MMGS_Get_triangleQuality(MMG5_pMesh mesh,MMG5_pSol met, int k);
 
int  MMGS_Set_scalarSol(MMG5_pSol met, double s,int pos);
int  MMGS_Set_scalarSols(MMG5_pSol met, double *s);
int MMGS_Set_vectorSol(MMG5_pSol met, double vx,double vy, double vz, int pos);
int MMGS_Set_vectorSols(MMG5_pSol met, double *sols);
int MMGS_Set_tensorSol(MMG5_pSol met, double m11,double m12, double m13,
                       double m22,double m23, double m33, int pos);
int MMGS_Set_tensorSols(MMG5_pSol met, double *sols);
  int  MMGS_Set_ithSol_inSolsAtVertices(MMG5_pSol sol,int i, double* s,int pos);
  int  MMGS_Set_ithSols_inSolsAtVertices(MMG5_pSol sol,int i, double* s);
 
/* check init */
int MMGS_Chk_meshData(MMG5_pMesh mesh, MMG5_pSol met);
 
int  MMGS_Set_iparameter(MMG5_pMesh mesh,MMG5_pSol sol, int iparam, int val);
int  MMGS_Set_dparameter(MMG5_pMesh mesh,MMG5_pSol sol, int dparam, double val);
int  MMGS_Set_localParameter(MMG5_pMesh mesh, MMG5_pSol sol, int typ, int ref,
                             double hmin, double hmax, double hausd);
 
int  MMGS_Get_meshSize(MMG5_pMesh mesh, int* np, int* nt, int* na);
int  MMGS_Get_solSize(MMG5_pMesh mesh, MMG5_pSol sol, int* typEntity, int* np,
                      int* typSol);
  int  MMGS_Get_solsAtVerticesSize(MMG5_pMesh mesh, MMG5_pSol* sol,int *nsols,
                                    int* nentities,int* typSol);
 
int  MMGS_Get_vertex(MMG5_pMesh mesh, double* c0, double* c1, double* c2, int* ref,
                     int* isCorner, int* isRequired);
 int  MMGS_GetByIdx_vertex(MMG5_pMesh mesh, double* c0, double* c1, double* c2, int* ref,
                           int* isCorner, int* isRequired,int idx);
 
int  MMGS_Get_vertices(MMG5_pMesh mesh, double* vertices, int* refs,
                        int* areCorners, int* areRequired);
int  MMGS_Get_triangle(MMG5_pMesh mesh, int* v0, int* v1, int* v2, int* ref,
                       int* isRequired);
int  MMGS_Get_triangles(MMG5_pMesh mesh, int* tria, int* refs,
                         int* areRequired);
int  MMGS_Get_edge(MMG5_pMesh mesh, int* e0, int* e1, int* ref,
                   int* isRidge, int* isRequired);
 
int  MMGS_Get_normalAtVertex(MMG5_pMesh mesh, int k, double *n0, double *n1, double *n2) ;
 
int  MMGS_Get_scalarSol(MMG5_pSol met, double* s);
int  MMGS_Get_scalarSols(MMG5_pSol met, double* s);
int MMGS_Get_vectorSol(MMG5_pSol met, double* vx, double* vy, double* vz);
int MMGS_Get_vectorSols(MMG5_pSol met, double* sols);
int MMGS_Get_tensorSol(MMG5_pSol met, double *m11,double *m12, double *m13,
                       double *m22,double *m23, double *m33);
int MMGS_Get_tensorSols(MMG5_pSol met, double *sols);
  int  MMGS_Get_ithSol_inSolsAtVertices(MMG5_pSol sol,int i, double* s,int pos);
  int  MMGS_Get_ithSols_inSolsAtVertices(MMG5_pSol sol,int i, double* s);
int MMGS_Get_iparameter(MMG5_pMesh mesh, int iparam);
 
/* input/output functions */
int  MMGS_loadMesh(MMG5_pMesh mesh, const char* filename);
 
int MMGS_loadVtpMesh(MMG5_pMesh mesh,MMG5_pSol sol,const char *filename);
int MMGS_loadVtpMesh_and_allData(MMG5_pMesh mesh,MMG5_pSol *sol,const char *filename);
int MMGS_loadVtuMesh(MMG5_pMesh mesh,MMG5_pSol sol,const char *filename);
int MMGS_loadVtuMesh_and_allData(MMG5_pMesh mesh,MMG5_pSol *sol,const char *filename);
int MMGS_loadVtkMesh(MMG5_pMesh mesh,MMG5_pSol sol,const char *filename);
int MMGS_loadVtkMesh_and_allData(MMG5_pMesh mesh,MMG5_pSol *sol,const char *filename);
int MMGS_loadMshMesh(MMG5_pMesh mesh,MMG5_pSol sol,const char *filename);
int MMGS_loadMshMesh_and_allData(MMG5_pMesh mesh,MMG5_pSol *sol,const char *filename);
int  MMGS_saveMesh(MMG5_pMesh mesh, const char *filename);
 int MMGS_saveMshMesh(MMG5_pMesh mesh,MMG5_pSol sol,const char *filename);
int MMGS_saveMshMesh_and_allData(MMG5_pMesh mesh,MMG5_pSol *sol,const char *filename);
  int MMGS_saveVtkMesh(MMG5_pMesh mesh,MMG5_pSol sol,const char *filename);
  int MMGS_saveVtkMesh_and_allData(MMG5_pMesh mesh,MMG5_pSol *sol,const char *filename);
  int MMGS_saveVtuMesh(MMG5_pMesh mesh,MMG5_pSol sol,const char *filename);
  int MMGS_saveVtuMesh_and_allData(MMG5_pMesh mesh,MMG5_pSol *sol,const char *filename);
  int MMGS_saveVtpMesh(MMG5_pMesh mesh,MMG5_pSol sol,const char *filename);
  int MMGS_saveVtpMesh_and_allData(MMG5_pMesh mesh,MMG5_pSol *sol,const char *filename);
int  MMGS_loadSol(MMG5_pMesh mesh,MMG5_pSol met, const char* filename);
int  MMGS_loadAllSols(MMG5_pMesh mesh,MMG5_pSol *sol, const char* filename);
int  MMGS_saveSol(MMG5_pMesh mesh, MMG5_pSol met, const char *filename);
int MMGS_saveAllSols(MMG5_pMesh  mesh,MMG5_pSol *sol ,const char *filename);
 
int MMGS_Free_allSols(MMG5_pMesh mesh,MMG5_pSol *sol);
 
/* deallocations */
int MMGS_Free_all(const int starter,...);
 
int MMGS_Free_structures(const int starter,...);
 
int MMGS_Free_names(const int starter,...);
 
/* library */
int  MMGS_mmgslib(MMG5_pMesh mesh, MMG5_pSol met);
 
int  MMGS_mmgsls(MMG5_pMesh mesh,  MMG5_pSol sol,MMG5_pSol met);
 
void  MMGS_setfunc(MMG5_pMesh mesh,MMG5_pSol met);
 
  int MMGS_Get_numberOfNonBdyEdges(MMG5_pMesh mesh, int* nb_edges);
 
  int MMGS_Get_nonBdyEdge(MMG5_pMesh mesh, int* e0, int* e1, int* ref, int idx);
 
 
/* Tools for the library */
int MMGS_doSol(MMG5_pMesh mesh,MMG5_pSol met);
 
int MMGS_Set_constantSize(MMG5_pMesh mesh,MMG5_pSol met);
 
int MMGS_usage(char *prog);
int  MMGS_parsar(int argc,char *argv[],MMG5_pMesh mesh,MMG5_pSol met,MMG5_pSol sol);
int MMGS_defaultValues(MMG5_pMesh mesh);
int MMGS_stockOptions(MMG5_pMesh mesh, MMG5_Info *info);
void MMGS_destockOptions(MMG5_pMesh mesh, MMG5_Info *info);
 
int MMGS_Get_adjaTri(MMG5_pMesh mesh, int kel, int listri[3]);
 
int MMGS_Get_adjaVerticesFast(MMG5_pMesh mesh, int ip,int start, int lispoi[MMGS_LMAX]);
 
int MMGS_Compute_eigenv(double m[6],double lambda[3],double vp[3][3]);
 
  void MMGS_Free_solutions(MMG5_pMesh mesh,MMG5_pSol sol);
 
#ifdef __cplusplus
}
#endif
 
#endif