solver_lin Function

private pure function solver_lin(A, b, method) result(x)

Arguments

Type IntentOptional Attributes Name
real(kind=rk), intent(in), dimension(:, :), contiguous :: A
real(kind=rk), intent(in), dimension(:), contiguous :: b
character(len=*), intent(in), optional :: method

Return Value real(kind=rk), dimension(max(1, size(A, 2)))


Calls

proc~~solver_lin~~CallsGraph proc~solver_lin forsolver::solver_lin proc~gels_rel forsolver::gels_rel proc~solver_lin->proc~gels_rel proc~gesv_rel forsolver::gesv_rel proc~solver_lin->proc~gesv_rel interface~gels external_interfaces_solver::gels proc~gels_rel->interface~gels interface~gesv external_interfaces_solver::gesv proc~gesv_rel->interface~gesv

Called by

proc~~solver_lin~~CalledByGraph proc~solver_lin forsolver::solver_lin interface~solve forsolver::solve interface~solve->proc~solver_lin proc~modified_newton_method_t1 forsolver::modified_newton_method_T1 proc~modified_newton_method_t1->interface~solve proc~modified_quasi_cs_newton_method_t1 forsolver::modified_quasi_cs_newton_method_T1 proc~modified_quasi_cs_newton_method_t1->interface~solve proc~modified_quasi_fd_newton_method_t1 forsolver::modified_quasi_fd_newton_method_T1 proc~modified_quasi_fd_newton_method_t1->interface~solve proc~newton_method_t1 forsolver::newton_method_T1 proc~newton_method_t1->interface~solve proc~quasi_cs_newton_method_t1 forsolver::quasi_cs_newton_method_T1 proc~quasi_cs_newton_method_t1->interface~solve proc~quasi_fd_newton_method_t1 forsolver::quasi_fd_newton_method_T1 proc~quasi_fd_newton_method_t1->interface~solve program~test_solver1 test_solver1 program~test_solver1->interface~solve program~test_solver15 test_solver15 program~test_solver15->interface~solve program~test_solver16 test_solver16 program~test_solver16->interface~solve program~test_solver17 test_solver17 program~test_solver17->interface~solve program~test_solver2 test_solver2 program~test_solver2->interface~solve proc~newton_complex_step_rel_t1 forsolver::nlsolver%newton_complex_step_rel_T1 proc~newton_complex_step_rel_t1->proc~modified_quasi_cs_newton_method_t1 proc~newton_complex_step_rel_t1->proc~quasi_cs_newton_method_t1 proc~newton_rel_t1 forsolver::nlsolver%newton_rel_T1 proc~newton_rel_t1->proc~modified_newton_method_t1 proc~newton_rel_t1->proc~modified_quasi_fd_newton_method_t1 proc~newton_rel_t1->proc~newton_method_t1 proc~newton_rel_t1->proc~quasi_fd_newton_method_t1 none~solve forsolver::nlsolver%solve none~solve->proc~newton_complex_step_rel_t1 none~solve->proc~newton_rel_t1 program~test_solver10 test_solver10 program~test_solver10->none~solve program~test_solver11 test_solver11 program~test_solver11->none~solve program~test_solver12 test_solver12 program~test_solver12->none~solve program~test_solver13 test_solver13 program~test_solver13->none~solve program~test_solver14 test_solver14 program~test_solver14->none~solve program~test_solver3 test_solver3 program~test_solver3->none~solve program~test_solver4 test_solver4 program~test_solver4->none~solve program~test_solver5 test_solver5 program~test_solver5->none~solve program~test_solver6 test_solver6 program~test_solver6->none~solve program~test_solver7 test_solver7 program~test_solver7->none~solve program~test_solver8 test_solver8 program~test_solver8->none~solve program~test_solver9 test_solver9 program~test_solver9->none~solve

Source Code

   pure function solver_lin(A, b, method) result(x)

      ! inputs
      real(rk),     dimension(:, :), contiguous, intent(in) :: A      ! input matrix A
      real(rk),     dimension(:),    contiguous, intent(in) :: b      ! right-hand side matrix b
      character(*), optional,                    intent(in) :: method

      ! outputs:
      real(rk), dimension(max(1, size(A, 2))) :: x      ! solution matrix x
      
      ! local variables
      integer :: info ! result info

      ! call solver
      if (present(method)) then
         select case (method)
          case ('gesv')
            call gesv_rel(A, b, x, info)
          case ('gels')
            call gels_rel(A, b, x, info)
         end select
      else
         if (size(A,1)==size(A,2)) then
            call gesv_rel(A, b, x, info)
         else
            call gels_rel(A, b, x, info)
         end if
      end if

   end function solver_lin