diff --git a/mlprec/impl/smoother/Makefile b/mlprec/impl/smoother/Makefile index f7b32ae9..d223fb92 100644 --- a/mlprec/impl/smoother/Makefile +++ b/mlprec/impl/smoother/Makefile @@ -12,6 +12,8 @@ mld_c_as_smoother_apply_vect.o \ mld_c_as_smoother_bld.o \ mld_c_as_smoother_check.o \ mld_c_as_smoother_clone.o \ +mld_c_as_smoother_clone_settings.o \ +mld_c_as_smoother_clear_data.o \ mld_c_as_smoother_cnv.o \ mld_c_as_smoother_csetc.o \ mld_c_as_smoother_cseti.o \ @@ -41,6 +43,8 @@ mld_c_jac_smoother_bld.o \ mld_c_jac_smoother_descr.o \ mld_c_jac_smoother_dmp.o \ mld_c_jac_smoother_clone.o \ +mld_c_jac_smoother_clone_settings.o \ +mld_c_jac_smoother_clear_data.o \ mld_c_jac_smoother_cnv.o \ mld_c_jac_smoother_csetc.o \ mld_c_jac_smoother_cseti.o \ @@ -53,6 +57,8 @@ mld_d_as_smoother_apply_vect.o \ mld_d_as_smoother_bld.o \ mld_d_as_smoother_check.o \ mld_d_as_smoother_clone.o \ +mld_d_as_smoother_clone_settings.o \ +mld_d_as_smoother_clear_data.o \ mld_d_as_smoother_cnv.o \ mld_d_as_smoother_csetc.o \ mld_d_as_smoother_cseti.o \ @@ -82,6 +88,8 @@ mld_d_jac_smoother_bld.o \ mld_d_jac_smoother_descr.o \ mld_d_jac_smoother_dmp.o \ mld_d_jac_smoother_clone.o \ +mld_d_jac_smoother_clone_settings.o \ +mld_d_jac_smoother_clear_data.o \ mld_d_jac_smoother_cnv.o \ mld_d_jac_smoother_csetc.o \ mld_d_jac_smoother_cseti.o \ @@ -94,6 +102,8 @@ mld_s_as_smoother_apply_vect.o \ mld_s_as_smoother_bld.o \ mld_s_as_smoother_check.o \ mld_s_as_smoother_clone.o \ +mld_s_as_smoother_clone_settings.o \ +mld_s_as_smoother_clear_data.o \ mld_s_as_smoother_cnv.o \ mld_s_as_smoother_csetc.o \ mld_s_as_smoother_cseti.o \ @@ -123,6 +133,8 @@ mld_s_jac_smoother_bld.o \ mld_s_jac_smoother_descr.o \ mld_s_jac_smoother_dmp.o \ mld_s_jac_smoother_clone.o \ +mld_s_jac_smoother_clone_settings.o \ +mld_s_jac_smoother_clear_data.o \ mld_s_jac_smoother_cnv.o \ mld_s_jac_smoother_csetc.o \ mld_s_jac_smoother_cseti.o \ @@ -135,6 +147,8 @@ mld_z_as_smoother_apply_vect.o \ mld_z_as_smoother_bld.o \ mld_z_as_smoother_check.o \ mld_z_as_smoother_clone.o \ +mld_z_as_smoother_clone_settings.o \ +mld_z_as_smoother_clear_data.o \ mld_z_as_smoother_cnv.o \ mld_z_as_smoother_csetc.o \ mld_z_as_smoother_cseti.o \ @@ -164,6 +178,8 @@ mld_z_jac_smoother_bld.o \ mld_z_jac_smoother_descr.o \ mld_z_jac_smoother_dmp.o \ mld_z_jac_smoother_clone.o \ +mld_z_jac_smoother_clone_settings.o \ +mld_z_jac_smoother_clear_data.o \ mld_z_jac_smoother_cnv.o \ mld_z_jac_smoother_csetc.o \ mld_z_jac_smoother_cseti.o \ diff --git a/mlprec/impl/smoother/mld_c_as_smoother_clear_data.f90 b/mlprec/impl/smoother/mld_c_as_smoother_clear_data.f90 new file mode 100644 index 00000000..ec8dfccc --- /dev/null +++ b/mlprec/impl/smoother/mld_c_as_smoother_clear_data.f90 @@ -0,0 +1,70 @@ +! +! +! 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_as_smoother_clear_data(sm,info) + + use psb_base_mod + use mld_c_as_smoother, mld_protect_name => mld_c_as_smoother_clear_data + Implicit None + ! Arguments + class(mld_c_as_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='c_as_smoother_clear_data' + + call psb_erractionsave(err_act) + + info = 0 + call sm%nd%free() + sm%nd_nnz_tot = 0 + call sm%desc_data%free(info) + if ((info==0).and.allocated(sm%sv)) then + call sm%sv%clear_data(info) + end if + if (info /= 0) then + info = psb_err_internal_error_ + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_c_as_smoother_clear_data diff --git a/mlprec/impl/smoother/mld_c_as_smoother_clone_settings.f90 b/mlprec/impl/smoother/mld_c_as_smoother_clone_settings.f90 new file mode 100644 index 00000000..67b3ade9 --- /dev/null +++ b/mlprec/impl/smoother/mld_c_as_smoother_clone_settings.f90 @@ -0,0 +1,95 @@ +! +! +! MLD2P4 version 2.2 +! MultiLevel Domain Decomposition Parallel Preconditioners Package +! asd 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_as_smoother_clone_settings(sm,smout,info) + + use psb_base_mod + use mld_c_as_smoother, mld_protect_name => mld_c_as_smoother_clone_settings + Implicit None + ! Arguments + class(mld_c_as_smoother_type), intent(inout) :: sm + class(mld_c_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='c_as_smoother_clone_settings' + + call psb_erractionsave(err_act) + + info = psb_success_ + + select type(smout) + class is(mld_c_as_smoother_type) + smout%novr = sm%novr + smout%restr = sm%restr + smout%prol = sm%prol + if (allocated(smout%sv)) then + if (.not.same_type_as(sm%sv,smout%sv)) then + call smout%sv%free(info) + if (info == 0) deallocate(smout%sv,stat=info) + end if + end if + if (info /= 0) then + info = psb_err_internal_error_ + else + if (allocated(smout%sv)) then + if (same_type_as(sm%sv,smout%sv)) then + call sm%sv%clone_settings(smout%sv,info) + else + info = psb_err_internal_error_ + end if + else + allocate(smout%sv,mold=sm%sv,stat=info) + if (info == 0) call sm%sv%clone_settings(smout%sv,info) + if (info /= 0) info = psb_err_internal_error_ + end if + end if + class default + info = psb_err_internal_error_ + end select + + if (info /= 0) then + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_c_as_smoother_clone_settings diff --git a/mlprec/impl/smoother/mld_c_jac_smoother_clear_data.f90 b/mlprec/impl/smoother/mld_c_jac_smoother_clear_data.f90 new file mode 100644 index 00000000..0dbf3d9e --- /dev/null +++ b/mlprec/impl/smoother/mld_c_jac_smoother_clear_data.f90 @@ -0,0 +1,70 @@ +! +! +! 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_jac_smoother_clear_data(sm,info) + + use psb_base_mod + use mld_c_jac_smoother, mld_protect_name => mld_c_jac_smoother_clear_data + Implicit None + ! Arguments + class(mld_c_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='c_jac_smoother_clear_data' + + call psb_erractionsave(err_act) + + info = 0 + call sm%nd%free() + sm%nd_nnz_tot = 0 + sm%pa => null() + if ((info==0).and.allocated(sm%sv)) then + call sm%sv%clear_data(info) + end if + if (info /= 0) then + info = psb_err_internal_error_ + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_c_jac_smoother_clear_data diff --git a/mlprec/impl/smoother/mld_c_jac_smoother_clone_settings.f90 b/mlprec/impl/smoother/mld_c_jac_smoother_clone_settings.f90 new file mode 100644 index 00000000..fb8c4567 --- /dev/null +++ b/mlprec/impl/smoother/mld_c_jac_smoother_clone_settings.f90 @@ -0,0 +1,101 @@ +! +! +! MLD2P4 version 2.2 +! MultiLevel Domain Decomposition Parallel Preconditioners Package +! asd 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_jac_smoother_clone_settings(sm,smout,info) + + use psb_base_mod + use mld_c_jac_smoother, mld_protect_name => mld_c_jac_smoother_clone_settings + Implicit None + ! Arguments + class(mld_c_jac_smoother_type), intent(inout) :: sm + class(mld_c_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='c_jac_smoother_clone_settings' + + call psb_erractionsave(err_act) + + info = psb_success_ + + select type(smout) + class is(mld_c_jac_smoother_type) + + smout%pa => null() + smout%nd_nnz_tot = 0 + smout%checkres = sm%checkres + smout%printres = sm%printres + smout%checkiter = sm%checkiter + smout%printiter = sm%printiter + smout%tol = sm%tol + + if (allocated(smout%sv)) then + if (.not.same_type_as(sm%sv,smout%sv)) then + call smout%sv%free(info) + if (info == 0) deallocate(smout%sv,stat=info) + end if + end if + if (info /= 0) then + info = psb_err_internal_error_ + else + if (allocated(smout%sv)) then + if (same_type_as(sm%sv,smout%sv)) then + call sm%sv%clone_settings(smout%sv,info) + else + info = psb_err_internal_error_ + end if + else + allocate(smout%sv,mold=sm%sv,stat=info) + if (info == 0) call sm%sv%clone_settings(smout%sv,info) + if (info /= 0) info = psb_err_internal_error_ + end if + end if + class default + info = psb_err_internal_error_ + end select + + if (info /= 0) then + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_c_jac_smoother_clone_settings diff --git a/mlprec/impl/smoother/mld_d_as_smoother_clear_data.f90 b/mlprec/impl/smoother/mld_d_as_smoother_clear_data.f90 new file mode 100644 index 00000000..898e6f8f --- /dev/null +++ b/mlprec/impl/smoother/mld_d_as_smoother_clear_data.f90 @@ -0,0 +1,70 @@ +! +! +! 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_as_smoother_clear_data(sm,info) + + use psb_base_mod + use mld_d_as_smoother, mld_protect_name => mld_d_as_smoother_clear_data + Implicit None + ! Arguments + class(mld_d_as_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='d_as_smoother_clear_data' + + call psb_erractionsave(err_act) + + info = 0 + call sm%nd%free() + sm%nd_nnz_tot = 0 + call sm%desc_data%free(info) + if ((info==0).and.allocated(sm%sv)) then + call sm%sv%clear_data(info) + end if + if (info /= 0) then + info = psb_err_internal_error_ + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_d_as_smoother_clear_data diff --git a/mlprec/impl/smoother/mld_d_as_smoother_clone_settings.f90 b/mlprec/impl/smoother/mld_d_as_smoother_clone_settings.f90 new file mode 100644 index 00000000..1f713a38 --- /dev/null +++ b/mlprec/impl/smoother/mld_d_as_smoother_clone_settings.f90 @@ -0,0 +1,95 @@ +! +! +! MLD2P4 version 2.2 +! MultiLevel Domain Decomposition Parallel Preconditioners Package +! asd 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_as_smoother_clone_settings(sm,smout,info) + + use psb_base_mod + use mld_d_as_smoother, mld_protect_name => mld_d_as_smoother_clone_settings + Implicit None + ! Arguments + class(mld_d_as_smoother_type), intent(inout) :: sm + class(mld_d_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='d_as_smoother_clone_settings' + + call psb_erractionsave(err_act) + + info = psb_success_ + + select type(smout) + class is(mld_d_as_smoother_type) + smout%novr = sm%novr + smout%restr = sm%restr + smout%prol = sm%prol + if (allocated(smout%sv)) then + if (.not.same_type_as(sm%sv,smout%sv)) then + call smout%sv%free(info) + if (info == 0) deallocate(smout%sv,stat=info) + end if + end if + if (info /= 0) then + info = psb_err_internal_error_ + else + if (allocated(smout%sv)) then + if (same_type_as(sm%sv,smout%sv)) then + call sm%sv%clone_settings(smout%sv,info) + else + info = psb_err_internal_error_ + end if + else + allocate(smout%sv,mold=sm%sv,stat=info) + if (info == 0) call sm%sv%clone_settings(smout%sv,info) + if (info /= 0) info = psb_err_internal_error_ + end if + end if + class default + info = psb_err_internal_error_ + end select + + if (info /= 0) then + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_d_as_smoother_clone_settings diff --git a/mlprec/impl/smoother/mld_d_jac_smoother_clear_data.f90 b/mlprec/impl/smoother/mld_d_jac_smoother_clear_data.f90 new file mode 100644 index 00000000..83e0dcc5 --- /dev/null +++ b/mlprec/impl/smoother/mld_d_jac_smoother_clear_data.f90 @@ -0,0 +1,70 @@ +! +! +! 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_jac_smoother_clear_data(sm,info) + + use psb_base_mod + use mld_d_jac_smoother, mld_protect_name => mld_d_jac_smoother_clear_data + Implicit None + ! Arguments + class(mld_d_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='d_jac_smoother_clear_data' + + call psb_erractionsave(err_act) + + info = 0 + call sm%nd%free() + sm%nd_nnz_tot = 0 + sm%pa => null() + if ((info==0).and.allocated(sm%sv)) then + call sm%sv%clear_data(info) + end if + if (info /= 0) then + info = psb_err_internal_error_ + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_d_jac_smoother_clear_data diff --git a/mlprec/impl/smoother/mld_d_jac_smoother_clone_settings.f90 b/mlprec/impl/smoother/mld_d_jac_smoother_clone_settings.f90 new file mode 100644 index 00000000..5591b832 --- /dev/null +++ b/mlprec/impl/smoother/mld_d_jac_smoother_clone_settings.f90 @@ -0,0 +1,101 @@ +! +! +! MLD2P4 version 2.2 +! MultiLevel Domain Decomposition Parallel Preconditioners Package +! asd 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_jac_smoother_clone_settings(sm,smout,info) + + use psb_base_mod + use mld_d_jac_smoother, mld_protect_name => mld_d_jac_smoother_clone_settings + Implicit None + ! Arguments + class(mld_d_jac_smoother_type), intent(inout) :: sm + class(mld_d_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='d_jac_smoother_clone_settings' + + call psb_erractionsave(err_act) + + info = psb_success_ + + select type(smout) + class is(mld_d_jac_smoother_type) + + smout%pa => null() + smout%nd_nnz_tot = 0 + smout%checkres = sm%checkres + smout%printres = sm%printres + smout%checkiter = sm%checkiter + smout%printiter = sm%printiter + smout%tol = sm%tol + + if (allocated(smout%sv)) then + if (.not.same_type_as(sm%sv,smout%sv)) then + call smout%sv%free(info) + if (info == 0) deallocate(smout%sv,stat=info) + end if + end if + if (info /= 0) then + info = psb_err_internal_error_ + else + if (allocated(smout%sv)) then + if (same_type_as(sm%sv,smout%sv)) then + call sm%sv%clone_settings(smout%sv,info) + else + info = psb_err_internal_error_ + end if + else + allocate(smout%sv,mold=sm%sv,stat=info) + if (info == 0) call sm%sv%clone_settings(smout%sv,info) + if (info /= 0) info = psb_err_internal_error_ + end if + end if + class default + info = psb_err_internal_error_ + end select + + if (info /= 0) then + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_d_jac_smoother_clone_settings diff --git a/mlprec/impl/smoother/mld_s_as_smoother_clear_data.f90 b/mlprec/impl/smoother/mld_s_as_smoother_clear_data.f90 new file mode 100644 index 00000000..fa16c549 --- /dev/null +++ b/mlprec/impl/smoother/mld_s_as_smoother_clear_data.f90 @@ -0,0 +1,70 @@ +! +! +! 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_as_smoother_clear_data(sm,info) + + use psb_base_mod + use mld_s_as_smoother, mld_protect_name => mld_s_as_smoother_clear_data + Implicit None + ! Arguments + class(mld_s_as_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='s_as_smoother_clear_data' + + call psb_erractionsave(err_act) + + info = 0 + call sm%nd%free() + sm%nd_nnz_tot = 0 + call sm%desc_data%free(info) + if ((info==0).and.allocated(sm%sv)) then + call sm%sv%clear_data(info) + end if + if (info /= 0) then + info = psb_err_internal_error_ + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_s_as_smoother_clear_data diff --git a/mlprec/impl/smoother/mld_s_as_smoother_clone_settings.f90 b/mlprec/impl/smoother/mld_s_as_smoother_clone_settings.f90 new file mode 100644 index 00000000..e65ca290 --- /dev/null +++ b/mlprec/impl/smoother/mld_s_as_smoother_clone_settings.f90 @@ -0,0 +1,95 @@ +! +! +! MLD2P4 version 2.2 +! MultiLevel Domain Decomposition Parallel Preconditioners Package +! asd 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_as_smoother_clone_settings(sm,smout,info) + + use psb_base_mod + use mld_s_as_smoother, mld_protect_name => mld_s_as_smoother_clone_settings + Implicit None + ! Arguments + class(mld_s_as_smoother_type), intent(inout) :: sm + class(mld_s_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='s_as_smoother_clone_settings' + + call psb_erractionsave(err_act) + + info = psb_success_ + + select type(smout) + class is(mld_s_as_smoother_type) + smout%novr = sm%novr + smout%restr = sm%restr + smout%prol = sm%prol + if (allocated(smout%sv)) then + if (.not.same_type_as(sm%sv,smout%sv)) then + call smout%sv%free(info) + if (info == 0) deallocate(smout%sv,stat=info) + end if + end if + if (info /= 0) then + info = psb_err_internal_error_ + else + if (allocated(smout%sv)) then + if (same_type_as(sm%sv,smout%sv)) then + call sm%sv%clone_settings(smout%sv,info) + else + info = psb_err_internal_error_ + end if + else + allocate(smout%sv,mold=sm%sv,stat=info) + if (info == 0) call sm%sv%clone_settings(smout%sv,info) + if (info /= 0) info = psb_err_internal_error_ + end if + end if + class default + info = psb_err_internal_error_ + end select + + if (info /= 0) then + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_s_as_smoother_clone_settings diff --git a/mlprec/impl/smoother/mld_s_jac_smoother_clear_data.f90 b/mlprec/impl/smoother/mld_s_jac_smoother_clear_data.f90 new file mode 100644 index 00000000..3838ef64 --- /dev/null +++ b/mlprec/impl/smoother/mld_s_jac_smoother_clear_data.f90 @@ -0,0 +1,70 @@ +! +! +! 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_jac_smoother_clear_data(sm,info) + + use psb_base_mod + use mld_s_jac_smoother, mld_protect_name => mld_s_jac_smoother_clear_data + Implicit None + ! Arguments + class(mld_s_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='s_jac_smoother_clear_data' + + call psb_erractionsave(err_act) + + info = 0 + call sm%nd%free() + sm%nd_nnz_tot = 0 + sm%pa => null() + if ((info==0).and.allocated(sm%sv)) then + call sm%sv%clear_data(info) + end if + if (info /= 0) then + info = psb_err_internal_error_ + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_s_jac_smoother_clear_data diff --git a/mlprec/impl/smoother/mld_s_jac_smoother_clone_settings.f90 b/mlprec/impl/smoother/mld_s_jac_smoother_clone_settings.f90 new file mode 100644 index 00000000..c16d6d6d --- /dev/null +++ b/mlprec/impl/smoother/mld_s_jac_smoother_clone_settings.f90 @@ -0,0 +1,101 @@ +! +! +! MLD2P4 version 2.2 +! MultiLevel Domain Decomposition Parallel Preconditioners Package +! asd 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_jac_smoother_clone_settings(sm,smout,info) + + use psb_base_mod + use mld_s_jac_smoother, mld_protect_name => mld_s_jac_smoother_clone_settings + Implicit None + ! Arguments + class(mld_s_jac_smoother_type), intent(inout) :: sm + class(mld_s_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='s_jac_smoother_clone_settings' + + call psb_erractionsave(err_act) + + info = psb_success_ + + select type(smout) + class is(mld_s_jac_smoother_type) + + smout%pa => null() + smout%nd_nnz_tot = 0 + smout%checkres = sm%checkres + smout%printres = sm%printres + smout%checkiter = sm%checkiter + smout%printiter = sm%printiter + smout%tol = sm%tol + + if (allocated(smout%sv)) then + if (.not.same_type_as(sm%sv,smout%sv)) then + call smout%sv%free(info) + if (info == 0) deallocate(smout%sv,stat=info) + end if + end if + if (info /= 0) then + info = psb_err_internal_error_ + else + if (allocated(smout%sv)) then + if (same_type_as(sm%sv,smout%sv)) then + call sm%sv%clone_settings(smout%sv,info) + else + info = psb_err_internal_error_ + end if + else + allocate(smout%sv,mold=sm%sv,stat=info) + if (info == 0) call sm%sv%clone_settings(smout%sv,info) + if (info /= 0) info = psb_err_internal_error_ + end if + end if + class default + info = psb_err_internal_error_ + end select + + if (info /= 0) then + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_s_jac_smoother_clone_settings diff --git a/mlprec/impl/smoother/mld_z_as_smoother_clear_data.f90 b/mlprec/impl/smoother/mld_z_as_smoother_clear_data.f90 new file mode 100644 index 00000000..f5055e9b --- /dev/null +++ b/mlprec/impl/smoother/mld_z_as_smoother_clear_data.f90 @@ -0,0 +1,70 @@ +! +! +! 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_as_smoother_clear_data(sm,info) + + use psb_base_mod + use mld_z_as_smoother, mld_protect_name => mld_z_as_smoother_clear_data + Implicit None + ! Arguments + class(mld_z_as_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='z_as_smoother_clear_data' + + call psb_erractionsave(err_act) + + info = 0 + call sm%nd%free() + sm%nd_nnz_tot = 0 + call sm%desc_data%free(info) + if ((info==0).and.allocated(sm%sv)) then + call sm%sv%clear_data(info) + end if + if (info /= 0) then + info = psb_err_internal_error_ + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_z_as_smoother_clear_data diff --git a/mlprec/impl/smoother/mld_z_as_smoother_clone_settings.f90 b/mlprec/impl/smoother/mld_z_as_smoother_clone_settings.f90 new file mode 100644 index 00000000..9a239b17 --- /dev/null +++ b/mlprec/impl/smoother/mld_z_as_smoother_clone_settings.f90 @@ -0,0 +1,95 @@ +! +! +! MLD2P4 version 2.2 +! MultiLevel Domain Decomposition Parallel Preconditioners Package +! asd 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_as_smoother_clone_settings(sm,smout,info) + + use psb_base_mod + use mld_z_as_smoother, mld_protect_name => mld_z_as_smoother_clone_settings + Implicit None + ! Arguments + class(mld_z_as_smoother_type), intent(inout) :: sm + class(mld_z_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='z_as_smoother_clone_settings' + + call psb_erractionsave(err_act) + + info = psb_success_ + + select type(smout) + class is(mld_z_as_smoother_type) + smout%novr = sm%novr + smout%restr = sm%restr + smout%prol = sm%prol + if (allocated(smout%sv)) then + if (.not.same_type_as(sm%sv,smout%sv)) then + call smout%sv%free(info) + if (info == 0) deallocate(smout%sv,stat=info) + end if + end if + if (info /= 0) then + info = psb_err_internal_error_ + else + if (allocated(smout%sv)) then + if (same_type_as(sm%sv,smout%sv)) then + call sm%sv%clone_settings(smout%sv,info) + else + info = psb_err_internal_error_ + end if + else + allocate(smout%sv,mold=sm%sv,stat=info) + if (info == 0) call sm%sv%clone_settings(smout%sv,info) + if (info /= 0) info = psb_err_internal_error_ + end if + end if + class default + info = psb_err_internal_error_ + end select + + if (info /= 0) then + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_z_as_smoother_clone_settings diff --git a/mlprec/impl/smoother/mld_z_jac_smoother_clear_data.f90 b/mlprec/impl/smoother/mld_z_jac_smoother_clear_data.f90 new file mode 100644 index 00000000..9ce9e070 --- /dev/null +++ b/mlprec/impl/smoother/mld_z_jac_smoother_clear_data.f90 @@ -0,0 +1,70 @@ +! +! +! 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_jac_smoother_clear_data(sm,info) + + use psb_base_mod + use mld_z_jac_smoother, mld_protect_name => mld_z_jac_smoother_clear_data + Implicit None + ! Arguments + class(mld_z_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='z_jac_smoother_clear_data' + + call psb_erractionsave(err_act) + + info = 0 + call sm%nd%free() + sm%nd_nnz_tot = 0 + sm%pa => null() + if ((info==0).and.allocated(sm%sv)) then + call sm%sv%clear_data(info) + end if + if (info /= 0) then + info = psb_err_internal_error_ + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_z_jac_smoother_clear_data diff --git a/mlprec/impl/smoother/mld_z_jac_smoother_clone_settings.f90 b/mlprec/impl/smoother/mld_z_jac_smoother_clone_settings.f90 new file mode 100644 index 00000000..416676a6 --- /dev/null +++ b/mlprec/impl/smoother/mld_z_jac_smoother_clone_settings.f90 @@ -0,0 +1,101 @@ +! +! +! MLD2P4 version 2.2 +! MultiLevel Domain Decomposition Parallel Preconditioners Package +! asd 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_jac_smoother_clone_settings(sm,smout,info) + + use psb_base_mod + use mld_z_jac_smoother, mld_protect_name => mld_z_jac_smoother_clone_settings + Implicit None + ! Arguments + class(mld_z_jac_smoother_type), intent(inout) :: sm + class(mld_z_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + integer(psb_ipk_) :: err_act + character(len=20) :: name='z_jac_smoother_clone_settings' + + call psb_erractionsave(err_act) + + info = psb_success_ + + select type(smout) + class is(mld_z_jac_smoother_type) + + smout%pa => null() + smout%nd_nnz_tot = 0 + smout%checkres = sm%checkres + smout%printres = sm%printres + smout%checkiter = sm%checkiter + smout%printiter = sm%printiter + smout%tol = sm%tol + + if (allocated(smout%sv)) then + if (.not.same_type_as(sm%sv,smout%sv)) then + call smout%sv%free(info) + if (info == 0) deallocate(smout%sv,stat=info) + end if + end if + if (info /= 0) then + info = psb_err_internal_error_ + else + if (allocated(smout%sv)) then + if (same_type_as(sm%sv,smout%sv)) then + call sm%sv%clone_settings(smout%sv,info) + else + info = psb_err_internal_error_ + end if + else + allocate(smout%sv,mold=sm%sv,stat=info) + if (info == 0) call sm%sv%clone_settings(smout%sv,info) + if (info /= 0) info = psb_err_internal_error_ + end if + end if + class default + info = psb_err_internal_error_ + end select + + if (info /= 0) then + call psb_errpush(info,name) + goto 9999 + end if + + call psb_erractionrestore(err_act) + return + +9999 call psb_error_handler(err_act) + + return +end subroutine mld_z_jac_smoother_clone_settings diff --git a/mlprec/mld_c_as_smoother.f90 b/mlprec/mld_c_as_smoother.f90 index c01250ae..dd7bec53 100644 --- a/mlprec/mld_c_as_smoother.f90 +++ b/mlprec/mld_c_as_smoother.f90 @@ -70,13 +70,15 @@ module mld_c_as_smoother integer(psb_ipk_) :: novr, restr, prol integer(psb_lpk_) :: nd_nnz_tot contains + procedure, pass(sm) :: apply_v => mld_c_as_smoother_apply_vect + procedure, pass(sm) :: apply_a => mld_c_as_smoother_apply procedure, pass(sm) :: check => mld_c_as_smoother_check procedure, pass(sm) :: dump => mld_c_as_smoother_dmp procedure, pass(sm) :: build => mld_c_as_smoother_bld procedure, pass(sm) :: cnv => mld_c_as_smoother_cnv procedure, pass(sm) :: clone => mld_c_as_smoother_clone - procedure, pass(sm) :: apply_v => mld_c_as_smoother_apply_vect - procedure, pass(sm) :: apply_a => mld_c_as_smoother_apply + procedure, pass(sm) :: clone_settings => mld_c_as_smoother_clone_settings + procedure, pass(sm) :: clear_data => mld_c_as_smoother_clear_data procedure, pass(sm) :: restr_a => mld_c_as_smoother_restr_a procedure, pass(sm) :: prol_a => mld_c_as_smoother_prol_a procedure, pass(sm) :: restr_v => mld_c_as_smoother_restr_v @@ -312,7 +314,36 @@ module mld_c_as_smoother integer(psb_ipk_), intent(out) :: info end subroutine mld_c_as_smoother_clone end interface + + + interface + subroutine mld_c_as_smoother_clone_settings(sm,smout,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_smoother_type, mld_c_as_smoother_type, psb_ipk_ + Implicit None + + ! Arguments + class(mld_c_as_smoother_type), intent(inout) :: sm + class(mld_c_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_c_as_smoother_clone_settings + end interface + + interface + subroutine mld_c_as_smoother_clear_data(sm,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_as_smoother_type, psb_ipk_ + Implicit None + + ! Arguments + class(mld_c_as_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_c_as_smoother_clear_data + end interface + contains function c_as_smoother_sizeof(sm) result(val) diff --git a/mlprec/mld_c_jac_smoother.f90 b/mlprec/mld_c_jac_smoother.f90 index 39b63d07..df30aed0 100644 --- a/mlprec/mld_c_jac_smoother.f90 +++ b/mlprec/mld_c_jac_smoother.f90 @@ -67,12 +67,14 @@ module mld_c_jac_smoother integer(psb_ipk_) :: printiter real(psb_dpk_) :: tol contains + procedure, pass(sm) :: apply_v => mld_c_jac_smoother_apply_vect + procedure, pass(sm) :: apply_a => mld_c_jac_smoother_apply procedure, pass(sm) :: dump => mld_c_jac_smoother_dmp procedure, pass(sm) :: build => mld_c_jac_smoother_bld procedure, pass(sm) :: cnv => mld_c_jac_smoother_cnv procedure, pass(sm) :: clone => mld_c_jac_smoother_clone - procedure, pass(sm) :: apply_v => mld_c_jac_smoother_apply_vect - procedure, pass(sm) :: apply_a => mld_c_jac_smoother_apply + procedure, pass(sm) :: clone_settings => mld_c_jac_smoother_clone_settings + procedure, pass(sm) :: clear_data => mld_c_jac_smoother_clear_data procedure, pass(sm) :: free => c_jac_smoother_free procedure, pass(sm) :: cseti => mld_c_jac_smoother_cseti procedure, pass(sm) :: csetc => mld_c_jac_smoother_csetc @@ -197,6 +199,25 @@ module mld_c_jac_smoother end subroutine mld_c_jac_smoother_clone end interface + interface + subroutine mld_c_jac_smoother_clone_settings(sm,smout,info) + import :: mld_c_jac_smoother_type, psb_spk_, & + & mld_c_base_smoother_type, psb_ipk_ + class(mld_c_jac_smoother_type), intent(inout) :: sm + class(mld_c_base_smoother_type), allocatable, intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_c_jac_smoother_clone_settings + end interface + + interface + subroutine mld_c_jac_smoother_clear_data(sm,info) + import :: mld_c_jac_smoother_type, psb_spk_, & + & mld_c_base_smoother_type, psb_ipk_ + class(mld_c_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_c_jac_smoother_clear_data + end interface + interface subroutine mld_c_jac_smoother_descr(sm,info,iout,coarse) import :: mld_c_jac_smoother_type, psb_ipk_ @@ -272,6 +293,25 @@ module mld_c_jac_smoother end subroutine mld_c_l1_jac_smoother_clone end interface + interface + subroutine mld_c_l1_jac_smoother_clone_settings(sm,smout,info) + import :: mld_c_l1_jac_smoother_type, & + & mld_c_base_smoother_type, psb_ipk_ + class(mld_c_l1_jac_smoother_type), intent(inout) :: sm + class(mld_c_base_smoother_type), allocatable, intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_c_l1_jac_smoother_clone_settings + end interface + + interface + subroutine mld_c_l1_jac_smoother_clear_data(sm,info) + import :: mld_c_l1_jac_smoother_type, & + & mld_c_base_smoother_type, psb_ipk_ + class(mld_c_l1_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_c_l1_jac_smoother_clear_data + end interface + interface subroutine mld_c_l1_jac_smoother_descr(sm,info,iout,coarse) import :: mld_c_l1_jac_smoother_type, psb_ipk_ diff --git a/mlprec/mld_d_as_smoother.f90 b/mlprec/mld_d_as_smoother.f90 index ff249ff4..c544b3c5 100644 --- a/mlprec/mld_d_as_smoother.f90 +++ b/mlprec/mld_d_as_smoother.f90 @@ -70,13 +70,15 @@ module mld_d_as_smoother integer(psb_ipk_) :: novr, restr, prol integer(psb_lpk_) :: nd_nnz_tot contains + procedure, pass(sm) :: apply_v => mld_d_as_smoother_apply_vect + procedure, pass(sm) :: apply_a => mld_d_as_smoother_apply procedure, pass(sm) :: check => mld_d_as_smoother_check procedure, pass(sm) :: dump => mld_d_as_smoother_dmp procedure, pass(sm) :: build => mld_d_as_smoother_bld procedure, pass(sm) :: cnv => mld_d_as_smoother_cnv procedure, pass(sm) :: clone => mld_d_as_smoother_clone - procedure, pass(sm) :: apply_v => mld_d_as_smoother_apply_vect - procedure, pass(sm) :: apply_a => mld_d_as_smoother_apply + procedure, pass(sm) :: clone_settings => mld_d_as_smoother_clone_settings + procedure, pass(sm) :: clear_data => mld_d_as_smoother_clear_data procedure, pass(sm) :: restr_a => mld_d_as_smoother_restr_a procedure, pass(sm) :: prol_a => mld_d_as_smoother_prol_a procedure, pass(sm) :: restr_v => mld_d_as_smoother_restr_v @@ -312,7 +314,36 @@ module mld_d_as_smoother integer(psb_ipk_), intent(out) :: info end subroutine mld_d_as_smoother_clone end interface + + + interface + subroutine mld_d_as_smoother_clone_settings(sm,smout,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_smoother_type, mld_d_as_smoother_type, psb_ipk_ + Implicit None + + ! Arguments + class(mld_d_as_smoother_type), intent(inout) :: sm + class(mld_d_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_d_as_smoother_clone_settings + end interface + + interface + subroutine mld_d_as_smoother_clear_data(sm,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_as_smoother_type, psb_ipk_ + Implicit None + + ! Arguments + class(mld_d_as_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_d_as_smoother_clear_data + end interface + contains function d_as_smoother_sizeof(sm) result(val) diff --git a/mlprec/mld_d_jac_smoother.f90 b/mlprec/mld_d_jac_smoother.f90 index 45871afb..a312333f 100644 --- a/mlprec/mld_d_jac_smoother.f90 +++ b/mlprec/mld_d_jac_smoother.f90 @@ -67,12 +67,14 @@ module mld_d_jac_smoother integer(psb_ipk_) :: printiter real(psb_dpk_) :: tol contains + procedure, pass(sm) :: apply_v => mld_d_jac_smoother_apply_vect + procedure, pass(sm) :: apply_a => mld_d_jac_smoother_apply procedure, pass(sm) :: dump => mld_d_jac_smoother_dmp procedure, pass(sm) :: build => mld_d_jac_smoother_bld procedure, pass(sm) :: cnv => mld_d_jac_smoother_cnv procedure, pass(sm) :: clone => mld_d_jac_smoother_clone - procedure, pass(sm) :: apply_v => mld_d_jac_smoother_apply_vect - procedure, pass(sm) :: apply_a => mld_d_jac_smoother_apply + procedure, pass(sm) :: clone_settings => mld_d_jac_smoother_clone_settings + procedure, pass(sm) :: clear_data => mld_d_jac_smoother_clear_data procedure, pass(sm) :: free => d_jac_smoother_free procedure, pass(sm) :: cseti => mld_d_jac_smoother_cseti procedure, pass(sm) :: csetc => mld_d_jac_smoother_csetc @@ -197,6 +199,25 @@ module mld_d_jac_smoother end subroutine mld_d_jac_smoother_clone end interface + interface + subroutine mld_d_jac_smoother_clone_settings(sm,smout,info) + import :: mld_d_jac_smoother_type, psb_dpk_, & + & mld_d_base_smoother_type, psb_ipk_ + class(mld_d_jac_smoother_type), intent(inout) :: sm + class(mld_d_base_smoother_type), allocatable, intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_d_jac_smoother_clone_settings + end interface + + interface + subroutine mld_d_jac_smoother_clear_data(sm,info) + import :: mld_d_jac_smoother_type, psb_dpk_, & + & mld_d_base_smoother_type, psb_ipk_ + class(mld_d_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_d_jac_smoother_clear_data + end interface + interface subroutine mld_d_jac_smoother_descr(sm,info,iout,coarse) import :: mld_d_jac_smoother_type, psb_ipk_ @@ -272,6 +293,25 @@ module mld_d_jac_smoother end subroutine mld_d_l1_jac_smoother_clone end interface + interface + subroutine mld_d_l1_jac_smoother_clone_settings(sm,smout,info) + import :: mld_d_l1_jac_smoother_type, & + & mld_d_base_smoother_type, psb_ipk_ + class(mld_d_l1_jac_smoother_type), intent(inout) :: sm + class(mld_d_base_smoother_type), allocatable, intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_d_l1_jac_smoother_clone_settings + end interface + + interface + subroutine mld_d_l1_jac_smoother_clear_data(sm,info) + import :: mld_d_l1_jac_smoother_type, & + & mld_d_base_smoother_type, psb_ipk_ + class(mld_d_l1_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_d_l1_jac_smoother_clear_data + end interface + interface subroutine mld_d_l1_jac_smoother_descr(sm,info,iout,coarse) import :: mld_d_l1_jac_smoother_type, psb_ipk_ diff --git a/mlprec/mld_s_as_smoother.f90 b/mlprec/mld_s_as_smoother.f90 index 65ebb371..0165f6bb 100644 --- a/mlprec/mld_s_as_smoother.f90 +++ b/mlprec/mld_s_as_smoother.f90 @@ -70,13 +70,15 @@ module mld_s_as_smoother integer(psb_ipk_) :: novr, restr, prol integer(psb_lpk_) :: nd_nnz_tot contains + procedure, pass(sm) :: apply_v => mld_s_as_smoother_apply_vect + procedure, pass(sm) :: apply_a => mld_s_as_smoother_apply procedure, pass(sm) :: check => mld_s_as_smoother_check procedure, pass(sm) :: dump => mld_s_as_smoother_dmp procedure, pass(sm) :: build => mld_s_as_smoother_bld procedure, pass(sm) :: cnv => mld_s_as_smoother_cnv procedure, pass(sm) :: clone => mld_s_as_smoother_clone - procedure, pass(sm) :: apply_v => mld_s_as_smoother_apply_vect - procedure, pass(sm) :: apply_a => mld_s_as_smoother_apply + procedure, pass(sm) :: clone_settings => mld_s_as_smoother_clone_settings + procedure, pass(sm) :: clear_data => mld_s_as_smoother_clear_data procedure, pass(sm) :: restr_a => mld_s_as_smoother_restr_a procedure, pass(sm) :: prol_a => mld_s_as_smoother_prol_a procedure, pass(sm) :: restr_v => mld_s_as_smoother_restr_v @@ -312,7 +314,36 @@ module mld_s_as_smoother integer(psb_ipk_), intent(out) :: info end subroutine mld_s_as_smoother_clone end interface + + + interface + subroutine mld_s_as_smoother_clone_settings(sm,smout,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_smoother_type, mld_s_as_smoother_type, psb_ipk_ + Implicit None + + ! Arguments + class(mld_s_as_smoother_type), intent(inout) :: sm + class(mld_s_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_s_as_smoother_clone_settings + end interface + + interface + subroutine mld_s_as_smoother_clear_data(sm,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_as_smoother_type, psb_ipk_ + Implicit None + + ! Arguments + class(mld_s_as_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_s_as_smoother_clear_data + end interface + contains function s_as_smoother_sizeof(sm) result(val) diff --git a/mlprec/mld_s_jac_smoother.f90 b/mlprec/mld_s_jac_smoother.f90 index bed7c140..278f997d 100644 --- a/mlprec/mld_s_jac_smoother.f90 +++ b/mlprec/mld_s_jac_smoother.f90 @@ -67,12 +67,14 @@ module mld_s_jac_smoother integer(psb_ipk_) :: printiter real(psb_dpk_) :: tol contains + procedure, pass(sm) :: apply_v => mld_s_jac_smoother_apply_vect + procedure, pass(sm) :: apply_a => mld_s_jac_smoother_apply procedure, pass(sm) :: dump => mld_s_jac_smoother_dmp procedure, pass(sm) :: build => mld_s_jac_smoother_bld procedure, pass(sm) :: cnv => mld_s_jac_smoother_cnv procedure, pass(sm) :: clone => mld_s_jac_smoother_clone - procedure, pass(sm) :: apply_v => mld_s_jac_smoother_apply_vect - procedure, pass(sm) :: apply_a => mld_s_jac_smoother_apply + procedure, pass(sm) :: clone_settings => mld_s_jac_smoother_clone_settings + procedure, pass(sm) :: clear_data => mld_s_jac_smoother_clear_data procedure, pass(sm) :: free => s_jac_smoother_free procedure, pass(sm) :: cseti => mld_s_jac_smoother_cseti procedure, pass(sm) :: csetc => mld_s_jac_smoother_csetc @@ -197,6 +199,25 @@ module mld_s_jac_smoother end subroutine mld_s_jac_smoother_clone end interface + interface + subroutine mld_s_jac_smoother_clone_settings(sm,smout,info) + import :: mld_s_jac_smoother_type, psb_spk_, & + & mld_s_base_smoother_type, psb_ipk_ + class(mld_s_jac_smoother_type), intent(inout) :: sm + class(mld_s_base_smoother_type), allocatable, intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_s_jac_smoother_clone_settings + end interface + + interface + subroutine mld_s_jac_smoother_clear_data(sm,info) + import :: mld_s_jac_smoother_type, psb_spk_, & + & mld_s_base_smoother_type, psb_ipk_ + class(mld_s_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_s_jac_smoother_clear_data + end interface + interface subroutine mld_s_jac_smoother_descr(sm,info,iout,coarse) import :: mld_s_jac_smoother_type, psb_ipk_ @@ -272,6 +293,25 @@ module mld_s_jac_smoother end subroutine mld_s_l1_jac_smoother_clone end interface + interface + subroutine mld_s_l1_jac_smoother_clone_settings(sm,smout,info) + import :: mld_s_l1_jac_smoother_type, & + & mld_s_base_smoother_type, psb_ipk_ + class(mld_s_l1_jac_smoother_type), intent(inout) :: sm + class(mld_s_base_smoother_type), allocatable, intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_s_l1_jac_smoother_clone_settings + end interface + + interface + subroutine mld_s_l1_jac_smoother_clear_data(sm,info) + import :: mld_s_l1_jac_smoother_type, & + & mld_s_base_smoother_type, psb_ipk_ + class(mld_s_l1_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_s_l1_jac_smoother_clear_data + end interface + interface subroutine mld_s_l1_jac_smoother_descr(sm,info,iout,coarse) import :: mld_s_l1_jac_smoother_type, psb_ipk_ diff --git a/mlprec/mld_z_as_smoother.f90 b/mlprec/mld_z_as_smoother.f90 index a32b430c..d2e7b8b1 100644 --- a/mlprec/mld_z_as_smoother.f90 +++ b/mlprec/mld_z_as_smoother.f90 @@ -70,13 +70,15 @@ module mld_z_as_smoother integer(psb_ipk_) :: novr, restr, prol integer(psb_lpk_) :: nd_nnz_tot contains + procedure, pass(sm) :: apply_v => mld_z_as_smoother_apply_vect + procedure, pass(sm) :: apply_a => mld_z_as_smoother_apply procedure, pass(sm) :: check => mld_z_as_smoother_check procedure, pass(sm) :: dump => mld_z_as_smoother_dmp procedure, pass(sm) :: build => mld_z_as_smoother_bld procedure, pass(sm) :: cnv => mld_z_as_smoother_cnv procedure, pass(sm) :: clone => mld_z_as_smoother_clone - procedure, pass(sm) :: apply_v => mld_z_as_smoother_apply_vect - procedure, pass(sm) :: apply_a => mld_z_as_smoother_apply + procedure, pass(sm) :: clone_settings => mld_z_as_smoother_clone_settings + procedure, pass(sm) :: clear_data => mld_z_as_smoother_clear_data procedure, pass(sm) :: restr_a => mld_z_as_smoother_restr_a procedure, pass(sm) :: prol_a => mld_z_as_smoother_prol_a procedure, pass(sm) :: restr_v => mld_z_as_smoother_restr_v @@ -312,7 +314,36 @@ module mld_z_as_smoother integer(psb_ipk_), intent(out) :: info end subroutine mld_z_as_smoother_clone end interface + + + interface + subroutine mld_z_as_smoother_clone_settings(sm,smout,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_smoother_type, mld_z_as_smoother_type, psb_ipk_ + Implicit None + + ! Arguments + class(mld_z_as_smoother_type), intent(inout) :: sm + class(mld_z_base_smoother_type), intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_z_as_smoother_clone_settings + end interface + + interface + subroutine mld_z_as_smoother_clear_data(sm,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_as_smoother_type, psb_ipk_ + Implicit None + + ! Arguments + class(mld_z_as_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_z_as_smoother_clear_data + end interface + contains function z_as_smoother_sizeof(sm) result(val) diff --git a/mlprec/mld_z_jac_smoother.f90 b/mlprec/mld_z_jac_smoother.f90 index 950bc077..f9eeda76 100644 --- a/mlprec/mld_z_jac_smoother.f90 +++ b/mlprec/mld_z_jac_smoother.f90 @@ -67,12 +67,14 @@ module mld_z_jac_smoother integer(psb_ipk_) :: printiter real(psb_dpk_) :: tol contains + procedure, pass(sm) :: apply_v => mld_z_jac_smoother_apply_vect + procedure, pass(sm) :: apply_a => mld_z_jac_smoother_apply procedure, pass(sm) :: dump => mld_z_jac_smoother_dmp procedure, pass(sm) :: build => mld_z_jac_smoother_bld procedure, pass(sm) :: cnv => mld_z_jac_smoother_cnv procedure, pass(sm) :: clone => mld_z_jac_smoother_clone - procedure, pass(sm) :: apply_v => mld_z_jac_smoother_apply_vect - procedure, pass(sm) :: apply_a => mld_z_jac_smoother_apply + procedure, pass(sm) :: clone_settings => mld_z_jac_smoother_clone_settings + procedure, pass(sm) :: clear_data => mld_z_jac_smoother_clear_data procedure, pass(sm) :: free => z_jac_smoother_free procedure, pass(sm) :: cseti => mld_z_jac_smoother_cseti procedure, pass(sm) :: csetc => mld_z_jac_smoother_csetc @@ -197,6 +199,25 @@ module mld_z_jac_smoother end subroutine mld_z_jac_smoother_clone end interface + interface + subroutine mld_z_jac_smoother_clone_settings(sm,smout,info) + import :: mld_z_jac_smoother_type, psb_dpk_, & + & mld_z_base_smoother_type, psb_ipk_ + class(mld_z_jac_smoother_type), intent(inout) :: sm + class(mld_z_base_smoother_type), allocatable, intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_z_jac_smoother_clone_settings + end interface + + interface + subroutine mld_z_jac_smoother_clear_data(sm,info) + import :: mld_z_jac_smoother_type, psb_dpk_, & + & mld_z_base_smoother_type, psb_ipk_ + class(mld_z_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_z_jac_smoother_clear_data + end interface + interface subroutine mld_z_jac_smoother_descr(sm,info,iout,coarse) import :: mld_z_jac_smoother_type, psb_ipk_ @@ -272,6 +293,25 @@ module mld_z_jac_smoother end subroutine mld_z_l1_jac_smoother_clone end interface + interface + subroutine mld_z_l1_jac_smoother_clone_settings(sm,smout,info) + import :: mld_z_l1_jac_smoother_type, & + & mld_z_base_smoother_type, psb_ipk_ + class(mld_z_l1_jac_smoother_type), intent(inout) :: sm + class(mld_z_base_smoother_type), allocatable, intent(inout) :: smout + integer(psb_ipk_), intent(out) :: info + end subroutine mld_z_l1_jac_smoother_clone_settings + end interface + + interface + subroutine mld_z_l1_jac_smoother_clear_data(sm,info) + import :: mld_z_l1_jac_smoother_type, & + & mld_z_base_smoother_type, psb_ipk_ + class(mld_z_l1_jac_smoother_type), intent(inout) :: sm + integer(psb_ipk_), intent(out) :: info + end subroutine mld_z_l1_jac_smoother_clear_data + end interface + interface subroutine mld_z_l1_jac_smoother_descr(sm,info,iout,coarse) import :: mld_z_l1_jac_smoother_type, psb_ipk_