fdurastante
/
psblas3
mirror of https://github.com/sfilippone/psblas3.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Fixed PRECAPLY name.

Browse Source
psblas3-type-indexed
Salvatore Filippone 18 years ago
parent c61e0af41e
commit 0bbc1835bc
8 changed files with 18 additions and 18 deletions
Show Stats Download Patch File Download Diff File

4
krylov/psb_dbicg.f90
Unescape Escape View File

@ -259,8 +259,8 @@ subroutine psb_dbicg(a,prec,b,x,eps,desc_a,info,&
itx = itx + 1 itx = itx + 1
if (debug) write(*,*) 'iteration: ',itx if (debug) write(*,*) 'iteration: ',itx
call psb_prc_aply(prec,r,z,desc_a,info,work=aux) call psb_precaply(prec,r,z,desc_a,info,work=aux)
call psb_prc_aply(prec,rt,zt,desc_a,info,trans='t',work=aux) call psb_precaply(prec,rt,zt,desc_a,info,trans='t',work=aux)
rho_old = rho rho_old = rho
rho = psb_gedot(rt,z,desc_a,info) rho = psb_gedot(rt,z,desc_a,info)

2
krylov/psb_dcg.f90
Unescape Escape View File

@ -206,7 +206,7 @@ Subroutine psb_dcg(a,prec,b,x,eps,desc_a,info,&
it = it + 1 it = it + 1
itx = itx + 1 itx = itx + 1
Call psb_prc_aply(prec,r,z,desc_a,info,work=aux) Call psb_precaply(prec,r,z,desc_a,info,work=aux)
rho_old = rho rho_old = rho
rho = psb_gedot(r,z,desc_a,info) rho = psb_gedot(r,z,desc_a,info)

4
krylov/psb_dcgs.f90
Unescape Escape View File

@ -259,7 +259,7 @@ Subroutine psb_dcgs(a,prec,b,x,eps,desc_a,info,&
End If End If
Call psb_prc_aply(prec,p,f,desc_a,info,work=aux) Call psb_precaply(prec,p,f,desc_a,info,work=aux)
Call psb_spmm(done,a,f,dzero,v,desc_a,info,& Call psb_spmm(done,a,f,dzero,v,desc_a,info,&
& work=aux) & work=aux)
@ -277,7 +277,7 @@ Subroutine psb_dcgs(a,prec,b,x,eps,desc_a,info,&
Call psb_geaxpby(done,uv,dzero,s,desc_a,info) Call psb_geaxpby(done,uv,dzero,s,desc_a,info)
Call psb_geaxpby(done,q,done,s,desc_a,info) Call psb_geaxpby(done,q,done,s,desc_a,info)
Call psb_prc_aply(prec,s,z,desc_a,info,work=aux) Call psb_precaply(prec,s,z,desc_a,info,work=aux)
Call psb_geaxpby(alpha,z,done,x,desc_a,info) Call psb_geaxpby(alpha,z,done,x,desc_a,info)

6
krylov/psb_dcgstab.f90
Unescape Escape View File

@ -285,7 +285,7 @@ Subroutine psb_dcgstab(a,prec,b,x,eps,desc_a,info,&
Call psb_geaxpby(done,r,beta,p,desc_a,info) Call psb_geaxpby(done,r,beta,p,desc_a,info)
End If End If
Call psb_prc_aply(prec,p,f,desc_a,info,work=aux) Call psb_precaply(prec,p,f,desc_a,info,work=aux)
Call psb_spmm(done,a,f,dzero,v,desc_a,info,& Call psb_spmm(done,a,f,dzero,v,desc_a,info,&
& work=aux) & work=aux)
@ -308,9 +308,9 @@ Subroutine psb_dcgstab(a,prec,b,x,eps,desc_a,info,&
goto 9999 goto 9999
end if end if
Call psb_prc_aply(prec,s,z,desc_a,info,work=aux) Call psb_precaply(prec,s,z,desc_a,info,work=aux)
if(info.ne.0) then if(info.ne.0) then
call psb_errpush(4010,name,a_err='psb_prc_aply') call psb_errpush(4010,name,a_err='psb_precaply')
goto 9999 goto 9999
end if end if

6
krylov/psb_dcgstabl.f90
Unescape Escape View File

@ -225,7 +225,7 @@ Subroutine psb_dcgstabl(a,prec,b,x,eps,desc_a,info,&
Call psb_geaxpby(done,b,dzero,r,desc_a,info) Call psb_geaxpby(done,b,dzero,r,desc_a,info)
Call psb_spmm(-done,a,x,done,r,desc_a,info,work=aux) Call psb_spmm(-done,a,x,done,r,desc_a,info,work=aux)
call psb_prc_aply(prec,r,desc_a,info) call psb_precaply(prec,r,desc_a,info)
Call psb_geaxpby(done,r,dzero,rt0,desc_a,info) Call psb_geaxpby(done,r,dzero,rt0,desc_a,info)
Call psb_geaxpby(done,r,dzero,rh(:,0),desc_a,info) Call psb_geaxpby(done,r,dzero,rh(:,0),desc_a,info)
@ -285,7 +285,7 @@ Subroutine psb_dcgstabl(a,prec,b,x,eps,desc_a,info,&
If (debug) Write(0,*) 'bicg part: ',rh(1,0),beta If (debug) Write(0,*) 'bicg part: ',rh(1,0),beta
Call psb_spmm(done,a,uh(:,j),dzero,uh(:,j+1),desc_a,info,work=aux) Call psb_spmm(done,a,uh(:,j),dzero,uh(:,j+1),desc_a,info,work=aux)
call psb_prc_aply(prec,uh(:,j+1),desc_a,info) call psb_precaply(prec,uh(:,j+1),desc_a,info)
gamma(j) = psb_gedot(uh(:,j+1),rt0,desc_a,info) gamma(j) = psb_gedot(uh(:,j+1),rt0,desc_a,info)
If (gamma(j)==dzero) Then If (gamma(j)==dzero) Then
@ -299,7 +299,7 @@ Subroutine psb_dcgstabl(a,prec,b,x,eps,desc_a,info,&
Call psb_geaxpby(alpha,uh(:,0),done,x,desc_a,info) Call psb_geaxpby(alpha,uh(:,0),done,x,desc_a,info)
Call psb_spmm(done,a,rh(:,j),dzero,rh(:,j+1),desc_a,info,work=aux) Call psb_spmm(done,a,rh(:,j),dzero,rh(:,j+1),desc_a,info,work=aux)
call psb_prc_aply(prec,rh(:,j+1),desc_a,info) call psb_precaply(prec,rh(:,j+1),desc_a,info)
Enddo Enddo

4
krylov/psb_dgmresr.f90
Unescape Escape View File

@ -214,7 +214,7 @@ Subroutine psb_dgmresr(a,prec,b,x,eps,desc_a,info,&
End If End If
Call psb_spmm(-done,a,x,done,v(:,1),desc_a,info,work=aux) Call psb_spmm(-done,a,x,done,v(:,1),desc_a,info,work=aux)
call psb_prc_aply(prec,v(:,1),desc_a,info) call psb_precaply(prec,v(:,1),desc_a,info)
rs(1) = psb_genrm2(v(:,1),desc_a,info) rs(1) = psb_genrm2(v(:,1),desc_a,info)
if (info.ne.0) Then if (info.ne.0) Then
info=4011 info=4011
@ -255,7 +255,7 @@ Subroutine psb_dgmresr(a,prec,b,x,eps,desc_a,info,&
itx = itx + 1 itx = itx + 1
Call psb_spmm(done,a,v(:,i),dzero,w,desc_a,info,work=aux) Call psb_spmm(done,a,v(:,i),dzero,w,desc_a,info,work=aux)
call psb_prc_aply(prec,w,desc_a,info) call psb_precaply(prec,w,desc_a,info)
do k = 1, i do k = 1, i
h(k,i) = psb_gedot(v(:,k),w,desc_a,info) h(k,i) = psb_gedot(v(:,k),w,desc_a,info)

4
krylov/psb_zcgs.f90
Unescape Escape View File

@ -262,7 +262,7 @@ Subroutine psb_zcgs(a,prec,b,x,eps,desc_a,info,&
End If End If
Call psb_prc_aply(prec,p,f,desc_a,info,work=aux) Call psb_precaply(prec,p,f,desc_a,info,work=aux)
Call psb_spmm(zone,a,f,zzero,v,desc_a,info,& Call psb_spmm(zone,a,f,zzero,v,desc_a,info,&
& work=aux) & work=aux)
@ -280,7 +280,7 @@ Subroutine psb_zcgs(a,prec,b,x,eps,desc_a,info,&
Call psb_geaxpby(zone,uv,zzero,s,desc_a,info) Call psb_geaxpby(zone,uv,zzero,s,desc_a,info)
Call psb_geaxpby(zone,q,zone,s,desc_a,info) Call psb_geaxpby(zone,q,zone,s,desc_a,info)
Call psb_prc_aply(prec,s,z,desc_a,info,work=aux) Call psb_precaply(prec,s,z,desc_a,info,work=aux)
Call psb_geaxpby(alpha,z,zone,x,desc_a,info) Call psb_geaxpby(alpha,z,zone,x,desc_a,info)

6
krylov/psb_zcgstab.f90
Unescape Escape View File

@ -274,7 +274,7 @@ Subroutine psb_zcgstab(a,prec,b,x,eps,desc_a,info,&
Call psb_geaxpby(zone,r,beta,p,desc_a,info) Call psb_geaxpby(zone,r,beta,p,desc_a,info)
End If End If
Call psb_prc_aply(prec,p,f,desc_a,info,work=aux) Call psb_precaply(prec,p,f,desc_a,info,work=aux)
Call psb_spmm(zone,a,f,zzero,v,desc_a,info,& Call psb_spmm(zone,a,f,zzero,v,desc_a,info,&
& work=aux) & work=aux)
@ -297,9 +297,9 @@ Subroutine psb_zcgstab(a,prec,b,x,eps,desc_a,info,&
goto 9999 goto 9999
end if end if
Call psb_prc_aply(prec,s,z,desc_a,info,work=aux) Call psb_precaply(prec,s,z,desc_a,info,work=aux)
if(info.ne.0) then if(info.ne.0) then
call psb_errpush(4010,name,a_err='psb_prc_aply') call psb_errpush(4010,name,a_err='psb_precaply')
goto 9999 goto 9999
end if end if

Write Preview
Loading…
Cancel
Save