diff --git a/base/modules/serial/psb_c_csc_mat_mod.f90 b/base/modules/serial/psb_c_csc_mat_mod.f90 index 300052f3..ad59f0f1 100644 --- a/base/modules/serial/psb_c_csc_mat_mod.f90 +++ b/base/modules/serial/psb_c_csc_mat_mod.f90 @@ -71,7 +71,6 @@ module psb_c_csc_mat_mod procedure, pass(a) :: inner_cssv => psb_c_csc_cssv procedure, pass(a) :: scals => psb_c_csc_scals procedure, pass(a) :: scalv => psb_c_csc_scal -! procedure, pass(a) :: scalpid => psb_c_csc_scalplusidentity procedure, pass(a) :: maxval => psb_c_csc_maxval procedure, pass(a) :: spnm1 => psb_c_csc_csnm1 procedure, pass(a) :: rowsum => psb_c_csc_rowsum @@ -128,7 +127,6 @@ module psb_c_csc_mat_mod procedure, pass(a) :: sizeof => lc_csc_sizeof procedure, pass(a) :: scals => psb_lc_csc_scals procedure, pass(a) :: scalv => psb_lc_csc_scal -! procedure, pass(a) :: scalpid => psb_lc_csc_scalplusidentity procedure, pass(a) :: maxval => psb_lc_csc_maxval procedure, pass(a) :: spnm1 => psb_lc_csc_csnm1 procedure, pass(a) :: rowsum => psb_lc_csc_rowsum @@ -563,18 +561,6 @@ module psb_c_csc_mat_mod end subroutine psb_c_csc_scals end interface - !> \memberof psb_c_csc_sparse_mat - !! \see psb_c_base_mat_mod::psb_c_base_scalplusidentity - ! interface - ! subroutine psb_c_csc_scalplusidentity(d,a,info) - ! import - ! class(psb_c_csc_sparse_mat), intent(inout) :: a - ! complex(psb_spk_), intent(in) :: d - ! integer(psb_ipk_), intent(out) :: info - ! end subroutine psb_c_csc_scalplusidentity - ! end interface - - ! ! lc ! @@ -926,17 +912,6 @@ module psb_c_csc_mat_mod end subroutine psb_lc_csc_scals end interface - !> \memberof psb_lc_csc_sparse_mat - !! \see psb_lc_base_mat_mod::psb_lc_base_scalplusidentity -! interface -! subroutine psb_lc_csc_scalplusidentity(d,a,info) -! import -! class(psb_lc_csc_sparse_mat), intent(inout) :: a -! complex(psb_spk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! end subroutine psb_lc_csc_scalplusidentity -! end interface - contains ! == =================================== diff --git a/base/modules/serial/psb_c_csr_mat_mod.f90 b/base/modules/serial/psb_c_csr_mat_mod.f90 index 4f7d267d..e9004b3c 100644 --- a/base/modules/serial/psb_c_csr_mat_mod.f90 +++ b/base/modules/serial/psb_c_csr_mat_mod.f90 @@ -73,7 +73,6 @@ module psb_c_csr_mat_mod procedure, pass(a) :: inner_cssv => psb_c_csr_cssv procedure, pass(a) :: scals => psb_c_csr_scals procedure, pass(a) :: scalv => psb_c_csr_scal -! procedure, pass(a) :: scalpid => psb_c_csr_scalplusidentity procedure, pass(a) :: maxval => psb_c_csr_maxval procedure, pass(a) :: spnmi => psb_c_csr_csnmi procedure, pass(a) :: rowsum => psb_c_csr_rowsum @@ -580,18 +579,6 @@ module psb_c_csr_mat_mod end subroutine psb_c_csr_scals end interface - !> \memberof psb_c_csr_sparse_mat - !! \see psb_c_base_mat_mod::psb_c_base_scalplusidentity - ! interface - ! subroutine psb_c_csr_scalplusidentity(d,a,info) - ! import - ! class(psb_c_csr_sparse_mat), intent(inout) :: a - ! complex(psb_spk_), intent(in) :: d - ! integer(psb_ipk_), intent(out) :: info - ! end subroutine psb_c_csr_scalplusidentity - ! end interface - - !> \namespace psb_base_mod \class psb_lc_csr_sparse_mat !! \extends psb_lc_base_mat_mod::psb_lc_base_sparse_mat !! diff --git a/base/modules/serial/psb_d_csc_mat_mod.f90 b/base/modules/serial/psb_d_csc_mat_mod.f90 index 5da0ccd3..42cefb72 100644 --- a/base/modules/serial/psb_d_csc_mat_mod.f90 +++ b/base/modules/serial/psb_d_csc_mat_mod.f90 @@ -71,7 +71,6 @@ module psb_d_csc_mat_mod procedure, pass(a) :: inner_cssv => psb_d_csc_cssv procedure, pass(a) :: scals => psb_d_csc_scals procedure, pass(a) :: scalv => psb_d_csc_scal -! procedure, pass(a) :: scalpid => psb_d_csc_scalplusidentity procedure, pass(a) :: maxval => psb_d_csc_maxval procedure, pass(a) :: spnm1 => psb_d_csc_csnm1 procedure, pass(a) :: rowsum => psb_d_csc_rowsum @@ -128,7 +127,6 @@ module psb_d_csc_mat_mod procedure, pass(a) :: sizeof => ld_csc_sizeof procedure, pass(a) :: scals => psb_ld_csc_scals procedure, pass(a) :: scalv => psb_ld_csc_scal -! procedure, pass(a) :: scalpid => psb_ld_csc_scalplusidentity procedure, pass(a) :: maxval => psb_ld_csc_maxval procedure, pass(a) :: spnm1 => psb_ld_csc_csnm1 procedure, pass(a) :: rowsum => psb_ld_csc_rowsum @@ -563,18 +561,6 @@ module psb_d_csc_mat_mod end subroutine psb_d_csc_scals end interface - !> \memberof psb_d_csc_sparse_mat - !! \see psb_d_base_mat_mod::psb_d_base_scalplusidentity - ! interface - ! subroutine psb_d_csc_scalplusidentity(d,a,info) - ! import - ! class(psb_d_csc_sparse_mat), intent(inout) :: a - ! real(psb_dpk_), intent(in) :: d - ! integer(psb_ipk_), intent(out) :: info - ! end subroutine psb_d_csc_scalplusidentity - ! end interface - - ! ! ld ! @@ -926,17 +912,6 @@ module psb_d_csc_mat_mod end subroutine psb_ld_csc_scals end interface - !> \memberof psb_ld_csc_sparse_mat - !! \see psb_ld_base_mat_mod::psb_ld_base_scalplusidentity -! interface -! subroutine psb_ld_csc_scalplusidentity(d,a,info) -! import -! class(psb_ld_csc_sparse_mat), intent(inout) :: a -! real(psb_dpk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! end subroutine psb_ld_csc_scalplusidentity -! end interface - contains ! == =================================== diff --git a/base/modules/serial/psb_d_csr_mat_mod.f90 b/base/modules/serial/psb_d_csr_mat_mod.f90 index 74dfe300..e6f926e8 100644 --- a/base/modules/serial/psb_d_csr_mat_mod.f90 +++ b/base/modules/serial/psb_d_csr_mat_mod.f90 @@ -73,7 +73,6 @@ module psb_d_csr_mat_mod procedure, pass(a) :: inner_cssv => psb_d_csr_cssv procedure, pass(a) :: scals => psb_d_csr_scals procedure, pass(a) :: scalv => psb_d_csr_scal -! procedure, pass(a) :: scalpid => psb_d_csr_scalplusidentity procedure, pass(a) :: maxval => psb_d_csr_maxval procedure, pass(a) :: spnmi => psb_d_csr_csnmi procedure, pass(a) :: rowsum => psb_d_csr_rowsum @@ -580,18 +579,6 @@ module psb_d_csr_mat_mod end subroutine psb_d_csr_scals end interface - !> \memberof psb_d_csr_sparse_mat - !! \see psb_d_base_mat_mod::psb_d_base_scalplusidentity - ! interface - ! subroutine psb_d_csr_scalplusidentity(d,a,info) - ! import - ! class(psb_d_csr_sparse_mat), intent(inout) :: a - ! real(psb_dpk_), intent(in) :: d - ! integer(psb_ipk_), intent(out) :: info - ! end subroutine psb_d_csr_scalplusidentity - ! end interface - - !> \namespace psb_base_mod \class psb_ld_csr_sparse_mat !! \extends psb_ld_base_mat_mod::psb_ld_base_sparse_mat !! diff --git a/base/modules/serial/psb_s_csc_mat_mod.f90 b/base/modules/serial/psb_s_csc_mat_mod.f90 index f513cb72..822659d9 100644 --- a/base/modules/serial/psb_s_csc_mat_mod.f90 +++ b/base/modules/serial/psb_s_csc_mat_mod.f90 @@ -71,7 +71,6 @@ module psb_s_csc_mat_mod procedure, pass(a) :: inner_cssv => psb_s_csc_cssv procedure, pass(a) :: scals => psb_s_csc_scals procedure, pass(a) :: scalv => psb_s_csc_scal -! procedure, pass(a) :: scalpid => psb_s_csc_scalplusidentity procedure, pass(a) :: maxval => psb_s_csc_maxval procedure, pass(a) :: spnm1 => psb_s_csc_csnm1 procedure, pass(a) :: rowsum => psb_s_csc_rowsum @@ -128,7 +127,6 @@ module psb_s_csc_mat_mod procedure, pass(a) :: sizeof => ls_csc_sizeof procedure, pass(a) :: scals => psb_ls_csc_scals procedure, pass(a) :: scalv => psb_ls_csc_scal -! procedure, pass(a) :: scalpid => psb_ls_csc_scalplusidentity procedure, pass(a) :: maxval => psb_ls_csc_maxval procedure, pass(a) :: spnm1 => psb_ls_csc_csnm1 procedure, pass(a) :: rowsum => psb_ls_csc_rowsum @@ -563,18 +561,6 @@ module psb_s_csc_mat_mod end subroutine psb_s_csc_scals end interface - !> \memberof psb_s_csc_sparse_mat - !! \see psb_s_base_mat_mod::psb_s_base_scalplusidentity - ! interface - ! subroutine psb_s_csc_scalplusidentity(d,a,info) - ! import - ! class(psb_s_csc_sparse_mat), intent(inout) :: a - ! real(psb_spk_), intent(in) :: d - ! integer(psb_ipk_), intent(out) :: info - ! end subroutine psb_s_csc_scalplusidentity - ! end interface - - ! ! ls ! @@ -926,17 +912,6 @@ module psb_s_csc_mat_mod end subroutine psb_ls_csc_scals end interface - !> \memberof psb_ls_csc_sparse_mat - !! \see psb_ls_base_mat_mod::psb_ls_base_scalplusidentity -! interface -! subroutine psb_ls_csc_scalplusidentity(d,a,info) -! import -! class(psb_ls_csc_sparse_mat), intent(inout) :: a -! real(psb_spk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! end subroutine psb_ls_csc_scalplusidentity -! end interface - contains ! == =================================== diff --git a/base/modules/serial/psb_s_csr_mat_mod.f90 b/base/modules/serial/psb_s_csr_mat_mod.f90 index 4a3cae6e..01d36829 100644 --- a/base/modules/serial/psb_s_csr_mat_mod.f90 +++ b/base/modules/serial/psb_s_csr_mat_mod.f90 @@ -73,7 +73,6 @@ module psb_s_csr_mat_mod procedure, pass(a) :: inner_cssv => psb_s_csr_cssv procedure, pass(a) :: scals => psb_s_csr_scals procedure, pass(a) :: scalv => psb_s_csr_scal -! procedure, pass(a) :: scalpid => psb_s_csr_scalplusidentity procedure, pass(a) :: maxval => psb_s_csr_maxval procedure, pass(a) :: spnmi => psb_s_csr_csnmi procedure, pass(a) :: rowsum => psb_s_csr_rowsum @@ -580,18 +579,6 @@ module psb_s_csr_mat_mod end subroutine psb_s_csr_scals end interface - !> \memberof psb_s_csr_sparse_mat - !! \see psb_s_base_mat_mod::psb_s_base_scalplusidentity - ! interface - ! subroutine psb_s_csr_scalplusidentity(d,a,info) - ! import - ! class(psb_s_csr_sparse_mat), intent(inout) :: a - ! real(psb_spk_), intent(in) :: d - ! integer(psb_ipk_), intent(out) :: info - ! end subroutine psb_s_csr_scalplusidentity - ! end interface - - !> \namespace psb_base_mod \class psb_ls_csr_sparse_mat !! \extends psb_ls_base_mat_mod::psb_ls_base_sparse_mat !! diff --git a/base/modules/serial/psb_z_csc_mat_mod.f90 b/base/modules/serial/psb_z_csc_mat_mod.f90 index 05c62cb8..b3b3e8fe 100644 --- a/base/modules/serial/psb_z_csc_mat_mod.f90 +++ b/base/modules/serial/psb_z_csc_mat_mod.f90 @@ -71,7 +71,6 @@ module psb_z_csc_mat_mod procedure, pass(a) :: inner_cssv => psb_z_csc_cssv procedure, pass(a) :: scals => psb_z_csc_scals procedure, pass(a) :: scalv => psb_z_csc_scal -! procedure, pass(a) :: scalpid => psb_z_csc_scalplusidentity procedure, pass(a) :: maxval => psb_z_csc_maxval procedure, pass(a) :: spnm1 => psb_z_csc_csnm1 procedure, pass(a) :: rowsum => psb_z_csc_rowsum @@ -128,7 +127,6 @@ module psb_z_csc_mat_mod procedure, pass(a) :: sizeof => lz_csc_sizeof procedure, pass(a) :: scals => psb_lz_csc_scals procedure, pass(a) :: scalv => psb_lz_csc_scal -! procedure, pass(a) :: scalpid => psb_lz_csc_scalplusidentity procedure, pass(a) :: maxval => psb_lz_csc_maxval procedure, pass(a) :: spnm1 => psb_lz_csc_csnm1 procedure, pass(a) :: rowsum => psb_lz_csc_rowsum @@ -563,18 +561,6 @@ module psb_z_csc_mat_mod end subroutine psb_z_csc_scals end interface - !> \memberof psb_z_csc_sparse_mat - !! \see psb_z_base_mat_mod::psb_z_base_scalplusidentity - ! interface - ! subroutine psb_z_csc_scalplusidentity(d,a,info) - ! import - ! class(psb_z_csc_sparse_mat), intent(inout) :: a - ! complex(psb_dpk_), intent(in) :: d - ! integer(psb_ipk_), intent(out) :: info - ! end subroutine psb_z_csc_scalplusidentity - ! end interface - - ! ! lz ! @@ -926,17 +912,6 @@ module psb_z_csc_mat_mod end subroutine psb_lz_csc_scals end interface - !> \memberof psb_lz_csc_sparse_mat - !! \see psb_lz_base_mat_mod::psb_lz_base_scalplusidentity -! interface -! subroutine psb_lz_csc_scalplusidentity(d,a,info) -! import -! class(psb_lz_csc_sparse_mat), intent(inout) :: a -! complex(psb_dpk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! end subroutine psb_lz_csc_scalplusidentity -! end interface - contains ! == =================================== diff --git a/base/modules/serial/psb_z_csr_mat_mod.f90 b/base/modules/serial/psb_z_csr_mat_mod.f90 index cf11e1dd..406ef6b4 100644 --- a/base/modules/serial/psb_z_csr_mat_mod.f90 +++ b/base/modules/serial/psb_z_csr_mat_mod.f90 @@ -73,7 +73,6 @@ module psb_z_csr_mat_mod procedure, pass(a) :: inner_cssv => psb_z_csr_cssv procedure, pass(a) :: scals => psb_z_csr_scals procedure, pass(a) :: scalv => psb_z_csr_scal -! procedure, pass(a) :: scalpid => psb_z_csr_scalplusidentity procedure, pass(a) :: maxval => psb_z_csr_maxval procedure, pass(a) :: spnmi => psb_z_csr_csnmi procedure, pass(a) :: rowsum => psb_z_csr_rowsum @@ -580,18 +579,6 @@ module psb_z_csr_mat_mod end subroutine psb_z_csr_scals end interface - !> \memberof psb_z_csr_sparse_mat - !! \see psb_z_base_mat_mod::psb_z_base_scalplusidentity - ! interface - ! subroutine psb_z_csr_scalplusidentity(d,a,info) - ! import - ! class(psb_z_csr_sparse_mat), intent(inout) :: a - ! complex(psb_dpk_), intent(in) :: d - ! integer(psb_ipk_), intent(out) :: info - ! end subroutine psb_z_csr_scalplusidentity - ! end interface - - !> \namespace psb_base_mod \class psb_lz_csr_sparse_mat !! \extends psb_lz_base_mat_mod::psb_lz_base_sparse_mat !! diff --git a/base/serial/impl/psb_c_csc_impl.f90 b/base/serial/impl/psb_c_csc_impl.f90 index 419ee70d..62f0b738 100644 --- a/base/serial/impl/psb_c_csc_impl.f90 +++ b/base/serial/impl/psb_c_csc_impl.f90 @@ -1485,49 +1485,6 @@ subroutine psb_c_csc_scals(d,a,info) end subroutine psb_c_csc_scals -! subroutine psb_c_csc_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_c_csc_mat_mod, psb_protect_name => psb_c_csc_scalplusidentity -! implicit none -! class(psb_c_csc_sparse_mat), intent(inout) :: a -! complex(psb_spk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_ipk_) :: err_act,mnm, i, j, k, m -! integer(psb_ipk_) :: ierr(5) -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%icp(i),a%icp(i+1)-1 -! j=a%ia(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + cone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! -! return -! -! end subroutine psb_c_csc_scalplusidentity - ! == =================================== ! @@ -3107,49 +3064,6 @@ subroutine psb_lc_csc_scals(d,a,info) end subroutine psb_lc_csc_scals -! subroutine psb_lc_csc_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_c_csc_mat_mod, psb_protect_name => psb_lc_csc_scalplusidentity -! implicit none -! class(psb_lc_csc_sparse_mat), intent(inout) :: a -! complex(psb_spk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_lpk_) :: mnm, i, j, k, m -! integer(psb_ipk_) :: err_act, ierr(5) -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%icp(i),a%icp(i+1)-1 -! j=a%ia(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + cone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! -! return -! -! end subroutine psb_lc_csc_scalplusidentity - function psb_lc_csc_maxval(a) result(res) use psb_error_mod use psb_c_csc_mat_mod, psb_protect_name => psb_lc_csc_maxval diff --git a/base/serial/impl/psb_c_csr_impl.f90 b/base/serial/impl/psb_c_csr_impl.f90 index 22f0850f..2c333b38 100644 --- a/base/serial/impl/psb_c_csr_impl.f90 +++ b/base/serial/impl/psb_c_csr_impl.f90 @@ -1677,50 +1677,6 @@ subroutine psb_c_csr_scals(d,a,info) end subroutine psb_c_csr_scals -! subroutine psb_c_csr_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_c_csr_mat_mod, psb_protect_name => psb_c_csr_scalplusidentity -! implicit none -! class(psb_c_csr_sparse_mat), intent(inout) :: a -! complex(psb_spk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_ipk_) :: err_act,mnm, i, j, k, m -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%irp(i),a%irp(i+1)-1 -! j=a%ja(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + cone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! return -! -! end subroutine psb_c_csr_scalplusidentity - - - - ! == =================================== ! ! diff --git a/base/serial/impl/psb_d_csc_impl.f90 b/base/serial/impl/psb_d_csc_impl.f90 index f55ce774..624332e5 100644 --- a/base/serial/impl/psb_d_csc_impl.f90 +++ b/base/serial/impl/psb_d_csc_impl.f90 @@ -1485,49 +1485,6 @@ subroutine psb_d_csc_scals(d,a,info) end subroutine psb_d_csc_scals -! subroutine psb_d_csc_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_d_csc_mat_mod, psb_protect_name => psb_d_csc_scalplusidentity -! implicit none -! class(psb_d_csc_sparse_mat), intent(inout) :: a -! real(psb_dpk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_ipk_) :: err_act,mnm, i, j, k, m -! integer(psb_ipk_) :: ierr(5) -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%icp(i),a%icp(i+1)-1 -! j=a%ia(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + done -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! -! return -! -! end subroutine psb_d_csc_scalplusidentity - ! == =================================== ! @@ -3107,49 +3064,6 @@ subroutine psb_ld_csc_scals(d,a,info) end subroutine psb_ld_csc_scals -! subroutine psb_ld_csc_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_d_csc_mat_mod, psb_protect_name => psb_ld_csc_scalplusidentity -! implicit none -! class(psb_ld_csc_sparse_mat), intent(inout) :: a -! real(psb_dpk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_lpk_) :: mnm, i, j, k, m -! integer(psb_ipk_) :: err_act, ierr(5) -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%icp(i),a%icp(i+1)-1 -! j=a%ia(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + done -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! -! return -! -! end subroutine psb_ld_csc_scalplusidentity - function psb_ld_csc_maxval(a) result(res) use psb_error_mod use psb_d_csc_mat_mod, psb_protect_name => psb_ld_csc_maxval diff --git a/base/serial/impl/psb_d_csr_impl.f90 b/base/serial/impl/psb_d_csr_impl.f90 index 70c2515e..1a99ae9a 100644 --- a/base/serial/impl/psb_d_csr_impl.f90 +++ b/base/serial/impl/psb_d_csr_impl.f90 @@ -1677,50 +1677,6 @@ subroutine psb_d_csr_scals(d,a,info) end subroutine psb_d_csr_scals -! subroutine psb_d_csr_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_d_csr_mat_mod, psb_protect_name => psb_d_csr_scalplusidentity -! implicit none -! class(psb_d_csr_sparse_mat), intent(inout) :: a -! real(psb_dpk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_ipk_) :: err_act,mnm, i, j, k, m -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%irp(i),a%irp(i+1)-1 -! j=a%ja(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + done -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! return -! -! end subroutine psb_d_csr_scalplusidentity - - - - ! == =================================== ! ! diff --git a/base/serial/impl/psb_s_csc_impl.f90 b/base/serial/impl/psb_s_csc_impl.f90 index fcd00755..40eeb40e 100644 --- a/base/serial/impl/psb_s_csc_impl.f90 +++ b/base/serial/impl/psb_s_csc_impl.f90 @@ -1485,49 +1485,6 @@ subroutine psb_s_csc_scals(d,a,info) end subroutine psb_s_csc_scals -! subroutine psb_s_csc_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_s_csc_mat_mod, psb_protect_name => psb_s_csc_scalplusidentity -! implicit none -! class(psb_s_csc_sparse_mat), intent(inout) :: a -! real(psb_spk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_ipk_) :: err_act,mnm, i, j, k, m -! integer(psb_ipk_) :: ierr(5) -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%icp(i),a%icp(i+1)-1 -! j=a%ia(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + sone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! -! return -! -! end subroutine psb_s_csc_scalplusidentity - ! == =================================== ! @@ -3107,49 +3064,6 @@ subroutine psb_ls_csc_scals(d,a,info) end subroutine psb_ls_csc_scals -! subroutine psb_ls_csc_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_s_csc_mat_mod, psb_protect_name => psb_ls_csc_scalplusidentity -! implicit none -! class(psb_ls_csc_sparse_mat), intent(inout) :: a -! real(psb_spk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_lpk_) :: mnm, i, j, k, m -! integer(psb_ipk_) :: err_act, ierr(5) -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%icp(i),a%icp(i+1)-1 -! j=a%ia(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + sone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! -! return -! -! end subroutine psb_ls_csc_scalplusidentity - function psb_ls_csc_maxval(a) result(res) use psb_error_mod use psb_s_csc_mat_mod, psb_protect_name => psb_ls_csc_maxval diff --git a/base/serial/impl/psb_s_csr_impl.f90 b/base/serial/impl/psb_s_csr_impl.f90 index 18c285ab..d60af4f4 100644 --- a/base/serial/impl/psb_s_csr_impl.f90 +++ b/base/serial/impl/psb_s_csr_impl.f90 @@ -1677,50 +1677,6 @@ subroutine psb_s_csr_scals(d,a,info) end subroutine psb_s_csr_scals -! subroutine psb_s_csr_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_s_csr_mat_mod, psb_protect_name => psb_s_csr_scalplusidentity -! implicit none -! class(psb_s_csr_sparse_mat), intent(inout) :: a -! real(psb_spk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_ipk_) :: err_act,mnm, i, j, k, m -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%irp(i),a%irp(i+1)-1 -! j=a%ja(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + sone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! return -! -! end subroutine psb_s_csr_scalplusidentity - - - - ! == =================================== ! ! diff --git a/base/serial/impl/psb_z_csc_impl.f90 b/base/serial/impl/psb_z_csc_impl.f90 index 715673ce..3cde422f 100644 --- a/base/serial/impl/psb_z_csc_impl.f90 +++ b/base/serial/impl/psb_z_csc_impl.f90 @@ -1485,49 +1485,6 @@ subroutine psb_z_csc_scals(d,a,info) end subroutine psb_z_csc_scals -! subroutine psb_z_csc_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_z_csc_mat_mod, psb_protect_name => psb_z_csc_scalplusidentity -! implicit none -! class(psb_z_csc_sparse_mat), intent(inout) :: a -! complex(psb_dpk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_ipk_) :: err_act,mnm, i, j, k, m -! integer(psb_ipk_) :: ierr(5) -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%icp(i),a%icp(i+1)-1 -! j=a%ia(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + zone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! -! return -! -! end subroutine psb_z_csc_scalplusidentity - ! == =================================== ! @@ -3107,49 +3064,6 @@ subroutine psb_lz_csc_scals(d,a,info) end subroutine psb_lz_csc_scals -! subroutine psb_lz_csc_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_z_csc_mat_mod, psb_protect_name => psb_lz_csc_scalplusidentity -! implicit none -! class(psb_lz_csc_sparse_mat), intent(inout) :: a -! complex(psb_dpk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_lpk_) :: mnm, i, j, k, m -! integer(psb_ipk_) :: err_act, ierr(5) -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%icp(i),a%icp(i+1)-1 -! j=a%ia(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + zone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! -! return -! -! end subroutine psb_lz_csc_scalplusidentity - function psb_lz_csc_maxval(a) result(res) use psb_error_mod use psb_z_csc_mat_mod, psb_protect_name => psb_lz_csc_maxval diff --git a/base/serial/impl/psb_z_csr_impl.f90 b/base/serial/impl/psb_z_csr_impl.f90 index 6071f563..43b40b3a 100644 --- a/base/serial/impl/psb_z_csr_impl.f90 +++ b/base/serial/impl/psb_z_csr_impl.f90 @@ -1677,50 +1677,6 @@ subroutine psb_z_csr_scals(d,a,info) end subroutine psb_z_csr_scals -! subroutine psb_z_csr_scalplusidentity(d,a,info) -! use psb_error_mod -! use psb_const_mod -! use psb_z_csr_mat_mod, psb_protect_name => psb_z_csr_scalplusidentity -! implicit none -! class(psb_z_csr_sparse_mat), intent(inout) :: a -! complex(psb_dpk_), intent(in) :: d -! integer(psb_ipk_), intent(out) :: info -! -! integer(psb_ipk_) :: err_act,mnm, i, j, k, m -! character(len=20) :: name='scalplusidentity' -! logical, parameter :: debug=.false. -! -! info = psb_success_ -! call psb_erractionsave(err_act) -! if (a%is_dev()) call a%sync() -! -! if (a%is_unit()) then -! call a%make_nonunit() -! end if -! -! mnm = min(a%get_nrows(),a%get_ncols()) -! do i=1,a%get_nzeros() -! a%val(i) = a%val(i) * d -! do k=a%irp(i),a%irp(i+1)-1 -! j=a%ja(k) -! if ((j == i) .and.(j <= mnm )) then -! a%val(k) = a%val(k) + zone -! endif -! enddo -! enddo -! call a%set_host() -! -! call psb_erractionrestore(err_act) -! return -! -! 9999 call psb_error_handler(err_act) -! return -! -! end subroutine psb_z_csr_scalplusidentity - - - - ! == =================================== ! !