main.F90

Go to the documentation of this file.

!> @author
!> Cecile Dobrzynski, Charles Dapogny, Pascal Frey and Algiane Froehly
!> @brief
!>  Example for using mmg3dlib (basic use)
 
PROGRAM main
 
  IMPLICIT NONE
 
!> Include here the mmg3d library hader file
! if the header file is in the "include" directory
! #include "libmmg3df.h"
 
! if the header file is in "include/mmg/mmg3d"
#include "mmg/mmg3d/libmmg3df.h"
 
  mmg5_data_ptr_t    :: mmgmesh
  mmg5_data_ptr_t    :: mmgsol
  INTEGER            :: ier,argc
  CHARACTER(len=300) :: exec_name,filename,fileout
 
  WRITE(*,*) "  -- TEST MMG3DLIB"
 
  argc =  command_argument_count();
  CALL get_command_argument(0, exec_name)
 
  IF ( argc /=2 ) THEN
     print*," Usage: ",trim(exec_name)," input_file_name output_filename"
     CALL exit(1);
  ENDIF
 
  ! Name and path of the mesh file
  CALL get_command_argument(1, filename)
  CALL get_command_argument(2, fileout)
 
  !> ------------------------------ STEP   I --------------------------
  !! 1) Initialisation of mesh and sol structures
  !! 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: your MMG5_pMesh (that store your mesh)
  !! MMG5_ARG_ppMet: next arg will be a pointer over a MMG5_pSol storing a metric
  !! mmgSol: your MMG5_pSol (that store your metric) */
 
  mmgmesh = 0
  mmgsol  = 0
 
  CALL mmg3d_init_mesh(mmg5_arg_start, &
       mmg5_arg_ppmesh,mmgmesh,mmg5_arg_ppmet,mmgsol, &
       mmg5_arg_end)
 
  !> 2) Build mesh in MMG5 format
  !! Two solutions: just use the MMG3D_loadMesh function that will read a .mesh(b)
  !!    file formatted or manually set your mesh using the MMG3D_Set* functions
 
  !> with MMG3D_loadMesh function
  CALL mmg3d_loadmesh(mmgmesh,trim(adjustl(filename)),&
       len(trim(adjustl(filename))),ier)
  IF ( ier == 0 )  CALL exit(102)
 
  !> 3) Build sol in MMG5 format
  !! Two solutions: just use the MMG3D_loadSol function that will read a .sol(b)
  !!    file formatted or manually set your sol using the MMG3D_Set* functions
 
  !> With MMG3D_loadSol function
  CALL mmg3d_loadsol(mmgmesh,mmgsol,trim(adjustl(filename)),&
       len(trim(adjustl(filename))),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(104)
  ENDIF
 
  !> 4) (not mandatory): check if the number of given entities match with mesh size
  CALL mmg3d_chk_meshdata(mmgmesh,mmgsol,ier)
  IF ( ier /= 1 ) CALL exit(105)
 
  !> ------------------------------ STEP  II --------------------------
  !! remesh function
  CALL mmg3d_mmg3dlib(mmgmesh,mmgsol,ier)
 
  IF ( ier == mmg5_strongfailure ) THEN
     print*,"BAD ENDING OF MMG3DLIB: UNABLE TO SAVE MESH"
     stop 2
  ELSE IF ( ier == mmg5_lowfailure ) THEN
     print*,"BAD ENDING OF MMG3DLIB"
  ENDIF
 
  !> ------------------------------ STEP III --------------------------
  !! get results
  !! Two solutions: just use the MMG3D_saveMesh/MMG3D_saveSol functions
  !!    that will write .mesh(b)/.sol formatted files or manually get your mesh/sol
  !!    using the MMG3D_getMesh/MMG3D_getSol functions
 
  !> 1) Automatically save the mesh
  CALL mmg3d_savemesh(mmgmesh,trim(adjustl(fileout)),&
       len(trim(adjustl(fileout))),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(106)
  ENDIF
 
  !> 2) Automatically save the solution
  CALL mmg3d_savesol(mmgmesh,mmgsol,trim(adjustl(fileout)),&
       len(trim(adjustl(fileout))),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(107)
  ENDIF
 
  !> 3) Free the MMG3D5 structures
  CALL mmg3d_free_all(mmg5_arg_start, &
       mmg5_arg_ppmesh,mmgmesh,mmg5_arg_ppmet,mmgsol, &
       mmg5_arg_end)
 
END PROGRAM main