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

Leading dimension of RS in call to dtrsm

Browse Source
psblas3-type-indexed
Salvatore Filippone 19 years ago
parent 008249e5e4
commit 6212203102
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File

4
src/methd/psb_dgmresr.f90
Unescape Escape View File

@ -300,7 +300,7 @@ Subroutine psb_dgmresr(a,prec,b,x,eps,desc_a,info,&
endif endif
if (rerr < eps ) then if (rerr < eps ) then
call dtrsm('l','u','n','n',i,1,done,rr,size(rr,1),rs,nl) call dtrsm('l','u','n','n',i,1,done,rr,size(rr,1),rs,size(rs))
if (debug) write(0,*) 'Rebuild x-> RS:',rs(21:nl) if (debug) write(0,*) 'Rebuild x-> RS:',rs(21:nl)
do k=1, i do k=1, i
call psb_geaxpby(rs(k),v(:,k),done,x,desc_a,info) call psb_geaxpby(rs(k),v(:,k),done,x,desc_a,info)
@ -314,7 +314,7 @@ Subroutine psb_dgmresr(a,prec,b,x,eps,desc_a,info,&
end Do inner end Do inner
if (debug) write(0,*) 'Before DTRSM :',rs(1:nl) if (debug) write(0,*) 'Before DTRSM :',rs(1:nl)
call dtrsm('l','u','n','n',nl,1,done,rr,size(rr,1),rs,nl) call dtrsm('l','u','n','n',nl,1,done,rr,size(rr,1),rs,size(rs))
if (debug) write(0,*) 'Rebuild x-> RS:',rs(21:nl) if (debug) write(0,*) 'Rebuild x-> RS:',rs(21:nl)
do k=1, nl do k=1, nl
call psb_geaxpby(rs(k),v(:,k),done,x,desc_a,info) call psb_geaxpby(rs(k),v(:,k),done,x,desc_a,info)

Write Preview
Loading…
Cancel
Save