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

Now base compiles cleanly.

Browse Source
ILmat
Salvatore Filippone 8 years ago
parent 4f4bde7ddc
commit 1bb802328b
32 changed files with 395 additions and 382 deletions
Show Stats Download Patch File Download Diff File

12
base/modules/Makefile
Unescape Escape View File

@ -212,12 +212,12 @@ serial/psb_c_csc_mat_mod.o serial/psb_c_csr_mat_mod.o: serial/psb_c_base_mat_mod
serial/psb_z_csc_mat_mod.o serial/psb_z_csr_mat_mod.o: serial/psb_z_base_mat_mod.o
serial/psb_mat_mod.o: serial/psb_vect_mod.o serial/psb_s_mat_mod.o serial/psb_d_mat_mod.o serial/psb_c_mat_mod.o serial/psb_z_mat_mod.o
serial/psb_serial_mod.o: serial/psb_s_serial_mod.o serial/psb_d_serial_mod.o serial/psb_c_serial_mod.o serial/psb_z_serial_mod.o
serial/psb_i_vect_mod.o: serial/psb_i_base_vect_mod.o
serial/psb_l_vect_mod.o: serial/psb_l_base_vect_mod.o serial/psb_i_vect_mod.o
serial/psb_s_vect_mod.o: serial/psb_s_base_vect_mod.o serial/psb_i_vect_mod.o
serial/psb_d_vect_mod.o: serial/psb_d_base_vect_mod.o serial/psb_i_vect_mod.o
serial/psb_c_vect_mod.o: serial/psb_c_base_vect_mod.o serial/psb_i_vect_mod.o
serial/psb_z_vect_mod.o: serial/psb_z_base_vect_mod.o serial/psb_i_vect_mod.o
serial/psb_i_vect_mod.o: serial/psb_i_base_vect_mod.o serial/psb_l_vect_mod.o
serial/psb_l_vect_mod.o: serial/psb_l_base_vect_mod.o
serial/psb_s_vect_mod.o: serial/psb_s_base_vect_mod.o serial/psb_l_vect_mod.o
serial/psb_d_vect_mod.o: serial/psb_d_base_vect_mod.o serial/psb_l_vect_mod.o
serial/psb_c_vect_mod.o: serial/psb_c_base_vect_mod.o serial/psb_l_vect_mod.o
serial/psb_z_vect_mod.o: serial/psb_z_base_vect_mod.o serial/psb_l_vect_mod.o
serial/psb_s_serial_mod.o serial/psb_d_serial_mod.o serial/psb_c_serial_mod.o serial/psb_z_serial_mod.o: serial/psb_mat_mod.o auxil/psb_string_mod.o auxil/psb_sort_mod.o auxil/psi_serial_mod.o
serial/psb_vect_mod.o: serial/psb_i_vect_mod.o serial/psb_l_vect_mod.o serial/psb_d_vect_mod.o serial/psb_s_vect_mod.o serial/psb_c_vect_mod.o serial/psb_z_vect_mod.o

4
base/modules/tools/psb_c_tools_a_mod.f90
Unescape Escape View File

