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

psblas3:

Browse Source 
base/modules/psb_c_psblas_mod.f90
 base/modules/psb_c_tools_mod.f90
 base/modules/psb_d_psblas_mod.f90
 base/modules/psb_d_tools_mod.f90
 base/modules/psb_s_psblas_mod.f90
 base/modules/psb_s_tools_mod.f90
 base/modules/psb_serial_mod.f90
 base/modules/psb_z_psblas_mod.f90
 base/modules/psb_z_tools_mod.f90
 base/modules/psi_mod.f90
 base/serial/impl/psb_c_csc_impl.f90
 base/serial/impl/psb_c_csr_impl.f90
 base/serial/impl/psb_d_csc_impl.f90
 base/serial/impl/psb_d_csr_impl.f90
 base/serial/impl/psb_s_csc_impl.f90
 base/serial/impl/psb_s_csr_impl.f90
 base/serial/impl/psb_z_csc_impl.f90
 base/serial/impl/psb_z_csr_impl.f90

Second batch of fixes for unused variables.
psblas3-type-indexed
Salvatore Filippone 14 years ago
parent f093a14c8a
commit a291147ae3
18 changed files with 284 additions and 284 deletions
Show Stats Download Patch File Download Diff File

54
base/modules/psb_c_psblas_mod.f90
Unescape Escape View File

