fdurastante/amg4psblas
1
0
Fork 0
mirror of https://github.com/sfilippone/amg4psblas.git synced 2026-10-06 22:55:12 +00:00
Code Issues Packages Projects Releases Wiki Activity

mld2p4-2:

Browse Source 
configure
 examples/fileread/data_input.f90
 examples/fileread/mld_cexample_1lev.f90
 examples/fileread/mld_cexample_ml.f90
 examples/fileread/mld_dexample_1lev.f90
 examples/fileread/mld_dexample_ml.f90
 examples/fileread/mld_sexample_1lev.f90
 examples/fileread/mld_sexample_ml.f90
 examples/fileread/mld_zexample_1lev.f90
 examples/fileread/mld_zexample_ml.f90
 examples/pdegen/data_input.f90
 examples/pdegen/mld_dexample_1lev.f90
 examples/pdegen/mld_dexample_ml.f90
 examples/pdegen/mld_sexample_1lev.f90
 examples/pdegen/mld_sexample_ml.f90
 mlprec/Makefile
 mlprec/mld_base_prec_type.f90
 mlprec/mld_c_as_smoother.f90
 mlprec/mld_c_diag_solver.f90
 mlprec/mld_c_id_solver.f90
 mlprec/mld_c_ilu_solver.f90
 mlprec/mld_c_inner_mod.f90
 mlprec/mld_c_jac_smoother.f90
 mlprec/mld_c_move_alloc_mod.f90
 mlprec/mld_c_prec_mod.f90
 mlprec/mld_c_prec_type.f90
 mlprec/mld_c_slu_solver.f90
 mlprec/mld_caggrmap_bld.f90
 mlprec/mld_caggrmat_asb.f90
 mlprec/mld_caggrmat_nosmth_asb.F90
 mlprec/mld_caggrmat_smth_asb.F90
 mlprec/mld_ccoarse_bld.f90
 mlprec/mld_cilu0_fact.f90
 mlprec/mld_ciluk_fact.f90
 mlprec/mld_cilut_fact.f90
 mlprec/mld_cmlprec_aply.f90
 mlprec/mld_cmlprec_bld.f90
 mlprec/mld_cprecaply.f90
 mlprec/mld_cprecbld.f90
 mlprec/mld_cprecinit.F90
 mlprec/mld_cprecset.F90
 mlprec/mld_cslu_bld.f90
 mlprec/mld_cslud_bld.f90
 mlprec/mld_csp_renum.f90
 mlprec/mld_cumf_bld.f90
 mlprec/mld_d_as_smoother.f90
 mlprec/mld_d_diag_solver.f90
 mlprec/mld_d_id_solver.f90
 mlprec/mld_d_ilu_solver.f90
 mlprec/mld_d_inner_mod.f90
 mlprec/mld_d_jac_smoother.f90
 mlprec/mld_d_move_alloc_mod.f90
 mlprec/mld_d_prec_mod.f90
 mlprec/mld_d_prec_type.f90
 mlprec/mld_d_slu_solver.f90
 mlprec/mld_d_sludist_solver.f90
 mlprec/mld_d_umf_solver.f90
 mlprec/mld_daggrmap_bld.f90
 mlprec/mld_daggrmat_asb.f90
 mlprec/mld_daggrmat_minnrg_asb.F90
 mlprec/mld_daggrmat_nosmth_asb.F90
 mlprec/mld_daggrmat_smth_asb.F90
 mlprec/mld_dcoarse_bld.f90
 mlprec/mld_dilu0_fact.f90
 mlprec/mld_diluk_fact.f90
 mlprec/mld_dilut_fact.f90
 mlprec/mld_dmlprec_aply.f90
 mlprec/mld_dmlprec_bld.f90
 mlprec/mld_dprecaply.f90
 mlprec/mld_dprecbld.f90
 mlprec/mld_dprecinit.F90
 mlprec/mld_dprecset.F90
 mlprec/mld_dslu_bld.f90
 mlprec/mld_dslud_bld.f90
 mlprec/mld_dsp_renum.f90
 mlprec/mld_s_as_smoother.f90
 mlprec/mld_s_diag_solver.f90
 mlprec/mld_s_id_solver.f90
 mlprec/mld_s_ilu_solver.f90
 mlprec/mld_s_inner_mod.f90
 mlprec/mld_s_jac_smoother.f90
 mlprec/mld_s_move_alloc_mod.f90
 mlprec/mld_s_prec_mod.f90
 mlprec/mld_s_prec_type.f90
 mlprec/mld_s_slu_solver.f90
 mlprec/mld_saggrmap_bld.f90
 mlprec/mld_saggrmat_asb.f90
 mlprec/mld_saggrmat_nosmth_asb.F90
 mlprec/mld_saggrmat_smth_asb.F90
 mlprec/mld_scoarse_bld.f90
 mlprec/mld_silu0_fact.f90
 mlprec/mld_siluk_fact.f90
 mlprec/mld_silut_fact.f90
 mlprec/mld_smlprec_aply.f90
 mlprec/mld_smlprec_bld.f90
 mlprec/mld_sprecaply.f90
 mlprec/mld_sprecbld.f90
 mlprec/mld_sprecinit.F90
 mlprec/mld_sprecset.F90
 mlprec/mld_sslu_bld.f90
 mlprec/mld_sslud_bld.f90
 mlprec/mld_ssp_renum.f90
 mlprec/mld_sumf_bld.f90
 mlprec/mld_z_as_smoother.f90
 mlprec/mld_z_diag_solver.f90
 mlprec/mld_z_id_solver.f90
 mlprec/mld_z_ilu_solver.f90
 mlprec/mld_z_inner_mod.f90
 mlprec/mld_z_jac_smoother.f90
 mlprec/mld_z_move_alloc_mod.f90
 mlprec/mld_z_prec_mod.f90
 mlprec/mld_z_prec_type.f90
 mlprec/mld_z_slu_solver.f90
 mlprec/mld_z_umf_solver.f90
 mlprec/mld_zaggrmap_bld.f90
 mlprec/mld_zaggrmat_asb.f90
 mlprec/mld_zaggrmat_nosmth_asb.F90
 mlprec/mld_zaggrmat_smth_asb.F90
 mlprec/mld_zas_aply.f90
 mlprec/mld_zas_bld.f90
 mlprec/mld_zbaseprec_aply.f90
 mlprec/mld_zbaseprec_bld.f90
 mlprec/mld_zcoarse_bld.f90
 mlprec/mld_zdiag_bld.f90
 mlprec/mld_zfact_bld.f90
 mlprec/mld_zilu0_fact.f90
 mlprec/mld_zilu_bld.f90
 mlprec/mld_ziluk_fact.f90
 mlprec/mld_zilut_fact.f90
 mlprec/mld_zmlprec_aply.f90
 mlprec/mld_zmlprec_bld.f90
 mlprec/mld_zprecaply.f90
 mlprec/mld_zprecbld.f90
 mlprec/mld_zprecinit.F90
 mlprec/mld_zprecset.F90
 mlprec/mld_zslu_bld.f90
 mlprec/mld_zslud_bld.f90
 mlprec/mld_zsp_renum.f90
 mlprec/mld_zumf_bld.f90
 tests/fileread/cf_sample.f90
 tests/fileread/data_input.f90
 tests/fileread/df_sample.f90
 tests/fileread/sf_sample.f90
 tests/fileread/zf_sample.f90
 tests/newslv/data_input.f90
 tests/newslv/mld_d_tlu_solver.f90
 tests/newslv/ppde.f90
 tests/newslv/spde.f90
 tests/pdegen/data_input.f90
 tests/pdegen/ppde.f90
 tests/pdegen/runs/ppde.inp
 tests/pdegen/spde.f90

