prec/psb_cilu_fct.f90
 prec/psb_dilu_fct.f90
 prec/psb_silu_fct.f90
 prec/psb_zilu_fct.f90

Workaround for ICE in gfortran. Have to figure out if the original
code was valid or not.
psblas3-type-indexed
Salvatore Filippone 15 years ago
parent 6ed696b359
commit 3c46ea898c

@ -174,7 +174,7 @@ contains
if ((mod(i,nrb) == 1).or.(nrb == 1)) then if ((mod(i,nrb) == 1).or.(nrb == 1)) then
irb = min(ma-i+1,nrb) irb = min(ma-i+1,nrb)
call a%a%csget(i,i+irb-1,trw,info) call aa%csget(i,i+irb-1,trw,info)
if(info /= psb_success_) then if(info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='a%csget' ch_err='a%csget'
@ -312,7 +312,7 @@ contains
if ((mod(i,nrb) == 1).or.(nrb == 1)) then if ((mod(i,nrb) == 1).or.(nrb == 1)) then
irb = min(ma-i+1,nrb) irb = min(ma-i+1,nrb)
call b%a%csget(i-ma,i-ma+irb-1,trw,info) call aa%csget(i-ma,i-ma+irb-1,trw,info)
nz = trw%get_nzeros() nz = trw%get_nzeros()
if(info /= psb_success_) then if(info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_

@ -177,7 +177,7 @@ contains
if ((mod(i,nrb) == 1).or.(nrb == 1)) then if ((mod(i,nrb) == 1).or.(nrb == 1)) then
irb = min(ma-i+1,nrb) irb = min(ma-i+1,nrb)
call a%a%csget(i,i+irb-1,trw,info) call aa%csget(i,i+irb-1,trw,info)
if(info /= psb_success_) then if(info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='a%csget' ch_err='a%csget'
@ -315,7 +315,7 @@ contains
if ((mod(i,nrb) == 1).or.(nrb == 1)) then if ((mod(i,nrb) == 1).or.(nrb == 1)) then
irb = min(ma-i+1,nrb) irb = min(ma-i+1,nrb)
call b%a%csget(i-ma,i-ma+irb-1,trw,info) call aa%csget(i-ma,i-ma+irb-1,trw,info)
nz = trw%get_nzeros() nz = trw%get_nzeros()
if(info /= psb_success_) then if(info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_

@ -176,7 +176,7 @@ contains
if ((mod(i,nrb) == 1).or.(nrb == 1)) then if ((mod(i,nrb) == 1).or.(nrb == 1)) then
irb = min(ma-i+1,nrb) irb = min(ma-i+1,nrb)
call a%a%csget(i,i+irb-1,trw,info) call aa%csget(i,i+irb-1,trw,info)
if(info /= psb_success_) then if(info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='a%csget' ch_err='a%csget'
@ -314,7 +314,7 @@ contains
if ((mod(i,nrb) == 1).or.(nrb == 1)) then if ((mod(i,nrb) == 1).or.(nrb == 1)) then
irb = min(ma-i+1,nrb) irb = min(ma-i+1,nrb)
call b%a%csget(i-ma,i-ma+irb-1,trw,info) call aa%csget(i-ma,i-ma+irb-1,trw,info)
nz = trw%get_nzeros() nz = trw%get_nzeros()
if(info /= psb_success_) then if(info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_

@ -174,7 +174,7 @@ contains
if ((mod(i,nrb) == 1).or.(nrb == 1)) then if ((mod(i,nrb) == 1).or.(nrb == 1)) then
irb = min(ma-i+1,nrb) irb = min(ma-i+1,nrb)
call a%a%csget(i,i+irb-1,trw,info) call aa%csget(i,i+irb-1,trw,info)
if(info /= psb_success_) then if(info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='a%csget' ch_err='a%csget'
@ -312,7 +312,7 @@ contains
if ((mod(i,nrb) == 1).or.(nrb == 1)) then if ((mod(i,nrb) == 1).or.(nrb == 1)) then
irb = min(ma-i+1,nrb) irb = min(ma-i+1,nrb)
call b%a%csget(i-ma,i-ma+irb-1,trw,info) call aa%csget(i-ma,i-ma+irb-1,trw,info)
nz = trw%get_nzeros() nz = trw%get_nzeros()
if(info /= psb_success_) then if(info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_

Loading…
Cancel
Save