main.cpp

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/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.
** =============================================================================
*/
 
#include <cassert>
#include <cstdio>
#include <cstdlib>
#include <csignal>
#include <string>
#include <cctype>
#include <cmath>
#include <cfloat>
 
using namespace std;
 
// if the "include/mmg" dir is in your include path
//#include "libmmg.h"
// if your include path do not contain the "mmg" subdirectories
#include "mmg/libmmg.h"
 
int main(int argc,char *argv[]) {
  MMG5_pMesh      mmgMesh;
  MMG5_pSol       mmgSol;
  int             ier;
  string          filename,filename_os,filename_o2d,filename_o3d;
  std::size_t     found;
 
  fprintf(stdout,"  -- TEST MMGLIB \n");
 
  if ( argc != 4 ) {
    printf(" Usage: %s 2d_filein 3d_filein fileout\n",argv[0]);
    return(1);
  }
 
  filename = argv[1];
 
  /* args of InitMesh:
   * MMG5_ARG_start: we start to give the args of a variadic func
   * MMG5_ARG_ppMesh: next arg will be a pointer over a MMG5_pMesh
   * &mmgMesh: pointer toward your MMG5_pMesh (that store your mesh)
   * MMG5_ARG_ppMet: next arg will be a pointer over a MMG5_pSol storing a metric
   * &mmgSol: pointer toward your MMG5_pSol (that store your metric) */
  mmgMesh = NULL;
  mmgSol  = NULL;
  MMG2D_Init_mesh(MMG5_ARG_start,
                  MMG5_ARG_ppMesh,&mmgMesh, MMG5_ARG_ppMet,&mmgSol,
                  MMG5_ARG_end);
 
  if ( MMG2D_loadMesh(mmgMesh,filename.c_str()) != 1 )  exit(EXIT_FAILURE);
  if ( MMG2D_loadSol(mmgMesh,mmgSol,filename.c_str()) != 1 )  exit(EXIT_FAILURE);
 
  ier = MMG2D_mmg2dlib(mmgMesh,mmgSol);
 
  if ( ier == MMG5_STRONGFAILURE ) {
    fprintf(stdout,"BAD ENDING OF MMG2DLIB: UNABLE TO SAVE MESH\n");
    return(ier);
  } else if ( ier == MMG5_LOWFAILURE )
    fprintf(stdout,"BAD ENDING OF MMG2DLIB\n");
 
  filename_o2d = argv[3];
 
  found = filename_o2d.find(".mesh");
  if (found==std::string::npos) {
    found = filename_o2d.find(".msh");
  }
 
  if (found!=std::string::npos) {
    filename_o2d.replace(found,1,"\0");
  }
  filename_o2d += ".2d";
  /*save result*/
  if ( MMG2D_saveMesh(mmgMesh,filename_o2d.c_str()) != 1 )  exit(EXIT_FAILURE);
  /*save metric*/
  if ( MMG2D_saveSol(mmgMesh,mmgSol,filename_o2d.c_str()) != 1 )  exit(EXIT_FAILURE);
 
  MMG2D_Free_all(MMG5_ARG_start,
                 MMG5_ARG_ppMesh,&mmgMesh, MMG5_ARG_ppMet,&mmgSol,
                 MMG5_ARG_end);
 
  /* Name and path of the mesh file */
  filename    = argv[2];
  filename_os = argv[3];
 
  found = filename_os.find(".mesh");
  if (found==std::string::npos) {
    found = filename_os.find(".msh");
  }
 
  if (found!=std::string::npos) {
    filename_os.replace(found,1,"\0");
  }
  filename_os += ".s";
 
  /* args of InitMesh
   * MMG5_ARG_start: we start to give the args of a variadic func
   * MMG5_ARG_ppMesh: next arg will be a pointer over a MMG5_pMesh
   * &mmgMesh: pointer toward your MMG5_pMesh (that store your mesh)
   * MMG5_ARG_ppMet: next arg will be a pointer over a MMG5_pSol storing a metric
   * &mmgSol: pointer toward your MMG5_pSol (that store your metric) */
 
  mmgMesh = NULL;
  mmgSol  = NULL;
  MMGS_Init_mesh(MMG5_ARG_start,
                 MMG5_ARG_ppMesh,&mmgMesh, MMG5_ARG_ppMet,&mmgSol,
                 MMG5_ARG_end);
 
  if ( MMGS_loadMesh(mmgMesh,filename.c_str()) != 1 )  exit(EXIT_FAILURE);
  if ( MMGS_loadSol(mmgMesh,mmgSol,filename.c_str()) != 1 )  exit(EXIT_FAILURE);
 
 
  ier = MMGS_mmgslib(mmgMesh,mmgSol);
 
  if ( ier == MMG5_STRONGFAILURE ) {
    fprintf(stdout,"BAD ENDING OF MMGSLIB: UNABLE TO SAVE MESH\n");
    return(ier);
  } else if ( ier == MMG5_LOWFAILURE )
    fprintf(stdout,"BAD ENDING OF MMGSLIB\n");
 
  if ( MMGS_saveMesh(mmgMesh,filename_os.c_str()) != 1 )  exit(EXIT_FAILURE);
 
  if ( MMGS_saveSol(mmgMesh,mmgSol,filename_os.c_str()) != 1 )  exit(EXIT_FAILURE);
 
  MMGS_Free_all(MMG5_ARG_start,
                MMG5_ARG_ppMesh,&mmgMesh, MMG5_ARG_ppMet,&mmgSol,
                MMG5_ARG_end);
 
  filename_o3d = argv[2];
 
  found = filename_o3d.find(".mesh");
  if (found==std::string::npos) {
    found = filename_o3d.find(".msh");
  }
 
  if (found!=std::string::npos) {
    filename_o3d.replace(found,1,"\0");
  }
  filename_o3d += ".3d";
 
  /* args of InitMesh:
   * MMG5_ARG_start: we start to give the args of a variadic func
   * MMG5_ARG_ppMesh: next arg will be a pointer over a MMG5_pMesh
   * &mmgMesh: pointer toward your MMG5_pMesh (that store your mesh)
   * MMG5_ARG_ppMet: next arg will be a pointer over a MMG5_pSol storing a metric
   * &mmgSol: pointer toward your MMG5_pSol (that store your metric) */
  mmgMesh = NULL;
  mmgSol  = NULL;
  MMG3D_Init_mesh(MMG5_ARG_start,
                  MMG5_ARG_ppMesh,&mmgMesh, MMG5_ARG_ppMet,&mmgSol,
                  MMG5_ARG_end);
 
  if ( MMG3D_loadMesh(mmgMesh,filename.c_str()) != 1 )  exit(EXIT_FAILURE);
 
  if ( MMG3D_loadSol(mmgMesh,mmgSol,filename.c_str()) != 1 )  exit(EXIT_FAILURE);
 
  ier = MMG3D_mmg3dlib(mmgMesh,mmgSol);
 
  if ( ier == MMG5_STRONGFAILURE ) {
    fprintf(stdout,"BAD ENDING OF MMG3DLIB: UNABLE TO SAVE MESH\n");
    return(ier);
  } else if ( ier == MMG5_LOWFAILURE )
    fprintf(stdout,"BAD ENDING OF MMG3DLIB\n");
 
  if ( MMG3D_saveMesh(mmgMesh,filename_o3d.c_str()) != 1 )  exit(EXIT_FAILURE);
 
  if ( MMG3D_saveSol(mmgMesh,mmgSol,filename_o3d.c_str()) != 1 )  exit(EXIT_FAILURE);
 
  MMG3D_Free_all(MMG5_ARG_start,
                 MMG5_ARG_ppMesh,&mmgMesh,MMG5_ARG_ppMet,&mmgSol,
                 MMG5_ARG_end);
 
  return(0);
}