solve_syst Subroutine

public subroutine solve_syst(mat, step)

Note

General hat subroutine that handles the resolution steps:

  • ini solver initialization
  • ana solver analyzis when it’s proposed by the solver
  • fac solver factorization
  • sol solver solution
  • fre solver memory release when it’s proposed by the solver
  • end solver end

Arguments

Type IntentOptional Attributes Name
type(MAT_SOLV), intent(inout) :: mat

high level system type
character(len=*), intent(in) :: step

‘ini’=initialize, ‘ana’=analyze, ‘fac’=factorize, ‘sol’=solve, ‘fre’=free memory, ‘end’=close solver

Calls

proc~~solve_syst~~CallsGraph proc~solve_syst solve_syst proc~analyse_solver analyse_solver proc~solve_syst->proc~analyse_solver proc~close_solver close_solver proc~solve_syst->proc~close_solver proc~factorize_solver factorize_solver proc~solve_syst->proc~factorize_solver proc~freefact_solver freefact_solver proc~solve_syst->proc~freefact_solver proc~init_solver init_solver proc~solve_syst->proc~init_solver proc~solution_solver solution_solver proc~solve_syst->proc~solution_solver dmumps dmumps proc~analyse_solver->dmumps ma48_analyse ma48_analyse proc~analyse_solver->ma48_analyse proc~prep_superlu prep_superlu proc~analyse_solver->proc~prep_superlu proc~s_umfpack_di_symbolic s_umfpack_di_symbolic proc~analyse_solver->proc~s_umfpack_di_symbolic proc~close_solver->dmumps infog infog proc~close_solver->infog ma48_finalize ma48_finalize proc~close_solver->ma48_finalize mpi_finalize mpi_finalize proc~close_solver->mpi_finalize proc~close_superlu close_superlu proc~close_solver->proc~close_superlu proc~umfpack_di_free_numeric umfpack_di_free_numeric proc~close_solver->proc~umfpack_di_free_numeric proc~umfpack_di_free_symbolic umfpack_di_free_symbolic proc~close_solver->proc~umfpack_di_free_symbolic proc~umfpack_di_report_info umfpack_di_report_info proc~close_solver->proc~umfpack_di_report_info proc~factorize_solver->dmumps ma48_factorize ma48_factorize proc~factorize_solver->ma48_factorize proc~fact_superlu fact_superlu proc~factorize_solver->proc~fact_superlu proc~s_umfpack_di_numeric s_umfpack_di_numeric proc~factorize_solver->proc~s_umfpack_di_numeric proc~factorize_solver->proc~umfpack_di_free_numeric proc~free_superlu free_superlu proc~freefact_solver->proc~free_superlu proc~freefact_solver->proc~umfpack_di_free_numeric proc~init_solver->dmumps icntl icntl proc~init_solver->icntl proc~init_solver->infog ma48_initialize ma48_initialize proc~init_solver->ma48_initialize proc~init_solver->mpi_finalize mpi_init mpi_init proc~init_solver->mpi_init proc~init_superlu init_superlu proc~init_solver->proc~init_superlu proc~umfpack_di_defaults umfpack_di_defaults proc~init_solver->proc~umfpack_di_defaults proc~umfpack_di_report_control umfpack_di_report_control proc~init_solver->proc~umfpack_di_report_control proc~solution_solver->proc~solve_syst proc~solution_solver->dmumps proc~solution_solver->infog ma48_solve ma48_solve proc~solution_solver->ma48_solve proc~solution_solver->mpi_finalize proc~s_umfpack_di_solve s_umfpack_di_solve proc~solution_solver->proc~s_umfpack_di_solve proc~solv_superlu solv_superlu proc~solution_solver->proc~solv_superlu rhs rhs proc~solution_solver->rhs interface~destroy_compcol_matrix Destroy_CompCol_Matrix proc~close_superlu->interface~destroy_compcol_matrix interface~destroy_dense_matrix Destroy_Dense_Matrix proc~close_superlu->interface~destroy_dense_matrix interface~dgssvx dgssvx proc~fact_superlu->interface~dgssvx interface~statfree StatFree proc~fact_superlu->interface~statfree interface~statinit StatInit proc~fact_superlu->interface~statinit interface~statprint StatPrint proc~fact_superlu->interface~statprint interface~set_default_options set_default_options proc~init_superlu->interface~set_default_options interface~dcreate_compcol_matrix dCreate_CompCol_Matrix proc~prep_superlu->interface~dcreate_compcol_matrix interface~dcreate_dense_matrix dCreate_Dense_Matrix proc~prep_superlu->interface~dcreate_dense_matrix proc~umfpack_di_numeric umfpack_di_numeric proc~s_umfpack_di_numeric->proc~umfpack_di_numeric proc~umfpack_di_solve umfpack_di_solve proc~s_umfpack_di_solve->proc~umfpack_di_solve proc~umfpack_di_symbolic umfpack_di_symbolic proc~s_umfpack_di_symbolic->proc~umfpack_di_symbolic proc~solv_superlu->proc~prep_superlu proc~solv_superlu->interface~destroy_compcol_matrix interface~destroy_supernode_matrix Destroy_SuperNode_Matrix proc~solv_superlu->interface~destroy_supernode_matrix proc~solv_superlu->interface~dgssvx proc~solv_superlu->interface~statfree proc~solv_superlu->interface~statinit proc~solv_superlu->interface~statprint interface~c_umfpack_di_defaults c_umfpack_di_defaults proc~umfpack_di_defaults->interface~c_umfpack_di_defaults interface~c_umfpack_di_free_numeric c_umfpack_di_free_numeric proc~umfpack_di_free_numeric->interface~c_umfpack_di_free_numeric interface~c_umfpack_di_free_symbolic c_umfpack_di_free_symbolic proc~umfpack_di_free_symbolic->interface~c_umfpack_di_free_symbolic interface~c_umfpack_di_report_control c_umfpack_di_report_control proc~umfpack_di_report_control->interface~c_umfpack_di_report_control interface~c_umfpack_di_report_info c_umfpack_di_report_info proc~umfpack_di_report_info->interface~c_umfpack_di_report_info interface~c_umfpack_di_numeric c_umfpack_di_numeric proc~umfpack_di_numeric->interface~c_umfpack_di_numeric interface~c_umfpack_di_solve c_umfpack_di_solve proc~umfpack_di_solve->interface~c_umfpack_di_solve interface~c_umfpack_di_symbolic c_umfpack_di_symbolic proc~umfpack_di_symbolic->interface~c_umfpack_di_symbolic

