Initialize alpha in CG.

6 years ago · aaf623bfc2
parent 9d70696e96
commit aaf623bfc2
4 changed files with 4 additions and 4 deletions
--- a/krylov/psb_ccg.F90
+++ b/krylov/psb_ccg.F90
@ -203,7 +203,7 @@ subroutine psb_ccg_vect(a,prec,b,x,eps,desc_a,info,&
    istebz = 0
  end if
  itx=0
-
+  alpha = czero

  restart: do 
 ! =   
--- a/krylov/psb_dcg.F90
+++ b/krylov/psb_dcg.F90
@ -211,7 +211,7 @@ subroutine psb_dcg_vect(a,prec,b,x,eps,desc_a,info,&
    end if
  end if
  itx=0
-
+  alpha = dzero

  restart: do 
 ! =   
--- a/krylov/psb_scg.F90
+++ b/krylov/psb_scg.F90
@ -211,7 +211,7 @@ subroutine psb_scg_vect(a,prec,b,x,eps,desc_a,info,&
    end if
  end if
  itx=0
-
+  alpha = szero

  restart: do 
 ! =   
--- a/krylov/psb_zcg.F90
+++ b/krylov/psb_zcg.F90
@ -203,7 +203,7 @@ subroutine psb_zcg_vect(a,prec,b,x,eps,desc_a,info,&
    istebz = 0
  end if
  itx=0
-
+  alpha = zzero

  restart: do 
 ! =