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

Take out precset implementation, leave only prec%set

Browse Source
newG2L
Salvatore Filippone 4 years ago
parent ed58551ffd
commit 08f3de57f0
5 changed files with 23 additions and 342 deletions
Show Stats Download Patch File Download Diff File

46
prec/impl/Makefile
Unescape Escape View File

@ -8,35 +8,35 @@ OBJS=psb_s_prec_type_impl.o psb_d_prec_type_impl.o \
psb_c_prec_type_impl.o psb_z_prec_type_impl.o \ psb_c_prec_type_impl.o psb_z_prec_type_impl.o \
psb_d_diagprec_impl.o psb_d_bjacprec_impl.o psb_d_nullprec_impl.o \ psb_d_diagprec_impl.o psb_d_bjacprec_impl.o psb_d_nullprec_impl.o \
psb_dilu_fct.o psb_d_ilu0_fact.o psb_d_iluk_fact.o psb_d_ilut_fact.o \ psb_dilu_fct.o psb_d_ilu0_fact.o psb_d_iluk_fact.o psb_d_ilut_fact.o \
psb_dprecbld.o psb_dprecset.o psb_dprecinit.o \ psb_dprecbld.o psb_dprecinit.o \
psb_s_diagprec_impl.o psb_s_bjacprec_impl.o psb_s_nullprec_impl.o \ psb_s_diagprec_impl.o psb_s_bjacprec_impl.o psb_s_nullprec_impl.o \
psb_silu_fct.o psb_s_ilu0_fact.o psb_s_iluk_fact.o psb_s_ilut_fact.o \ psb_silu_fct.o psb_s_ilu0_fact.o psb_s_iluk_fact.o psb_s_ilut_fact.o \
psb_sprecbld.o psb_sprecset.o psb_sprecinit.o \ psb_sprecbld.o psb_sprecinit.o \
psb_c_diagprec_impl.o psb_c_bjacprec_impl.o psb_c_nullprec_impl.o \ psb_c_diagprec_impl.o psb_c_bjacprec_impl.o psb_c_nullprec_impl.o \
psb_cilu_fct.o psb_c_ilu0_fact.o psb_c_iluk_fact.o psb_c_ilut_fact.o \ psb_cilu_fct.o psb_c_ilu0_fact.o psb_c_iluk_fact.o psb_c_ilut_fact.o \
psb_cprecbld.o psb_cprecset.o psb_cprecinit.o \ psb_cprecbld.o psb_cprecinit.o \
psb_z_diagprec_impl.o psb_z_bjacprec_impl.o psb_z_nullprec_impl.o \ psb_z_diagprec_impl.o psb_z_bjacprec_impl.o psb_z_nullprec_impl.o \
psb_zilu_fct.o psb_z_ilu0_fact.o psb_z_iluk_fact.o psb_z_ilut_fact.o \ psb_zilu_fct.o psb_z_ilu0_fact.o psb_z_iluk_fact.o psb_z_ilut_fact.o \
psb_zprecbld.o psb_zprecset.o psb_zprecinit.o \ psb_zprecbld.o psb_zprecinit.o \
psb_c_sparsify.o psb_d_sparsify.o psb_s_sparsify.o psb_z_sparsify.o \ psb_c_sparsify.o psb_d_sparsify.o psb_s_sparsify.o psb_z_sparsify.o \
psb_crwclip.o psb_drwclip.o psb_srwclip.o psb_zrwclip.o \ psb_crwclip.o psb_drwclip.o psb_srwclip.o psb_zrwclip.o \
psb_c_sp_drop.o psb_d_sp_drop.o psb_s_sp_drop.o psb_z_sp_drop.o \ psb_c_sp_drop.o psb_d_sp_drop.o psb_s_sp_drop.o psb_z_sp_drop.o \
psb_dsparse_biconjg_llk_noth.o psb_dsparse_biconjg_llk.o \ psb_dsparse_biconjg_llk_noth.o psb_dsparse_biconjg_llk.o \
psb_dsparse_biconjg_mlk.o psb_dsparse_biconjg_s_ft_llk.o \ psb_dsparse_biconjg_mlk.o psb_dsparse_biconjg_s_ft_llk.o \
psb_dsparse_biconjg_s_llk.o \ psb_dsparse_biconjg_s_llk.o \
psb_csparse_biconjg_llk_noth.o psb_csparse_biconjg_llk.o \ psb_csparse_biconjg_llk_noth.o psb_csparse_biconjg_llk.o \
psb_csparse_biconjg_mlk.o psb_csparse_biconjg_s_ft_llk.o \ psb_csparse_biconjg_mlk.o psb_csparse_biconjg_s_ft_llk.o \
psb_csparse_biconjg_s_llk.o \ psb_csparse_biconjg_s_llk.o \
psb_zsparse_biconjg_llk_noth.o psb_zsparse_biconjg_llk.o \ psb_zsparse_biconjg_llk_noth.o psb_zsparse_biconjg_llk.o \
psb_zsparse_biconjg_mlk.o psb_zsparse_biconjg_s_ft_llk.o \ psb_zsparse_biconjg_mlk.o psb_zsparse_biconjg_s_ft_llk.o \
psb_zsparse_biconjg_s_llk.o \ psb_zsparse_biconjg_s_llk.o \
psb_ssparse_biconjg_llk_noth.o psb_ssparse_biconjg_llk.o \ psb_ssparse_biconjg_llk_noth.o psb_ssparse_biconjg_llk.o \
psb_ssparse_biconjg_mlk.o psb_ssparse_biconjg_s_ft_llk.o \ psb_ssparse_biconjg_mlk.o psb_ssparse_biconjg_s_ft_llk.o \
psb_ssparse_biconjg_s_llk.o \ psb_ssparse_biconjg_s_llk.o \
psb_d_ainv_bld.o psb_c_ainv_bld.o psb_s_ainv_bld.o \ psb_d_ainv_bld.o psb_c_ainv_bld.o psb_s_ainv_bld.o \
psb_z_ainv_bld.o \ psb_z_ainv_bld.o \
psb_c_invt_fact.o psb_d_invt_fact.o psb_s_invt_fact.o psb_z_invt_fact.o\ psb_c_invt_fact.o psb_d_invt_fact.o psb_s_invt_fact.o psb_z_invt_fact.o\
psb_c_invk_fact.o psb_d_invk_fact.o psb_s_invk_fact.o psb_z_invk_fact.o psb_c_invk_fact.o psb_d_invk_fact.o psb_s_invk_fact.o psb_z_invk_fact.o
LIBNAME=$(PRECLIBNAME) LIBNAME=$(PRECLIBNAME)
COBJS= COBJS=