@ -96,7 +96,7 @@ Module psb_c_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
complex(psb_spk_),intent(inout) :: x(:,:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -108,7 +108,7 @@ Module psb_c_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
complex(psb_spk_),intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl

25
base/modules/tools/psb_c_tools_mod.f90
Unescape Escape View File

@ -30,9 +30,10 @@
!
!
Module psb_c_tools_mod
use psb_desc_mod, only : psb_desc_type, psb_spk_, psb_ipk_
use psb_c_vect_mod, only : psb_c_base_vect_type, psb_c_vect_type, psb_i_vect_type
use psb_desc_mod, only : psb_desc_type, psb_spk_, psb_ipk_, psb_lpk_
use psb_c_vect_mod, only : psb_c_base_vect_type, psb_c_vect_type
use psb_c_mat_mod, only : psb_cspmat_type, psb_c_base_sparse_mat
use psb_l_vect_mod, only : psb_l_vect_type
use psb_c_multivect_mod, only : psb_c_base_multivect_type, psb_c_multivect_type
interface psb_geall
@ -126,7 +127,7 @@ Module psb_c_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_c_vect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -138,7 +139,7 @@ Module psb_c_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_c_vect_type), intent(inout) :: x
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_c_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -150,7 +151,7 @@ Module psb_c_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_c_vect_type), intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -162,7 +163,7 @@ Module psb_c_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_c_multivect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -240,9 +241,10 @@ Module psb_c_tools_mod
implicit none
type(psb_desc_type), intent(inout) :: desc_a
type(psb_cspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
complex(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild
logical, intent(in), optional :: local
end subroutine psb_cspins
@ -253,7 +255,7 @@ Module psb_c_tools_mod
type(psb_desc_type), intent(inout) :: desc_a
type(psb_cspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz
type(psb_i_vect_type), intent(inout) :: ia,ja
type(psb_l_vect_type), intent(inout) :: ia,ja
type(psb_c_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild
@ -265,9 +267,10 @@ Module psb_c_tools_mod
type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac
type(psb_cspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
complex(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
end subroutine psb_cspins_2desc
end interface

4
base/modules/tools/psb_d_tools_a_mod.f90
Unescape Escape View File

@ -96,7 +96,7 @@ Module psb_d_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
real(psb_dpk_),intent(inout) :: x(:,:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -108,7 +108,7 @@ Module psb_d_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
real(psb_dpk_),intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl

25
base/modules/tools/psb_d_tools_mod.f90
Unescape Escape View File

@ -30,9 +30,10 @@
!
!
Module psb_d_tools_mod
use psb_desc_mod, only : psb_desc_type, psb_dpk_, psb_ipk_
use psb_d_vect_mod, only : psb_d_base_vect_type, psb_d_vect_type, psb_i_vect_type
use psb_desc_mod, only : psb_desc_type, psb_dpk_, psb_ipk_, psb_lpk_
use psb_d_vect_mod, only : psb_d_base_vect_type, psb_d_vect_type
use psb_d_mat_mod, only : psb_dspmat_type, psb_d_base_sparse_mat
use psb_l_vect_mod, only : psb_l_vect_type
use psb_d_multivect_mod, only : psb_d_base_multivect_type, psb_d_multivect_type
interface psb_geall
@ -126,7 +127,7 @@ Module psb_d_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_d_vect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -138,7 +139,7 @@ Module psb_d_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_d_vect_type), intent(inout) :: x
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_d_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -150,7 +151,7 @@ Module psb_d_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_d_vect_type), intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -162,7 +163,7 @@ Module psb_d_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_d_multivect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -240,9 +241,10 @@ Module psb_d_tools_mod
implicit none
type(psb_desc_type), intent(inout) :: desc_a
type(psb_dspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
real(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild
logical, intent(in), optional :: local
end subroutine psb_dspins
@ -253,7 +255,7 @@ Module psb_d_tools_mod
type(psb_desc_type), intent(inout) :: desc_a
type(psb_dspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz
type(psb_i_vect_type), intent(inout) :: ia,ja
type(psb_l_vect_type), intent(inout) :: ia,ja
type(psb_d_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild
@ -265,9 +267,10 @@ Module psb_d_tools_mod
type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac
type(psb_dspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
real(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
end subroutine psb_dspins_2desc
end interface

4
base/modules/tools/psb_e_tools_a_mod.f90
Unescape Escape View File

@ -96,7 +96,7 @@ Module psb_e_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
integer(psb_epk_),intent(inout) :: x(:,:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_epk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -108,7 +108,7 @@ Module psb_e_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
integer(psb_epk_),intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_epk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl

9
base/modules/tools/psb_i_tools_mod.f90
Unescape Escape View File

@ -32,6 +32,7 @@
Module psb_i_tools_mod
use psb_desc_mod, only : psb_desc_type, psb_ipk_, psb_lpk_, psb_success_
use psb_i_vect_mod, only : psb_i_base_vect_type, psb_i_vect_type
use psb_l_vect_mod, only : psb_l_vect_type
use psb_i_multivect_mod, only : psb_i_base_multivect_type, psb_i_multivect_type
interface psb_geall
@ -125,7 +126,7 @@ Module psb_i_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_i_vect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_ipk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -137,7 +138,7 @@ Module psb_i_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_i_vect_type), intent(inout) :: x
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_i_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -149,7 +150,7 @@ Module psb_i_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_i_vect_type), intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_ipk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -161,7 +162,7 @@ Module psb_i_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_i_multivect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_ipk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl

11
base/modules/tools/psb_l_tools_mod.f90
Unescape Escape View File

@ -31,7 +31,8 @@
!
Module psb_l_tools_mod
use psb_desc_mod, only : psb_desc_type, psb_ipk_, psb_lpk_, psb_success_
use psb_l_vect_mod, only : psb_l_base_vect_type, psb_l_vect_type, psb_i_vect_type
use psb_l_vect_mod, only : psb_l_base_vect_type, psb_l_vect_type
! use psb_i_vect_mod, only : psb_i_vect_type
use psb_l_multivect_mod, only : psb_l_base_multivect_type, psb_l_multivect_type
interface psb_geall
@ -125,7 +126,7 @@ Module psb_l_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_l_vect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -137,7 +138,7 @@ Module psb_l_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_l_vect_type), intent(inout) :: x
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -149,7 +150,7 @@ Module psb_l_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_l_vect_type), intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -161,7 +162,7 @@ Module psb_l_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_l_multivect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl

4
base/modules/tools/psb_m_tools_a_mod.f90
Unescape Escape View File

@ -96,7 +96,7 @@ Module psb_m_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
integer(psb_mpk_),intent(inout) :: x(:,:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_mpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -108,7 +108,7 @@ Module psb_m_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
integer(psb_mpk_),intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_mpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl

4
base/modules/tools/psb_s_tools_a_mod.f90
Unescape Escape View File

@ -96,7 +96,7 @@ Module psb_s_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
real(psb_spk_),intent(inout) :: x(:,:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -108,7 +108,7 @@ Module psb_s_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
real(psb_spk_),intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl

25
base/modules/tools/psb_s_tools_mod.f90
Unescape Escape View File

@ -30,9 +30,10 @@
!
!
Module psb_s_tools_mod
use psb_desc_mod, only : psb_desc_type, psb_spk_, psb_ipk_
use psb_s_vect_mod, only : psb_s_base_vect_type, psb_s_vect_type, psb_i_vect_type
use psb_desc_mod, only : psb_desc_type, psb_spk_, psb_ipk_, psb_lpk_
use psb_s_vect_mod, only : psb_s_base_vect_type, psb_s_vect_type
use psb_s_mat_mod, only : psb_sspmat_type, psb_s_base_sparse_mat
use psb_l_vect_mod, only : psb_l_vect_type
use psb_s_multivect_mod, only : psb_s_base_multivect_type, psb_s_multivect_type
interface psb_geall
@ -126,7 +127,7 @@ Module psb_s_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_s_vect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -138,7 +139,7 @@ Module psb_s_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_s_vect_type), intent(inout) :: x
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_s_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -150,7 +151,7 @@ Module psb_s_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_s_vect_type), intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -162,7 +163,7 @@ Module psb_s_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_s_multivect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -240,9 +241,10 @@ Module psb_s_tools_mod
implicit none
type(psb_desc_type), intent(inout) :: desc_a
type(psb_sspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
real(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild
logical, intent(in), optional :: local
end subroutine psb_sspins
@ -253,7 +255,7 @@ Module psb_s_tools_mod
type(psb_desc_type), intent(inout) :: desc_a
type(psb_sspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz
type(psb_i_vect_type), intent(inout) :: ia,ja
type(psb_l_vect_type), intent(inout) :: ia,ja
type(psb_s_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild
@ -265,9 +267,10 @@ Module psb_s_tools_mod
type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac
type(psb_sspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
real(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
end subroutine psb_sspins_2desc
end interface

4
base/modules/tools/psb_z_tools_a_mod.f90
Unescape Escape View File

@ -96,7 +96,7 @@ Module psb_z_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
complex(psb_dpk_),intent(inout) :: x(:,:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -108,7 +108,7 @@ Module psb_z_tools_a_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
complex(psb_dpk_),intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl

25
base/modules/tools/psb_z_tools_mod.f90
Unescape Escape View File

@ -30,9 +30,10 @@
!
!
Module psb_z_tools_mod
use psb_desc_mod, only : psb_desc_type, psb_dpk_, psb_ipk_
use psb_z_vect_mod, only : psb_z_base_vect_type, psb_z_vect_type, psb_i_vect_type
use psb_desc_mod, only : psb_desc_type, psb_dpk_, psb_ipk_, psb_lpk_
use psb_z_vect_mod, only : psb_z_base_vect_type, psb_z_vect_type
use psb_z_mat_mod, only : psb_zspmat_type, psb_z_base_sparse_mat
use psb_l_vect_mod, only : psb_l_vect_type
use psb_z_multivect_mod, only : psb_z_base_multivect_type, psb_z_multivect_type
interface psb_geall
@ -126,7 +127,7 @@ Module psb_z_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_z_vect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -138,7 +139,7 @@ Module psb_z_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_z_vect_type), intent(inout) :: x
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_z_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -150,7 +151,7 @@ Module psb_z_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_z_vect_type), intent(inout) :: x(:)
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -162,7 +163,7 @@ Module psb_z_tools_mod
integer(psb_ipk_), intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a
type(psb_z_multivect_type), intent(inout) :: x
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:,:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), optional, intent(in) :: dupl
@ -240,9 +241,10 @@ Module psb_z_tools_mod
implicit none
type(psb_desc_type), intent(inout) :: desc_a
type(psb_zspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
complex(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild
logical, intent(in), optional :: local
end subroutine psb_zspins
@ -253,7 +255,7 @@ Module psb_z_tools_mod
type(psb_desc_type), intent(inout) :: desc_a
type(psb_zspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz
type(psb_i_vect_type), intent(inout) :: ia,ja
type(psb_l_vect_type), intent(inout) :: ia,ja
type(psb_z_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild
@ -265,9 +267,10 @@ Module psb_z_tools_mod
type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac
type(psb_zspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
complex(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
end subroutine psb_zspins_2desc
end interface

2
base/tools/psb_cd_switch_ovl_indxmap.f90
Unescape Escape View File

@ -45,7 +45,7 @@ Subroutine psb_cd_switch_ovl_indxmap(desc,info)
integer(psb_ipk_), intent(out) :: info
! .. Local Scalars ..
integer(psb_ipk_) :: i, j, np, me, mglob, ictxt, n_row, n_col
integer(psb_ipk_) :: i, j, np, me, ictxt, n_row, n_col
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: err_act

2
base/tools/psb_cdall.f90
Unescape Escape View File

@ -54,7 +54,7 @@ subroutine psb_cdall(ictxt, desc, info,mg,ng,parts,vg,vl,flag,nl,repl, globalche
character(len=20) :: name
integer(psb_ipk_) :: err_act, flag_, i, me, np, nlp, nnv, lr
integer(psb_lpk_) :: n_
integer(psb_ipk_), allocatable :: itmpsz(:)
integer(psb_lpk_), allocatable :: itmpsz(:)
integer(psb_mpk_) :: iictxt

3
base/tools/psb_cdrep.f90
Unescape Escape View File

@ -115,7 +115,8 @@ subroutine psb_cdrep(m, ictxt, desc, info)
!locals
integer(psb_ipk_) :: i,np,me,err,err_act
integer(psb_lpk_) :: n
integer(psb_ipk_) :: int_err(5),exch(2), thalo(1), tovr(1), text(1)
integer(psb_ipk_) :: int_err(5),exch(2)
integer(psb_lpk_) :: thalo(1), tovr(1), text(1)
integer(psb_ipk_) :: debug_level, debug_unit
integer(psb_mpk_) :: iictxt
character(len=20) :: name

31
base/tools/psb_cins.f90
Unescape Escape View File

@ -56,7 +56,7 @@ subroutine psb_cins_vect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:)
type(psb_c_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -66,7 +66,8 @@ subroutine psb_cins_vect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -170,7 +171,7 @@ subroutine psb_cins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_c_vect_type), intent(inout) :: val
type(psb_c_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -180,7 +181,8 @@ subroutine psb_cins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
complex(psb_spk_), allocatable :: lval(:)
@ -246,15 +248,14 @@ subroutine psb_cins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
local_ = .false.
endif
if (irw%is_dev()) call irw%sync()
if (local_) then
call x%ins(m,irw,val,dupl_,info)
irl(1:m) = irw%v%v(1:m)
else
irl = irw%get_vect()
lval = val%get_vect()
call desc_a%indxmap%g2lip(irl(1:m),info,owned=.true.)
call x%ins(m,irl,lval,dupl_,info)
call desc_a%indxmap%g2l(irw%v%v(1:m),irl(1:m),info,owned=.true.)
end if
call x%ins(m,irl,lval,dupl_,info)
if (info /= 0) then
call psb_errpush(info,name)
goto 9999
@ -280,7 +281,7 @@ subroutine psb_cins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:,:)
type(psb_c_vect_type), intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -290,7 +291,8 @@ subroutine psb_cins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5), n
& loc_rows,loc_cols,err_act, int_err(5), n
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -401,7 +403,7 @@ subroutine psb_cins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:,:)
type(psb_c_multivect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -411,7 +413,8 @@ subroutine psb_cins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

10
base/tools/psb_cins_a.f90
Unescape Escape View File

@ -57,7 +57,7 @@ subroutine psb_cinsvi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:)
complex(psb_spk_),intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -67,7 +67,8 @@ subroutine psb_cinsvi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -239,7 +240,7 @@ subroutine psb_cinsi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_spk_), intent(in) :: val(:,:)
complex(psb_spk_),intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a
@ -249,7 +250,8 @@ subroutine psb_cinsi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,loc_row,j,n,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np,me,dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

92
base/tools/psb_cspins.f90
Unescape Escape View File

@ -56,10 +56,11 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild,local)
!....parameters...
type(psb_desc_type), intent(inout) :: desc_a
type(psb_cspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
complex(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
!locals.....
integer(psb_ipk_) :: nrow, err_act, ncol, spstate
@ -159,31 +160,25 @@ subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild,local)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (local_) then
call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
ila(1:nz) = ia(1:nz)
jla(1:nz) = ja(1:nz)
else
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info)
if (info == 0) call desc_a%indxmap%g2l(ja(1:nz),jla(1:nz),info)
call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
end if
call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
else
info = psb_err_invalid_cd_state_
@ -210,9 +205,10 @@ subroutine psb_cspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac
type(psb_cspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
complex(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
!locals.....
integer(psb_ipk_) :: nrow, err_act, ncol, spstate
@ -327,7 +323,7 @@ subroutine psb_cspins_v(nz,ia,ja,val,a,desc_a,info,rebuild,local)
type(psb_desc_type), intent(inout) :: desc_a
type(psb_cspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz
type(psb_i_vect_type), intent(inout) :: ia,ja
type(psb_l_vect_type), intent(inout) :: ia,ja
type(psb_c_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
@ -433,33 +429,29 @@ subroutine psb_cspins_v(nz,ia,ja,val,a,desc_a,info,rebuild,local)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (ia%is_dev()) call ia%sync()
if (ja%is_dev()) call ja%sync()
if (val%is_dev()) call val%sync()
if (local_) then
call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
ila(1:nz) = ia%v%v(1:nz)
jla(1:nz) = ja%v%v(1:nz)
else
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (ia%is_dev()) call ia%sync()
if (ja%is_dev()) call ja%sync()
if (val%is_dev()) call val%sync()
call desc_a%indxmap%g2l(ia%v%v(1:nz),ila(1:nz),info)
if (info == 0) call desc_a%indxmap%g2l(ja%v%v(1:nz),jla(1:nz),info)
if (info == 0) call a%csput(nz,ila,jla,val%v%v,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
end if
if (info == 0) call a%csput(nz,ila,jla,val%v%v,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
else
info = psb_err_invalid_cd_state_

31
base/tools/psb_dins.f90
Unescape Escape View File

@ -56,7 +56,7 @@ subroutine psb_dins_vect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:)
type(psb_d_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -66,7 +66,8 @@ subroutine psb_dins_vect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -170,7 +171,7 @@ subroutine psb_dins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_d_vect_type), intent(inout) :: val
type(psb_d_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -180,7 +181,8 @@ subroutine psb_dins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
real(psb_dpk_), allocatable :: lval(:)
@ -246,15 +248,14 @@ subroutine psb_dins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
local_ = .false.
endif
if (irw%is_dev()) call irw%sync()
if (local_) then
call x%ins(m,irw,val,dupl_,info)
irl(1:m) = irw%v%v(1:m)
else
irl = irw%get_vect()
lval = val%get_vect()
call desc_a%indxmap%g2lip(irl(1:m),info,owned=.true.)
call x%ins(m,irl,lval,dupl_,info)
call desc_a%indxmap%g2l(irw%v%v(1:m),irl(1:m),info,owned=.true.)
end if
call x%ins(m,irl,lval,dupl_,info)
if (info /= 0) then
call psb_errpush(info,name)
goto 9999
@ -280,7 +281,7 @@ subroutine psb_dins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:,:)
type(psb_d_vect_type), intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -290,7 +291,8 @@ subroutine psb_dins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5), n
& loc_rows,loc_cols,err_act, int_err(5), n
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -401,7 +403,7 @@ subroutine psb_dins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:,:)
type(psb_d_multivect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -411,7 +413,8 @@ subroutine psb_dins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

10
base/tools/psb_dins_a.f90
Unescape Escape View File

@ -57,7 +57,7 @@ subroutine psb_dinsvi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:)
real(psb_dpk_),intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -67,7 +67,8 @@ subroutine psb_dinsvi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -239,7 +240,7 @@ subroutine psb_dinsi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_dpk_), intent(in) :: val(:,:)
real(psb_dpk_),intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a
@ -249,7 +250,8 @@ subroutine psb_dinsi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,loc_row,j,n,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np,me,dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

92
base/tools/psb_dspins.f90
Unescape Escape View File

@ -56,10 +56,11 @@ subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild,local)
!....parameters...
type(psb_desc_type), intent(inout) :: desc_a
type(psb_dspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
real(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
!locals.....
integer(psb_ipk_) :: nrow, err_act, ncol, spstate
@ -159,31 +160,25 @@ subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild,local)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (local_) then
call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
ila(1:nz) = ia(1:nz)
jla(1:nz) = ja(1:nz)
else
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info)
if (info == 0) call desc_a%indxmap%g2l(ja(1:nz),jla(1:nz),info)
call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
end if
call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
else
info = psb_err_invalid_cd_state_
@ -210,9 +205,10 @@ subroutine psb_dspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac
type(psb_dspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
real(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
!locals.....
integer(psb_ipk_) :: nrow, err_act, ncol, spstate
@ -327,7 +323,7 @@ subroutine psb_dspins_v(nz,ia,ja,val,a,desc_a,info,rebuild,local)
type(psb_desc_type), intent(inout) :: desc_a
type(psb_dspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz
type(psb_i_vect_type), intent(inout) :: ia,ja
type(psb_l_vect_type), intent(inout) :: ia,ja
type(psb_d_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
@ -433,33 +429,29 @@ subroutine psb_dspins_v(nz,ia,ja,val,a,desc_a,info,rebuild,local)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (ia%is_dev()) call ia%sync()
if (ja%is_dev()) call ja%sync()
if (val%is_dev()) call val%sync()
if (local_) then
call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
ila(1:nz) = ia%v%v(1:nz)
jla(1:nz) = ja%v%v(1:nz)
else
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (ia%is_dev()) call ia%sync()
if (ja%is_dev()) call ja%sync()
if (val%is_dev()) call val%sync()
call desc_a%indxmap%g2l(ia%v%v(1:nz),ila(1:nz),info)
if (info == 0) call desc_a%indxmap%g2l(ja%v%v(1:nz),jla(1:nz),info)
if (info == 0) call a%csput(nz,ila,jla,val%v%v,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
end if
if (info == 0) call a%csput(nz,ila,jla,val%v%v,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
else
info = psb_err_invalid_cd_state_

10
base/tools/psb_eins_a.f90
Unescape Escape View File

@ -57,7 +57,7 @@ subroutine psb_einsvi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_epk_), intent(in) :: val(:)
integer(psb_epk_),intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -67,7 +67,8 @@ subroutine psb_einsvi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -239,7 +240,7 @@ subroutine psb_einsi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_epk_), intent(in) :: val(:,:)
integer(psb_epk_),intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a
@ -249,7 +250,8 @@ subroutine psb_einsi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,loc_row,j,n,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np,me,dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

31
base/tools/psb_iins.f90
Unescape Escape View File

@ -56,7 +56,7 @@ subroutine psb_iins_vect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_ipk_), intent(in) :: val(:)
type(psb_i_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -66,7 +66,8 @@ subroutine psb_iins_vect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -170,7 +171,7 @@ subroutine psb_iins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_i_vect_type), intent(inout) :: val
type(psb_i_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -180,7 +181,8 @@ subroutine psb_iins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
integer(psb_ipk_), allocatable :: lval(:)
@ -246,15 +248,14 @@ subroutine psb_iins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
local_ = .false.
endif
if (irw%is_dev()) call irw%sync()
if (local_) then
call x%ins(m,irw,val,dupl_,info)
irl(1:m) = irw%v%v(1:m)
else
irl = irw%get_vect()
lval = val%get_vect()
call desc_a%indxmap%g2lip(irl(1:m),info,owned=.true.)
call x%ins(m,irl,lval,dupl_,info)
call desc_a%indxmap%g2l(irw%v%v(1:m),irl(1:m),info,owned=.true.)
end if
call x%ins(m,irl,lval,dupl_,info)
if (info /= 0) then
call psb_errpush(info,name)
goto 9999
@ -280,7 +281,7 @@ subroutine psb_iins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_ipk_), intent(in) :: val(:,:)
type(psb_i_vect_type), intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -290,7 +291,8 @@ subroutine psb_iins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5), n
& loc_rows,loc_cols,err_act, int_err(5), n
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -401,7 +403,7 @@ subroutine psb_iins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_ipk_), intent(in) :: val(:,:)
type(psb_i_multivect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -411,7 +413,8 @@ subroutine psb_iins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

31
base/tools/psb_lins.f90
Unescape Escape View File

@ -56,7 +56,7 @@ subroutine psb_lins_vect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: val(:)
type(psb_l_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -66,7 +66,8 @@ subroutine psb_lins_vect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -170,7 +171,7 @@ subroutine psb_lins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: val
type(psb_l_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -180,7 +181,8 @@ subroutine psb_lins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
integer(psb_lpk_), allocatable :: lval(:)
@ -246,15 +248,14 @@ subroutine psb_lins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
local_ = .false.
endif
if (irw%is_dev()) call irw%sync()
if (local_) then
call x%ins(m,irw,val,dupl_,info)
irl(1:m) = irw%v%v(1:m)
else
irl = irw%get_vect()
lval = val%get_vect()
call desc_a%indxmap%g2lip(irl(1:m),info,owned=.true.)
call x%ins(m,irl,lval,dupl_,info)
call desc_a%indxmap%g2l(irw%v%v(1:m),irl(1:m),info,owned=.true.)
end if
call x%ins(m,irl,lval,dupl_,info)
if (info /= 0) then
call psb_errpush(info,name)
goto 9999
@ -280,7 +281,7 @@ subroutine psb_lins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: val(:,:)
type(psb_l_vect_type), intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -290,7 +291,8 @@ subroutine psb_lins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5), n
& loc_rows,loc_cols,err_act, int_err(5), n
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -401,7 +403,7 @@ subroutine psb_lins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: val(:,:)
type(psb_l_multivect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -411,7 +413,8 @@ subroutine psb_lins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

10
base/tools/psb_mins_a.f90
Unescape Escape View File

@ -57,7 +57,7 @@ subroutine psb_minsvi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_mpk_), intent(in) :: val(:)
integer(psb_mpk_),intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -67,7 +67,8 @@ subroutine psb_minsvi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -239,7 +240,7 @@ subroutine psb_minsi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
integer(psb_mpk_), intent(in) :: val(:,:)
integer(psb_mpk_),intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a
@ -249,7 +250,8 @@ subroutine psb_minsi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,loc_row,j,n,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np,me,dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

31
base/tools/psb_sins.f90
Unescape Escape View File

@ -56,7 +56,7 @@ subroutine psb_sins_vect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:)
type(psb_s_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -66,7 +66,8 @@ subroutine psb_sins_vect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -170,7 +171,7 @@ subroutine psb_sins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_s_vect_type), intent(inout) :: val
type(psb_s_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -180,7 +181,8 @@ subroutine psb_sins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
real(psb_spk_), allocatable :: lval(:)
@ -246,15 +248,14 @@ subroutine psb_sins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
local_ = .false.
endif
if (irw%is_dev()) call irw%sync()
if (local_) then
call x%ins(m,irw,val,dupl_,info)
irl(1:m) = irw%v%v(1:m)
else
irl = irw%get_vect()
lval = val%get_vect()
call desc_a%indxmap%g2lip(irl(1:m),info,owned=.true.)
call x%ins(m,irl,lval,dupl_,info)
call desc_a%indxmap%g2l(irw%v%v(1:m),irl(1:m),info,owned=.true.)
end if
call x%ins(m,irl,lval,dupl_,info)
if (info /= 0) then
call psb_errpush(info,name)
goto 9999
@ -280,7 +281,7 @@ subroutine psb_sins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:,:)
type(psb_s_vect_type), intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -290,7 +291,8 @@ subroutine psb_sins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5), n
& loc_rows,loc_cols,err_act, int_err(5), n
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -401,7 +403,7 @@ subroutine psb_sins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:,:)
type(psb_s_multivect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -411,7 +413,8 @@ subroutine psb_sins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

10
base/tools/psb_sins_a.f90
Unescape Escape View File

@ -57,7 +57,7 @@ subroutine psb_sinsvi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:)
real(psb_spk_),intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -67,7 +67,8 @@ subroutine psb_sinsvi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -239,7 +240,7 @@ subroutine psb_sinsi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
real(psb_spk_), intent(in) :: val(:,:)
real(psb_spk_),intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a
@ -249,7 +250,8 @@ subroutine psb_sinsi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,loc_row,j,n,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np,me,dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

92
base/tools/psb_sspins.f90
Unescape Escape View File

@ -56,10 +56,11 @@ subroutine psb_sspins(nz,ia,ja,val,a,desc_a,info,rebuild,local)
!....parameters...
type(psb_desc_type), intent(inout) :: desc_a
type(psb_sspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
real(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
!locals.....
integer(psb_ipk_) :: nrow, err_act, ncol, spstate
@ -159,31 +160,25 @@ subroutine psb_sspins(nz,ia,ja,val,a,desc_a,info,rebuild,local)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (local_) then
call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
ila(1:nz) = ia(1:nz)
jla(1:nz) = ja(1:nz)
else
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info)
if (info == 0) call desc_a%indxmap%g2l(ja(1:nz),jla(1:nz),info)
call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
end if
call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
else
info = psb_err_invalid_cd_state_
@ -210,9 +205,10 @@ subroutine psb_sspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac
type(psb_sspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
real(psb_spk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
!locals.....
integer(psb_ipk_) :: nrow, err_act, ncol, spstate
@ -327,7 +323,7 @@ subroutine psb_sspins_v(nz,ia,ja,val,a,desc_a,info,rebuild,local)
type(psb_desc_type), intent(inout) :: desc_a
type(psb_sspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz
type(psb_i_vect_type), intent(inout) :: ia,ja
type(psb_l_vect_type), intent(inout) :: ia,ja
type(psb_s_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
@ -433,33 +429,29 @@ subroutine psb_sspins_v(nz,ia,ja,val,a,desc_a,info,rebuild,local)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (ia%is_dev()) call ia%sync()
if (ja%is_dev()) call ja%sync()
if (val%is_dev()) call val%sync()
if (local_) then
call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
ila(1:nz) = ia%v%v(1:nz)
jla(1:nz) = ja%v%v(1:nz)
else
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (ia%is_dev()) call ia%sync()
if (ja%is_dev()) call ja%sync()
if (val%is_dev()) call val%sync()
call desc_a%indxmap%g2l(ia%v%v(1:nz),ila(1:nz),info)
if (info == 0) call desc_a%indxmap%g2l(ja%v%v(1:nz),jla(1:nz),info)
if (info == 0) call a%csput(nz,ila,jla,val%v%v,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
end if
if (info == 0) call a%csput(nz,ila,jla,val%v%v,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
else
info = psb_err_invalid_cd_state_

31
base/tools/psb_zins.f90
Unescape Escape View File

@ -56,7 +56,7 @@ subroutine psb_zins_vect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:)
type(psb_z_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -66,7 +66,8 @@ subroutine psb_zins_vect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -170,7 +171,7 @@ subroutine psb_zins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
type(psb_i_vect_type), intent(inout) :: irw
type(psb_l_vect_type), intent(inout) :: irw
type(psb_z_vect_type), intent(inout) :: val
type(psb_z_vect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -180,7 +181,8 @@ subroutine psb_zins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
complex(psb_dpk_), allocatable :: lval(:)
@ -246,15 +248,14 @@ subroutine psb_zins_vect_v(m, irw, val, x, desc_a, info, dupl,local)
local_ = .false.
endif
if (irw%is_dev()) call irw%sync()
if (local_) then
call x%ins(m,irw,val,dupl_,info)
irl(1:m) = irw%v%v(1:m)
else
irl = irw%get_vect()
lval = val%get_vect()
call desc_a%indxmap%g2lip(irl(1:m),info,owned=.true.)
call x%ins(m,irl,lval,dupl_,info)
call desc_a%indxmap%g2l(irw%v%v(1:m),irl(1:m),info,owned=.true.)
end if
call x%ins(m,irl,lval,dupl_,info)
if (info /= 0) then
call psb_errpush(info,name)
goto 9999
@ -280,7 +281,7 @@ subroutine psb_zins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:,:)
type(psb_z_vect_type), intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -290,7 +291,8 @@ subroutine psb_zins_vect_r2(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5), n
& loc_rows,loc_cols,err_act, int_err(5), n
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -401,7 +403,7 @@ subroutine psb_zins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:,:)
type(psb_z_multivect_type), intent(inout) :: x
type(psb_desc_type), intent(in) :: desc_a
@ -411,7 +413,8 @@ subroutine psb_zins_multivect(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

10
base/tools/psb_zins_a.f90
Unescape Escape View File

@ -57,7 +57,7 @@ subroutine psb_zinsvi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:)
complex(psb_dpk_),intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a
@ -67,7 +67,8 @@ subroutine psb_zinsvi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np, me, dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_
@ -239,7 +240,7 @@ subroutine psb_zinsi(m, irw, val, x, desc_a, info, dupl,local)
!....parameters...
integer(psb_ipk_), intent(in) :: m
integer(psb_ipk_), intent(in) :: irw(:)
integer(psb_lpk_), intent(in) :: irw(:)
complex(psb_dpk_), intent(in) :: val(:,:)
complex(psb_dpk_),intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a
@ -249,7 +250,8 @@ subroutine psb_zinsi(m, irw, val, x, desc_a, info, dupl,local)
!locals.....
integer(psb_ipk_) :: ictxt,i,loc_row,j,n,&
& loc_rows,loc_cols,mglob,err_act, int_err(5)
& loc_rows,loc_cols,err_act, int_err(5)
integer(psb_lpk_) :: mglob
integer(psb_ipk_) :: np,me,dupl_
integer(psb_ipk_), allocatable :: irl(:)
logical :: local_

92
base/tools/psb_zspins.f90
Unescape Escape View File

@ -56,10 +56,11 @@ subroutine psb_zspins(nz,ia,ja,val,a,desc_a,info,rebuild,local)
!....parameters...
type(psb_desc_type), intent(inout) :: desc_a
type(psb_zspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
complex(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
!locals.....
integer(psb_ipk_) :: nrow, err_act, ncol, spstate
@ -159,31 +160,25 @@ subroutine psb_zspins(nz,ia,ja,val,a,desc_a,info,rebuild,local)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (local_) then
call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
ila(1:nz) = ia(1:nz)
jla(1:nz) = ja(1:nz)
else
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
call desc_a%indxmap%g2l(ia(1:nz),ila(1:nz),info)
if (info == 0) call desc_a%indxmap%g2l(ja(1:nz),jla(1:nz),info)
call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
end if
call a%csput(nz,ila,jla,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
else
info = psb_err_invalid_cd_state_
@ -210,9 +205,10 @@ subroutine psb_zspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac
type(psb_zspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz,ia(:),ja(:)
integer(psb_ipk_), intent(in) :: nz
integer(psb_lpk_), intent(in) :: ia(:),ja(:)
complex(psb_dpk_), intent(in) :: val(:)
integer(psb_ipk_), intent(out) :: info
integer(psb_ipk_), intent(out) :: info
!locals.....
integer(psb_ipk_) :: nrow, err_act, ncol, spstate
@ -327,7 +323,7 @@ subroutine psb_zspins_v(nz,ia,ja,val,a,desc_a,info,rebuild,local)
type(psb_desc_type), intent(inout) :: desc_a
type(psb_zspmat_type), intent(inout) :: a
integer(psb_ipk_), intent(in) :: nz
type(psb_i_vect_type), intent(inout) :: ia,ja
type(psb_l_vect_type), intent(inout) :: ia,ja
type(psb_z_vect_type), intent(inout) :: val
integer(psb_ipk_), intent(out) :: info
logical, intent(in), optional :: rebuild, local
@ -433,33 +429,29 @@ subroutine psb_zspins_v(nz,ia,ja,val,a,desc_a,info,rebuild,local)
nrow = desc_a%get_local_rows()
ncol = desc_a%get_local_cols()
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (ia%is_dev()) call ia%sync()
if (ja%is_dev()) call ja%sync()
if (val%is_dev()) call val%sync()
if (local_) then
call a%csput(nz,ia,ja,val,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
ila(1:nz) = ia%v%v(1:nz)
jla(1:nz) = ja%v%v(1:nz)
else
allocate(ila(nz),jla(nz),stat=info)
if (info /= psb_success_) then
ierr(1) = info
call psb_errpush(psb_err_from_subroutine_ai_,name,&
& a_err='allocate',i_err=ierr)
goto 9999
end if
if (ia%is_dev()) call ia%sync()
if (ja%is_dev()) call ja%sync()
if (val%is_dev()) call val%sync()
call desc_a%indxmap%g2l(ia%v%v(1:nz),ila(1:nz),info)
if (info == 0) call desc_a%indxmap%g2l(ja%v%v(1:nz),jla(1:nz),info)
if (info == 0) call a%csput(nz,ila,jla,val%v%v,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
end if
if (info == 0) call a%csput(nz,ila,jla,val%v%v,ione,nrow,ione,ncol,info)
if (info /= psb_success_) then
info=psb_err_from_subroutine_
call psb_errpush(info,name,a_err='a%csput')
goto 9999
end if
else
info = psb_err_invalid_cd_state_

Write Preview
Loading…
Cancel
Save