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

Fix SHFT implementation, step 2

Browse Source
maint-3.8.1
Salvatore Filippone 1 year ago
parent 250a6300ba
commit 25e9183e50
12 changed files with 16 additions and 16 deletions
Show Stats Download Patch File Download Diff File

4
prec/impl/psb_c_ilu0_fact.f90
Unescape Escape View File

@ -390,14 +390,14 @@ contains
! into lval/d(i)/uval ! into lval/d(i)/uval
! !
call ilu_copyin(i,ma,a,i,ione,m,l1,lja,lval,& call ilu_copyin(i,ma,a,i,ione,m,l1,lja,lval,&
& d(i),l2,uja,uval,ktrw,trw,upd,shft) & d(i),l2,uja,uval,ktrw,trw,upd,shft_)
else else
! !
! Copy the i-th local row of the matrix, stored in b ! Copy the i-th local row of the matrix, stored in b
! (as (i-ma)-th row), into lval/d(i)/uval ! (as (i-ma)-th row), into lval/d(i)/uval
! !
call ilu_copyin(i-ma,mb,b,i,ione,m,l1,lja,lval,& call ilu_copyin(i-ma,mb,b,i,ione,m,l1,lja,lval,&
& d(i),l2,uja,uval,ktrw,trw,upd,shft) & d(i),l2,uja,uval,ktrw,trw,upd,shft_)
endif endif
lirp(i+1) = l1 + 1 lirp(i+1) = l1 + 1

2
prec/impl/psb_c_iluk_fact.f90
Unescape Escape View File

@ -191,7 +191,7 @@ subroutine psb_ciluk_fact(fill_in,ialg,a,l,u,d,info,blck,shft)
! Compute the ILU(k) or the MILU(k) factorization, depending on ialg ! Compute the ILU(k) or the MILU(k) factorization, depending on ialg
! !
call psb_ciluk_factint(fill_in,ialg,a,blck_,& call psb_ciluk_factint(fill_in,ialg,a,blck_,&
& d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,shft) & d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,shft_)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='psb_ciluk_factint' ch_err='psb_ciluk_factint'

2
prec/impl/psb_c_ilut_fact.f90
Unescape Escape View File

@ -213,7 +213,7 @@ subroutine psb_cilut_fact(fill_in,thres,a,l,u,d,info,blck,iscale,shft)
! Compute the ILU(k,t) factorization ! Compute the ILU(k,t) factorization
! !
call psb_cilut_factint(fill_in,thres,a,blck_,& call psb_cilut_factint(fill_in,thres,a,blck_,&
& d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,scale,shft) & d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,scale,shft_)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='psb_cilut_factint' ch_err='psb_cilut_factint'

4
prec/impl/psb_d_ilu0_fact.f90
Unescape Escape View File

@ -390,14 +390,14 @@ contains
! into lval/d(i)/uval ! into lval/d(i)/uval
! !
call ilu_copyin(i,ma,a,i,ione,m,l1,lja,lval,& call ilu_copyin(i,ma,a,i,ione,m,l1,lja,lval,&
& d(i),l2,uja,uval,ktrw,trw,upd,shft) & d(i),l2,uja,uval,ktrw,trw,upd,shft_)
else else
! !
! Copy the i-th local row of the matrix, stored in b ! Copy the i-th local row of the matrix, stored in b
! (as (i-ma)-th row), into lval/d(i)/uval ! (as (i-ma)-th row), into lval/d(i)/uval
! !
call ilu_copyin(i-ma,mb,b,i,ione,m,l1,lja,lval,& call ilu_copyin(i-ma,mb,b,i,ione,m,l1,lja,lval,&
& d(i),l2,uja,uval,ktrw,trw,upd,shft) & d(i),l2,uja,uval,ktrw,trw,upd,shft_)
endif endif
lirp(i+1) = l1 + 1 lirp(i+1) = l1 + 1

2
prec/impl/psb_d_iluk_fact.f90
Unescape Escape View File

@ -191,7 +191,7 @@ subroutine psb_diluk_fact(fill_in,ialg,a,l,u,d,info,blck,shft)
! Compute the ILU(k) or the MILU(k) factorization, depending on ialg ! Compute the ILU(k) or the MILU(k) factorization, depending on ialg
! !
call psb_diluk_factint(fill_in,ialg,a,blck_,& call psb_diluk_factint(fill_in,ialg,a,blck_,&
& d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,shft) & d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,shft_)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='psb_diluk_factint' ch_err='psb_diluk_factint'

2
prec/impl/psb_d_ilut_fact.f90
Unescape Escape View File

@ -213,7 +213,7 @@ subroutine psb_dilut_fact(fill_in,thres,a,l,u,d,info,blck,iscale,shft)
! Compute the ILU(k,t) factorization ! Compute the ILU(k,t) factorization
! !
call psb_dilut_factint(fill_in,thres,a,blck_,& call psb_dilut_factint(fill_in,thres,a,blck_,&
& d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,scale,shft) & d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,scale,shft_)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='psb_dilut_factint' ch_err='psb_dilut_factint'