80
prec/impl/psb_cprecset.f90
Unescape Escape View File

@ -1,80 +0,0 @@
!
! Parallel Sparse BLAS version 3.5
! (C) Copyright 2006-2018
! Salvatore Filippone
! Alfredo Buttari
!
! Redistribution and use in source and binary forms, with or without
! modification, are permitted provided that the following conditions
! are met:
! 1. Redistributions of source code must retain the above copyright
! notice, this list of conditions and the following disclaimer.
! 2. Redistributions in binary form must reproduce the above copyright
! notice, this list of conditions, and the following disclaimer in the
! documentation and/or other materials provided with the distribution.
! 3. The name of the PSBLAS group or the names of its contributors may
! not be used to endorse or promote products derived from this
! software without specific written permission.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
! ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
! TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
! PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
! BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
! CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
! SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
! CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
! ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
!
subroutine psb_cprecseti(p,what,val,info)
use psb_base_mod
use psb_prec_mod, psb_protect_name => psb_cprecseti
implicit none
type(psb_cprec_type), intent(inout) :: p
integer(psb_ipk_) :: what, val
integer(psb_ipk_), intent(out) :: info
character(len=20) :: name='precset'
info = psb_success_
if (.not.allocated(p%prec)) then
info = 1124
call psb_errpush(info,name,a_err="preconditioner")
return
!!$ goto 9999
end if
call p%prec%precset(what,val,info)
return
end subroutine psb_cprecseti
subroutine psb_cprecsetr(p,what,val,info)
use psb_base_mod
use psb_prec_mod, psb_protect_name => psb_cprecsetr
implicit none
type(psb_cprec_type), intent(inout) :: p
integer(psb_ipk_) :: what
real(psb_spk_) :: val
integer(psb_ipk_), intent(out) :: info
character(len=20) :: name='precset'
info = psb_success_
if (.not.allocated(p%prec)) then
info = 1124
call psb_errpush(info,name,a_err="preconditioner")
return
!!$ goto 9999
end if
call p%prec%precset(what,val,info)
return
end subroutine psb_cprecsetr

80
prec/impl/psb_dprecset.f90
Unescape Escape View File

@ -1,80 +0,0 @@
!
! Parallel Sparse BLAS version 3.5
! (C) Copyright 2006-2018
! Salvatore Filippone
! Alfredo Buttari
!
! Redistribution and use in source and binary forms, with or without
! modification, are permitted provided that the following conditions
! are met:
! 1. Redistributions of source code must retain the above copyright
! notice, this list of conditions and the following disclaimer.
! 2. Redistributions in binary form must reproduce the above copyright
! notice, this list of conditions, and the following disclaimer in the
! documentation and/or other materials provided with the distribution.
! 3. The name of the PSBLAS group or the names of its contributors may
! not be used to endorse or promote products derived from this
! software without specific written permission.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
! ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
! TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
! PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
! BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
! CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
! SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
! CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
! ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
!
subroutine psb_dprecseti(p,what,val,info)
use psb_base_mod
use psb_prec_mod, psb_protect_name => psb_dprecseti
implicit none
type(psb_dprec_type), intent(inout) :: p
integer(psb_ipk_) :: what, val
integer(psb_ipk_), intent(out) :: info
character(len=20) :: name='precset'
info = psb_success_
if (.not.allocated(p%prec)) then
info = 1124
call psb_errpush(info,name,a_err="preconditioner")
return
!!$ goto 9999
end if
call p%prec%precset(what,val,info)
return
end subroutine psb_dprecseti
subroutine psb_dprecsetr(p,what,val,info)
use psb_base_mod
use psb_prec_mod, psb_protect_name => psb_dprecsetr
implicit none
type(psb_dprec_type), intent(inout) :: p
integer(psb_ipk_) :: what
real(psb_dpk_) :: val
integer(psb_ipk_), intent(out) :: info
character(len=20) :: name='precset'
info = psb_success_
if (.not.allocated(p%prec)) then
info = 1124
call psb_errpush(info,name,a_err="preconditioner")
return
!!$ goto 9999
end if
call p%prec%precset(what,val,info)
return
end subroutine psb_dprecsetr

