Telluride Solver Routine

call SOLVER (Solution, RHS, Y_EQ_AX_PROJECTION)

[snip]

SUBROUTINE SOLVER (Solution, RHS, Y_EQ_AX)

  [snip]

  if (TRIM(JT%solver) == 'gmres') then

     ! GMRES solver.
     call JT_GMRES (JT%status, RHS, JT%Control, &
                    Solution, ...,              &
                    User_y_eq_Ax=Y_EQ_AX,       &
                    User_Precond=PRECOND)

  else if (TRIM(JT%solver) == 'cg') then

     ! CG solver.
     call JT_CG (JT%status, RHS, JT%Control, &
                 Solution, ...,              &
                 User_y_eq_Ax=Y_EQ_AX,       &
                 User_Precond=PRECOND)
  end if

  [snip]