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

Added clear_data and clone_settings to GS solver

Browse Source
pizdaint-runs
Salvatore Filippone 5 years ago
parent 07a2a97b13
commit 9ba74e89cd
17 changed files with 664 additions and 4 deletions
Show Stats Download Patch File Download Diff File

8
mlprec/impl/solver/Makefile
Unescape Escape View File

@ -34,6 +34,8 @@ mld_c_gs_solver_cnv.o \
mld_c_gs_solver_dmp.o \ mld_c_gs_solver_dmp.o \
mld_c_gs_solver_apply.o \ mld_c_gs_solver_apply.o \
mld_c_gs_solver_apply_vect.o \ mld_c_gs_solver_apply_vect.o \
mld_c_gs_solver_clear_data.o \
mld_c_gs_solver_clone_settings.o \
mld_c_bwgs_solver_bld.o \ mld_c_bwgs_solver_bld.o \
mld_c_bwgs_solver_apply.o \ mld_c_bwgs_solver_apply.o \
mld_c_bwgs_solver_apply_vect.o \ mld_c_bwgs_solver_apply_vect.o \
@ -82,6 +84,8 @@ mld_d_gs_solver_cnv.o \
mld_d_gs_solver_dmp.o \ mld_d_gs_solver_dmp.o \
mld_d_gs_solver_apply.o \ mld_d_gs_solver_apply.o \
mld_d_gs_solver_apply_vect.o \ mld_d_gs_solver_apply_vect.o \
mld_d_gs_solver_clear_data.o \
mld_d_gs_solver_clone_settings.o \
mld_d_id_solver_apply.o \ mld_d_id_solver_apply.o \
mld_d_id_solver_apply_vect.o \ mld_d_id_solver_apply_vect.o \
mld_d_id_solver_clone.o \ mld_d_id_solver_clone.o \
@ -124,6 +128,8 @@ mld_s_gs_solver_cnv.o \
mld_s_gs_solver_dmp.o \ mld_s_gs_solver_dmp.o \
mld_s_gs_solver_apply.o \ mld_s_gs_solver_apply.o \
mld_s_gs_solver_apply_vect.o \ mld_s_gs_solver_apply_vect.o \
mld_s_gs_solver_clear_data.o \
mld_s_gs_solver_clone_settings.o \
mld_s_bwgs_solver_bld.o \ mld_s_bwgs_solver_bld.o \
mld_s_bwgs_solver_apply.o \ mld_s_bwgs_solver_apply.o \
mld_s_bwgs_solver_apply_vect.o \ mld_s_bwgs_solver_apply_vect.o \
@ -169,6 +175,8 @@ mld_z_gs_solver_cnv.o \
mld_z_gs_solver_dmp.o \ mld_z_gs_solver_dmp.o \
mld_z_gs_solver_apply.o \ mld_z_gs_solver_apply.o \
mld_z_gs_solver_apply_vect.o \ mld_z_gs_solver_apply_vect.o \
mld_z_gs_solver_clear_data.o \
mld_z_gs_solver_clone_settings.o \
mld_z_bwgs_solver_bld.o \ mld_z_bwgs_solver_bld.o \
mld_z_bwgs_solver_apply.o \ mld_z_bwgs_solver_apply.o \
mld_z_bwgs_solver_apply_vect.o \ mld_z_bwgs_solver_apply_vect.o \

2
mlprec/impl/solver/mld_c_base_solver_clone_settings.f90
Unescape Escape View File