79
prec/impl/psb_sprecset.f90
Unescape Escape View File

@ -1,79 +0,0 @@
!
! Parallel Sparse BLAS version 3.5
! (C) Copyright 2006-2018
! Salvatore Filippone
! Alfredo Buttari
!
! Redistribution and use in source and binary forms, with or without
! modification, are permitted provided that the following conditions
! are met:
! 1. Redistributions of source code must retain the above copyright
! notice, this list of conditions and the following disclaimer.
! 2. Redistributions in binary form must reproduce the above copyright
! notice, this list of conditions, and the following disclaimer in the
! documentation and/or other materials provided with the distribution.
! 3. The name of the PSBLAS group or the names of its contributors may
! not be used to endorse or promote products derived from this
! software without specific written permission.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
! ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
! TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
! PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
! BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
! CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
! SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
! CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
! ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
!
subroutine psb_sprecseti(p,what,val,info)
use psb_base_mod
use psb_prec_mod, psb_protect_name => psb_sprecseti
implicit none
type(psb_sprec_type), intent(inout) :: p
integer(psb_ipk_) :: what, val
integer(psb_ipk_), intent(out) :: info
character(len=20) :: name='precset'
info = psb_success_
if (.not.allocated(p%prec)) then
info = 1124
call psb_errpush(info,name,a_err="preconditioner")
return
!!$ goto 9999
end if
call p%prec%precset(what,val,info)
return
end subroutine psb_sprecseti
subroutine psb_sprecsetr(p,what,val,info)
use psb_base_mod
use psb_prec_mod, psb_protect_name => psb_sprecsetr
implicit none
type(psb_sprec_type), intent(inout) :: p
integer(psb_ipk_) :: what
real(psb_spk_) :: val
integer(psb_ipk_), intent(out) :: info
character(len=20) :: name='precset'
info = psb_success_
if (.not.allocated(p%prec)) then
info = 1124
call psb_errpush(info,name,a_err="preconditioner")
return
!!$ goto 9999
end if
call p%prec%precset(what,val,info)
return
end subroutine psb_sprecsetr

80
prec/impl/psb_zprecset.f90
Unescape Escape View File

@ -1,80 +0,0 @@
!
! Parallel Sparse BLAS version 3.5
! (C) Copyright 2006-2018
! Salvatore Filippone
! Alfredo Buttari
!
! Redistribution and use in source and binary forms, with or without
! modification, are permitted provided that the following conditions
! are met:
! 1. Redistributions of source code must retain the above copyright
! notice, this list of conditions and the following disclaimer.
! 2. Redistributions in binary form must reproduce the above copyright
! notice, this list of conditions, and the following disclaimer in the
! documentation and/or other materials provided with the distribution.
! 3. The name of the PSBLAS group or the names of its contributors may
! not be used to endorse or promote products derived from this
! software without specific written permission.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
! ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
! TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
! PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
! BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
! CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
! SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
! INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
! CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
! ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
! POSSIBILITY OF SUCH DAMAGE.
!
!
subroutine psb_zprecseti(p,what,val,info)
use psb_base_mod
use psb_prec_mod, psb_protect_name => psb_zprecseti
implicit none
type(psb_zprec_type), intent(inout) :: p
integer(psb_ipk_) :: what, val
integer(psb_ipk_), intent(out) :: info
character(len=20) :: name='precset'
info = psb_success_
if (.not.allocated(p%prec)) then
info = 1124
call psb_errpush(info,name,a_err="preconditioner")
return
!!$ goto 9999
end if
call p%prec%precset(what,val,info)
return
end subroutine psb_zprecseti
subroutine psb_zprecsetr(p,what,val,info)
use psb_base_mod
use psb_prec_mod, psb_protect_name => psb_zprecsetr
implicit none
type(psb_zprec_type), intent(inout) :: p
integer(psb_ipk_) :: what
real(psb_dpk_) :: val
integer(psb_ipk_), intent(out) :: info
character(len=20) :: name='precset'
info = psb_success_
if (.not.allocated(p%prec)) then
info = 1124
call psb_errpush(info,name,a_err="preconditioner")
return
!!$ goto 9999
end if
call p%prec%precset(what,val,info)
return
end subroutine psb_zprecsetr
Write Preview
Loading…
Cancel
Save