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

Fix dump of TPROL and POLY smoother

Browse Source
l1-and-0-aggr
sfilippone 2 months ago
parent 68a9cceaa0
commit b704d50df1
6 changed files with 10 additions and 22 deletions
Show Stats Download Patch File Download Diff File

6
amgprec/impl/level/amg_c_base_onelev_dump.f90
Unescape Escape View File

@ -127,8 +127,7 @@ subroutine amg_c_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
ivr = lv%linmap%p_desc_U%get_global_indices(owned=.false.) ivr = lv%linmap%p_desc_U%get_global_indices(owned=.false.)
write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx' write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx'
! !
! This is not implemented yet. call lv%tprol%print(fname,head=head,ivr=ivr)
!call lv%tprol%print(fname,head=head,ivr=ivr)
end if end if
end if end if
else else
@ -151,8 +150,7 @@ subroutine amg_c_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
if (tprol_) then if (tprol_) then
write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx' write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx'
! !
! This is not implemented yet. call lv%tprol%print(fname,head=head)
!call lv%tprol%print(fname,head=head)
end if end if
end if end if
end if end if

6
amgprec/impl/level/amg_d_base_onelev_dump.f90
Unescape Escape View File

@ -127,8 +127,7 @@ subroutine amg_d_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
ivr = lv%linmap%p_desc_U%get_global_indices(owned=.false.) ivr = lv%linmap%p_desc_U%get_global_indices(owned=.false.)
write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx' write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx'
! !
! This is not implemented yet. call lv%tprol%print(fname,head=head,ivr=ivr)
!call lv%tprol%print(fname,head=head,ivr=ivr)
end if end if
end if end if
else else
@ -151,8 +150,7 @@ subroutine amg_d_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
if (tprol_) then if (tprol_) then
write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx' write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx'
! !
! This is not implemented yet. call lv%tprol%print(fname,head=head)
!call lv%tprol%print(fname,head=head)
end if end if
end if end if
end if end if

6
amgprec/impl/level/amg_s_base_onelev_dump.f90
Unescape Escape View File

@ -127,8 +127,7 @@ subroutine amg_s_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
ivr = lv%linmap%p_desc_U%get_global_indices(owned=.false.) ivr = lv%linmap%p_desc_U%get_global_indices(owned=.false.)
write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx' write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx'
! !
! This is not implemented yet. call lv%tprol%print(fname,head=head,ivr=ivr)
!call lv%tprol%print(fname,head=head,ivr=ivr)
end if end if
end if end if
else else
@ -151,8 +150,7 @@ subroutine amg_s_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
if (tprol_) then if (tprol_) then
write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx' write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx'
! !
! This is not implemented yet. call lv%tprol%print(fname,head=head)
!call lv%tprol%print(fname,head=head)
end if end if
end if end if
end if end if

6
amgprec/impl/level/amg_z_base_onelev_dump.f90
Unescape Escape View File

@ -127,8 +127,7 @@ subroutine amg_z_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
ivr = lv%linmap%p_desc_U%get_global_indices(owned=.false.) ivr = lv%linmap%p_desc_U%get_global_indices(owned=.false.)
write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx' write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx'
! !
! This is not implemented yet. call lv%tprol%print(fname,head=head,ivr=ivr)
!call lv%tprol%print(fname,head=head,ivr=ivr)
end if end if
end if end if
else else
@ -151,8 +150,7 @@ subroutine amg_z_base_onelev_dump(lv,level,info,prefix,head,ac,rp,&
if (tprol_) then if (tprol_) then
write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx' write(fname(lname+1:),'(a,i3.3,a)')'_l',level,'_tprol.mtx'
! !
! This is not implemented yet. call lv%tprol%print(fname,head=head)
!call lv%tprol%print(fname,head=head)
end if end if
end if end if
end if end if

4
amgprec/impl/smoother/amg_d_poly_smoother_dmp.f90
Unescape Escape View File

@ -77,11 +77,9 @@ subroutine amg_d_poly_smoother_dmp(sm,desc,level,info,prefix,head,smoother,solve
end if end if
lname = len_trim(prefix_) lname = len_trim(prefix_)
fname = trim(prefix_) fname = trim(prefix_)
write(fname(lname+1:lname+5),'(a,i3.3)') '_poly',iam write(fname(lname+1:lname+8),'(a,i3.3)') '_poly',iam
lname = lname + 8 lname = lname + 8
! to be completed ! to be completed
! At base level do nothing for the smoother ! At base level do nothing for the smoother
if (allocated(sm%sv)) & if (allocated(sm%sv)) &

4
amgprec/impl/smoother/amg_s_poly_smoother_dmp.f90
Unescape Escape View File

@ -77,11 +77,9 @@ subroutine amg_s_poly_smoother_dmp(sm,desc,level,info,prefix,head,smoother,solve
end if end if
lname = len_trim(prefix_) lname = len_trim(prefix_)
fname = trim(prefix_) fname = trim(prefix_)
write(fname(lname+1:lname+5),'(a,i3.3)') '_poly',iam write(fname(lname+1:lname+8),'(a,i3.3)') '_poly',iam
lname = lname + 8 lname = lname + 8
! to be completed ! to be completed
! At base level do nothing for the smoother ! At base level do nothing for the smoother
if (allocated(sm%sv)) & if (allocated(sm%sv)) &

Write Preview
Loading…
Cancel
Save