main.F90

Go to the documentation of this file.

 
 
PROGRAM main
 
  IMPLICIT NONE
 
! if the "include/mmg" dir is in your include path
!#include "libmmgf.h"
 
! if your include path do not contain the "mmg/mmg" subdirectories
#include "mmg/libmmgf.h"
 
  mmg5_data_ptr_t    :: mmgmesh
  mmg5_data_ptr_t    :: mmgsol
  INTEGER            :: ier,argc
  CHARACTER(len=350) :: tdfile,tdfileout, bdfile, bdfileout,exec_name,fileout
 
  WRITE(*,*) "  -- TEST MMGLIB"
 
  argc =  command_argument_count();
  CALL get_command_argument(0, exec_name)
 
  IF ( argc /=3 ) THEN
     print*,argc," Usage: ",trim(exec_name),&
          " 2d_file_name 3d_file_name output_file_name"
     CALL exit(1);
  ENDIF
 
  ! Name and path of the mesh file
  CALL get_command_argument(1, bdfile)
  CALL get_command_argument(2, tdfile)
  CALL get_command_argument(3, fileout)
 
 
 
  mmgmesh = 0
  mmgsol  = 0
 
  CALL mmg2d_init_mesh(mmg5_arg_start, &
       mmg5_arg_ppmesh,mmgmesh,mmg5_arg_ppmet,mmgsol, &
       mmg5_arg_end)
 
 
  CALL mmg2d_loadmesh(mmgmesh,trim(adjustl(bdfile)),&
       len(trim(adjustl(bdfile))),ier)
  IF ( ier /= 1 )  CALL exit(102)
 
  CALL mmg2d_loadsol(mmgmesh,mmgsol,trim(adjustl(bdfile)),len(trim(adjustl(bdfile))),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(104)
  ENDIF
 
 
  CALL mmg2d_mmg2dlib(mmgmesh,mmgsol,ier)
 
  IF ( ier == mmg5_strongfailure ) THEN
     print*,"BAD ENDING OF MMG2DLIB: UNABLE TO SAVE MESH"
     stop 2
  ELSE IF ( ier == mmg5_lowfailure ) THEN
     print*,"BAD ENDING OF MMG2DLIB"
  ENDIF
 
 
  CALL mmg2d_savemesh(mmgmesh,trim(adjustl(fileout)) // ".2d.mesh",&
       len(trim(adjustl(fileout)) // ".2d.mesh"),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(106)
  ENDIF
 
 
  !! b) save the solution in a file named bdfileout
  CALL mmg2d_savesol(mmgmesh,mmgsol,trim(adjustl(fileout)) // ".2d.mesh",&
       len(trim(adjustl(fileout)) // ".2d.mesh"),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(107)
  ENDIF
 
  CALL mmg2d_free_all(mmg5_arg_start, &
       mmg5_arg_ppmesh,mmgmesh,mmg5_arg_ppmet,mmgsol, &
       mmg5_arg_end)
 
 
  !! ------------------------------ 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 mmgs_init_mesh(mmg5_arg_start, &
       mmg5_arg_ppmesh,mmgmesh,mmg5_arg_ppmet,mmgsol, &
       mmg5_arg_end)
 
 
  CALL mmgs_loadmesh(mmgmesh,trim(adjustl(tdfile)),&
       len(trim(adjustl(tdfile))),ier);
  IF ( ier /= 1 )  CALL exit(102)
 
 
  CALL mmgs_loadsol(mmgmesh,mmgsol,trim(adjustl(tdfile)),&
       len(trim(adjustl(tdfile))),ier);
  IF ( ier /= 1 ) THEN
     CALL exit(104)
  ENDIF
 
  CALL mmgs_chk_meshdata(mmgmesh,mmgsol,ier)
  IF ( ier /=1 ) CALL exit(105)
 
  CALL mmgs_mmgslib(mmgmesh,mmgsol,ier)
 
  IF ( ier == mmg5_strongfailure ) THEN
     print*,"BAD ENDING OF MMGSLIB: UNABLE TO SAVE MESH"
     stop 2
  ELSE IF ( ier == mmg5_lowfailure ) THEN
     print*,"BAD ENDING OF MMGSLIB"
  ENDIF
 
 
  CALL mmgs_savemesh(mmgmesh,trim(adjustl(fileout)) // ".s.mesh",&
       len(trim(adjustl(fileout)) // ".s.mesh"),ier);
  IF ( ier /= 1 ) THEN
     CALL exit(108)
  ENDIF
 
  CALL mmgs_savesol(mmgmesh,mmgsol,trim(adjustl(fileout))//".s.sol",&
       len(trim(adjustl(fileout))//".s.sol"),ier);
  IF ( ier /= 1 ) THEN
     CALL exit(107)
  ENDIF
 
  CALL mmgs_free_all(mmg5_arg_start, &
       mmg5_arg_ppmesh,mmgmesh,mmg5_arg_ppmet,mmgsol, &
       mmg5_arg_end)
 
 
  !! ------------------------------ 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)
 
 
  CALL mmg3d_loadmesh(mmgmesh,trim(adjustl(tdfile)),&
       len(trim(adjustl(tdfile))),ier)
  IF ( ier /= 1 )  CALL exit(102)
 
 
  CALL mmg3d_loadsol(mmgmesh,mmgsol,trim(adjustl(tdfile)),len(trim(adjustl(tdfile))),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(104)
  ENDIF
 
  CALL mmg3d_chk_meshdata(mmgmesh,mmgsol,ier)
  IF ( ier /= 1 ) CALL exit(105)
 
  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
 
 
  CALL mmg3d_savemesh(mmgmesh,trim(adjustl(fileout)) // ".3d.mesh",&
       len(trim(adjustl(fileout)) // ".3d.mesh"),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(108)
  ENDIF
 
  CALL mmg3d_savesol(mmgmesh,mmgsol,trim(adjustl(fileout)) // ".3d.sol",&
       len(trim(adjustl(fileout)) // ".3d.sol"),ier)
  IF ( ier /= 1 ) THEN
     CALL exit(107)
  ENDIF
 
 
  CALL mmg3d_free_all(mmg5_arg_start, &
       mmg5_arg_ppmesh,mmgmesh,mmg5_arg_ppmet,mmgsol, &
       mmg5_arg_end)
 
END PROGRAM main