main.c

Go to the documentation of this file.

 
#include <assert.h>
#include <stdio.h>
#include <stdlib.h>
#include <signal.h>
#include <string.h>
#include <ctype.h>
#include <math.h>
#include <float.h>

// if the header file is in the "include" directory
// #include "libmmg3d.h"
// if the header file is in "include/mmg/mmg3d"
#include "mmg/mmg3d/libmmg3d.h"
 
int main(int argc,char *argv[]) {
  MMG5_pMesh      mmgMesh;
  MMG5_pSol       mmgMetric, mmgDisp;
  int             ier;
  char            *inname, *outname;
 
  fprintf(stdout, "  -- TEST MMG3DMOV \n");
 
  if ( argc != 3 ) {
    printf(" Usage: %s filein fileout \n", argv[0]);
    return(1);
  }
 
  /* Name and path of the mesh files */
  inname = (char *) calloc(strlen(argv[1]) + 1, sizeof(char));
  if ( inname == NULL ) {
    perror("  ## Memory problem: calloc");
    exit(EXIT_FAILURE);
  }
  strcpy(inname,argv[1]);
 
  outname = (char *) calloc(strlen(argv[2]) + 1, sizeof(char));
  if ( outname == NULL ) {
    perror("  ## Memory problem: calloc");
    exit(EXIT_FAILURE);
  }
  strcpy(outname,argv[2]);

  /* args of InitMesh:
   * MMG5_ARG_start: we start to give the args of a variadic function
   * MMG5_ARG_ppMesh: next arg will be a pointer to an MMG5_pMesh
   * &mmgMesh: pointer to your MMG5_pMesh (that stores your mesh)
   * MMG5_ARG_ppMet: next arg will be a pointer to MMG5_pSol storing a metric
   * &mmgMetric: pointer to an MMG5_pSol that stores your metric field
   * MMG5_ARG_ppDisp: next arg will be a pointer to MMG5_pSol storing a displacement
   * &mmgDisp: pointer to an MMG5_pSol that stores your displacement field
   * In this example the metric field is not used but it must be provided
   * to MMG3D_mmg3dmov() and therefore it must be initialized; it cannot
   * remain a NULL pointer. */
  mmgMesh = NULL;
  mmgMetric = NULL;
  mmgDisp = NULL;
  MMG3D_Init_mesh(MMG5_ARG_start,
                  MMG5_ARG_ppMesh, &mmgMesh,
                  MMG5_ARG_ppMet, &mmgMetric,
                  MMG5_ARG_ppDisp, &mmgDisp,
                  MMG5_ARG_end);


  if ( MMG3D_loadMesh(mmgMesh, inname) != 1 )  exit(EXIT_FAILURE);


  /* Ask for lagrangian motion (mode 1) */
  if ( MMG3D_Set_iparameter(mmgMesh, mmgDisp, MMG3D_IPARAM_lag, 1) != 1 )
    exit(EXIT_FAILURE);

  if ( MMG3D_loadSol(mmgMesh, mmgDisp, inname) != 1 )
    exit(EXIT_FAILURE);

  if ( MMG3D_Chk_meshData(mmgMesh, mmgDisp) != 1 ) exit(EXIT_FAILURE);

 

 
  /* debug mode ON (default value = OFF) */
  if ( MMG3D_Set_iparameter(mmgMesh, mmgDisp, MMG3D_IPARAM_debug, 1) != 1 )
    exit(EXIT_FAILURE);

  ier = MMG3D_mmg3dmov(mmgMesh, mmgMetric, mmgDisp);
 
  if ( ier == MMG5_STRONGFAILURE ) {
    fprintf(stdout,"BAD ENDING OF MMG3DMOV: UNABLE TO SAVE MESH\n");
    return(ier);
  } else if ( ier == MMG5_LOWFAILURE )
    fprintf(stdout,"BAD ENDING OF MMG3DMOV\n");
 
  /* (Not mandatory) Automatically save the mesh */
  if ( MMG3D_saveMesh(mmgMesh, outname) != 1 )
    exit(EXIT_FAILURE);
 
  /* 9) free the MMG3D5 structures */
  MMG3D_Free_all(MMG5_ARG_start,
                 MMG5_ARG_ppMesh, &mmgMesh,
                 MMG5_ARG_ppMet, &mmgMetric,
                 MMG5_ARG_ppDisp, &mmgDisp,
                 MMG5_ARG_end);
  free(inname);
  inname = NULL;
 
  free(outname);
  outname = NULL;
 
  return(ier);
}