choral.F90

Go to the documentation of this file.

module choral

  use basic_tools
  use choral_env
  use real_type

  use io
  use funclib
  use abstract_interfaces
  use csr_mod
  use graph_mod
  use algebra_lin

  use mumps_mod
  use fe_mod
  use quad_mod
  use mesh_mod
  use eigen_mod
  use precond_mod
  use ionicmodel_mod
  use stimulation_mod

  use fespace_mod
  use mesh_interfaces
  use quadmesh_mod
  use mesh_tools
  use mesh_mmg

  use ode_def
  use ode_problem_mod
  use ode_solution_mod
  use ode_output_mod
  use krylov_mod
  use integral
  use diffusion
  use fespacexk_mod

  use ode_opsplt_mod
  use newton_mod
  use cardiomodel_mod
  use elasticity

  use ode_solver_mod

  implicit none
  private


  !! !!!!!!!!!!!!!!!!!!!!!!
  !!
  !! DERIVED TYPES
  !!
  public :: mesh, csr, fespace, krylov, ionicmodel, cardiomodel
  public :: quadmesh, precond, graph, newton, fespacexk
  public :: ode_problem, ode_solution, ode_output, ode_solver
  public :: mumps, eigen, stimulation

  !! !!!!!!!!!!!!!!!!!!!!!!
  !!
  !! GENERIC ROUTINES
  !!
  public :: re, equal
  public :: clear, set, assemble, valid
  public :: print, write, read, solve
  public :: allocmem, freemem
  public :: matvecprod, add, transpose
  public :: l2_product, l2_dist, l2_dist_grad


  !! !!!!!!!!!!!!!!!!!!!!!!
  !!
  !! ROUTINES
  !!
  public :: crossprod_3d
  public :: choral_init
  public :: clock, quit, warning, inttostring
  public :: one_r3, emetric, vector_field_e_x
  public :: f_c0, f_c1, f_c2, f_c3, f_c4, f_c5
  public :: graph_prod
  public :: scale, copy
  public :: add_submatrix, gmsh_addview
  public :: eval_scal_fe_at_points
  public :: define_interfaces, closest_node
  public :: diffusion_massmat, diffusion_stiffmat
  public :: diffusion_massmat_vect, diffusion_mixed_divmat
  public :: diffusion_matrix_pattern
  public :: diffusion_neumann_rhs, diffusion_dirichlet
  public :: l2_dist_vect, l2_dist_grad_proj, integ
  public :: interp_scal_func, interp_scal_fe
  public :: maxedgelength, minedgelength
  public :: elasticity_matrix_pattern    
  public :: elasticity_massmat,  elasticity_neumann_rhs         
  public :: elasticity_stiffmat, elasticity_dirichlet         
  public :: s_prefactor, name_ode_problem, name_ode_method
  public :: initialcond
  public :: scale_conductivities, monodomain_assemble
  public :: check_ode_method, name_ode_solver_type
  public :: name_ode_method_opsplt
  public :: order_ode_method, set_ode_solver_output
  public :: celerity_l2_dist, celerity_l2_norm
  public :: mesh_analyse
  public :: remesh_levelset
  public :: interp_vect_func, extract_component



  !! !!!!!!!!!!!!!!!!!!!!!!
  !!
  !! GLOBAL VARIABLES
  !!
  public :: choral_dir, home_dir, gmsh_dir, app_dir 
  public :: sp, dp, rp, real_tol, pi 
  public :: fe_name, quad_name

  !! !!!!!!!!!!!!!!!!!!!!!!
  !!
  !! ABSTRACT INTERFACES
  !!
  public :: r3tor3, r3tor, rtor, rntornxl


contains

  subroutine choral_init(verb)

    use choral_variables
    use fe_mod
    use cell_mod
    use quad_mod
    use mumps_mod
    !$ use OMP_LIB
    
    integer, intent(in), optional :: verb

    integer :: rang
    logical :: b

    if (present(verb)) then
       choral_verb = verb
    end if

    if (choral_verb>0) then
       write(*,*)"choral          : choral_init    = &
            &INITIALISATION "

       write(*,*)"  PREPROCESSOR OPTIONS"

#if RPC==4
       write(*,*)"    REAL NUMBER PRECISION        = SIMPLE"

#elif RPC==8
       write(*,*)"    REAL NUMBER PRECISION        = DOUBLE"

#elif RPC==12
       write(*,*)"    REAL NUMBER PRECISION        = TRIPLE"

#elif RPC==16
       write(*,*)"    REAL NUMBER PRECISION        = QUADRUPLE"

#else
       write(*,*)"    REAL NUMBER PRECISION        = DOUBLE"
       call warning("VARIABLE 'RPC' reset", 1)
#endif

#if DBG>0
       write(*,*)"    DEBUG MODE                   = ON"
#else
       write(*,*)"    DEBUG MODE                   = OFF" 
#endif

#ifdef WOMP
       write(*,*)"    OPENMP                       = ON"
       !$OMP PARALLEL 
       rang = omp_get_num_threads()
       !$OMP END PARALLEL
       write(*,*)"      number of threads          =", int(rang,1)
       if (choral_verb>2) then 
          !$OMP PARALLEL PRIVATE(rang)
          rang = omp_get_thread_num()
          write(*,*)"    Thread number       ",&
               & int(rang,1), "   = acitve"
          !$OMP END PARALLEL       
       end if
#else
       write(*,*)"    OPENMP                       = OFF"       
#endif

#ifdef WMMG
       write(*,*)"    MMG                          = ON"
#endif

#ifdef WARPACK
       write(*,*)"    ARPACK                       = ON"
#endif

#ifdef WMUMPS
       write(*,*)"    MUMPS                        = ON"
#endif

#ifdef NOWARN
       write(*,*)"    NO WARNING ALLOWED           = ON"
#else
       write(*,*)"    NO WARNING ALLOWED           = OFF" 
#endif

    end if


    if (choral_verb>1) then   
       write(*,*)"  GLOBAL VARIABLES"
       write(*,*)"    REAL PRECISION (RP), KIND RP =", int(rp,1)
       write(*,*)"    VERBOSITY,       CHORAL_VERB =", int(choral_verb,1)
    end if
    
    b = ( choral_verb > 1 )

    call cell_init(b)
    call fe_init(b)
    call quad_init(b)
    call mumps_init()

    if (choral_verb>0) write(*,*)

  end subroutine choral_init



end module choral