From cee05fc0218c35cb86c1158f867869f0af4364a6 Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Fri, 10 Aug 2018 12:56:08 +0100 Subject: [PATCH] Allocate work vectors with scratch option. --- krylov/psb_ccg.F90 | 2 +- krylov/psb_dcg.F90 | 2 +- krylov/psb_scg.F90 | 2 +- krylov/psb_zcg.F90 | 2 +- 4 files changed, 4 insertions(+), 4 deletions(-) diff --git a/krylov/psb_ccg.F90 b/krylov/psb_ccg.F90 index 5ac253e7..bd16e657 100644 --- a/krylov/psb_ccg.F90 +++ b/krylov/psb_ccg.F90 @@ -172,7 +172,7 @@ subroutine psb_ccg_vect(a,prec,b,x,eps,desc_a,info,& naux=4*n_col allocate(aux(naux), stat=info) if (info == psb_success_) call psb_geall(wwrk,desc_a,info,n=5_psb_ipk_) - if (info == psb_success_) call psb_geasb(wwrk,desc_a,info,mold=x%v) + if (info == psb_success_) call psb_geasb(wwrk,desc_a,info,mold=x%v,scratch=.true.) if (info /= psb_success_) then info=psb_err_from_subroutine_non_ call psb_errpush(info,name) diff --git a/krylov/psb_dcg.F90 b/krylov/psb_dcg.F90 index e8d08f03..74d41d69 100644 --- a/krylov/psb_dcg.F90 +++ b/krylov/psb_dcg.F90 @@ -172,7 +172,7 @@ subroutine psb_dcg_vect(a,prec,b,x,eps,desc_a,info,& naux=4*n_col allocate(aux(naux), stat=info) if (info == psb_success_) call psb_geall(wwrk,desc_a,info,n=5_psb_ipk_) - if (info == psb_success_) call psb_geasb(wwrk,desc_a,info,mold=x%v) + if (info == psb_success_) call psb_geasb(wwrk,desc_a,info,mold=x%v,scratch=.true.) if (info /= psb_success_) then info=psb_err_from_subroutine_non_ call psb_errpush(info,name) diff --git a/krylov/psb_scg.F90 b/krylov/psb_scg.F90 index fad4538c..92241791 100644 --- a/krylov/psb_scg.F90 +++ b/krylov/psb_scg.F90 @@ -172,7 +172,7 @@ subroutine psb_scg_vect(a,prec,b,x,eps,desc_a,info,& naux=4*n_col allocate(aux(naux), stat=info) if (info == psb_success_) call psb_geall(wwrk,desc_a,info,n=5_psb_ipk_) - if (info == psb_success_) call psb_geasb(wwrk,desc_a,info,mold=x%v) + if (info == psb_success_) call psb_geasb(wwrk,desc_a,info,mold=x%v,scratch=.true.) if (info /= psb_success_) then info=psb_err_from_subroutine_non_ call psb_errpush(info,name) diff --git a/krylov/psb_zcg.F90 b/krylov/psb_zcg.F90 index f1d63cc7..691a0f08 100644 --- a/krylov/psb_zcg.F90 +++ b/krylov/psb_zcg.F90 @@ -172,7 +172,7 @@ subroutine psb_zcg_vect(a,prec,b,x,eps,desc_a,info,& naux=4*n_col allocate(aux(naux), stat=info) if (info == psb_success_) call psb_geall(wwrk,desc_a,info,n=5_psb_ipk_) - if (info == psb_success_) call psb_geasb(wwrk,desc_a,info,mold=x%v) + if (info == psb_success_) call psb_geasb(wwrk,desc_a,info,mold=x%v,scratch=.true.) if (info /= psb_success_) then info=psb_err_from_subroutine_non_ call psb_errpush(info,name)