4
prec/impl/psb_s_ilu0_fact.f90
Unescape Escape View File

@ -390,14 +390,14 @@ contains
! into lval/d(i)/uval ! into lval/d(i)/uval
! !
call ilu_copyin(i,ma,a,i,ione,m,l1,lja,lval,& call ilu_copyin(i,ma,a,i,ione,m,l1,lja,lval,&
& d(i),l2,uja,uval,ktrw,trw,upd,shft) & d(i),l2,uja,uval,ktrw,trw,upd,shft_)
else else
! !
! Copy the i-th local row of the matrix, stored in b ! Copy the i-th local row of the matrix, stored in b
! (as (i-ma)-th row), into lval/d(i)/uval ! (as (i-ma)-th row), into lval/d(i)/uval
! !
call ilu_copyin(i-ma,mb,b,i,ione,m,l1,lja,lval,& call ilu_copyin(i-ma,mb,b,i,ione,m,l1,lja,lval,&
& d(i),l2,uja,uval,ktrw,trw,upd,shft) & d(i),l2,uja,uval,ktrw,trw,upd,shft_)
endif endif
lirp(i+1) = l1 + 1 lirp(i+1) = l1 + 1

2
prec/impl/psb_s_iluk_fact.f90
Unescape Escape View File

@ -191,7 +191,7 @@ subroutine psb_siluk_fact(fill_in,ialg,a,l,u,d,info,blck,shft)
! Compute the ILU(k) or the MILU(k) factorization, depending on ialg ! Compute the ILU(k) or the MILU(k) factorization, depending on ialg
! !
call psb_siluk_factint(fill_in,ialg,a,blck_,& call psb_siluk_factint(fill_in,ialg,a,blck_,&
& d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,shft) & d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,shft_)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='psb_siluk_factint' ch_err='psb_siluk_factint'

2
prec/impl/psb_s_ilut_fact.f90
Unescape Escape View File

@ -213,7 +213,7 @@ subroutine psb_silut_fact(fill_in,thres,a,l,u,d,info,blck,iscale,shft)
! Compute the ILU(k,t) factorization ! Compute the ILU(k,t) factorization
! !
call psb_silut_factint(fill_in,thres,a,blck_,& call psb_silut_factint(fill_in,thres,a,blck_,&
& d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,scale,shft) & d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,scale,shft_)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='psb_silut_factint' ch_err='psb_silut_factint'

4
prec/impl/psb_z_ilu0_fact.f90
Unescape Escape View File

@ -390,14 +390,14 @@ contains
! into lval/d(i)/uval ! into lval/d(i)/uval
! !
call ilu_copyin(i,ma,a,i,ione,m,l1,lja,lval,& call ilu_copyin(i,ma,a,i,ione,m,l1,lja,lval,&
& d(i),l2,uja,uval,ktrw,trw,upd,shft) & d(i),l2,uja,uval,ktrw,trw,upd,shft_)
else else
! !
! Copy the i-th local row of the matrix, stored in b ! Copy the i-th local row of the matrix, stored in b
! (as (i-ma)-th row), into lval/d(i)/uval ! (as (i-ma)-th row), into lval/d(i)/uval
! !
call ilu_copyin(i-ma,mb,b,i,ione,m,l1,lja,lval,& call ilu_copyin(i-ma,mb,b,i,ione,m,l1,lja,lval,&
& d(i),l2,uja,uval,ktrw,trw,upd,shft) & d(i),l2,uja,uval,ktrw,trw,upd,shft_)
endif endif
lirp(i+1) = l1 + 1 lirp(i+1) = l1 + 1

2
prec/impl/psb_z_iluk_fact.f90
Unescape Escape View File

@ -191,7 +191,7 @@ subroutine psb_ziluk_fact(fill_in,ialg,a,l,u,d,info,blck,shft)
! Compute the ILU(k) or the MILU(k) factorization, depending on ialg ! Compute the ILU(k) or the MILU(k) factorization, depending on ialg
! !
call psb_ziluk_factint(fill_in,ialg,a,blck_,& call psb_ziluk_factint(fill_in,ialg,a,blck_,&
& d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,shft) & d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,shft_)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='psb_ziluk_factint' ch_err='psb_ziluk_factint'

2
prec/impl/psb_z_ilut_fact.f90
Unescape Escape View File

@ -213,7 +213,7 @@ subroutine psb_zilut_fact(fill_in,thres,a,l,u,d,info,blck,iscale,shft)
! Compute the ILU(k,t) factorization ! Compute the ILU(k,t) factorization
! !
call psb_zilut_factint(fill_in,thres,a,blck_,& call psb_zilut_factint(fill_in,thres,a,blck_,&
& d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,scale,shft) & d,ll%val,ll%ja,ll%irp,uu%val,uu%ja,uu%irp,l1,l2,info,scale,shft_)
if (info /= psb_success_) then if (info /= psb_success_) then
info=psb_err_from_subroutine_ info=psb_err_from_subroutine_
ch_err='psb_zilut_factint' ch_err='psb_zilut_factint'

Write Preview
Loading…
Cancel
Save