Called by

proc~~solve_syst~~CalledByGraph proc~solve_syst solve_syst proc~solution_solver solution_solver proc~solve_syst->proc~solution_solver proc~solution_solver->proc~solve_syst program~test_solvers test_solvers program~test_solvers->proc~solve_syst

Source Code

   subroutine solve_syst(mat, step)
   implicit none
   type(MAT_SOLV), intent(inout) :: mat   !! *high level system type*
   character(len=*), intent(in)  :: step  !! *'ini'=initialize, 'ana'=analyze, 'fac'=factorize, 'sol'=solve, 'fre'=free memory, 'end'=close solver*

      if ( index(step, 'ini')/=0 ) then ; call init_solver(     mat) ; mat%ini = .true. ; return ; endif
      if ( index(step, 'ana')/=0 ) then ; call analyse_solver(  mat) ; mat%ana = .true. ; return ; endif
      if ( index(step, 'fac')/=0 ) then ; call factorize_solver(mat) ; mat%fac = .true. ; return ; endif
      if ( index(step, 'sol')/=0 ) then ; call solution_solver( mat) ; mat%sol = .true. ; return ; endif
      if ( index(step, 'fre')/=0 ) then ; call freefact_solver( mat) ; mat%fre = .true. ; return ; endif
      if ( index(step, 'end')/=0 ) then ; call close_solver( mat)    ; mat%end = .true. ; return ; endif
      stop 'Bad step chosen in SOLVE_SYST'

   return
   endsubroutine solve_syst