Define CG as recursive.

7 years ago · bf3212e793
parent 9d70696e96
commit bf3212e793
12 changed files with 16 additions and 16 deletions
--- a/krylov/psb_ccg.F90
+++ b/krylov/psb_ccg.F90
@ -95,7 +95,7 @@
 !                                         estimate of) residual. 
 ! 
 !
-subroutine psb_ccg_vect(a,prec,b,x,eps,desc_a,info,&
+recursive subroutine psb_ccg_vect(a,prec,b,x,eps,desc_a,info,&
     & itmax,iter,err,itrace,istop,cond)
  use psb_base_mod
  use psb_prec_mod
--- a/krylov/psb_cfcg.F90
+++ b/krylov/psb_cfcg.F90
@ -101,9 +101,9 @@
 !                                            stopped when  |r| <= eps * |b|
 !                                         where r is the (preconditioned, recursive
 !                                         estimate of) residual. 
-! 
 !
-subroutine psb_cfcg_vect(a,prec,b,x,eps,desc_a,info,&
+!
+recursive subroutine psb_cfcg_vect(a,prec,b,x,eps,desc_a,info,&
     & itmax,iter,err,itrace,istop,cond)
  use psb_base_mod
  use psb_prec_mod
--- a/krylov/psb_ckrylov.f90
+++ b/krylov/psb_ckrylov.f90
@ -129,7 +129,7 @@ Subroutine psb_ckrylov_vect(method,a,prec,b,x,eps,desc_a,info,&
      integer(psb_ipk_), Optional, Intent(out)       :: iter
      Real(psb_spk_), Optional, Intent(out) :: err
    end subroutine psb_ckryl_rest_vect
-    Subroutine psb_ckryl_cond_vect(a,prec,b,x,eps,desc_a,info,&
+    recursive Subroutine psb_ckryl_cond_vect(a,prec,b,x,eps,desc_a,info,&
         &itmax,iter,err, itrace,istop,cond)
      import :: psb_ipk_, psb_spk_, psb_desc_type, &
           & psb_cspmat_type, psb_cprec_type, psb_c_vect_type
--- a/krylov/psb_dcg.F90
+++ b/krylov/psb_dcg.F90
@ -95,7 +95,7 @@
 !                                         estimate of) residual. 
 ! 
 !
-subroutine psb_dcg_vect(a,prec,b,x,eps,desc_a,info,&
+recursive subroutine psb_dcg_vect(a,prec,b,x,eps,desc_a,info,&
     & itmax,iter,err,itrace,istop,cond)
  use psb_base_mod
  use psb_prec_mod
--- a/krylov/psb_dfcg.F90
+++ b/krylov/psb_dfcg.F90
@ -101,9 +101,9 @@
 !                                            stopped when  |r| <= eps * |b|
 !                                         where r is the (preconditioned, recursive
 !                                         estimate of) residual. 
-! 
 !
-subroutine psb_dfcg_vect(a,prec,b,x,eps,desc_a,info,&
+!
+recursive subroutine psb_dfcg_vect(a,prec,b,x,eps,desc_a,info,&
     & itmax,iter,err,itrace,istop,cond)
  use psb_base_mod
  use psb_prec_mod
--- a/krylov/psb_dkrylov.f90
+++ b/krylov/psb_dkrylov.f90
@ -129,7 +129,7 @@ Subroutine psb_dkrylov_vect(method,a,prec,b,x,eps,desc_a,info,&
      integer(psb_ipk_), Optional, Intent(out)       :: iter
      Real(psb_dpk_), Optional, Intent(out) :: err
    end subroutine psb_dkryl_rest_vect
-    Subroutine psb_dkryl_cond_vect(a,prec,b,x,eps,desc_a,info,&
+    recursive Subroutine psb_dkryl_cond_vect(a,prec,b,x,eps,desc_a,info,&
         &itmax,iter,err, itrace,istop,cond)
      import :: psb_ipk_, psb_dpk_, psb_desc_type, &
           & psb_dspmat_type, psb_dprec_type, psb_d_vect_type
--- a/krylov/psb_scg.F90
+++ b/krylov/psb_scg.F90
@ -95,7 +95,7 @@
 !                                         estimate of) residual. 
 ! 
 !
-subroutine psb_scg_vect(a,prec,b,x,eps,desc_a,info,&
+recursive subroutine psb_scg_vect(a,prec,b,x,eps,desc_a,info,&
     & itmax,iter,err,itrace,istop,cond)
  use psb_base_mod
  use psb_prec_mod
--- a/krylov/psb_sfcg.F90
+++ b/krylov/psb_sfcg.F90
@ -101,9 +101,9 @@
 !                                            stopped when  |r| <= eps * |b|
 !                                         where r is the (preconditioned, recursive
 !                                         estimate of) residual. 
-! 
 !
-subroutine psb_sfcg_vect(a,prec,b,x,eps,desc_a,info,&
+!
+recursive subroutine psb_sfcg_vect(a,prec,b,x,eps,desc_a,info,&
     & itmax,iter,err,itrace,istop,cond)
  use psb_base_mod
  use psb_prec_mod
--- a/krylov/psb_skrylov.f90
+++ b/krylov/psb_skrylov.f90
@ -129,7 +129,7 @@ Subroutine psb_skrylov_vect(method,a,prec,b,x,eps,desc_a,info,&
      integer(psb_ipk_), Optional, Intent(out)       :: iter
      Real(psb_spk_), Optional, Intent(out) :: err
    end subroutine psb_skryl_rest_vect
-    Subroutine psb_skryl_cond_vect(a,prec,b,x,eps,desc_a,info,&
+    recursive Subroutine psb_skryl_cond_vect(a,prec,b,x,eps,desc_a,info,&
         &itmax,iter,err, itrace,istop,cond)
      import :: psb_ipk_, psb_spk_, psb_desc_type, &
           & psb_sspmat_type, psb_sprec_type, psb_s_vect_type
--- a/krylov/psb_zcg.F90
+++ b/krylov/psb_zcg.F90
@ -95,7 +95,7 @@
 !                                         estimate of) residual. 
 ! 
 !
-subroutine psb_zcg_vect(a,prec,b,x,eps,desc_a,info,&
+recursive subroutine psb_zcg_vect(a,prec,b,x,eps,desc_a,info,&
     & itmax,iter,err,itrace,istop,cond)
  use psb_base_mod
  use psb_prec_mod
--- a/krylov/psb_zfcg.F90
+++ b/krylov/psb_zfcg.F90
@ -101,9 +101,9 @@
 !                                            stopped when  |r| <= eps * |b|
 !                                         where r is the (preconditioned, recursive
 !                                         estimate of) residual. 
-! 
 !
-subroutine psb_zfcg_vect(a,prec,b,x,eps,desc_a,info,&
+!
+recursive subroutine psb_zfcg_vect(a,prec,b,x,eps,desc_a,info,&
     & itmax,iter,err,itrace,istop,cond)
  use psb_base_mod
  use psb_prec_mod
--- a/krylov/psb_zkrylov.f90
+++ b/krylov/psb_zkrylov.f90
@ -129,7 +129,7 @@ Subroutine psb_zkrylov_vect(method,a,prec,b,x,eps,desc_a,info,&
      integer(psb_ipk_), Optional, Intent(out)       :: iter
      Real(psb_dpk_), Optional, Intent(out) :: err
    end subroutine psb_zkryl_rest_vect
-    Subroutine psb_zkryl_cond_vect(a,prec,b,x,eps,desc_a,info,&
+    recursive Subroutine psb_zkryl_cond_vect(a,prec,b,x,eps,desc_a,info,&
         &itmax,iter,err, itrace,istop,cond)
      import :: psb_ipk_, psb_dpk_, psb_desc_type, &
           & psb_zspmat_type, psb_zprec_type, psb_z_vect_type