@ -33,14 +33,14 @@ module psb_c_psblas_mod
interface psb_gedot interface psb_gedot
function psb_cdotv(x, y, desc_a,info) function psb_cdotv(x, y, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_) :: psb_cdotv complex(psb_spk_) :: psb_cdotv
complex(psb_spk_), intent(in) :: x(:), y(:) complex(psb_spk_), intent(in) :: x(:), y(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end function psb_cdotv end function psb_cdotv
function psb_cdot(x, y, desc_a, info, jx, jy) function psb_cdot(x, y, desc_a, info, jx, jy)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_) :: psb_cdot complex(psb_spk_) :: psb_cdot
complex(psb_spk_), intent(in) :: x(:,:), y(:,:) complex(psb_spk_), intent(in) :: x(:,:), y(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -52,14 +52,14 @@ module psb_c_psblas_mod
interface psb_gedots interface psb_gedots
subroutine psb_cdotvs(res,x, y, desc_a, info) subroutine psb_cdotvs(res,x, y, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_), intent(out) :: res complex(psb_spk_), intent(out) :: res
complex(psb_spk_), intent(in) :: x(:), y(:) complex(psb_spk_), intent(in) :: x(:), y(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_cdotvs end subroutine psb_cdotvs
subroutine psb_cmdots(res,x, y, desc_a,info) subroutine psb_cmdots(res,x, y, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_), intent(out) :: res(:) complex(psb_spk_), intent(out) :: res(:)
complex(psb_spk_), intent(in) :: x(:,:), y(:,:) complex(psb_spk_), intent(in) :: x(:,:), y(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -70,7 +70,7 @@ module psb_c_psblas_mod
interface psb_geaxpby interface psb_geaxpby
subroutine psb_caxpbyv(alpha, x, beta, y,& subroutine psb_caxpbyv(alpha, x, beta, y,&
& desc_a, info) & desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_), intent (in) :: x(:) complex(psb_spk_), intent (in) :: x(:)
complex(psb_spk_), intent (inout) :: y(:) complex(psb_spk_), intent (inout) :: y(:)
complex(psb_spk_), intent (in) :: alpha, beta complex(psb_spk_), intent (in) :: alpha, beta
@ -79,7 +79,7 @@ module psb_c_psblas_mod
end subroutine psb_caxpbyv end subroutine psb_caxpbyv
subroutine psb_caxpby(alpha, x, beta, y,& subroutine psb_caxpby(alpha, x, beta, y,&
& desc_a, info, n, jx, jy) & desc_a, info, n, jx, jy)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_), intent (in) :: x(:,:) complex(psb_spk_), intent (in) :: x(:,:)
complex(psb_spk_), intent (inout) :: y(:,:) complex(psb_spk_), intent (inout) :: y(:,:)
complex(psb_spk_), intent (in) :: alpha, beta complex(psb_spk_), intent (in) :: alpha, beta
@ -91,7 +91,7 @@ module psb_c_psblas_mod
interface psb_geamax interface psb_geamax
function psb_camax(x, desc_a, info, jx) function psb_camax(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_camax real(psb_spk_) psb_camax
complex(psb_spk_), intent (in) :: x(:,:) complex(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -99,7 +99,7 @@ module psb_c_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_camax end function psb_camax
function psb_camaxv(x, desc_a,info) function psb_camaxv(x, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_camaxv real(psb_spk_) psb_camaxv
complex(psb_spk_), intent (in) :: x(:) complex(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -109,14 +109,14 @@ module psb_c_psblas_mod
interface psb_geamaxs interface psb_geamaxs
subroutine psb_camaxvs(res,x,desc_a,info) subroutine psb_camaxvs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res real(psb_spk_), intent (out) :: res
complex(psb_spk_), intent (in) :: x(:) complex(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_camaxvs end subroutine psb_camaxvs
subroutine psb_cmamaxs(res,x,desc_a,info,jx) subroutine psb_cmamaxs(res,x,desc_a,info,jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res(:) real(psb_spk_), intent (out) :: res(:)
complex(psb_spk_), intent (in) :: x(:,:) complex(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -127,7 +127,7 @@ module psb_c_psblas_mod
interface psb_geasum interface psb_geasum
function psb_casum(x, desc_a, info, jx) function psb_casum(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_casum real(psb_spk_) psb_casum
complex(psb_spk_), intent (in) :: x(:,:) complex(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -135,7 +135,7 @@ module psb_c_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_casum end function psb_casum
function psb_casumv(x, desc_a, info) function psb_casumv(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_casumv real(psb_spk_) psb_casumv
complex(psb_spk_), intent (in) :: x(:) complex(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -145,14 +145,14 @@ module psb_c_psblas_mod
interface psb_geasums interface psb_geasums
subroutine psb_casumvs(res,x,desc_a,info) subroutine psb_casumvs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res real(psb_spk_), intent (out) :: res
complex(psb_spk_), intent (in) :: x(:) complex(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_casumvs end subroutine psb_casumvs
subroutine psb_cmasum(res,x,desc_a,info) subroutine psb_cmasum(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res(:) real(psb_spk_), intent (out) :: res(:)
complex(psb_spk_), intent (in) :: x(:,:) complex(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -163,7 +163,7 @@ module psb_c_psblas_mod
interface psb_genrm2 interface psb_genrm2
function psb_cnrm2(x, desc_a, info, jx) function psb_cnrm2(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_cnrm2 real(psb_spk_) psb_cnrm2
complex(psb_spk_), intent (in) :: x(:,:) complex(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -171,7 +171,7 @@ module psb_c_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_cnrm2 end function psb_cnrm2
function psb_cnrm2v(x, desc_a, info) function psb_cnrm2v(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_cnrm2v real(psb_spk_) psb_cnrm2v
complex(psb_spk_), intent (in) :: x(:) complex(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -181,30 +181,30 @@ module psb_c_psblas_mod
interface psb_genrm2s interface psb_genrm2s
subroutine psb_cnrm2vs(res,x,desc_a,info) subroutine psb_cnrm2vs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res real(psb_spk_), intent (out) :: res
complex(psb_spk_), intent (in) :: x(:) complex(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_cnrm2vs end subroutine psb_cnrm2vs
end interface end interface
interface psb_spnrmi interface psb_spnrmi
function psb_cnrmi(a, desc_a,info) function psb_cnrmi(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
real(psb_spk_) :: psb_cnrmi real(psb_spk_) :: psb_cnrmi
type(psb_cspmat_type), intent (in) :: a type(psb_cspmat_type), intent (in) :: a
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end function psb_cnrmi end function psb_cnrmi
end interface end interface
interface psb_spmm interface psb_spmm
subroutine psb_cspmm(alpha, a, x, beta, y, desc_a, info,& subroutine psb_cspmm(alpha, a, x, beta, y, desc_a, info,&
&trans, k, jx, jy,work,doswap) &trans, k, jx, jy,work,doswap)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_cspmat_type), intent(in) :: a type(psb_cspmat_type), intent(in) :: a
complex(psb_spk_), intent(inout) :: x(:,:) complex(psb_spk_), intent(inout) :: x(:,:)
@ -219,7 +219,7 @@ module psb_c_psblas_mod
end subroutine psb_cspmm end subroutine psb_cspmm
subroutine psb_cspmv(alpha, a, x, beta, y,& subroutine psb_cspmv(alpha, a, x, beta, y,&
& desc_a, info, trans, work,doswap) & desc_a, info, trans, work,doswap)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_cspmat_type), intent(in) :: a type(psb_cspmat_type), intent(in) :: a
complex(psb_spk_), intent(inout) :: x(:) complex(psb_spk_), intent(inout) :: x(:)
@ -237,7 +237,7 @@ module psb_c_psblas_mod
subroutine psb_cspsm(alpha, t, x, beta, y,& subroutine psb_cspsm(alpha, t, x, beta, y,&
& desc_a, info, trans, scale, choice,& & desc_a, info, trans, scale, choice,&
& diag, n, jx, jy, work) & diag, n, jx, jy, work)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_cspmat_type), intent(in) :: t type(psb_cspmat_type), intent(in) :: t
complex(psb_spk_), intent(in) :: x(:,:) complex(psb_spk_), intent(in) :: x(:,:)
@ -254,7 +254,7 @@ module psb_c_psblas_mod
subroutine psb_cspsv(alpha, t, x, beta, y,& subroutine psb_cspsv(alpha, t, x, beta, y,&
& desc_a, info, trans, scale, choice,& & desc_a, info, trans, scale, choice,&
& diag, work) & diag, work)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_cspmat_type), intent(in) :: t type(psb_cspmat_type), intent(in) :: t
complex(psb_spk_), intent(in) :: x(:) complex(psb_spk_), intent(in) :: x(:)

32
base/modules/psb_c_tools_mod.f90
Unescape Escape View File

@ -34,7 +34,7 @@ Module psb_c_tools_mod
interface psb_geall interface psb_geall
subroutine psb_calloc(x, desc_a, info, n, lb) subroutine psb_calloc(x, desc_a, info, n, lb)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
implicit none implicit none
complex(psb_spk_), allocatable, intent(out) :: x(:,:) complex(psb_spk_), allocatable, intent(out) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -42,7 +42,7 @@ Module psb_c_tools_mod
integer, optional, intent(in) :: n, lb integer, optional, intent(in) :: n, lb
end subroutine psb_calloc end subroutine psb_calloc
subroutine psb_callocv(x, desc_a,info,n) subroutine psb_callocv(x, desc_a,info,n)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_), allocatable, intent(out) :: x(:) complex(psb_spk_), allocatable, intent(out) :: x(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
@ -53,13 +53,13 @@ Module psb_c_tools_mod
interface psb_geasb interface psb_geasb
subroutine psb_casb(x, desc_a, info) subroutine psb_casb(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
complex(psb_spk_), allocatable, intent(inout) :: x(:,:) complex(psb_spk_), allocatable, intent(inout) :: x(:,:)
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_casb end subroutine psb_casb
subroutine psb_casbv(x, desc_a, info) subroutine psb_casbv(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
complex(psb_spk_), allocatable, intent(inout) :: x(:) complex(psb_spk_), allocatable, intent(inout) :: x(:)
integer, intent(out) :: info integer, intent(out) :: info
@ -69,7 +69,7 @@ Module psb_c_tools_mod
interface psb_sphalo interface psb_sphalo
Subroutine psb_csphalo(a,desc_a,blk,info,rowcnv,colcnv,& Subroutine psb_csphalo(a,desc_a,blk,info,rowcnv,colcnv,&
& rowscale,colscale,outfmt,data) & rowscale,colscale,outfmt,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
Type(psb_cspmat_type),Intent(in) :: a Type(psb_cspmat_type),Intent(in) :: a
Type(psb_cspmat_type),Intent(inout) :: blk Type(psb_cspmat_type),Intent(inout) :: blk
@ -83,13 +83,13 @@ Module psb_c_tools_mod
interface psb_gefree interface psb_gefree
subroutine psb_cfree(x, desc_a, info) subroutine psb_cfree(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_),allocatable, intent(inout) :: x(:,:) complex(psb_spk_),allocatable, intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_cfree end subroutine psb_cfree
subroutine psb_cfreev(x, desc_a, info) subroutine psb_cfreev(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
complex(psb_spk_),allocatable, intent(inout) :: x(:) complex(psb_spk_),allocatable, intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
@ -99,7 +99,7 @@ Module psb_c_tools_mod
interface psb_geins interface psb_geins
subroutine psb_cinsi(m,irw,val, x, desc_a,info,dupl) subroutine psb_cinsi(m,irw,val, x, desc_a,info,dupl)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: m integer, intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
complex(psb_spk_),intent(inout) :: x(:,:) complex(psb_spk_),intent(inout) :: x(:,:)
@ -109,7 +109,7 @@ Module psb_c_tools_mod
integer, optional, intent(in) :: dupl integer, optional, intent(in) :: dupl
end subroutine psb_cinsi end subroutine psb_cinsi
subroutine psb_cinsvi(m, irw,val, x,desc_a,info,dupl) subroutine psb_cinsvi(m, irw,val, x,desc_a,info,dupl)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: m integer, intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
complex(psb_spk_),intent(inout) :: x(:) complex(psb_spk_),intent(inout) :: x(:)
@ -123,7 +123,7 @@ Module psb_c_tools_mod
interface psb_cdbldext interface psb_cdbldext
Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info,extype) Subroutine psb_ccdbldext(a,desc_a,novr,desc_ov,info,extype)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
integer, intent(in) :: novr integer, intent(in) :: novr
Type(psb_cspmat_type), Intent(in) :: a Type(psb_cspmat_type), Intent(in) :: a
@ -136,7 +136,7 @@ Module psb_c_tools_mod
interface psb_spall interface psb_spall
subroutine psb_cspalloc(a, desc_a, info, nnz) subroutine psb_cspalloc(a, desc_a, info, nnz)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_cspmat_type), intent(inout) :: a type(psb_cspmat_type), intent(inout) :: a
@ -147,7 +147,7 @@ Module psb_c_tools_mod
interface psb_spasb interface psb_spasb
subroutine psb_cspasb(a,desc_a, info, afmt, upd, dupl,mold) subroutine psb_cspasb(a,desc_a, info, afmt, upd, dupl,mold)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type, psb_c_base_sparse_mat use psb_mat_mod, only : psb_cspmat_type, psb_c_base_sparse_mat
type(psb_cspmat_type), intent (inout) :: a type(psb_cspmat_type), intent (inout) :: a
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -160,7 +160,7 @@ Module psb_c_tools_mod
interface psb_spfree interface psb_spfree
subroutine psb_cspfree(a, desc_a,info) subroutine psb_cspfree(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_cspmat_type), intent(inout) ::a type(psb_cspmat_type), intent(inout) ::a
@ -171,7 +171,7 @@ Module psb_c_tools_mod
interface psb_spins interface psb_spins
subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild) subroutine psb_cspins(nz,ia,ja,val,a,desc_a,info,rebuild)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_desc_type), intent(inout) :: desc_a type(psb_desc_type), intent(inout) :: desc_a
type(psb_cspmat_type), intent(inout) :: a type(psb_cspmat_type), intent(inout) :: a
@ -181,7 +181,7 @@ Module psb_c_tools_mod
logical, intent(in), optional :: rebuild logical, intent(in), optional :: rebuild
end subroutine psb_cspins end subroutine psb_cspins
subroutine psb_cspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info) subroutine psb_cspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_desc_type), intent(in) :: desc_ar type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac type(psb_desc_type), intent(inout) :: desc_ac
@ -195,7 +195,7 @@ Module psb_c_tools_mod
interface psb_sprn interface psb_sprn
subroutine psb_csprn(a, desc_a,info,clear) subroutine psb_csprn(a, desc_a,info,clear)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_cspmat_type use psb_mat_mod, only : psb_cspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_cspmat_type), intent(inout) :: a type(psb_cspmat_type), intent(inout) :: a

46
base/modules/psb_d_psblas_mod.f90
Unescape Escape View File

@ -33,14 +33,14 @@ module psb_d_psblas_mod
interface psb_gedot interface psb_gedot
function psb_ddotv(x, y, desc_a,info) function psb_ddotv(x, y, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) :: psb_ddotv real(psb_dpk_) :: psb_ddotv
real(psb_dpk_), intent(in) :: x(:), y(:) real(psb_dpk_), intent(in) :: x(:), y(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end function psb_ddotv end function psb_ddotv
function psb_ddot(x, y, desc_a, info, jx, jy) function psb_ddot(x, y, desc_a, info, jx, jy)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) :: psb_ddot real(psb_dpk_) :: psb_ddot
real(psb_dpk_), intent(in) :: x(:,:), y(:,:) real(psb_dpk_), intent(in) :: x(:,:), y(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -52,14 +52,14 @@ module psb_d_psblas_mod
interface psb_gedots interface psb_gedots
subroutine psb_ddotvs(res,x, y, desc_a, info) subroutine psb_ddotvs(res,x, y, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent(out) :: res real(psb_dpk_), intent(out) :: res
real(psb_dpk_), intent(in) :: x(:), y(:) real(psb_dpk_), intent(in) :: x(:), y(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_ddotvs end subroutine psb_ddotvs
subroutine psb_dmdots(res,x, y, desc_a,info) subroutine psb_dmdots(res,x, y, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent(out) :: res(:) real(psb_dpk_), intent(out) :: res(:)
real(psb_dpk_), intent(in) :: x(:,:), y(:,:) real(psb_dpk_), intent(in) :: x(:,:), y(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -70,7 +70,7 @@ module psb_d_psblas_mod
interface psb_geaxpby interface psb_geaxpby
subroutine psb_daxpbyv(alpha, x, beta, y,& subroutine psb_daxpbyv(alpha, x, beta, y,&
& desc_a, info) & desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (in) :: x(:) real(psb_dpk_), intent (in) :: x(:)
real(psb_dpk_), intent (inout) :: y(:) real(psb_dpk_), intent (inout) :: y(:)
real(psb_dpk_), intent (in) :: alpha, beta real(psb_dpk_), intent (in) :: alpha, beta
@ -79,7 +79,7 @@ module psb_d_psblas_mod
end subroutine psb_daxpbyv end subroutine psb_daxpbyv
subroutine psb_daxpby(alpha, x, beta, y,& subroutine psb_daxpby(alpha, x, beta, y,&
& desc_a, info, n, jx, jy) & desc_a, info, n, jx, jy)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (in) :: x(:,:) real(psb_dpk_), intent (in) :: x(:,:)
real(psb_dpk_), intent (inout) :: y(:,:) real(psb_dpk_), intent (inout) :: y(:,:)
real(psb_dpk_), intent (in) :: alpha, beta real(psb_dpk_), intent (in) :: alpha, beta
@ -91,7 +91,7 @@ module psb_d_psblas_mod
interface psb_geamax interface psb_geamax
function psb_damax(x, desc_a, info, jx) function psb_damax(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_damax real(psb_dpk_) psb_damax
real(psb_dpk_), intent (in) :: x(:,:) real(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -99,7 +99,7 @@ module psb_d_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_damax end function psb_damax
function psb_damaxv(x, desc_a,info) function psb_damaxv(x, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_damaxv real(psb_dpk_) psb_damaxv
real(psb_dpk_), intent (in) :: x(:) real(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -109,14 +109,14 @@ module psb_d_psblas_mod
interface psb_geamaxs interface psb_geamaxs
subroutine psb_damaxvs(res,x,desc_a,info) subroutine psb_damaxvs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res real(psb_dpk_), intent (out) :: res
real(psb_dpk_), intent (in) :: x(:) real(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_damaxvs end subroutine psb_damaxvs
subroutine psb_dmamaxs(res,x,desc_a,info,jx) subroutine psb_dmamaxs(res,x,desc_a,info,jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res(:) real(psb_dpk_), intent (out) :: res(:)
real(psb_dpk_), intent (in) :: x(:,:) real(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -127,7 +127,7 @@ module psb_d_psblas_mod
interface psb_geasum interface psb_geasum
function psb_dasum(x, desc_a, info, jx) function psb_dasum(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_dasum real(psb_dpk_) psb_dasum
real(psb_dpk_), intent (in) :: x(:,:) real(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -135,7 +135,7 @@ module psb_d_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_dasum end function psb_dasum
function psb_dasumv(x, desc_a, info) function psb_dasumv(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_dasumv real(psb_dpk_) psb_dasumv
real(psb_dpk_), intent (in) :: x(:) real(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -145,14 +145,14 @@ module psb_d_psblas_mod
interface psb_geasums interface psb_geasums
subroutine psb_dasumvs(res,x,desc_a,info) subroutine psb_dasumvs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res real(psb_dpk_), intent (out) :: res
real(psb_dpk_), intent (in) :: x(:) real(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_dasumvs end subroutine psb_dasumvs
subroutine psb_dmasum(res,x,desc_a,info) subroutine psb_dmasum(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res(:) real(psb_dpk_), intent (out) :: res(:)
real(psb_dpk_), intent (in) :: x(:,:) real(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -163,7 +163,7 @@ module psb_d_psblas_mod
interface psb_genrm2 interface psb_genrm2
function psb_dnrm2(x, desc_a, info, jx) function psb_dnrm2(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_dnrm2 real(psb_dpk_) psb_dnrm2
real(psb_dpk_), intent (in) :: x(:,:) real(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -171,7 +171,7 @@ module psb_d_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_dnrm2 end function psb_dnrm2
function psb_dnrm2v(x, desc_a, info) function psb_dnrm2v(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_dnrm2v real(psb_dpk_) psb_dnrm2v
real(psb_dpk_), intent (in) :: x(:) real(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -181,7 +181,7 @@ module psb_d_psblas_mod
interface psb_genrm2s interface psb_genrm2s
subroutine psb_dnrm2vs(res,x,desc_a,info) subroutine psb_dnrm2vs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res real(psb_dpk_), intent (out) :: res
real(psb_dpk_), intent (in) :: x(:) real(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -192,7 +192,7 @@ module psb_d_psblas_mod
interface psb_spnrmi interface psb_spnrmi
function psb_dnrmi(a, desc_a,info) function psb_dnrmi(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
real(psb_dpk_) :: psb_dnrmi real(psb_dpk_) :: psb_dnrmi
type(psb_dspmat_type), intent (in) :: a type(psb_dspmat_type), intent (in) :: a
@ -203,7 +203,7 @@ module psb_d_psblas_mod
interface psb_spnrm1 interface psb_spnrm1
function psb_dspnrm1(a, desc_a,info) function psb_dspnrm1(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
real(psb_dpk_) :: psb_dspnrm1 real(psb_dpk_) :: psb_dspnrm1
type(psb_dspmat_type), intent (in) :: a type(psb_dspmat_type), intent (in) :: a
@ -215,7 +215,7 @@ module psb_d_psblas_mod
interface psb_spmm interface psb_spmm
subroutine psb_dspmm(alpha, a, x, beta, y, desc_a, info,& subroutine psb_dspmm(alpha, a, x, beta, y, desc_a, info,&
&trans, k, jx, jy,work,doswap) &trans, k, jx, jy,work,doswap)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_dspmat_type), intent(in) :: a type(psb_dspmat_type), intent(in) :: a
real(psb_dpk_), intent(inout) :: x(:,:) real(psb_dpk_), intent(inout) :: x(:,:)
@ -230,7 +230,7 @@ module psb_d_psblas_mod
end subroutine psb_dspmm end subroutine psb_dspmm
subroutine psb_dspmv(alpha, a, x, beta, y,& subroutine psb_dspmv(alpha, a, x, beta, y,&
& desc_a, info, trans, work,doswap) & desc_a, info, trans, work,doswap)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_dspmat_type), intent(in) :: a type(psb_dspmat_type), intent(in) :: a
real(psb_dpk_), intent(inout) :: x(:) real(psb_dpk_), intent(inout) :: x(:)
@ -248,7 +248,7 @@ module psb_d_psblas_mod
subroutine psb_dspsm(alpha, t, x, beta, y,& subroutine psb_dspsm(alpha, t, x, beta, y,&
& desc_a, info, trans, scale, choice,& & desc_a, info, trans, scale, choice,&
& diag, n, jx, jy, work) & diag, n, jx, jy, work)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_dspmat_type), intent(in) :: t type(psb_dspmat_type), intent(in) :: t
real(psb_dpk_), intent(in) :: x(:,:) real(psb_dpk_), intent(in) :: x(:,:)
@ -265,7 +265,7 @@ module psb_d_psblas_mod
subroutine psb_dspsv(alpha, t, x, beta, y,& subroutine psb_dspsv(alpha, t, x, beta, y,&
& desc_a, info, trans, scale, choice,& & desc_a, info, trans, scale, choice,&
& diag, work) & diag, work)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_dspmat_type), intent(in) :: t type(psb_dspmat_type), intent(in) :: t
real(psb_dpk_), intent(in) :: x(:) real(psb_dpk_), intent(in) :: x(:)

32
base/modules/psb_d_tools_mod.f90
Unescape Escape View File

@ -34,7 +34,7 @@ Module psb_d_tools_mod
interface psb_geall interface psb_geall
subroutine psb_dalloc(x, desc_a, info, n, lb) subroutine psb_dalloc(x, desc_a, info, n, lb)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
implicit none implicit none
real(psb_dpk_), allocatable, intent(out) :: x(:,:) real(psb_dpk_), allocatable, intent(out) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -42,7 +42,7 @@ Module psb_d_tools_mod
integer, optional, intent(in) :: n, lb integer, optional, intent(in) :: n, lb
end subroutine psb_dalloc end subroutine psb_dalloc
subroutine psb_dallocv(x, desc_a,info,n) subroutine psb_dallocv(x, desc_a,info,n)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), allocatable, intent(out) :: x(:) real(psb_dpk_), allocatable, intent(out) :: x(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer,intent(out) :: info integer,intent(out) :: info
@ -53,13 +53,13 @@ Module psb_d_tools_mod
interface psb_geasb interface psb_geasb
subroutine psb_dasb(x, desc_a, info) subroutine psb_dasb(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
real(psb_dpk_), allocatable, intent(inout) :: x(:,:) real(psb_dpk_), allocatable, intent(inout) :: x(:,:)
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_dasb end subroutine psb_dasb
subroutine psb_dasbv(x, desc_a, info) subroutine psb_dasbv(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
real(psb_dpk_), allocatable, intent(inout) :: x(:) real(psb_dpk_), allocatable, intent(inout) :: x(:)
integer, intent(out) :: info integer, intent(out) :: info
@ -69,7 +69,7 @@ Module psb_d_tools_mod
interface psb_sphalo interface psb_sphalo
Subroutine psb_dsphalo(a,desc_a,blk,info,rowcnv,colcnv,& Subroutine psb_dsphalo(a,desc_a,blk,info,rowcnv,colcnv,&
& rowscale,colscale,outfmt,data) & rowscale,colscale,outfmt,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
Type(psb_dspmat_type),Intent(in) :: a Type(psb_dspmat_type),Intent(in) :: a
Type(psb_dspmat_type),Intent(inout) :: blk Type(psb_dspmat_type),Intent(inout) :: blk
@ -83,13 +83,13 @@ Module psb_d_tools_mod
interface psb_gefree interface psb_gefree
subroutine psb_dfree(x, desc_a, info) subroutine psb_dfree(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_),allocatable, intent(inout) :: x(:,:) real(psb_dpk_),allocatable, intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_dfree end subroutine psb_dfree
subroutine psb_dfreev(x, desc_a, info) subroutine psb_dfreev(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_),allocatable, intent(inout) :: x(:) real(psb_dpk_),allocatable, intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
@ -98,7 +98,7 @@ Module psb_d_tools_mod
interface psb_geins interface psb_geins
subroutine psb_dinsi(m,irw,val, x,desc_a,info,dupl) subroutine psb_dinsi(m,irw,val, x,desc_a,info,dupl)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: m integer, intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
real(psb_dpk_),intent(inout) :: x(:,:) real(psb_dpk_),intent(inout) :: x(:,:)
@ -108,7 +108,7 @@ Module psb_d_tools_mod
integer, optional, intent(in) :: dupl integer, optional, intent(in) :: dupl
end subroutine psb_dinsi end subroutine psb_dinsi
subroutine psb_dinsvi(m,irw,val,x,desc_a,info,dupl) subroutine psb_dinsvi(m,irw,val,x,desc_a,info,dupl)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: m integer, intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
real(psb_dpk_),intent(inout) :: x(:) real(psb_dpk_),intent(inout) :: x(:)
@ -122,7 +122,7 @@ Module psb_d_tools_mod
interface psb_cdbldext interface psb_cdbldext
Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info,extype) Subroutine psb_dcdbldext(a,desc_a,novr,desc_ov,info,extype)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
Use psb_mat_mod, only : psb_dspmat_type Use psb_mat_mod, only : psb_dspmat_type
integer, intent(in) :: novr integer, intent(in) :: novr
Type(psb_dspmat_type), Intent(in) :: a Type(psb_dspmat_type), Intent(in) :: a
@ -135,7 +135,7 @@ Module psb_d_tools_mod
interface psb_spall interface psb_spall
subroutine psb_dspalloc(a, desc_a, info, nnz) subroutine psb_dspalloc(a, desc_a, info, nnz)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_dspmat_type), intent(inout) :: a type(psb_dspmat_type), intent(inout) :: a
@ -146,7 +146,7 @@ Module psb_d_tools_mod
interface psb_spasb interface psb_spasb
subroutine psb_dspasb(a,desc_a, info, afmt, upd, dupl,mold) subroutine psb_dspasb(a,desc_a, info, afmt, upd, dupl,mold)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type, psb_d_base_sparse_mat use psb_mat_mod, only : psb_dspmat_type, psb_d_base_sparse_mat
type(psb_dspmat_type), intent (inout) :: a type(psb_dspmat_type), intent (inout) :: a
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -159,7 +159,7 @@ Module psb_d_tools_mod
interface psb_spfree interface psb_spfree
subroutine psb_dspfree(a, desc_a,info) subroutine psb_dspfree(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_dspmat_type), intent(inout) :: a type(psb_dspmat_type), intent(inout) :: a
@ -170,7 +170,7 @@ Module psb_d_tools_mod
interface psb_spins interface psb_spins
subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild) subroutine psb_dspins(nz,ia,ja,val,a,desc_a,info,rebuild)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_desc_type), intent(inout) :: desc_a type(psb_desc_type), intent(inout) :: desc_a
type(psb_dspmat_type), intent(inout) :: a type(psb_dspmat_type), intent(inout) :: a
@ -180,7 +180,7 @@ Module psb_d_tools_mod
logical, intent(in), optional :: rebuild logical, intent(in), optional :: rebuild
end subroutine psb_dspins end subroutine psb_dspins
subroutine psb_dspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info) subroutine psb_dspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_dspmat_type), intent(inout) :: a type(psb_dspmat_type), intent(inout) :: a
type(psb_desc_type), intent(in) :: desc_ar type(psb_desc_type), intent(in) :: desc_ar
@ -194,7 +194,7 @@ Module psb_d_tools_mod
interface psb_sprn interface psb_sprn
subroutine psb_dsprn(a, desc_a,info,clear) subroutine psb_dsprn(a, desc_a,info,clear)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_dspmat_type use psb_mat_mod, only : psb_dspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_dspmat_type), intent(inout) :: a type(psb_dspmat_type), intent(inout) :: a

44
base/modules/psb_s_psblas_mod.f90
Unescape Escape View File

@ -33,14 +33,14 @@ module psb_s_psblas_mod
interface psb_gedot interface psb_gedot
function psb_sdotv(x, y, desc_a,info) function psb_sdotv(x, y, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) :: psb_sdotv real(psb_spk_) :: psb_sdotv
real(psb_spk_), intent(in) :: x(:), y(:) real(psb_spk_), intent(in) :: x(:), y(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end function psb_sdotv end function psb_sdotv
function psb_sdot(x, y, desc_a, info, jx, jy) function psb_sdot(x, y, desc_a, info, jx, jy)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) :: psb_sdot real(psb_spk_) :: psb_sdot
real(psb_spk_), intent(in) :: x(:,:), y(:,:) real(psb_spk_), intent(in) :: x(:,:), y(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -52,14 +52,14 @@ module psb_s_psblas_mod
interface psb_gedots interface psb_gedots
subroutine psb_sdotvs(res,x, y, desc_a, info) subroutine psb_sdotvs(res,x, y, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent(out) :: res real(psb_spk_), intent(out) :: res
real(psb_spk_), intent(in) :: x(:), y(:) real(psb_spk_), intent(in) :: x(:), y(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_sdotvs end subroutine psb_sdotvs
subroutine psb_smdots(res,x, y, desc_a,info) subroutine psb_smdots(res,x, y, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent(out) :: res(:) real(psb_spk_), intent(out) :: res(:)
real(psb_spk_), intent(in) :: x(:,:), y(:,:) real(psb_spk_), intent(in) :: x(:,:), y(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -70,7 +70,7 @@ module psb_s_psblas_mod
interface psb_geaxpby interface psb_geaxpby
subroutine psb_saxpbyv(alpha, x, beta, y,& subroutine psb_saxpbyv(alpha, x, beta, y,&
& desc_a, info) & desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (in) :: x(:) real(psb_spk_), intent (in) :: x(:)
real(psb_spk_), intent (inout) :: y(:) real(psb_spk_), intent (inout) :: y(:)
real(psb_spk_), intent (in) :: alpha, beta real(psb_spk_), intent (in) :: alpha, beta
@ -79,7 +79,7 @@ module psb_s_psblas_mod
end subroutine psb_saxpbyv end subroutine psb_saxpbyv
subroutine psb_saxpby(alpha, x, beta, y,& subroutine psb_saxpby(alpha, x, beta, y,&
& desc_a, info, n, jx, jy) & desc_a, info, n, jx, jy)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (in) :: x(:,:) real(psb_spk_), intent (in) :: x(:,:)
real(psb_spk_), intent (inout) :: y(:,:) real(psb_spk_), intent (inout) :: y(:,:)
real(psb_spk_), intent (in) :: alpha, beta real(psb_spk_), intent (in) :: alpha, beta
@ -91,7 +91,7 @@ module psb_s_psblas_mod
interface psb_geamax interface psb_geamax
function psb_samax(x, desc_a, info, jx) function psb_samax(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_samax real(psb_spk_) psb_samax
real(psb_spk_), intent (in) :: x(:,:) real(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -99,7 +99,7 @@ module psb_s_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_samax end function psb_samax
function psb_samaxv(x, desc_a,info) function psb_samaxv(x, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_samaxv real(psb_spk_) psb_samaxv
real(psb_spk_), intent (in) :: x(:) real(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -109,14 +109,14 @@ module psb_s_psblas_mod
interface psb_geamaxs interface psb_geamaxs
subroutine psb_samaxvs(res,x,desc_a,info) subroutine psb_samaxvs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res real(psb_spk_), intent (out) :: res
real(psb_spk_), intent (in) :: x(:) real(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_samaxvs end subroutine psb_samaxvs
subroutine psb_smamaxs(res,x,desc_a,info,jx) subroutine psb_smamaxs(res,x,desc_a,info,jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res(:) real(psb_spk_), intent (out) :: res(:)
real(psb_spk_), intent (in) :: x(:,:) real(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -127,7 +127,7 @@ module psb_s_psblas_mod
interface psb_geasum interface psb_geasum
function psb_sasum(x, desc_a, info, jx) function psb_sasum(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_sasum real(psb_spk_) psb_sasum
real(psb_spk_), intent (in) :: x(:,:) real(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -135,7 +135,7 @@ module psb_s_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_sasum end function psb_sasum
function psb_sasumv(x, desc_a, info) function psb_sasumv(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_sasumv real(psb_spk_) psb_sasumv
real(psb_spk_), intent (in) :: x(:) real(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -145,14 +145,14 @@ module psb_s_psblas_mod
interface psb_geasums interface psb_geasums
subroutine psb_sasumvs(res,x,desc_a,info) subroutine psb_sasumvs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res real(psb_spk_), intent (out) :: res
real(psb_spk_), intent (in) :: x(:) real(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_sasumvs end subroutine psb_sasumvs
subroutine psb_smasum(res,x,desc_a,info) subroutine psb_smasum(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res(:) real(psb_spk_), intent (out) :: res(:)
real(psb_spk_), intent (in) :: x(:,:) real(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -163,7 +163,7 @@ module psb_s_psblas_mod
interface psb_genrm2 interface psb_genrm2
function psb_snrm2(x, desc_a, info, jx) function psb_snrm2(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_snrm2 real(psb_spk_) psb_snrm2
real(psb_spk_), intent (in) :: x(:,:) real(psb_spk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -171,7 +171,7 @@ module psb_s_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_snrm2 end function psb_snrm2
function psb_snrm2v(x, desc_a, info) function psb_snrm2v(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_) psb_snrm2v real(psb_spk_) psb_snrm2v
real(psb_spk_), intent (in) :: x(:) real(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -181,7 +181,7 @@ module psb_s_psblas_mod
interface psb_genrm2s interface psb_genrm2s
subroutine psb_snrm2vs(res,x,desc_a,info) subroutine psb_snrm2vs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), intent (out) :: res real(psb_spk_), intent (out) :: res
real(psb_spk_), intent (in) :: x(:) real(psb_spk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -192,7 +192,7 @@ module psb_s_psblas_mod
interface psb_spnrmi interface psb_spnrmi
function psb_snrmi(a, desc_a,info) function psb_snrmi(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
real(psb_spk_) :: psb_snrmi real(psb_spk_) :: psb_snrmi
type(psb_sspmat_type), intent (in) :: a type(psb_sspmat_type), intent (in) :: a
@ -204,7 +204,7 @@ module psb_s_psblas_mod
interface psb_spmm interface psb_spmm
subroutine psb_sspmm(alpha, a, x, beta, y, desc_a, info,& subroutine psb_sspmm(alpha, a, x, beta, y, desc_a, info,&
&trans, k, jx, jy,work,doswap) &trans, k, jx, jy,work,doswap)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_sspmat_type), intent(in) :: a type(psb_sspmat_type), intent(in) :: a
real(psb_spk_), intent(inout) :: x(:,:) real(psb_spk_), intent(inout) :: x(:,:)
@ -219,7 +219,7 @@ module psb_s_psblas_mod
end subroutine psb_sspmm end subroutine psb_sspmm
subroutine psb_sspmv(alpha, a, x, beta, y,& subroutine psb_sspmv(alpha, a, x, beta, y,&
& desc_a, info, trans, work,doswap) & desc_a, info, trans, work,doswap)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_sspmat_type), intent(in) :: a type(psb_sspmat_type), intent(in) :: a
real(psb_spk_), intent(inout) :: x(:) real(psb_spk_), intent(inout) :: x(:)
@ -237,7 +237,7 @@ module psb_s_psblas_mod
subroutine psb_sspsm(alpha, t, x, beta, y,& subroutine psb_sspsm(alpha, t, x, beta, y,&
& desc_a, info, trans, scale, choice,& & desc_a, info, trans, scale, choice,&
& diag, n, jx, jy, work) & diag, n, jx, jy, work)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_sspmat_type), intent(in) :: t type(psb_sspmat_type), intent(in) :: t
real(psb_spk_), intent(in) :: x(:,:) real(psb_spk_), intent(in) :: x(:,:)
@ -254,7 +254,7 @@ module psb_s_psblas_mod
subroutine psb_sspsv(alpha, t, x, beta, y,& subroutine psb_sspsv(alpha, t, x, beta, y,&
& desc_a, info, trans, scale, choice,& & desc_a, info, trans, scale, choice,&
& diag, work) & diag, work)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_sspmat_type), intent(in) :: t type(psb_sspmat_type), intent(in) :: t
real(psb_spk_), intent(in) :: x(:) real(psb_spk_), intent(in) :: x(:)

32
base/modules/psb_s_tools_mod.f90
Unescape Escape View File

@ -34,7 +34,7 @@ Module psb_s_tools_mod
interface psb_geall interface psb_geall
subroutine psb_salloc(x, desc_a, info, n, lb) subroutine psb_salloc(x, desc_a, info, n, lb)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
implicit none implicit none
real(psb_spk_), allocatable, intent(out) :: x(:,:) real(psb_spk_), allocatable, intent(out) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -42,7 +42,7 @@ Module psb_s_tools_mod
integer, optional, intent(in) :: n, lb integer, optional, intent(in) :: n, lb
end subroutine psb_salloc end subroutine psb_salloc
subroutine psb_sallocv(x, desc_a,info,n) subroutine psb_sallocv(x, desc_a,info,n)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_), allocatable, intent(out) :: x(:) real(psb_spk_), allocatable, intent(out) :: x(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer,intent(out) :: info integer,intent(out) :: info
@ -53,13 +53,13 @@ Module psb_s_tools_mod
interface psb_geasb interface psb_geasb
subroutine psb_sasb(x, desc_a, info) subroutine psb_sasb(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
real(psb_spk_), allocatable, intent(inout) :: x(:,:) real(psb_spk_), allocatable, intent(inout) :: x(:,:)
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_sasb end subroutine psb_sasb
subroutine psb_sasbv(x, desc_a, info) subroutine psb_sasbv(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
real(psb_spk_), allocatable, intent(inout) :: x(:) real(psb_spk_), allocatable, intent(inout) :: x(:)
integer, intent(out) :: info integer, intent(out) :: info
@ -69,7 +69,7 @@ Module psb_s_tools_mod
interface psb_sphalo interface psb_sphalo
Subroutine psb_ssphalo(a,desc_a,blk,info,rowcnv,colcnv,& Subroutine psb_ssphalo(a,desc_a,blk,info,rowcnv,colcnv,&
& rowscale,colscale,outfmt,data) & rowscale,colscale,outfmt,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
Type(psb_sspmat_type),Intent(in) :: a Type(psb_sspmat_type),Intent(in) :: a
Type(psb_sspmat_type),Intent(inout) :: blk Type(psb_sspmat_type),Intent(inout) :: blk
@ -84,13 +84,13 @@ Module psb_s_tools_mod
interface psb_gefree interface psb_gefree
subroutine psb_sfree(x, desc_a, info) subroutine psb_sfree(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_),allocatable, intent(inout) :: x(:,:) real(psb_spk_),allocatable, intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_sfree end subroutine psb_sfree
subroutine psb_sfreev(x, desc_a, info) subroutine psb_sfreev(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
real(psb_spk_),allocatable, intent(inout) :: x(:) real(psb_spk_),allocatable, intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
@ -100,7 +100,7 @@ Module psb_s_tools_mod
interface psb_geins interface psb_geins
subroutine psb_sinsi(m,irw,val, x,desc_a,info,dupl) subroutine psb_sinsi(m,irw,val, x,desc_a,info,dupl)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: m integer, intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
real(psb_spk_),intent(inout) :: x(:,:) real(psb_spk_),intent(inout) :: x(:,:)
@ -110,7 +110,7 @@ Module psb_s_tools_mod
integer, optional, intent(in) :: dupl integer, optional, intent(in) :: dupl
end subroutine psb_sinsi end subroutine psb_sinsi
subroutine psb_sinsvi(m,irw,val,x,desc_a,info,dupl) subroutine psb_sinsvi(m,irw,val,x,desc_a,info,dupl)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: m integer, intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
real(psb_spk_),intent(inout) :: x(:) real(psb_spk_),intent(inout) :: x(:)
@ -123,7 +123,7 @@ Module psb_s_tools_mod
interface psb_cdbldext interface psb_cdbldext
Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info,extype) Subroutine psb_scdbldext(a,desc_a,novr,desc_ov,info,extype)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
Use psb_mat_mod, only : psb_sspmat_type Use psb_mat_mod, only : psb_sspmat_type
integer, intent(in) :: novr integer, intent(in) :: novr
Type(psb_sspmat_type), Intent(in) :: a Type(psb_sspmat_type), Intent(in) :: a
@ -136,7 +136,7 @@ Module psb_s_tools_mod
interface psb_spall interface psb_spall
subroutine psb_sspalloc(a, desc_a, info, nnz) subroutine psb_sspalloc(a, desc_a, info, nnz)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_sspmat_type), intent(inout) :: a type(psb_sspmat_type), intent(inout) :: a
@ -147,7 +147,7 @@ Module psb_s_tools_mod
interface psb_spasb interface psb_spasb
subroutine psb_sspasb(a,desc_a, info, afmt, upd, dupl, mold) subroutine psb_sspasb(a,desc_a, info, afmt, upd, dupl, mold)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_s_base_sparse_mat, psb_sspmat_type use psb_mat_mod, only : psb_s_base_sparse_mat, psb_sspmat_type
type(psb_sspmat_type), intent (inout) :: a type(psb_sspmat_type), intent (inout) :: a
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -160,7 +160,7 @@ Module psb_s_tools_mod
interface psb_spfree interface psb_spfree
subroutine psb_sspfree(a, desc_a,info) subroutine psb_sspfree(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_sspmat_type), intent(inout) :: a type(psb_sspmat_type), intent(inout) :: a
@ -171,7 +171,7 @@ Module psb_s_tools_mod
interface psb_spins interface psb_spins
subroutine psb_sspins(nz,ia,ja,val,a,desc_a,info,rebuild) subroutine psb_sspins(nz,ia,ja,val,a,desc_a,info,rebuild)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_desc_type), intent(inout) :: desc_a type(psb_desc_type), intent(inout) :: desc_a
type(psb_sspmat_type), intent(inout) :: a type(psb_sspmat_type), intent(inout) :: a
@ -181,7 +181,7 @@ Module psb_s_tools_mod
logical, intent(in), optional :: rebuild logical, intent(in), optional :: rebuild
end subroutine psb_sspins end subroutine psb_sspins
subroutine psb_sspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info) subroutine psb_sspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_desc_type), intent(in) :: desc_ar type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac type(psb_desc_type), intent(inout) :: desc_ac
@ -195,7 +195,7 @@ Module psb_s_tools_mod
interface psb_sprn interface psb_sprn
subroutine psb_ssprn(a, desc_a,info,clear) subroutine psb_ssprn(a, desc_a,info,clear)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
use psb_mat_mod, only : psb_sspmat_type use psb_mat_mod, only : psb_sspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_sspmat_type), intent(inout) :: a type(psb_sspmat_type), intent(inout) :: a

32
base/modules/psb_serial_mod.f90
Unescape Escape View File

@ -229,97 +229,97 @@ module psb_serial_mod
interface psb_geprt interface psb_geprt
subroutine psb_sgeprtn2(fname,a,head) subroutine psb_sgeprtn2(fname,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_spk_
character(len=*), intent(in) :: fname character(len=*), intent(in) :: fname
real(psb_spk_), intent(in) :: a(:,:) real(psb_spk_), intent(in) :: a(:,:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_sgeprtn2 end subroutine psb_sgeprtn2
subroutine psb_sgeprtn1(fname,a,head) subroutine psb_sgeprtn1(fname,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_spk_
character(len=*), intent(in) :: fname character(len=*), intent(in) :: fname
real(psb_spk_), intent(in) :: a(:) real(psb_spk_), intent(in) :: a(:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_sgeprtn1 end subroutine psb_sgeprtn1
subroutine psb_sgeprt2(iout,a,head) subroutine psb_sgeprt2(iout,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_spk_
integer, intent(in) :: iout integer, intent(in) :: iout
real(psb_spk_), intent(in) :: a(:,:) real(psb_spk_), intent(in) :: a(:,:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_sgeprt2 end subroutine psb_sgeprt2
subroutine psb_sgeprt1(iout,a,head) subroutine psb_sgeprt1(iout,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_spk_
integer, intent(in) :: iout integer, intent(in) :: iout
real(psb_spk_), intent(in) :: a(:) real(psb_spk_), intent(in) :: a(:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_sgeprt1 end subroutine psb_sgeprt1
subroutine psb_dgeprtn2(fname,a,head) subroutine psb_dgeprtn2(fname,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_dpk_
character(len=*), intent(in) :: fname character(len=*), intent(in) :: fname
real(psb_dpk_), intent(in) :: a(:,:) real(psb_dpk_), intent(in) :: a(:,:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_dgeprtn2 end subroutine psb_dgeprtn2
subroutine psb_dgeprtn1(fname,a,head) subroutine psb_dgeprtn1(fname,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_dpk_
character(len=*), intent(in) :: fname character(len=*), intent(in) :: fname
real(psb_dpk_), intent(in) :: a(:) real(psb_dpk_), intent(in) :: a(:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_dgeprtn1 end subroutine psb_dgeprtn1
subroutine psb_dgeprt2(iout,a,head) subroutine psb_dgeprt2(iout,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_dpk_
integer, intent(in) :: iout integer, intent(in) :: iout
real(psb_dpk_), intent(in) :: a(:,:) real(psb_dpk_), intent(in) :: a(:,:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_dgeprt2 end subroutine psb_dgeprt2
subroutine psb_dgeprt1(iout,a,head) subroutine psb_dgeprt1(iout,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_dpk_
integer, intent(in) :: iout integer, intent(in) :: iout
real(psb_dpk_), intent(in) :: a(:) real(psb_dpk_), intent(in) :: a(:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_dgeprt1 end subroutine psb_dgeprt1
subroutine psb_cgeprtn2(fname,a,head) subroutine psb_cgeprtn2(fname,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_spk_
character(len=*), intent(in) :: fname character(len=*), intent(in) :: fname
complex(psb_spk_), intent(in) :: a(:,:) complex(psb_spk_), intent(in) :: a(:,:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_cgeprtn2 end subroutine psb_cgeprtn2
subroutine psb_cgeprtn1(fname,a,head) subroutine psb_cgeprtn1(fname,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_spk_
character(len=*), intent(in) :: fname character(len=*), intent(in) :: fname
complex(psb_spk_), intent(in) :: a(:) complex(psb_spk_), intent(in) :: a(:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_cgeprtn1 end subroutine psb_cgeprtn1
subroutine psb_cgeprt2(iout,a,head) subroutine psb_cgeprt2(iout,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_spk_
integer, intent(in) :: iout integer, intent(in) :: iout
complex(psb_spk_), intent(in) :: a(:,:) complex(psb_spk_), intent(in) :: a(:,:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_cgeprt2 end subroutine psb_cgeprt2
subroutine psb_cgeprt1(iout,a,head) subroutine psb_cgeprt1(iout,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_spk_
integer, intent(in) :: iout integer, intent(in) :: iout
complex(psb_spk_), intent(in) :: a(:) complex(psb_spk_), intent(in) :: a(:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_cgeprt1 end subroutine psb_cgeprt1
subroutine psb_zgeprtn2(fname,a,head) subroutine psb_zgeprtn2(fname,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_dpk_
character(len=*), intent(in) :: fname character(len=*), intent(in) :: fname
complex(psb_dpk_), intent(in) :: a(:,:) complex(psb_dpk_), intent(in) :: a(:,:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_zgeprtn2 end subroutine psb_zgeprtn2
subroutine psb_zgeprtn1(fname,a,head) subroutine psb_zgeprtn1(fname,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_dpk_
character(len=*), intent(in) :: fname character(len=*), intent(in) :: fname
complex(psb_dpk_), intent(in) :: a(:) complex(psb_dpk_), intent(in) :: a(:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_zgeprtn1 end subroutine psb_zgeprtn1
subroutine psb_zgeprt2(iout,a,head) subroutine psb_zgeprt2(iout,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_dpk_
integer, intent(in) :: iout integer, intent(in) :: iout
complex(psb_dpk_), intent(in) :: a(:,:) complex(psb_dpk_), intent(in) :: a(:,:)
character(len=*), optional :: head character(len=*), optional :: head
end subroutine psb_zgeprt2 end subroutine psb_zgeprt2
subroutine psb_zgeprt1(iout,a,head) subroutine psb_zgeprt1(iout,a,head)
use psb_const_mod, only : psb_spk_, psb_dpk_ use psb_const_mod, only : psb_dpk_
integer, intent(in) :: iout integer, intent(in) :: iout
complex(psb_dpk_), intent(in) :: a(:) complex(psb_dpk_), intent(in) :: a(:)
character(len=*), optional :: head character(len=*), optional :: head

44
base/modules/psb_z_psblas_mod.f90
Unescape Escape View File

@ -33,14 +33,14 @@ module psb_z_psblas_mod
interface psb_gedot interface psb_gedot
function psb_zdotv(x, y, desc_a,info) function psb_zdotv(x, y, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_) :: psb_zdotv complex(psb_dpk_) :: psb_zdotv
complex(psb_dpk_), intent(in) :: x(:), y(:) complex(psb_dpk_), intent(in) :: x(:), y(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end function psb_zdotv end function psb_zdotv
function psb_zdot(x, y, desc_a, info, jx, jy) function psb_zdot(x, y, desc_a, info, jx, jy)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_) :: psb_zdot complex(psb_dpk_) :: psb_zdot
complex(psb_dpk_), intent(in) :: x(:,:), y(:,:) complex(psb_dpk_), intent(in) :: x(:,:), y(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -52,14 +52,14 @@ module psb_z_psblas_mod
interface psb_gedots interface psb_gedots
subroutine psb_zdotvs(res,x, y, desc_a, info) subroutine psb_zdotvs(res,x, y, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_), intent(out) :: res complex(psb_dpk_), intent(out) :: res
complex(psb_dpk_), intent(in) :: x(:), y(:) complex(psb_dpk_), intent(in) :: x(:), y(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_zdotvs end subroutine psb_zdotvs
subroutine psb_zmdots(res,x, y, desc_a,info) subroutine psb_zmdots(res,x, y, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_), intent(out) :: res(:) complex(psb_dpk_), intent(out) :: res(:)
complex(psb_dpk_), intent(in) :: x(:,:), y(:,:) complex(psb_dpk_), intent(in) :: x(:,:), y(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -70,7 +70,7 @@ module psb_z_psblas_mod
interface psb_geaxpby interface psb_geaxpby
subroutine psb_zaxpbyv(alpha, x, beta, y,& subroutine psb_zaxpbyv(alpha, x, beta, y,&
& desc_a, info) & desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_), intent (in) :: x(:) complex(psb_dpk_), intent (in) :: x(:)
complex(psb_dpk_), intent (inout) :: y(:) complex(psb_dpk_), intent (inout) :: y(:)
complex(psb_dpk_), intent (in) :: alpha, beta complex(psb_dpk_), intent (in) :: alpha, beta
@ -79,7 +79,7 @@ module psb_z_psblas_mod
end subroutine psb_zaxpbyv end subroutine psb_zaxpbyv
subroutine psb_zaxpby(alpha, x, beta, y,& subroutine psb_zaxpby(alpha, x, beta, y,&
& desc_a, info, n, jx, jy) & desc_a, info, n, jx, jy)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_), intent (in) :: x(:,:) complex(psb_dpk_), intent (in) :: x(:,:)
complex(psb_dpk_), intent (inout) :: y(:,:) complex(psb_dpk_), intent (inout) :: y(:,:)
complex(psb_dpk_), intent (in) :: alpha, beta complex(psb_dpk_), intent (in) :: alpha, beta
@ -91,7 +91,7 @@ module psb_z_psblas_mod
interface psb_geamax interface psb_geamax
function psb_zamax(x, desc_a, info, jx) function psb_zamax(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_zamax real(psb_dpk_) psb_zamax
complex(psb_dpk_), intent (in) :: x(:,:) complex(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -99,7 +99,7 @@ module psb_z_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_zamax end function psb_zamax
function psb_zamaxv(x, desc_a,info) function psb_zamaxv(x, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_zamaxv real(psb_dpk_) psb_zamaxv
complex(psb_dpk_), intent (in) :: x(:) complex(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -109,14 +109,14 @@ module psb_z_psblas_mod
interface psb_geamaxs interface psb_geamaxs
subroutine psb_zamaxvs(res,x,desc_a,info) subroutine psb_zamaxvs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res real(psb_dpk_), intent (out) :: res
complex(psb_dpk_), intent (in) :: x(:) complex(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_zamaxvs end subroutine psb_zamaxvs
subroutine psb_zmamaxs(res,x,desc_a,info,jx) subroutine psb_zmamaxs(res,x,desc_a,info,jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res(:) real(psb_dpk_), intent (out) :: res(:)
complex(psb_dpk_), intent (in) :: x(:,:) complex(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -127,7 +127,7 @@ module psb_z_psblas_mod
interface psb_geasum interface psb_geasum
function psb_zasum(x, desc_a, info, jx) function psb_zasum(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_zasum real(psb_dpk_) psb_zasum
complex(psb_dpk_), intent (in) :: x(:,:) complex(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -135,7 +135,7 @@ module psb_z_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_zasum end function psb_zasum
function psb_zasumv(x, desc_a, info) function psb_zasumv(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_zasumv real(psb_dpk_) psb_zasumv
complex(psb_dpk_), intent (in) :: x(:) complex(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -145,14 +145,14 @@ module psb_z_psblas_mod
interface psb_geasums interface psb_geasums
subroutine psb_zasumvs(res,x,desc_a,info) subroutine psb_zasumvs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res real(psb_dpk_), intent (out) :: res
complex(psb_dpk_), intent (in) :: x(:) complex(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_zasumvs end subroutine psb_zasumvs
subroutine psb_zmasum(res,x,desc_a,info) subroutine psb_zmasum(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res(:) real(psb_dpk_), intent (out) :: res(:)
complex(psb_dpk_), intent (in) :: x(:,:) complex(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -163,7 +163,7 @@ module psb_z_psblas_mod
interface psb_genrm2 interface psb_genrm2
function psb_znrm2(x, desc_a, info, jx) function psb_znrm2(x, desc_a, info, jx)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_znrm2 real(psb_dpk_) psb_znrm2
complex(psb_dpk_), intent (in) :: x(:,:) complex(psb_dpk_), intent (in) :: x(:,:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -171,7 +171,7 @@ module psb_z_psblas_mod
integer, intent(out) :: info integer, intent(out) :: info
end function psb_znrm2 end function psb_znrm2
function psb_znrm2v(x, desc_a, info) function psb_znrm2v(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_) psb_znrm2v real(psb_dpk_) psb_znrm2v
complex(psb_dpk_), intent (in) :: x(:) complex(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -181,7 +181,7 @@ module psb_z_psblas_mod
interface psb_genrm2s interface psb_genrm2s
subroutine psb_znrm2vs(res,x,desc_a,info) subroutine psb_znrm2vs(res,x,desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
real(psb_dpk_), intent (out) :: res real(psb_dpk_), intent (out) :: res
complex(psb_dpk_), intent (in) :: x(:) complex(psb_dpk_), intent (in) :: x(:)
type(psb_desc_type), intent (in) :: desc_a type(psb_desc_type), intent (in) :: desc_a
@ -192,7 +192,7 @@ module psb_z_psblas_mod
interface psb_spnrmi interface psb_spnrmi
function psb_znrmi(a, desc_a,info) function psb_znrmi(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
real(psb_dpk_) :: psb_znrmi real(psb_dpk_) :: psb_znrmi
type(psb_zspmat_type), intent (in) :: a type(psb_zspmat_type), intent (in) :: a
@ -204,7 +204,7 @@ module psb_z_psblas_mod
interface psb_spmm interface psb_spmm
subroutine psb_zspmm(alpha, a, x, beta, y, desc_a, info,& subroutine psb_zspmm(alpha, a, x, beta, y, desc_a, info,&
&trans, k, jx, jy,work,doswap) &trans, k, jx, jy,work,doswap)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_zspmat_type), intent(in) :: a type(psb_zspmat_type), intent(in) :: a
complex(psb_dpk_), intent(inout) :: x(:,:) complex(psb_dpk_), intent(inout) :: x(:,:)
@ -219,7 +219,7 @@ module psb_z_psblas_mod
end subroutine psb_zspmm end subroutine psb_zspmm
subroutine psb_zspmv(alpha, a, x, beta, y,& subroutine psb_zspmv(alpha, a, x, beta, y,&
& desc_a, info, trans, work,doswap) & desc_a, info, trans, work,doswap)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_zspmat_type), intent(in) :: a type(psb_zspmat_type), intent(in) :: a
complex(psb_dpk_), intent(inout) :: x(:) complex(psb_dpk_), intent(inout) :: x(:)
@ -237,7 +237,7 @@ module psb_z_psblas_mod
subroutine psb_zspsm(alpha, t, x, beta, y,& subroutine psb_zspsm(alpha, t, x, beta, y,&
& desc_a, info, trans, scale, choice,& & desc_a, info, trans, scale, choice,&
& diag, n, jx, jy, work) & diag, n, jx, jy, work)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_zspmat_type), intent(in) :: t type(psb_zspmat_type), intent(in) :: t
complex(psb_dpk_), intent(in) :: x(:,:) complex(psb_dpk_), intent(in) :: x(:,:)
@ -254,7 +254,7 @@ module psb_z_psblas_mod
subroutine psb_zspsv(alpha, t, x, beta, y,& subroutine psb_zspsv(alpha, t, x, beta, y,&
& desc_a, info, trans, scale, choice,& & desc_a, info, trans, scale, choice,&
& diag, work) & diag, work)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_zspmat_type), intent(in) :: t type(psb_zspmat_type), intent(in) :: t
complex(psb_dpk_), intent(in) :: x(:) complex(psb_dpk_), intent(in) :: x(:)

32
base/modules/psb_z_tools_mod.f90
Unescape Escape View File

@ -33,7 +33,7 @@ Module psb_z_tools_mod
interface psb_geall interface psb_geall
subroutine psb_zalloc(x, desc_a, info, n, lb) subroutine psb_zalloc(x, desc_a, info, n, lb)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
implicit none implicit none
complex(psb_dpk_), allocatable, intent(out) :: x(:,:) complex(psb_dpk_), allocatable, intent(out) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -41,7 +41,7 @@ Module psb_z_tools_mod
integer, optional, intent(in) :: n, lb integer, optional, intent(in) :: n, lb
end subroutine psb_zalloc end subroutine psb_zalloc
subroutine psb_zallocv(x, desc_a,info,n) subroutine psb_zallocv(x, desc_a,info,n)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_), allocatable, intent(out) :: x(:) complex(psb_dpk_), allocatable, intent(out) :: x(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
@ -52,13 +52,13 @@ Module psb_z_tools_mod
interface psb_geasb interface psb_geasb
subroutine psb_zasb(x, desc_a, info) subroutine psb_zasb(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
complex(psb_dpk_), allocatable, intent(inout) :: x(:,:) complex(psb_dpk_), allocatable, intent(inout) :: x(:,:)
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_zasb end subroutine psb_zasb
subroutine psb_zasbv(x, desc_a, info) subroutine psb_zasbv(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
complex(psb_dpk_), allocatable, intent(inout) :: x(:) complex(psb_dpk_), allocatable, intent(inout) :: x(:)
integer, intent(out) :: info integer, intent(out) :: info
@ -68,7 +68,7 @@ Module psb_z_tools_mod
interface psb_sphalo interface psb_sphalo
Subroutine psb_zsphalo(a,desc_a,blk,info,rowcnv,colcnv,& Subroutine psb_zsphalo(a,desc_a,blk,info,rowcnv,colcnv,&
& rowscale,colscale,outfmt,data) & rowscale,colscale,outfmt,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
Type(psb_zspmat_type),Intent(in) :: a Type(psb_zspmat_type),Intent(in) :: a
Type(psb_zspmat_type),Intent(inout) :: blk Type(psb_zspmat_type),Intent(inout) :: blk
@ -82,13 +82,13 @@ Module psb_z_tools_mod
interface psb_gefree interface psb_gefree
subroutine psb_zfree(x, desc_a, info) subroutine psb_zfree(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_),allocatable, intent(inout) :: x(:,:) complex(psb_dpk_),allocatable, intent(inout) :: x(:,:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psb_zfree end subroutine psb_zfree
subroutine psb_zfreev(x, desc_a, info) subroutine psb_zfreev(x, desc_a, info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
complex(psb_dpk_),allocatable, intent(inout) :: x(:) complex(psb_dpk_),allocatable, intent(inout) :: x(:)
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info integer, intent(out) :: info
@ -98,7 +98,7 @@ Module psb_z_tools_mod
interface psb_geins interface psb_geins
subroutine psb_zinsi(m,irw,val, x, desc_a,info,dupl) subroutine psb_zinsi(m,irw,val, x, desc_a,info,dupl)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: m integer, intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
complex(psb_dpk_),intent(inout) :: x(:,:) complex(psb_dpk_),intent(inout) :: x(:,:)
@ -108,7 +108,7 @@ Module psb_z_tools_mod
integer, optional, intent(in) :: dupl integer, optional, intent(in) :: dupl
end subroutine psb_zinsi end subroutine psb_zinsi
subroutine psb_zinsvi(m, irw,val, x,desc_a,info,dupl) subroutine psb_zinsvi(m, irw,val, x,desc_a,info,dupl)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: m integer, intent(in) :: m
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
complex(psb_dpk_),intent(inout) :: x(:) complex(psb_dpk_),intent(inout) :: x(:)
@ -123,7 +123,7 @@ Module psb_z_tools_mod
interface psb_cdbldext interface psb_cdbldext
Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info,extype) Subroutine psb_zcdbldext(a,desc_a,novr,desc_ov,info,extype)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
Use psb_mat_mod, only : psb_zspmat_type Use psb_mat_mod, only : psb_zspmat_type
integer, intent(in) :: novr integer, intent(in) :: novr
Type(psb_zspmat_type), Intent(in) :: a Type(psb_zspmat_type), Intent(in) :: a
@ -136,7 +136,7 @@ Module psb_z_tools_mod
interface psb_spall interface psb_spall
subroutine psb_zspalloc(a, desc_a, info, nnz) subroutine psb_zspalloc(a, desc_a, info, nnz)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_zspmat_type), intent(inout) :: a type(psb_zspmat_type), intent(inout) :: a
@ -147,7 +147,7 @@ Module psb_z_tools_mod
interface psb_spasb interface psb_spasb
subroutine psb_zspasb(a,desc_a, info, afmt, upd, dupl,mold) subroutine psb_zspasb(a,desc_a, info, afmt, upd, dupl,mold)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type, psb_z_base_sparse_mat use psb_mat_mod, only : psb_zspmat_type, psb_z_base_sparse_mat
type(psb_zspmat_type), intent (inout) :: a type(psb_zspmat_type), intent (inout) :: a
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
@ -160,7 +160,7 @@ Module psb_z_tools_mod
interface psb_spfree interface psb_spfree
subroutine psb_zspfree(a, desc_a,info) subroutine psb_zspfree(a, desc_a,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_zspmat_type), intent(inout) ::a type(psb_zspmat_type), intent(inout) ::a
@ -171,7 +171,7 @@ Module psb_z_tools_mod
interface psb_spins interface psb_spins
subroutine psb_zspins(nz,ia,ja,val,a,desc_a,info,rebuild) subroutine psb_zspins(nz,ia,ja,val,a,desc_a,info,rebuild)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_desc_type), intent(inout) :: desc_a type(psb_desc_type), intent(inout) :: desc_a
type(psb_zspmat_type), intent(inout) :: a type(psb_zspmat_type), intent(inout) :: a
@ -181,7 +181,7 @@ Module psb_z_tools_mod
logical, intent(in), optional :: rebuild logical, intent(in), optional :: rebuild
end subroutine psb_zspins end subroutine psb_zspins
subroutine psb_zspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info) subroutine psb_zspins_2desc(nz,ia,ja,val,a,desc_ar,desc_ac,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_desc_type), intent(in) :: desc_ar type(psb_desc_type), intent(in) :: desc_ar
type(psb_desc_type), intent(inout) :: desc_ac type(psb_desc_type), intent(inout) :: desc_ac
@ -195,7 +195,7 @@ Module psb_z_tools_mod
interface psb_sprn interface psb_sprn
subroutine psb_zsprn(a, desc_a,info,clear) subroutine psb_zsprn(a, desc_a,info,clear)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
use psb_mat_mod, only : psb_zspmat_type use psb_mat_mod, only : psb_zspmat_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
type(psb_zspmat_type), intent(inout) :: a type(psb_zspmat_type), intent(inout) :: a

60
base/modules/psi_mod.f90
Unescape Escape View File

@ -57,7 +57,7 @@ module psi_mod
interface interface
subroutine psi_crea_index(desc_a,index_in,index_out,glob_idx,nxch,nsnd,nrcv,info) subroutine psi_crea_index(desc_a,index_in,index_out,glob_idx,nxch,nsnd,nrcv,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
type(psb_desc_type), intent(in) :: desc_a type(psb_desc_type), intent(in) :: desc_a
integer, intent(out) :: info,nxch,nsnd,nrcv integer, intent(out) :: info,nxch,nsnd,nrcv
integer, intent(in) :: index_in(:) integer, intent(in) :: index_in(:)
@ -77,7 +77,7 @@ module psi_mod
interface interface
subroutine psi_desc_index(desc,index_in,dep_list,& subroutine psi_desc_index(desc,index_in,dep_list,&
& length_dl,nsnd,nrcv,desc_index,isglob_in,info) & length_dl,nsnd,nrcv,desc_index,isglob_in,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
type(psb_desc_type) :: desc type(psb_desc_type) :: desc
integer :: index_in(:),dep_list(:) integer :: index_in(:),dep_list(:)
integer,allocatable :: desc_index(:) integer,allocatable :: desc_index(:)
@ -101,7 +101,7 @@ module psi_mod
interface psi_swapdata interface psi_swapdata
subroutine psi_sswapdatam(flag,n,beta,y,desc_a,work,info,data) subroutine psi_sswapdatam(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
real(psb_spk_) :: y(:,:), beta real(psb_spk_) :: y(:,:), beta
@ -110,7 +110,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_sswapdatam end subroutine psi_sswapdatam
subroutine psi_sswapdatav(flag,beta,y,desc_a,work,info,data) subroutine psi_sswapdatav(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
real(psb_spk_) :: y(:), beta real(psb_spk_) :: y(:), beta
@ -135,7 +135,7 @@ module psi_mod
integer, intent(in) :: idx(:),totxch,totsnd,totrcv integer, intent(in) :: idx(:),totxch,totsnd,totrcv
end subroutine psi_sswapidxv end subroutine psi_sswapidxv
subroutine psi_dswapdatam(flag,n,beta,y,desc_a,work,info,data) subroutine psi_dswapdatam(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
real(psb_dpk_) :: y(:,:), beta real(psb_dpk_) :: y(:,:), beta
@ -144,7 +144,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_dswapdatam end subroutine psi_dswapdatam
subroutine psi_dswapdatav(flag,beta,y,desc_a,work,info,data) subroutine psi_dswapdatav(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
real(psb_dpk_) :: y(:), beta real(psb_dpk_) :: y(:), beta
@ -169,7 +169,7 @@ module psi_mod
integer, intent(in) :: idx(:),totxch,totsnd,totrcv integer, intent(in) :: idx(:),totxch,totsnd,totrcv
end subroutine psi_dswapidxv end subroutine psi_dswapidxv
subroutine psi_iswapdatam(flag,n,beta,y,desc_a,work,info,data) subroutine psi_iswapdatam(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
integer :: y(:,:), beta integer :: y(:,:), beta
@ -178,7 +178,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_iswapdatam end subroutine psi_iswapdatam
subroutine psi_iswapdatav(flag,beta,y,desc_a,work,info,data) subroutine psi_iswapdatav(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
integer :: y(:), beta integer :: y(:), beta
@ -203,7 +203,7 @@ module psi_mod
integer, intent(in) :: idx(:),totxch,totsnd,totrcv integer, intent(in) :: idx(:),totxch,totsnd,totrcv
end subroutine psi_iswapidxv end subroutine psi_iswapidxv
subroutine psi_cswapdatam(flag,n,beta,y,desc_a,work,info,data) subroutine psi_cswapdatam(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
complex(psb_spk_) :: y(:,:), beta complex(psb_spk_) :: y(:,:), beta
@ -212,7 +212,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_cswapdatam end subroutine psi_cswapdatam
subroutine psi_cswapdatav(flag,beta,y,desc_a,work,info,data) subroutine psi_cswapdatav(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
complex(psb_spk_) :: y(:), beta complex(psb_spk_) :: y(:), beta
@ -237,7 +237,7 @@ module psi_mod
integer, intent(in) :: idx(:),totxch,totsnd,totrcv integer, intent(in) :: idx(:),totxch,totsnd,totrcv
end subroutine psi_cswapidxv end subroutine psi_cswapidxv
subroutine psi_zswapdatam(flag,n,beta,y,desc_a,work,info,data) subroutine psi_zswapdatam(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
complex(psb_dpk_) :: y(:,:), beta complex(psb_dpk_) :: y(:,:), beta
@ -246,7 +246,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_zswapdatam end subroutine psi_zswapdatam
subroutine psi_zswapdatav(flag,beta,y,desc_a,work,info,data) subroutine psi_zswapdatav(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
complex(psb_dpk_) :: y(:), beta complex(psb_dpk_) :: y(:), beta
@ -275,7 +275,7 @@ module psi_mod
interface psi_swaptran interface psi_swaptran
subroutine psi_sswaptranm(flag,n,beta,y,desc_a,work,info,data) subroutine psi_sswaptranm(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
real(psb_spk_) :: y(:,:), beta real(psb_spk_) :: y(:,:), beta
@ -284,7 +284,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_sswaptranm end subroutine psi_sswaptranm
subroutine psi_sswaptranv(flag,beta,y,desc_a,work,info,data) subroutine psi_sswaptranv(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
real(psb_spk_) :: y(:), beta real(psb_spk_) :: y(:), beta
@ -309,7 +309,7 @@ module psi_mod
integer, intent(in) :: idx(:),totxch,totsnd,totrcv integer, intent(in) :: idx(:),totxch,totsnd,totrcv
end subroutine psi_stranidxv end subroutine psi_stranidxv
subroutine psi_dswaptranm(flag,n,beta,y,desc_a,work,info,data) subroutine psi_dswaptranm(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
real(psb_dpk_) :: y(:,:), beta real(psb_dpk_) :: y(:,:), beta
@ -318,7 +318,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_dswaptranm end subroutine psi_dswaptranm
subroutine psi_dswaptranv(flag,beta,y,desc_a,work,info,data) subroutine psi_dswaptranv(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
real(psb_dpk_) :: y(:), beta real(psb_dpk_) :: y(:), beta
@ -343,7 +343,7 @@ module psi_mod
integer, intent(in) :: idx(:),totxch,totsnd,totrcv integer, intent(in) :: idx(:),totxch,totsnd,totrcv
end subroutine psi_dtranidxv end subroutine psi_dtranidxv
subroutine psi_iswaptranm(flag,n,beta,y,desc_a,work,info,data) subroutine psi_iswaptranm(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
integer :: y(:,:), beta integer :: y(:,:), beta
@ -352,7 +352,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_iswaptranm end subroutine psi_iswaptranm
subroutine psi_iswaptranv(flag,beta,y,desc_a,work,info,data) subroutine psi_iswaptranv(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
integer :: y(:), beta integer :: y(:), beta
@ -377,7 +377,7 @@ module psi_mod
integer, intent(in) :: idx(:),totxch,totsnd,totrcv integer, intent(in) :: idx(:),totxch,totsnd,totrcv
end subroutine psi_itranidxv end subroutine psi_itranidxv
subroutine psi_cswaptranm(flag,n,beta,y,desc_a,work,info,data) subroutine psi_cswaptranm(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
complex(psb_spk_) :: y(:,:), beta complex(psb_spk_) :: y(:,:), beta
@ -386,7 +386,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_cswaptranm end subroutine psi_cswaptranm
subroutine psi_cswaptranv(flag,beta,y,desc_a,work,info,data) subroutine psi_cswaptranv(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_spk_
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
complex(psb_spk_) :: y(:), beta complex(psb_spk_) :: y(:), beta
@ -411,7 +411,7 @@ module psi_mod
integer, intent(in) :: idx(:),totxch,totsnd,totrcv integer, intent(in) :: idx(:),totxch,totsnd,totrcv
end subroutine psi_ctranidxv end subroutine psi_ctranidxv
subroutine psi_zswaptranm(flag,n,beta,y,desc_a,work,info,data) subroutine psi_zswaptranm(flag,n,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: flag, n integer, intent(in) :: flag, n
integer, intent(out) :: info integer, intent(out) :: info
complex(psb_dpk_) :: y(:,:), beta complex(psb_dpk_) :: y(:,:), beta
@ -420,7 +420,7 @@ module psi_mod
integer, optional :: data integer, optional :: data
end subroutine psi_zswaptranm end subroutine psi_zswaptranm
subroutine psi_zswaptranv(flag,beta,y,desc_a,work,info,data) subroutine psi_zswaptranv(flag,beta,y,desc_a,work,info,data)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type, psb_dpk_
integer, intent(in) :: flag integer, intent(in) :: flag
integer, intent(out) :: info integer, intent(out) :: info
complex(psb_dpk_) :: y(:), beta complex(psb_dpk_) :: y(:), beta
@ -457,7 +457,7 @@ module psi_mod
interface psi_fnd_owner interface psi_fnd_owner
subroutine psi_fnd_owner(nv,idx,iprc,desc,info) subroutine psi_fnd_owner(nv,idx,iprc,desc,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: nv integer, intent(in) :: nv
integer, intent(in) :: idx(:) integer, intent(in) :: idx(:)
integer, allocatable, intent(out) :: iprc(:) integer, allocatable, intent(out) :: iprc(:)
@ -468,7 +468,7 @@ module psi_mod
interface psi_ldsc_pre_halo interface psi_ldsc_pre_halo
subroutine psi_ldsc_pre_halo(desc,ext_hv,info) subroutine psi_ldsc_pre_halo(desc,ext_hv,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
type(psb_desc_type), intent(inout) :: desc type(psb_desc_type), intent(inout) :: desc
logical, intent(in) :: ext_hv logical, intent(in) :: ext_hv
integer, intent(out) :: info integer, intent(out) :: info
@ -477,7 +477,7 @@ module psi_mod
interface psi_bld_tmphalo interface psi_bld_tmphalo
subroutine psi_bld_tmphalo(desc,info) subroutine psi_bld_tmphalo(desc,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
type(psb_desc_type), intent(inout) :: desc type(psb_desc_type), intent(inout) :: desc
integer, intent(out) :: info integer, intent(out) :: info
end subroutine psi_bld_tmphalo end subroutine psi_bld_tmphalo
@ -486,7 +486,7 @@ module psi_mod
interface psi_bld_tmpovrl interface psi_bld_tmpovrl
subroutine psi_bld_tmpovrl(iv,desc,info) subroutine psi_bld_tmpovrl(iv,desc,info)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: iv(:) integer, intent(in) :: iv(:)
type(psb_desc_type), intent(inout) :: desc type(psb_desc_type), intent(inout) :: desc
integer, intent(out) :: info integer, intent(out) :: info
@ -534,7 +534,7 @@ module psi_mod
interface psi_idx_ins_cnv interface psi_idx_ins_cnv
subroutine psi_idx_ins_cnv1(nv,idxin,desc,info,mask) subroutine psi_idx_ins_cnv1(nv,idxin,desc,info,mask)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: nv integer, intent(in) :: nv
integer, intent(inout) :: idxin(:) integer, intent(inout) :: idxin(:)
type(psb_desc_type), intent(inout) :: desc type(psb_desc_type), intent(inout) :: desc
@ -542,7 +542,7 @@ module psi_mod
logical, intent(in), optional :: mask(:) logical, intent(in), optional :: mask(:)
end subroutine psi_idx_ins_cnv1 end subroutine psi_idx_ins_cnv1
subroutine psi_idx_ins_cnv2(nv,idxin,idxout,desc,info,mask) subroutine psi_idx_ins_cnv2(nv,idxin,idxout,desc,info,mask)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: nv, idxin(:) integer, intent(in) :: nv, idxin(:)
integer, intent(out) :: idxout(:) integer, intent(out) :: idxout(:)
type(psb_desc_type), intent(inout) :: desc type(psb_desc_type), intent(inout) :: desc
@ -550,7 +550,7 @@ module psi_mod
logical, intent(in), optional :: mask(:) logical, intent(in), optional :: mask(:)
end subroutine psi_idx_ins_cnv2 end subroutine psi_idx_ins_cnv2
subroutine psi_idx_ins_cnvs2(idxin,idxout,desc,info,mask) subroutine psi_idx_ins_cnvs2(idxin,idxout,desc,info,mask)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(in) :: idxin integer, intent(in) :: idxin
integer, intent(out) :: idxout integer, intent(out) :: idxout
type(psb_desc_type), intent(inout) :: desc type(psb_desc_type), intent(inout) :: desc
@ -558,7 +558,7 @@ module psi_mod
logical, intent(in), optional :: mask logical, intent(in), optional :: mask
end subroutine psi_idx_ins_cnvs2 end subroutine psi_idx_ins_cnvs2
subroutine psi_idx_ins_cnvs1(idxin,desc,info,mask) subroutine psi_idx_ins_cnvs1(idxin,desc,info,mask)
use psb_descriptor_type, only : psb_desc_type, psb_spk_, psb_dpk_ use psb_descriptor_type, only : psb_desc_type
integer, intent(inout) :: idxin integer, intent(inout) :: idxin
type(psb_desc_type), intent(inout) :: desc type(psb_desc_type), intent(inout) :: desc
integer, intent(out) :: info integer, intent(out) :: info

22
base/serial/impl/psb_c_csc_impl.f90
Unescape Escape View File

@ -862,7 +862,7 @@ subroutine psb_c_csc_cssv(alpha,a,x,beta,y,info,trans)
character, optional, intent(in) :: trans character, optional, intent(in) :: trans
character :: trans_ character :: trans_
integer :: i,j,k,m,n, nnz, ir, jc integer :: i,j,k,m, nnz, ir, jc
complex(psb_spk_) :: acc complex(psb_spk_) :: acc
complex(psb_spk_), allocatable :: tmp(:) complex(psb_spk_), allocatable :: tmp(:)
logical :: tra, ctra logical :: tra, ctra
@ -896,7 +896,7 @@ subroutine psb_c_csc_cssv(alpha,a,x,beta,y,info,trans)
if (size(x,1)<m) then if (size(x,1)<m) then
info = 36 info = 36
call psb_errpush(info,name,i_err=(/3,n,0,0,0/)) call psb_errpush(info,name,i_err=(/3,m,0,0,0/))
goto 9999 goto 9999
end if end if
@ -2258,7 +2258,7 @@ subroutine psb_c_cp_csc_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2285,7 +2285,7 @@ subroutine psb_c_cp_csc_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2328,7 +2328,7 @@ subroutine psb_c_mv_csc_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2372,7 +2372,7 @@ subroutine psb_c_mv_csc_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc, icl Integer :: nza, nr, i,j,irw, err_act, nc, icl
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2417,7 +2417,7 @@ subroutine psb_c_mv_csc_from_coo(a,b,info)
if (i >= icl) exit inner if (i >= icl) exit inner
if (i > nc) then if (i > nc) then
write(debug_unit,*) trim(name),& write(debug_unit,*) trim(name),&
& 'Strange situation: i>nr ',i,nc,j,nza,icl,idl & 'Strange situation: i>nr ',i,nc,j,nza,icl
exit outer exit outer
end if end if
a%icp(i+1) = a%icp(i) a%icp(i+1) = a%icp(i)
@ -2465,7 +2465,7 @@ subroutine psb_c_mv_csc_to_fmt(a,b,info)
!locals !locals
type(psb_c_coo_sparse_mat) :: tmp type(psb_c_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2505,7 +2505,7 @@ subroutine psb_c_cp_csc_to_fmt(a,b,info)
!locals !locals
type(psb_c_coo_sparse_mat) :: tmp type(psb_c_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2545,7 +2545,7 @@ subroutine psb_c_mv_csc_from_fmt(a,b,info)
!locals !locals
type(psb_c_coo_sparse_mat) :: tmp type(psb_c_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2586,7 +2586,7 @@ subroutine psb_c_cp_csc_from_fmt(a,b,info)
!locals !locals
type(psb_c_coo_sparse_mat) :: tmp type(psb_c_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name

18
base/serial/impl/psb_c_csr_impl.f90
Unescape Escape View File

@ -2422,7 +2422,7 @@ subroutine psb_c_cp_csr_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2449,7 +2449,7 @@ subroutine psb_c_cp_csr_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2491,7 +2491,7 @@ subroutine psb_c_mv_csr_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2536,7 +2536,7 @@ subroutine psb_c_mv_csr_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2581,7 +2581,7 @@ subroutine psb_c_mv_csr_from_coo(a,b,info)
if (i >= irw) exit inner if (i >= irw) exit inner
if (i>nr) then if (i>nr) then
write(debug_unit,*) trim(name),& write(debug_unit,*) trim(name),&
& 'Strange situation: i>nr ',i,nr,j,nza,irw,idl & 'Strange situation: i>nr ',i,nr,j,nza,irw
exit outer exit outer
end if end if
a%irp(i+1) = a%irp(i) a%irp(i+1) = a%irp(i)
@ -2628,7 +2628,7 @@ subroutine psb_c_mv_csr_to_fmt(a,b,info)
!locals !locals
type(psb_c_coo_sparse_mat) :: tmp type(psb_c_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2668,7 +2668,7 @@ subroutine psb_c_cp_csr_to_fmt(a,b,info)
!locals !locals
type(psb_c_coo_sparse_mat) :: tmp type(psb_c_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2707,7 +2707,7 @@ subroutine psb_c_mv_csr_from_fmt(a,b,info)
!locals !locals
type(psb_c_coo_sparse_mat) :: tmp type(psb_c_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2748,7 +2748,7 @@ subroutine psb_c_cp_csr_from_fmt(a,b,info)
!locals !locals
type(psb_c_coo_sparse_mat) :: tmp type(psb_c_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nz, nr, i,j,irw, idl,err_act, nc Integer :: nz, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name

22
base/serial/impl/psb_d_csc_impl.f90
Unescape Escape View File

@ -585,7 +585,7 @@ subroutine psb_d_csc_cssv(alpha,a,x,beta,y,info,trans)
character, optional, intent(in) :: trans character, optional, intent(in) :: trans
character :: trans_ character :: trans_
integer :: i,j,k,m,n, nnz, ir, jc integer :: i,j,k,m, nnz, ir, jc
real(psb_dpk_) :: acc real(psb_dpk_) :: acc
real(psb_dpk_), allocatable :: tmp(:) real(psb_dpk_), allocatable :: tmp(:)
logical :: tra logical :: tra
@ -617,7 +617,7 @@ subroutine psb_d_csc_cssv(alpha,a,x,beta,y,info,trans)
if (size(x,1)<m) then if (size(x,1)<m) then
info = 36 info = 36
call psb_errpush(info,name,i_err=(/3,n,0,0,0/)) call psb_errpush(info,name,i_err=(/3,m,0,0,0/))
goto 9999 goto 9999
end if end if
@ -2130,7 +2130,7 @@ subroutine psb_d_cp_csc_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2157,7 +2157,7 @@ subroutine psb_d_cp_csc_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2200,7 +2200,7 @@ subroutine psb_d_mv_csc_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2244,7 +2244,7 @@ subroutine psb_d_mv_csc_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc, icl Integer :: nza, nr, i,j,irw, err_act, nc, icl
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2289,7 +2289,7 @@ subroutine psb_d_mv_csc_from_coo(a,b,info)
if (i >= icl) exit inner if (i >= icl) exit inner
if (i > nc) then if (i > nc) then
write(debug_unit,*) trim(name),& write(debug_unit,*) trim(name),&
& 'Strange situation: i>nr ',i,nc,j,nza,icl,idl & 'Strange situation: i>nr ',i,nc,j,nza,icl
exit outer exit outer
end if end if
a%icp(i+1) = a%icp(i) a%icp(i+1) = a%icp(i)
@ -2337,7 +2337,7 @@ subroutine psb_d_mv_csc_to_fmt(a,b,info)
!locals !locals
type(psb_d_coo_sparse_mat) :: tmp type(psb_d_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2377,7 +2377,7 @@ subroutine psb_d_cp_csc_to_fmt(a,b,info)
!locals !locals
type(psb_d_coo_sparse_mat) :: tmp type(psb_d_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2417,7 +2417,7 @@ subroutine psb_d_mv_csc_from_fmt(a,b,info)
!locals !locals
type(psb_d_coo_sparse_mat) :: tmp type(psb_d_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2458,7 +2458,7 @@ subroutine psb_d_cp_csc_from_fmt(a,b,info)
!locals !locals
type(psb_d_coo_sparse_mat) :: tmp type(psb_d_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name

18
base/serial/impl/psb_d_csr_impl.f90
Unescape Escape View File

@ -2471,7 +2471,7 @@ subroutine psb_d_cp_csr_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2498,7 +2498,7 @@ subroutine psb_d_cp_csr_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2540,7 +2540,7 @@ subroutine psb_d_mv_csr_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2585,7 +2585,7 @@ subroutine psb_d_mv_csr_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2629,7 +2629,7 @@ subroutine psb_d_mv_csr_from_coo(a,b,info)
if (i >= irw) exit inner if (i >= irw) exit inner
if (i>nr) then if (i>nr) then
write(debug_unit,*) trim(name),& write(debug_unit,*) trim(name),&
& 'Strange situation: i>nr ',i,nr,j,nza,irw,idl & 'Strange situation: i>nr ',i,nr,j,nza,irw
exit outer exit outer
end if end if
a%irp(i+1) = a%irp(i) a%irp(i+1) = a%irp(i)
@ -2676,7 +2676,7 @@ subroutine psb_d_mv_csr_to_fmt(a,b,info)
!locals !locals
type(psb_d_coo_sparse_mat) :: tmp type(psb_d_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2716,7 +2716,7 @@ subroutine psb_d_cp_csr_to_fmt(a,b,info)
!locals !locals
type(psb_d_coo_sparse_mat) :: tmp type(psb_d_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2755,7 +2755,7 @@ subroutine psb_d_mv_csr_from_fmt(a,b,info)
!locals !locals
type(psb_d_coo_sparse_mat) :: tmp type(psb_d_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2796,7 +2796,7 @@ subroutine psb_d_cp_csr_from_fmt(a,b,info)
!locals !locals
type(psb_d_coo_sparse_mat) :: tmp type(psb_d_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nz, nr, i,j,irw, idl,err_act, nc Integer :: nz, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name

22
base/serial/impl/psb_s_csc_impl.f90
Unescape Escape View File

@ -585,7 +585,7 @@ subroutine psb_s_csc_cssv(alpha,a,x,beta,y,info,trans)
character, optional, intent(in) :: trans character, optional, intent(in) :: trans
character :: trans_ character :: trans_
integer :: i,j,k,m,n, nnz, ir, jc integer :: i,j,k,m, nnz, ir, jc
real(psb_spk_) :: acc real(psb_spk_) :: acc
real(psb_spk_), allocatable :: tmp(:) real(psb_spk_), allocatable :: tmp(:)
logical :: tra logical :: tra
@ -617,7 +617,7 @@ subroutine psb_s_csc_cssv(alpha,a,x,beta,y,info,trans)
if (size(x,1)<m) then if (size(x,1)<m) then
info = 36 info = 36
call psb_errpush(info,name,i_err=(/3,n,0,0,0/)) call psb_errpush(info,name,i_err=(/3,m,0,0,0/))
goto 9999 goto 9999
end if end if
@ -1895,7 +1895,7 @@ subroutine psb_s_cp_csc_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -1922,7 +1922,7 @@ subroutine psb_s_cp_csc_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -1965,7 +1965,7 @@ subroutine psb_s_mv_csc_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2009,7 +2009,7 @@ subroutine psb_s_mv_csc_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc, icl Integer :: nza, nr, i,j,irw, err_act, nc, icl
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2054,7 +2054,7 @@ subroutine psb_s_mv_csc_from_coo(a,b,info)
if (i >= icl) exit inner if (i >= icl) exit inner
if (i > nc) then if (i > nc) then
write(debug_unit,*) trim(name),& write(debug_unit,*) trim(name),&
& 'Strange situation: i>nr ',i,nc,j,nza,icl,idl & 'Strange situation: i>nr ',i,nc,j,nza,icl
exit outer exit outer
end if end if
a%icp(i+1) = a%icp(i) a%icp(i+1) = a%icp(i)
@ -2102,7 +2102,7 @@ subroutine psb_s_mv_csc_to_fmt(a,b,info)
!locals !locals
type(psb_s_coo_sparse_mat) :: tmp type(psb_s_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2142,7 +2142,7 @@ subroutine psb_s_cp_csc_to_fmt(a,b,info)
!locals !locals
type(psb_s_coo_sparse_mat) :: tmp type(psb_s_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2182,7 +2182,7 @@ subroutine psb_s_mv_csc_from_fmt(a,b,info)
!locals !locals
type(psb_s_coo_sparse_mat) :: tmp type(psb_s_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2223,7 +2223,7 @@ subroutine psb_s_cp_csc_from_fmt(a,b,info)
!locals !locals
type(psb_s_coo_sparse_mat) :: tmp type(psb_s_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name

18
base/serial/impl/psb_s_csr_impl.f90
Unescape Escape View File

@ -2229,7 +2229,7 @@ subroutine psb_s_cp_csr_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2256,7 +2256,7 @@ subroutine psb_s_cp_csr_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2298,7 +2298,7 @@ subroutine psb_s_mv_csr_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2343,7 +2343,7 @@ subroutine psb_s_mv_csr_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2388,7 +2388,7 @@ subroutine psb_s_mv_csr_from_coo(a,b,info)
if (i >= irw) exit inner if (i >= irw) exit inner
if (i>nr) then if (i>nr) then
write(debug_unit,*) trim(name),& write(debug_unit,*) trim(name),&
& 'Strange situation: i>nr ',i,nr,j,nza,irw,idl & 'Strange situation: i>nr ',i,nr,j,nza,irw
exit outer exit outer
end if end if
a%irp(i+1) = a%irp(i) a%irp(i+1) = a%irp(i)
@ -2435,7 +2435,7 @@ subroutine psb_s_mv_csr_to_fmt(a,b,info)
!locals !locals
type(psb_s_coo_sparse_mat) :: tmp type(psb_s_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2475,7 +2475,7 @@ subroutine psb_s_cp_csr_to_fmt(a,b,info)
!locals !locals
type(psb_s_coo_sparse_mat) :: tmp type(psb_s_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2514,7 +2514,7 @@ subroutine psb_s_mv_csr_from_fmt(a,b,info)
!locals !locals
type(psb_s_coo_sparse_mat) :: tmp type(psb_s_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2555,7 +2555,7 @@ subroutine psb_s_cp_csr_from_fmt(a,b,info)
!locals !locals
type(psb_s_coo_sparse_mat) :: tmp type(psb_s_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nz, nr, i,j,irw, idl,err_act, nc Integer :: nz, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name

22
base/serial/impl/psb_z_csc_impl.f90
Unescape Escape View File

@ -863,7 +863,7 @@ subroutine psb_z_csc_cssv(alpha,a,x,beta,y,info,trans)
character, optional, intent(in) :: trans character, optional, intent(in) :: trans
character :: trans_ character :: trans_
integer :: i,j,k,m,n, nnz, ir, jc integer :: i,j,k,m, nnz, ir, jc
complex(psb_dpk_) :: acc complex(psb_dpk_) :: acc
complex(psb_dpk_), allocatable :: tmp(:) complex(psb_dpk_), allocatable :: tmp(:)
logical :: tra, ctra logical :: tra, ctra
@ -897,7 +897,7 @@ subroutine psb_z_csc_cssv(alpha,a,x,beta,y,info,trans)
if (size(x,1)<m) then if (size(x,1)<m) then
info = 36 info = 36
call psb_errpush(info,name,i_err=(/3,n,0,0,0/)) call psb_errpush(info,name,i_err=(/3,m,0,0,0/))
goto 9999 goto 9999
end if end if
@ -2258,7 +2258,7 @@ subroutine psb_z_cp_csc_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2285,7 +2285,7 @@ subroutine psb_z_cp_csc_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2328,7 +2328,7 @@ subroutine psb_z_mv_csc_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2372,7 +2372,7 @@ subroutine psb_z_mv_csc_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc, icl Integer :: nza, nr, i,j,irw, err_act, nc, icl
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2417,7 +2417,7 @@ subroutine psb_z_mv_csc_from_coo(a,b,info)
if (i >= icl) exit inner if (i >= icl) exit inner
if (i > nc) then if (i > nc) then
write(debug_unit,*) trim(name),& write(debug_unit,*) trim(name),&
& 'Strange situation: i>nr ',i,nc,j,nza,icl,idl & 'Strange situation: i>nr ',i,nc,j,nza,icl
exit outer exit outer
end if end if
a%icp(i+1) = a%icp(i) a%icp(i+1) = a%icp(i)
@ -2465,7 +2465,7 @@ subroutine psb_z_mv_csc_to_fmt(a,b,info)
!locals !locals
type(psb_z_coo_sparse_mat) :: tmp type(psb_z_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2505,7 +2505,7 @@ subroutine psb_z_cp_csc_to_fmt(a,b,info)
!locals !locals
type(psb_z_coo_sparse_mat) :: tmp type(psb_z_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2545,7 +2545,7 @@ subroutine psb_z_mv_csc_from_fmt(a,b,info)
!locals !locals
type(psb_z_coo_sparse_mat) :: tmp type(psb_z_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2586,7 +2586,7 @@ subroutine psb_z_cp_csc_from_fmt(a,b,info)
!locals !locals
type(psb_z_coo_sparse_mat) :: tmp type(psb_z_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name

18
base/serial/impl/psb_z_csr_impl.f90
Unescape Escape View File

@ -2421,7 +2421,7 @@ subroutine psb_z_cp_csr_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2448,7 +2448,7 @@ subroutine psb_z_cp_csr_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2490,7 +2490,7 @@ subroutine psb_z_mv_csr_to_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, nc,i,j,irw, idl,err_act Integer :: nza, nr, nc,i,j,irw, err_act
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2535,7 +2535,7 @@ subroutine psb_z_mv_csr_from_coo(a,b,info)
integer, allocatable :: itemp(:) integer, allocatable :: itemp(:)
!locals !locals
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2580,7 +2580,7 @@ subroutine psb_z_mv_csr_from_coo(a,b,info)
if (i >= irw) exit inner if (i >= irw) exit inner
if (i>nr) then if (i>nr) then
write(debug_unit,*) trim(name),& write(debug_unit,*) trim(name),&
& 'Strange situation: i>nr ',i,nr,j,nza,irw,idl & 'Strange situation: i>nr ',i,nr,j,nza,irw
exit outer exit outer
end if end if
a%irp(i+1) = a%irp(i) a%irp(i+1) = a%irp(i)
@ -2627,7 +2627,7 @@ subroutine psb_z_mv_csr_to_fmt(a,b,info)
!locals !locals
type(psb_z_coo_sparse_mat) :: tmp type(psb_z_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2667,7 +2667,7 @@ subroutine psb_z_cp_csr_to_fmt(a,b,info)
!locals !locals
type(psb_z_coo_sparse_mat) :: tmp type(psb_z_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2706,7 +2706,7 @@ subroutine psb_z_mv_csr_from_fmt(a,b,info)
!locals !locals
type(psb_z_coo_sparse_mat) :: tmp type(psb_z_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nza, nr, i,j,irw, idl,err_act, nc Integer :: nza, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name
@ -2747,7 +2747,7 @@ subroutine psb_z_cp_csr_from_fmt(a,b,info)
!locals !locals
type(psb_z_coo_sparse_mat) :: tmp type(psb_z_coo_sparse_mat) :: tmp
logical :: rwshr_ logical :: rwshr_
Integer :: nz, nr, i,j,irw, idl,err_act, nc Integer :: nz, nr, i,j,irw, err_act, nc
Integer, Parameter :: maxtry=8 Integer, Parameter :: maxtry=8
integer :: debug_level, debug_unit integer :: debug_level, debug_unit
character(len=20) :: name character(len=20) :: name

Write Preview
Loading…
Cancel
Save