@ -45,7 +45,7 @@ subroutine mld_c_base_solver_clone_settings(sv,svout,info)
class(mld_c_base_solver_type), intent(inout) :: svout class(mld_c_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act integer(psb_ipk_) :: err_act
character(len=20) :: name='c_base_solver_clone' character(len=20) :: name='c_base_solver_clone_settings'
call psb_erractionsave(err_act) call psb_erractionsave(err_act)

65
mlprec/impl/solver/mld_c_gs_solver_clear_data.f90
Unescape Escape View File

@ -0,0 +1,65 @@
!
!
! MLD2P4 version 2.2
! MultiLevel Domain Decomposition Parallel Preconditioners Package
! based on PSBLAS (Parallel Sparse BLAS version 3.5)
!
! (C) Copyright 2008-2018
!
! Salvatore Filippone
! Pasqua D'Ambra
! Daniela di Serafino
!
! 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 MLD2P4 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 MLD2P4 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 mld_c_gs_solver_clear_data(sv,info)
use psb_base_mod
use mld_c_gs_solver, mld_protect_name => mld_c_gs_solver_clear_data
Implicit None
! Arguments
class(mld_c_gs_solver_type), intent(inout) :: sv
integer(psb_ipk_), intent(out) :: info
! Local variables
integer(psb_ipk_) :: err_act
info=psb_success_
call psb_erractionsave(err_act)
call sv%l%free()
call sv%u%free()
if (info /= 0) goto 9999
call psb_erractionrestore(err_act)
return
9999 call psb_error_handler(err_act)
return
end subroutine mld_c_gs_solver_clear_data

69
mlprec/impl/solver/mld_c_gs_solver_clone_settings.f90
Unescape Escape View File

@ -0,0 +1,69 @@
!
!
! MLD2P4 version 2.2
! MultiLevel Domain Decomposition Parallel Preconditioners Package
! based on PSBLAS (Parallel Sparse BLAS version 3.5)
!
! (C) Copyright 2008-2018
!
! Salvatore Filippone
! Pasqua D'Ambra
! Daniela di Serafino
!
! 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 MLD2P4 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 MLD2P4 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 mld_c_gs_solver_clone_settings(sv,svout,info)
use psb_base_mod
use mld_c_gs_solver, mld_protect_name => mld_c_gs_solver_clone_settings
Implicit None
! Arguments
class(mld_c_gs_solver_type), intent(inout) :: sv
class(mld_c_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act
character(len=20) :: name='c_gs_solver_clone_settings'
call psb_erractionsave(err_act)
select type(svout)
class is(mld_c_gs_solver_type)
svout%sweeps = sv%sweeps
svout%eps = sv%eps
class default
info = psb_err_internal_error_
call psb_errpush(info,name)
goto 9999
end select
call psb_erractionrestore(err_act)
return
9999 call psb_error_handler(err_act)
return
end subroutine mld_c_gs_solver_clone_settings

2
mlprec/impl/solver/mld_d_base_solver_clone_settings.f90
Unescape Escape View File

@ -45,7 +45,7 @@ subroutine mld_d_base_solver_clone_settings(sv,svout,info)
class(mld_d_base_solver_type), intent(inout) :: svout class(mld_d_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act integer(psb_ipk_) :: err_act
character(len=20) :: name='d_base_solver_clone' character(len=20) :: name='d_base_solver_clone_settings'
call psb_erractionsave(err_act) call psb_erractionsave(err_act)

65
mlprec/impl/solver/mld_d_gs_solver_clear_data.f90
Unescape Escape View File

@ -0,0 +1,65 @@
!
!
! MLD2P4 version 2.2
! MultiLevel Domain Decomposition Parallel Preconditioners Package
! based on PSBLAS (Parallel Sparse BLAS version 3.5)
!
! (C) Copyright 2008-2018
!
! Salvatore Filippone
! Pasqua D'Ambra
! Daniela di Serafino
!
! 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 MLD2P4 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 MLD2P4 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 mld_d_gs_solver_clear_data(sv,info)
use psb_base_mod
use mld_d_gs_solver, mld_protect_name => mld_d_gs_solver_clear_data
Implicit None
! Arguments
class(mld_d_gs_solver_type), intent(inout) :: sv
integer(psb_ipk_), intent(out) :: info
! Local variables
integer(psb_ipk_) :: err_act
info=psb_success_
call psb_erractionsave(err_act)
call sv%l%free()
call sv%u%free()
if (info /= 0) goto 9999
call psb_erractionrestore(err_act)
return
9999 call psb_error_handler(err_act)
return
end subroutine mld_d_gs_solver_clear_data

69
mlprec/impl/solver/mld_d_gs_solver_clone_settings.f90
Unescape Escape View File

@ -0,0 +1,69 @@
!
!
! MLD2P4 version 2.2
! MultiLevel Domain Decomposition Parallel Preconditioners Package
! based on PSBLAS (Parallel Sparse BLAS version 3.5)
!
! (C) Copyright 2008-2018
!
! Salvatore Filippone
! Pasqua D'Ambra
! Daniela di Serafino
!
! 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 MLD2P4 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 MLD2P4 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 mld_d_gs_solver_clone_settings(sv,svout,info)
use psb_base_mod
use mld_d_gs_solver, mld_protect_name => mld_d_gs_solver_clone_settings
Implicit None
! Arguments
class(mld_d_gs_solver_type), intent(inout) :: sv
class(mld_d_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act
character(len=20) :: name='d_gs_solver_clone_settings'
call psb_erractionsave(err_act)
select type(svout)
class is(mld_d_gs_solver_type)
svout%sweeps = sv%sweeps
svout%eps = sv%eps
class default
info = psb_err_internal_error_
call psb_errpush(info,name)
goto 9999
end select
call psb_erractionrestore(err_act)
return
9999 call psb_error_handler(err_act)
return
end subroutine mld_d_gs_solver_clone_settings

2
mlprec/impl/solver/mld_s_base_solver_clone_settings.f90
Unescape Escape View File

@ -45,7 +45,7 @@ subroutine mld_s_base_solver_clone_settings(sv,svout,info)
class(mld_s_base_solver_type), intent(inout) :: svout class(mld_s_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act integer(psb_ipk_) :: err_act
character(len=20) :: name='s_base_solver_clone' character(len=20) :: name='s_base_solver_clone_settings'
call psb_erractionsave(err_act) call psb_erractionsave(err_act)

65
mlprec/impl/solver/mld_s_gs_solver_clear_data.f90
Unescape Escape View File

@ -0,0 +1,65 @@
!
!
! MLD2P4 version 2.2
! MultiLevel Domain Decomposition Parallel Preconditioners Package
! based on PSBLAS (Parallel Sparse BLAS version 3.5)
!
! (C) Copyright 2008-2018
!
! Salvatore Filippone
! Pasqua D'Ambra
! Daniela di Serafino
!
! 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 MLD2P4 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 MLD2P4 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 mld_s_gs_solver_clear_data(sv,info)
use psb_base_mod
use mld_s_gs_solver, mld_protect_name => mld_s_gs_solver_clear_data
Implicit None
! Arguments
class(mld_s_gs_solver_type), intent(inout) :: sv
integer(psb_ipk_), intent(out) :: info
! Local variables
integer(psb_ipk_) :: err_act
info=psb_success_
call psb_erractionsave(err_act)
call sv%l%free()
call sv%u%free()
if (info /= 0) goto 9999
call psb_erractionrestore(err_act)
return
9999 call psb_error_handler(err_act)
return
end subroutine mld_s_gs_solver_clear_data

69
mlprec/impl/solver/mld_s_gs_solver_clone_settings.f90
Unescape Escape View File

@ -0,0 +1,69 @@
!
!
! MLD2P4 version 2.2
! MultiLevel Domain Decomposition Parallel Preconditioners Package
! based on PSBLAS (Parallel Sparse BLAS version 3.5)
!
! (C) Copyright 2008-2018
!
! Salvatore Filippone
! Pasqua D'Ambra
! Daniela di Serafino
!
! 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 MLD2P4 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 MLD2P4 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 mld_s_gs_solver_clone_settings(sv,svout,info)
use psb_base_mod
use mld_s_gs_solver, mld_protect_name => mld_s_gs_solver_clone_settings
Implicit None
! Arguments
class(mld_s_gs_solver_type), intent(inout) :: sv
class(mld_s_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act
character(len=20) :: name='s_gs_solver_clone_settings'
call psb_erractionsave(err_act)
select type(svout)
class is(mld_s_gs_solver_type)
svout%sweeps = sv%sweeps
svout%eps = sv%eps
class default
info = psb_err_internal_error_
call psb_errpush(info,name)
goto 9999
end select
call psb_erractionrestore(err_act)
return
9999 call psb_error_handler(err_act)
return
end subroutine mld_s_gs_solver_clone_settings

2
mlprec/impl/solver/mld_z_base_solver_clone_settings.f90
Unescape Escape View File

@ -45,7 +45,7 @@ subroutine mld_z_base_solver_clone_settings(sv,svout,info)
class(mld_z_base_solver_type), intent(inout) :: svout class(mld_z_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act integer(psb_ipk_) :: err_act
character(len=20) :: name='z_base_solver_clone' character(len=20) :: name='z_base_solver_clone_settings'
call psb_erractionsave(err_act) call psb_erractionsave(err_act)

65
mlprec/impl/solver/mld_z_gs_solver_clear_data.f90
Unescape Escape View File

@ -0,0 +1,65 @@
!
!
! MLD2P4 version 2.2
! MultiLevel Domain Decomposition Parallel Preconditioners Package
! based on PSBLAS (Parallel Sparse BLAS version 3.5)
!
! (C) Copyright 2008-2018
!
! Salvatore Filippone
! Pasqua D'Ambra
! Daniela di Serafino
!
! 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 MLD2P4 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 MLD2P4 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 mld_z_gs_solver_clear_data(sv,info)
use psb_base_mod
use mld_z_gs_solver, mld_protect_name => mld_z_gs_solver_clear_data
Implicit None
! Arguments
class(mld_z_gs_solver_type), intent(inout) :: sv
integer(psb_ipk_), intent(out) :: info
! Local variables
integer(psb_ipk_) :: err_act
info=psb_success_
call psb_erractionsave(err_act)
call sv%l%free()
call sv%u%free()
if (info /= 0) goto 9999
call psb_erractionrestore(err_act)
return
9999 call psb_error_handler(err_act)
return
end subroutine mld_z_gs_solver_clear_data

69
mlprec/impl/solver/mld_z_gs_solver_clone_settings.f90
Unescape Escape View File

@ -0,0 +1,69 @@
!
!
! MLD2P4 version 2.2
! MultiLevel Domain Decomposition Parallel Preconditioners Package
! based on PSBLAS (Parallel Sparse BLAS version 3.5)
!
! (C) Copyright 2008-2018
!
! Salvatore Filippone
! Pasqua D'Ambra
! Daniela di Serafino
!
! 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 MLD2P4 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 MLD2P4 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 mld_z_gs_solver_clone_settings(sv,svout,info)
use psb_base_mod
use mld_z_gs_solver, mld_protect_name => mld_z_gs_solver_clone_settings
Implicit None
! Arguments
class(mld_z_gs_solver_type), intent(inout) :: sv
class(mld_z_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_) :: err_act
character(len=20) :: name='z_gs_solver_clone_settings'
call psb_erractionsave(err_act)
select type(svout)
class is(mld_z_gs_solver_type)
svout%sweeps = sv%sweeps
svout%eps = sv%eps
class default
info = psb_err_internal_error_
call psb_errpush(info,name)
goto 9999
end select
call psb_erractionrestore(err_act)
return
9999 call psb_error_handler(err_act)
return
end subroutine mld_z_gs_solver_clone_settings

29
mlprec/mld_c_gs_solver.f90
Unescape Escape View File

@ -63,6 +63,8 @@ module mld_c_gs_solver
procedure, pass(sv) :: dump => mld_c_gs_solver_dmp procedure, pass(sv) :: dump => mld_c_gs_solver_dmp
procedure, pass(sv) :: check => c_gs_solver_check procedure, pass(sv) :: check => c_gs_solver_check
procedure, pass(sv) :: clone => mld_c_gs_solver_clone procedure, pass(sv) :: clone => mld_c_gs_solver_clone
procedure, pass(sv) :: clone_settings => mld_c_gs_solver_clone_settings
procedure, pass(sv) :: clear_data => mld_c_gs_solver_clear_data
procedure, pass(sv) :: build => mld_c_gs_solver_bld procedure, pass(sv) :: build => mld_c_gs_solver_bld
procedure, pass(sv) :: cnv => mld_c_gs_solver_cnv procedure, pass(sv) :: cnv => mld_c_gs_solver_cnv
procedure, pass(sv) :: apply_v => mld_c_gs_solver_apply_vect procedure, pass(sv) :: apply_v => mld_c_gs_solver_apply_vect
@ -247,6 +249,33 @@ module mld_c_gs_solver
end subroutine mld_c_gs_solver_clone end subroutine mld_c_gs_solver_clone
end interface end interface
interface
subroutine mld_c_gs_solver_clone_settings(sv,svout,info)
import :: psb_desc_type, psb_cspmat_type, psb_c_base_sparse_mat, &
& psb_c_vect_type, psb_c_base_vect_type, psb_spk_, &
& mld_c_base_solver_type, mld_c_gs_solver_type, psb_ipk_
Implicit None
! Arguments
class(mld_c_gs_solver_type), intent(inout) :: sv
class(mld_c_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info
end subroutine mld_c_gs_solver_clone_settings
end interface
interface
subroutine mld_c_gs_solver_clear_data(sv,info)
import :: psb_desc_type, psb_cspmat_type, psb_c_base_sparse_mat, &
& psb_c_vect_type, psb_c_base_vect_type, psb_spk_, &
& mld_c_gs_solver_type, psb_ipk_
Implicit None
! Arguments
class(mld_c_gs_solver_type), intent(inout) :: sv
integer(psb_ipk_), intent(out) :: info
end subroutine mld_c_gs_solver_clear_data
end interface
contains contains
subroutine c_gs_solver_default(sv) subroutine c_gs_solver_default(sv)

29
mlprec/mld_d_gs_solver.f90
Unescape Escape View File

@ -63,6 +63,8 @@ module mld_d_gs_solver
procedure, pass(sv) :: dump => mld_d_gs_solver_dmp procedure, pass(sv) :: dump => mld_d_gs_solver_dmp
procedure, pass(sv) :: check => d_gs_solver_check procedure, pass(sv) :: check => d_gs_solver_check
procedure, pass(sv) :: clone => mld_d_gs_solver_clone procedure, pass(sv) :: clone => mld_d_gs_solver_clone
procedure, pass(sv) :: clone_settings => mld_d_gs_solver_clone_settings
procedure, pass(sv) :: clear_data => mld_d_gs_solver_clear_data
procedure, pass(sv) :: build => mld_d_gs_solver_bld procedure, pass(sv) :: build => mld_d_gs_solver_bld
procedure, pass(sv) :: cnv => mld_d_gs_solver_cnv procedure, pass(sv) :: cnv => mld_d_gs_solver_cnv
procedure, pass(sv) :: apply_v => mld_d_gs_solver_apply_vect procedure, pass(sv) :: apply_v => mld_d_gs_solver_apply_vect
@ -247,6 +249,33 @@ module mld_d_gs_solver
end subroutine mld_d_gs_solver_clone end subroutine mld_d_gs_solver_clone
end interface end interface
interface
subroutine mld_d_gs_solver_clone_settings(sv,svout,info)
import :: psb_desc_type, psb_dspmat_type, psb_d_base_sparse_mat, &
& psb_d_vect_type, psb_d_base_vect_type, psb_dpk_, &
& mld_d_base_solver_type, mld_d_gs_solver_type, psb_ipk_
Implicit None
! Arguments
class(mld_d_gs_solver_type), intent(inout) :: sv
class(mld_d_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info
end subroutine mld_d_gs_solver_clone_settings
end interface
interface
subroutine mld_d_gs_solver_clear_data(sv,info)
import :: psb_desc_type, psb_dspmat_type, psb_d_base_sparse_mat, &
& psb_d_vect_type, psb_d_base_vect_type, psb_dpk_, &
& mld_d_gs_solver_type, psb_ipk_
Implicit None
! Arguments
class(mld_d_gs_solver_type), intent(inout) :: sv
integer(psb_ipk_), intent(out) :: info
end subroutine mld_d_gs_solver_clear_data
end interface
contains contains
subroutine d_gs_solver_default(sv) subroutine d_gs_solver_default(sv)

29
mlprec/mld_s_gs_solver.f90
Unescape Escape View File

@ -63,6 +63,8 @@ module mld_s_gs_solver
procedure, pass(sv) :: dump => mld_s_gs_solver_dmp procedure, pass(sv) :: dump => mld_s_gs_solver_dmp
procedure, pass(sv) :: check => s_gs_solver_check procedure, pass(sv) :: check => s_gs_solver_check
procedure, pass(sv) :: clone => mld_s_gs_solver_clone procedure, pass(sv) :: clone => mld_s_gs_solver_clone
procedure, pass(sv) :: clone_settings => mld_s_gs_solver_clone_settings
procedure, pass(sv) :: clear_data => mld_s_gs_solver_clear_data
procedure, pass(sv) :: build => mld_s_gs_solver_bld procedure, pass(sv) :: build => mld_s_gs_solver_bld
procedure, pass(sv) :: cnv => mld_s_gs_solver_cnv procedure, pass(sv) :: cnv => mld_s_gs_solver_cnv
procedure, pass(sv) :: apply_v => mld_s_gs_solver_apply_vect procedure, pass(sv) :: apply_v => mld_s_gs_solver_apply_vect
@ -247,6 +249,33 @@ module mld_s_gs_solver
end subroutine mld_s_gs_solver_clone end subroutine mld_s_gs_solver_clone
end interface end interface
interface
subroutine mld_s_gs_solver_clone_settings(sv,svout,info)
import :: psb_desc_type, psb_sspmat_type, psb_s_base_sparse_mat, &
& psb_s_vect_type, psb_s_base_vect_type, psb_spk_, &
& mld_s_base_solver_type, mld_s_gs_solver_type, psb_ipk_
Implicit None
! Arguments
class(mld_s_gs_solver_type), intent(inout) :: sv
class(mld_s_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info
end subroutine mld_s_gs_solver_clone_settings
end interface
interface
subroutine mld_s_gs_solver_clear_data(sv,info)
import :: psb_desc_type, psb_sspmat_type, psb_s_base_sparse_mat, &
& psb_s_vect_type, psb_s_base_vect_type, psb_spk_, &
& mld_s_gs_solver_type, psb_ipk_
Implicit None
! Arguments
class(mld_s_gs_solver_type), intent(inout) :: sv
integer(psb_ipk_), intent(out) :: info
end subroutine mld_s_gs_solver_clear_data
end interface
contains contains
subroutine s_gs_solver_default(sv) subroutine s_gs_solver_default(sv)

29
mlprec/mld_z_gs_solver.f90
Unescape Escape View File

@ -63,6 +63,8 @@ module mld_z_gs_solver
procedure, pass(sv) :: dump => mld_z_gs_solver_dmp procedure, pass(sv) :: dump => mld_z_gs_solver_dmp
procedure, pass(sv) :: check => z_gs_solver_check procedure, pass(sv) :: check => z_gs_solver_check
procedure, pass(sv) :: clone => mld_z_gs_solver_clone procedure, pass(sv) :: clone => mld_z_gs_solver_clone
procedure, pass(sv) :: clone_settings => mld_z_gs_solver_clone_settings
procedure, pass(sv) :: clear_data => mld_z_gs_solver_clear_data
procedure, pass(sv) :: build => mld_z_gs_solver_bld procedure, pass(sv) :: build => mld_z_gs_solver_bld
procedure, pass(sv) :: cnv => mld_z_gs_solver_cnv procedure, pass(sv) :: cnv => mld_z_gs_solver_cnv
procedure, pass(sv) :: apply_v => mld_z_gs_solver_apply_vect procedure, pass(sv) :: apply_v => mld_z_gs_solver_apply_vect
@ -247,6 +249,33 @@ module mld_z_gs_solver
end subroutine mld_z_gs_solver_clone end subroutine mld_z_gs_solver_clone
end interface end interface
interface
subroutine mld_z_gs_solver_clone_settings(sv,svout,info)
import :: psb_desc_type, psb_zspmat_type, psb_z_base_sparse_mat, &
& psb_z_vect_type, psb_z_base_vect_type, psb_dpk_, &
& mld_z_base_solver_type, mld_z_gs_solver_type, psb_ipk_
Implicit None
! Arguments
class(mld_z_gs_solver_type), intent(inout) :: sv
class(mld_z_base_solver_type), intent(inout) :: svout
integer(psb_ipk_), intent(out) :: info
end subroutine mld_z_gs_solver_clone_settings
end interface
interface
subroutine mld_z_gs_solver_clear_data(sv,info)
import :: psb_desc_type, psb_zspmat_type, psb_z_base_sparse_mat, &
& psb_z_vect_type, psb_z_base_vect_type, psb_dpk_, &
& mld_z_gs_solver_type, psb_ipk_
Implicit None
! Arguments
class(mld_z_gs_solver_type), intent(inout) :: sv
integer(psb_ipk_), intent(out) :: info
end subroutine mld_z_gs_solver_clear_data
end interface
contains contains
subroutine z_gs_solver_default(sv) subroutine z_gs_solver_default(sv)

Write Preview
Loading…
Cancel
Save