Adapted psb_base_mod.
This commit is contained in:
Salvatore Filippone
2011-03-21 17:08:37 +00:00
parent 53d6992dea
commit 5fc55bfc4c
152 changed files with 875 additions and 861 deletions
Download Patch File Download Diff File Expand all files Collapse all files
tests/newslv/data_input.f90
+4 -4
View File
@@ -73,7 +73,7 @@ contains
end subroutine read_int
subroutine read_single(val,file,marker)
use psb_sparse_mod
use psb_base_mod
real(psb_spk_), intent(out) :: val
integer, intent(in) :: file
character(len=1), optional, intent(in) :: marker
@@ -83,7 +83,7 @@ contains
end subroutine read_single
subroutine read_double(val,file,marker)
use psb_sparse_mod
use psb_base_mod
real(psb_dpk_), intent(out) :: val
integer, intent(in) :: file
character(len=1), optional, intent(in) :: marker
@@ -131,7 +131,7 @@ contains
read(charbuf(1:idx-1),*) val
end subroutine string_read_int
subroutine string_read_single(val,file,marker)
use psb_sparse_mod
use psb_base_mod
real(psb_spk_), intent(out) :: val
character(len=*), intent(in) :: file
character(len=1), optional, intent(in) :: marker
@@ -150,7 +150,7 @@ contains
read(charbuf(1:idx-1),*) val
end subroutine string_read_single
subroutine string_read_double(val,file,marker)
use psb_sparse_mod
use psb_base_mod
real(psb_dpk_), intent(out) :: val
character(len=*), intent(in) :: file
character(len=1), optional, intent(in) :: marker
tests/newslv/mld_d_tlu_solver.f90
+14 -14
View File
@@ -75,7 +75,7 @@ module mld_d_tlu_solver
interface mld_ilu0_fact
subroutine mld_dilu0_fact(ialg,a,l,u,d,info,blck,upd)
use psb_sparse_mod, only : psb_dspmat_type, psb_dpk_
use psb_base_mod, only : psb_dspmat_type, psb_dpk_
integer, intent(in) :: ialg
integer, intent(out) :: info
type(psb_dspmat_type),intent(in) :: a
@@ -88,7 +88,7 @@ module mld_d_tlu_solver
interface mld_iluk_fact
subroutine mld_diluk_fact(fill_in,ialg,a,l,u,d,info,blck)
use psb_sparse_mod, only : psb_dspmat_type, psb_dpk_
use psb_base_mod, only : psb_dspmat_type, psb_dpk_
integer, intent(in) :: fill_in,ialg
integer, intent(out) :: info
type(psb_dspmat_type),intent(in) :: a
@@ -100,7 +100,7 @@ module mld_d_tlu_solver
interface mld_ilut_fact
subroutine mld_dilut_fact(fill_in,thres,a,l,u,d,info,blck)
use psb_sparse_mod, only : psb_dspmat_type, psb_dpk_
use psb_base_mod, only : psb_dspmat_type, psb_dpk_
integer, intent(in) :: fill_in
real(psb_dpk_), intent(in) :: thres
integer, intent(out) :: info
@@ -124,7 +124,7 @@ contains
subroutine d_tlu_solver_default(sv)
use psb_sparse_mod
use psb_base_mod
Implicit None
@@ -140,7 +140,7 @@ contains
subroutine d_tlu_solver_check(sv,info)
use psb_sparse_mod
use psb_base_mod
Implicit None
@@ -181,7 +181,7 @@ contains
subroutine d_tlu_solver_apply(alpha,sv,x,beta,y,desc_data,trans,work,info)
use psb_sparse_mod
use psb_base_mod
type(psb_desc_type), intent(in) :: desc_data
class(mld_d_tlu_solver_type), intent(in) :: sv
real(psb_dpk_),intent(in) :: x(:)
@@ -287,7 +287,7 @@ contains
subroutine d_tlu_solver_bld(a,desc_a,sv,upd,info,b)
use psb_sparse_mod
use psb_base_mod
Implicit None
@@ -459,7 +459,7 @@ contains
subroutine d_tlu_solver_seti(sv,what,val,info)
use psb_sparse_mod
use psb_base_mod
Implicit None
@@ -498,7 +498,7 @@ contains
subroutine d_tlu_solver_setc(sv,what,val,info)
use psb_sparse_mod
use psb_base_mod
Implicit None
@@ -536,7 +536,7 @@ contains
subroutine d_tlu_solver_setr(sv,what,val,info)
use psb_sparse_mod
use psb_base_mod
Implicit None
@@ -574,7 +574,7 @@ contains
subroutine d_tlu_solver_free(sv,info)
use psb_sparse_mod
use psb_base_mod
Implicit None
@@ -612,7 +612,7 @@ contains
subroutine d_tlu_solver_descr(sv,info,iout,coarse)
use psb_sparse_mod
use psb_base_mod
Implicit None
@@ -660,7 +660,7 @@ contains
end subroutine d_tlu_solver_descr
function d_tlu_solver_sizeof(sv) result(val)
use psb_sparse_mod
use psb_base_mod
implicit none
! Arguments
class(mld_d_tlu_solver_type), intent(in) :: sv
@@ -676,7 +676,7 @@ contains
end function d_tlu_solver_sizeof
subroutine d_tlu_solver_dmp(sv,ictxt,level,info,prefix,head,solver)
use psb_sparse_mod
use psb_base_mod
implicit none
class(mld_d_tlu_solver_type), intent(in) :: sv
integer, intent(in) :: ictxt,level
tests/newslv/ppde.f90
+9 -9
View File
@@ -69,7 +69,7 @@
! Note that if a1=a2=a3=a4=0., the PDE is the well-known Laplace equation.
!
program ppde
use psb_sparse_mod
use psb_base_mod
use mld_prec_mod
use psb_krylov_mod
use psb_util_mod
@@ -429,7 +429,7 @@ contains
!
! Note that if a1=a2=a3=a4=0., the PDE is the well-known Laplace equation.
!
use psb_sparse_mod
use psb_base_mod
implicit none
integer :: idim
integer, parameter :: nb=20
@@ -706,43 +706,43 @@ end program ppde
! functions parametrizing the differential equation
!
function a1(x,y,z)
use psb_sparse_mod, only : psb_dpk_
use psb_base_mod, only : psb_dpk_
real(psb_dpk_) :: a1
real(psb_dpk_) :: x,y,z
a1=1.d0
end function a1
function a2(x,y,z)
use psb_sparse_mod, only : psb_dpk_
use psb_base_mod, only : psb_dpk_
real(psb_dpk_) :: a2
real(psb_dpk_) :: x,y,z
a2=2.d1*y
end function a2
function a3(x,y,z)
use psb_sparse_mod, only : psb_dpk_
use psb_base_mod, only : psb_dpk_
real(psb_dpk_) :: a3
real(psb_dpk_) :: x,y,z
a3=1.d0
end function a3
function a4(x,y,z)
use psb_sparse_mod, only : psb_dpk_
use psb_base_mod, only : psb_dpk_
real(psb_dpk_) :: a4
real(psb_dpk_) :: x,y,z
a4=1.d0
end function a4
function b1(x,y,z)
use psb_sparse_mod, only : psb_dpk_
use psb_base_mod, only : psb_dpk_
real(psb_dpk_) :: b1
real(psb_dpk_) :: x,y,z
b1=1.d0
end function b1
function b2(x,y,z)
use psb_sparse_mod, only : psb_dpk_
use psb_base_mod, only : psb_dpk_
real(psb_dpk_) :: b2
real(psb_dpk_) :: x,y,z
b2=1.d0
end function b2
function b3(x,y,z)
use psb_sparse_mod, only : psb_dpk_
use psb_base_mod, only : psb_dpk_
real(psb_dpk_) :: b3
real(psb_dpk_) :: x,y,z
b3=1.d0
tests/newslv/spde.f90
+9 -9
View File
@@ -69,7 +69,7 @@
! Note that if a1=a2=a3=a4=0., the PDE is the well-known Laplace equation.
!
program spde
use psb_sparse_mod
use psb_base_mod
use mld_prec_mod
use psb_krylov_mod
use psb_util_mod
@@ -423,7 +423,7 @@ contains
!
! Note that if a1=a2=a3=a4=0., the PDE is the well-known Laplace equation.
!
use psb_sparse_mod
use psb_base_mod
implicit none
integer :: idim
integer, parameter :: nb=20
@@ -700,43 +700,43 @@ end program spde
! functions parametrizing the differential equation
!
function a1(x,y,z)
use psb_sparse_mod, only : psb_spk_
use psb_base_mod, only : psb_spk_
real(psb_spk_) :: a1
real(psb_spk_) :: x,y,z
a1=1.e0
end function a1
function a2(x,y,z)
use psb_sparse_mod, only : psb_spk_
use psb_base_mod, only : psb_spk_
real(psb_spk_) :: a2
real(psb_spk_) :: x,y,z
a2=2.e1*y
end function a2
function a3(x,y,z)
use psb_sparse_mod, only : psb_spk_
use psb_base_mod, only : psb_spk_
real(psb_spk_) :: a3
real(psb_spk_) :: x,y,z
a3=1.e0
end function a3
function a4(x,y,z)
use psb_sparse_mod, only : psb_spk_
use psb_base_mod, only : psb_spk_
real(psb_spk_) :: a4
real(psb_spk_) :: x,y,z
a4=1.e0
end function a4
function b1(x,y,z)
use psb_sparse_mod, only : psb_spk_
use psb_base_mod, only : psb_spk_
real(psb_spk_) :: b1
real(psb_spk_) :: x,y,z
b1=1.e0
end function b1
function b2(x,y,z)
use psb_sparse_mod, only : psb_spk_
use psb_base_mod, only : psb_spk_
real(psb_spk_) :: b2
real(psb_spk_) :: x,y,z
b2=1.e0
end function b2
function b3(x,y,z)
use psb_sparse_mod, only : psb_spk_
use psb_base_mod, only : psb_spk_
real(psb_spk_) :: b3
real(psb_spk_) :: x,y,z
b3=1.e0