Fix setelem definition and impl.

base/modules/tools/psb_c_tools_mod.F90

@@ -454,6 +454,17 @@ Module psb_c_tools_mod
     end function
   end interface
 
+  interface psb_setelem
+    subroutine psb_c_setelem(index,val,x,desc_a,info)
+      import
+      type(psb_c_vect_type), intent(inout) :: x
+      integer(psb_lpk_), intent(in)          :: index
+      type(psb_desc_type), intent(inout)     :: desc_a
+      integer(psb_ipk_), intent(out)         :: info
+      complex(psb_spk_)                        ::val
+    end subroutine psb_c_setelem
+  end interface
+
   interface psb_remap
     subroutine psb_c_remap(np_remap, desc_in, a_in, &
          & ipd, isrc, nrsrc, naggr, desc_out, a_out, info)

base/modules/tools/psb_d_tools_mod.F90

@@ -454,6 +454,17 @@ Module psb_d_tools_mod
     end function
   end interface
 
+  interface psb_setelem
+    subroutine psb_d_setelem(index,val,x,desc_a,info)
+      import
+      type(psb_d_vect_type), intent(inout) :: x
+      integer(psb_lpk_), intent(in)          :: index
+      type(psb_desc_type), intent(inout)     :: desc_a
+      integer(psb_ipk_), intent(out)         :: info
+      real(psb_dpk_)                        ::val
+    end subroutine psb_d_setelem
+  end interface
+
   interface psb_remap
     subroutine psb_d_remap(np_remap, desc_in, a_in, &
          & ipd, isrc, nrsrc, naggr, desc_out, a_out, info)

base/modules/tools/psb_s_tools_mod.F90

@@ -454,6 +454,17 @@ Module psb_s_tools_mod
     end function
   end interface
 
+  interface psb_setelem
+    subroutine psb_s_setelem(index,val,x,desc_a,info)
+      import
+      type(psb_s_vect_type), intent(inout) :: x
+      integer(psb_lpk_), intent(in)          :: index
+      type(psb_desc_type), intent(inout)     :: desc_a
+      integer(psb_ipk_), intent(out)         :: info
+      real(psb_spk_)                        ::val
+    end subroutine psb_s_setelem
+  end interface
+
   interface psb_remap
     subroutine psb_s_remap(np_remap, desc_in, a_in, &
          & ipd, isrc, nrsrc, naggr, desc_out, a_out, info)

base/modules/tools/psb_z_tools_mod.F90

@@ -454,6 +454,17 @@ Module psb_z_tools_mod
     end function
   end interface
 
+  interface psb_setelem
+    subroutine psb_z_setelem(index,val,x,desc_a,info)
+      import
+      type(psb_z_vect_type), intent(inout) :: x
+      integer(psb_lpk_), intent(in)          :: index
+      type(psb_desc_type), intent(inout)     :: desc_a
+      integer(psb_ipk_), intent(out)         :: info
+      complex(psb_dpk_)                        ::val
+    end subroutine psb_z_setelem
+  end interface
+
   interface psb_remap
     subroutine psb_z_remap(np_remap, desc_in, a_in, &
          & ipd, isrc, nrsrc, naggr, desc_out, a_out, info)

base/tools/Makefile

@@ -28,7 +28,8 @@ FOBJS = psb_cdall.o psb_cdals.o psb_cdalv.o psb_cd_inloc.o psb_cdins.o psb_cdprt
 	psb_s_map.o psb_d_map.o psb_c_map.o psb_z_map.o \
 	psb_s_par_csr_spspmm.o psb_d_par_csr_spspmm.o psb_c_par_csr_spspmm.o psb_z_par_csr_spspmm.o \
 	psb_s_glob_transpose.o psb_d_glob_transpose.o psb_c_glob_transpose.o psb_z_glob_transpose.o \
-	psb_cgetelem.o psb_dgetelem.o psb_sgetelem.o psb_zgetelem.o   
+	psb_cgetelem.o psb_dgetelem.o psb_sgetelem.o psb_zgetelem.o   \
+	psb_csetelem.o psb_dsetelem.o psb_ssetelem.o psb_zsetelem.o   
 
 MPFOBJS = psb_icdasb.o psb_ssphalo.o  psb_dsphalo.o psb_csphalo.o psb_zsphalo.o  \
 	psb_dcdbldext.o psb_zcdbldext.o psb_scdbldext.o psb_ccdbldext.o \

base/tools/psb_csetelem.f90

@@ -0,0 +1,112 @@
+!
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone
+!        Alfredo Buttari
+!
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!
+! Function: psb_c_setelem
+!    Set  entries into  a dense  vector. Note: the row indices in index
+!    are assumed to be in global numbering and are converted on the fly.
+!    Row indices not belonging to the current process have to be in the halo,
+!    othewise failure is ensured.
+!
+! Arguments:
+!    x       - type(psb_c_vect_type) The source vector
+!    desc_a  - type(psb_desc_type).    The communication descriptor.
+!    index   - integer. Row index of x of the value to extract
+!    iam     - integer. Index of the process requesting the value
+!    info    - integer.                       return code
+
+
+subroutine psb_c_setelem(index,val,x,desc_a,info) 
+  use psb_base_mod, psb_protect_name => psb_c_setelem
+  use psi_mod
+  implicit none
+
+  type(psb_c_vect_type), intent(inout) :: x
+  type(psb_desc_type), intent(inout)     :: desc_a
+  integer(psb_lpk_), intent(in)          :: index
+  integer(psb_ipk_), intent(out)         :: info
+  complex(psb_spk_)                        :: val
+
+  !locals
+  integer(psb_ipk_) :: localindex(1)
+  type(psb_ctxt_type) :: ctxt
+  integer(psb_ipk_) :: np, me, err_act
+  integer(psb_lpk_) :: gindex(1)
+  integer(psb_lpk_), allocatable :: myidx(:),mylocal(:)
+  character(len=20) :: name
+  logical, parameter :: debug = .false.
+
+  info = 0
+  
+  gindex(1) = index
+  
+  if (psb_errstatus_fatal()) return
+  info=psb_success_
+  call psb_erractionsave(err_act)
+  name = 'psb_c_getelem'
+
+  if (.not.desc_a%is_ok()) then
+    info = psb_err_invalid_cd_state_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+
+  ctxt = desc_a%get_context()
+
+  call psb_info(ctxt, me, np)
+  if (np == -1) then
+    info = psb_err_context_error_
+    call psb_errpush(info,name)
+    goto 9999
+  endif
+
+  call desc_a%indxmap%g2l(gindex,localindex,info,owned=.false.)
+  if(debug.and.(localindex(1) < 1)) then
+    write(*,*)"Process ",me," owns ",desc_a%get_local_rows()," rows"," Global index is ",gindex,"Local index is ",localindex
+    myidx = desc_a%get_global_indices(owned=.false.)
+    mylocal = desc_a%get_global_indices(owned=.true.)
+    write(*,*)"My (local+halo) indexes are: ",myidx
+    write(*,*)"My (local) indexes are: ",mylocal
+  end if
+  if ( localindex(1) < 1) then
+    info = psb_err_internal_error_
+    call psb_errpush(info,name,a_err="Index not in the HALO")
+    goto 9999
+  else
+    call  x%set_entry(localindex(1),val)
+  end if
+  call psb_erractionrestore(err_act)
+  return
+
+9999 call psb_error_handler(ctxt,err_act)
+
+  return
+
+end subroutine psb_c_setelem
+

base/tools/psb_dsetelem.f90

@@ -0,0 +1,112 @@
+!
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone
+!        Alfredo Buttari
+!
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!
+! Function: psb_d_setelem
+!    Set  entries into  a dense  vector. Note: the row indices in index
+!    are assumed to be in global numbering and are converted on the fly.
+!    Row indices not belonging to the current process have to be in the halo,
+!    othewise failure is ensured.
+!
+! Arguments:
+!    x       - type(psb_d_vect_type) The source vector
+!    desc_a  - type(psb_desc_type).    The communication descriptor.
+!    index   - integer. Row index of x of the value to extract
+!    iam     - integer. Index of the process requesting the value
+!    info    - integer.                       return code
+
+
+subroutine psb_d_setelem(index,val,x,desc_a,info) 
+  use psb_base_mod, psb_protect_name => psb_d_setelem
+  use psi_mod
+  implicit none
+
+  type(psb_d_vect_type), intent(inout) :: x
+  type(psb_desc_type), intent(inout)     :: desc_a
+  integer(psb_lpk_), intent(in)          :: index
+  integer(psb_ipk_), intent(out)         :: info
+  real(psb_dpk_)                        :: val
+
+  !locals
+  integer(psb_ipk_) :: localindex(1)
+  type(psb_ctxt_type) :: ctxt
+  integer(psb_ipk_) :: np, me, err_act
+  integer(psb_lpk_) :: gindex(1)
+  integer(psb_lpk_), allocatable :: myidx(:),mylocal(:)
+  character(len=20) :: name
+  logical, parameter :: debug = .false.
+
+  info = 0
+  
+  gindex(1) = index
+  
+  if (psb_errstatus_fatal()) return
+  info=psb_success_
+  call psb_erractionsave(err_act)
+  name = 'psb_d_getelem'
+
+  if (.not.desc_a%is_ok()) then
+    info = psb_err_invalid_cd_state_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+
+  ctxt = desc_a%get_context()
+
+  call psb_info(ctxt, me, np)
+  if (np == -1) then
+    info = psb_err_context_error_
+    call psb_errpush(info,name)
+    goto 9999
+  endif
+
+  call desc_a%indxmap%g2l(gindex,localindex,info,owned=.false.)
+  if(debug.and.(localindex(1) < 1)) then
+    write(*,*)"Process ",me," owns ",desc_a%get_local_rows()," rows"," Global index is ",gindex,"Local index is ",localindex
+    myidx = desc_a%get_global_indices(owned=.false.)
+    mylocal = desc_a%get_global_indices(owned=.true.)
+    write(*,*)"My (local+halo) indexes are: ",myidx
+    write(*,*)"My (local) indexes are: ",mylocal
+  end if
+  if ( localindex(1) < 1) then
+    info = psb_err_internal_error_
+    call psb_errpush(info,name,a_err="Index not in the HALO")
+    goto 9999
+  else
+    call  x%set_entry(localindex(1),val)
+  end if
+  call psb_erractionrestore(err_act)
+  return
+
+9999 call psb_error_handler(ctxt,err_act)
+
+  return
+
+end subroutine psb_d_setelem
+

base/tools/psb_ssetelem.f90

@@ -0,0 +1,112 @@
+!
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone
+!        Alfredo Buttari
+!
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!
+! Function: psb_s_setelem
+!    Set  entries into  a dense  vector. Note: the row indices in index
+!    are assumed to be in global numbering and are converted on the fly.
+!    Row indices not belonging to the current process have to be in the halo,
+!    othewise failure is ensured.
+!
+! Arguments:
+!    x       - type(psb_s_vect_type) The source vector
+!    desc_a  - type(psb_desc_type).    The communication descriptor.
+!    index   - integer. Row index of x of the value to extract
+!    iam     - integer. Index of the process requesting the value
+!    info    - integer.                       return code
+
+
+subroutine psb_s_setelem(index,val,x,desc_a,info) 
+  use psb_base_mod, psb_protect_name => psb_s_setelem
+  use psi_mod
+  implicit none
+
+  type(psb_s_vect_type), intent(inout) :: x
+  type(psb_desc_type), intent(inout)     :: desc_a
+  integer(psb_lpk_), intent(in)          :: index
+  integer(psb_ipk_), intent(out)         :: info
+  real(psb_spk_)                        :: val
+
+  !locals
+  integer(psb_ipk_) :: localindex(1)
+  type(psb_ctxt_type) :: ctxt
+  integer(psb_ipk_) :: np, me, err_act
+  integer(psb_lpk_) :: gindex(1)
+  integer(psb_lpk_), allocatable :: myidx(:),mylocal(:)
+  character(len=20) :: name
+  logical, parameter :: debug = .false.
+
+  info = 0
+  
+  gindex(1) = index
+  
+  if (psb_errstatus_fatal()) return
+  info=psb_success_
+  call psb_erractionsave(err_act)
+  name = 'psb_s_getelem'
+
+  if (.not.desc_a%is_ok()) then
+    info = psb_err_invalid_cd_state_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+
+  ctxt = desc_a%get_context()
+
+  call psb_info(ctxt, me, np)
+  if (np == -1) then
+    info = psb_err_context_error_
+    call psb_errpush(info,name)
+    goto 9999
+  endif
+
+  call desc_a%indxmap%g2l(gindex,localindex,info,owned=.false.)
+  if(debug.and.(localindex(1) < 1)) then
+    write(*,*)"Process ",me," owns ",desc_a%get_local_rows()," rows"," Global index is ",gindex,"Local index is ",localindex
+    myidx = desc_a%get_global_indices(owned=.false.)
+    mylocal = desc_a%get_global_indices(owned=.true.)
+    write(*,*)"My (local+halo) indexes are: ",myidx
+    write(*,*)"My (local) indexes are: ",mylocal
+  end if
+  if ( localindex(1) < 1) then
+    info = psb_err_internal_error_
+    call psb_errpush(info,name,a_err="Index not in the HALO")
+    goto 9999
+  else
+    call  x%set_entry(localindex(1),val)
+  end if
+  call psb_erractionrestore(err_act)
+  return
+
+9999 call psb_error_handler(ctxt,err_act)
+
+  return
+
+end subroutine psb_s_setelem
+

base/tools/psb_zsetelem.f90

@@ -0,0 +1,112 @@
+!
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone
+!        Alfredo Buttari
+!
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!
+! Function: psb_z_setelem
+!    Set  entries into  a dense  vector. Note: the row indices in index
+!    are assumed to be in global numbering and are converted on the fly.
+!    Row indices not belonging to the current process have to be in the halo,
+!    othewise failure is ensured.
+!
+! Arguments:
+!    x       - type(psb_z_vect_type) The source vector
+!    desc_a  - type(psb_desc_type).    The communication descriptor.
+!    index   - integer. Row index of x of the value to extract
+!    iam     - integer. Index of the process requesting the value
+!    info    - integer.                       return code
+
+
+subroutine psb_z_setelem(index,val,x,desc_a,info) 
+  use psb_base_mod, psb_protect_name => psb_z_setelem
+  use psi_mod
+  implicit none
+
+  type(psb_z_vect_type), intent(inout) :: x
+  type(psb_desc_type), intent(inout)     :: desc_a
+  integer(psb_lpk_), intent(in)          :: index
+  integer(psb_ipk_), intent(out)         :: info
+  complex(psb_dpk_)                        :: val
+
+  !locals
+  integer(psb_ipk_) :: localindex(1)
+  type(psb_ctxt_type) :: ctxt
+  integer(psb_ipk_) :: np, me, err_act
+  integer(psb_lpk_) :: gindex(1)
+  integer(psb_lpk_), allocatable :: myidx(:),mylocal(:)
+  character(len=20) :: name
+  logical, parameter :: debug = .false.
+
+  info = 0
+  
+  gindex(1) = index
+  
+  if (psb_errstatus_fatal()) return
+  info=psb_success_
+  call psb_erractionsave(err_act)
+  name = 'psb_z_getelem'
+
+  if (.not.desc_a%is_ok()) then
+    info = psb_err_invalid_cd_state_
+    call psb_errpush(info,name)
+    goto 9999
+  end if
+
+  ctxt = desc_a%get_context()
+
+  call psb_info(ctxt, me, np)
+  if (np == -1) then
+    info = psb_err_context_error_
+    call psb_errpush(info,name)
+    goto 9999
+  endif
+
+  call desc_a%indxmap%g2l(gindex,localindex,info,owned=.false.)
+  if(debug.and.(localindex(1) < 1)) then
+    write(*,*)"Process ",me," owns ",desc_a%get_local_rows()," rows"," Global index is ",gindex,"Local index is ",localindex
+    myidx = desc_a%get_global_indices(owned=.false.)
+    mylocal = desc_a%get_global_indices(owned=.true.)
+    write(*,*)"My (local+halo) indexes are: ",myidx
+    write(*,*)"My (local) indexes are: ",mylocal
+  end if
+  if ( localindex(1) < 1) then
+    info = psb_err_internal_error_
+    call psb_errpush(info,name,a_err="Index not in the HALO")
+    goto 9999
+  else
+    call  x%set_entry(localindex(1),val)
+  end if
+  call psb_erractionrestore(err_act)
+  return
+
+9999 call psb_error_handler(ctxt,err_act)
+
+  return
+
+end subroutine psb_z_setelem
+

cbind/base/psb_c_tools_cbind_mod.F90

@@ -692,7 +692,7 @@ contains
 
     ixb = psb_c_get_index_base()
     if (ixb == 1) then
-      call  psb_setelem(index,val,xp,descp,info)
+      call psb_setelem(index,val,xp,descp,info)
     else
       call psb_setelem(index+(1-ixb),val,xp,descp,info)
     end if

cbind/base/psb_d_tools_cbind_mod.F90

@@ -702,7 +702,7 @@ contains
 
     ixb = psb_c_get_index_base()
     if (ixb == 1) then
-      call  psb_setelem(index,val,xp,descp,info)
+      call psb_setelem(index,val,xp,descp,info)
     else
       call psb_setelem(index+(1-ixb),val,xp,descp,info)
     end if

cbind/base/psb_s_tools_cbind_mod.F90

@@ -702,7 +702,7 @@ contains
 
     ixb = psb_c_get_index_base()
     if (ixb == 1) then
-      call  psb_setelem(index,val,xp,descp,info)
+      call psb_setelem(index,val,xp,descp,info)
     else
       call psb_setelem(index+(1-ixb),val,xp,descp,info)
     end if

cbind/base/psb_z_tools_cbind_mod.F90

@@ -692,7 +692,7 @@ contains
 
     ixb = psb_c_get_index_base()
     if (ixb == 1) then
-      call  psb_setelem(index,val,xp,descp,info)
+      call psb_setelem(index,val,xp,descp,info)
     else
       call psb_setelem(index+(1-ixb),val,xp,descp,info)
     end if

linsolve/psb_c_linsolve_mod.f90

@@ -0,0 +1,81 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+!
+! File: psb_linsolve_mod.f90
+!  Interfaces for linear solvers.
+!
+Module psb_c_linsolve_mod
+
+  use psb_const_mod
+  public 
+
+  interface psb_krylov
+    Subroutine psb_ckrylov_vect(method,a,prec,b,x,eps,desc_a,info,&
+         & itmax,iter,err,itrace,irst,istop,cond,s1,s2)
+      use psb_base_mod, only  : psb_ipk_, psb_desc_type, psb_cspmat_type, &
+           & psb_spk_, psb_c_vect_type
+      use psb_prec_mod, only : psb_cprec_type
+      character(len=*)                      :: method
+      Type(psb_cspmat_type), Intent(in)     :: a
+      Type(psb_desc_type), Intent(in)       :: desc_a
+      class(psb_cprec_type), intent(inout)  :: prec 
+      type(psb_c_vect_type), Intent(inout)  :: b
+      type(psb_c_vect_type), Intent(inout)  :: x
+      Real(psb_spk_), Intent(in)             :: eps
+      integer(psb_ipk_), intent(out)            :: info
+      integer(psb_ipk_), Optional, Intent(in)   :: itmax, itrace, irst,istop
+      integer(psb_ipk_), Optional, Intent(out)  :: iter
+      Real(psb_spk_), Optional, Intent(out) :: err,cond
+      type(psb_c_vect_type), optional  :: s1,s2
+    end Subroutine psb_ckrylov_vect
+  end interface
+
+  interface psb_richardson
+    Subroutine psb_crichardson_vect(a,prec,b,x,eps,desc_a,info,&
+         & itmax,iter,err,itrace,istop)
+      use psb_base_mod, only  : psb_ipk_, psb_desc_type, psb_cspmat_type, &
+           & psb_spk_, psb_c_vect_type
+      use psb_prec_mod, only : psb_cprec_type
+      Type(psb_cspmat_type), Intent(in)     :: a
+      Type(psb_desc_type), Intent(in)       :: desc_a
+      class(psb_cprec_type), intent(inout)  :: prec 
+      type(psb_c_vect_type), Intent(inout)  :: b
+      type(psb_c_vect_type), Intent(inout)  :: x
+      Real(psb_spk_), Intent(in)            :: eps
+      integer(psb_ipk_), intent(out)            :: info
+      integer(psb_ipk_), Optional, Intent(in)   :: itmax, itrace, istop
+      integer(psb_ipk_), Optional, Intent(out)  :: iter
+      Real(psb_spk_), Optional, Intent(out) :: err
+    end Subroutine psb_crichardson_vect
+  end interface
+  
+end module psb_c_linsolve_mod

linsolve/psb_d_linsolve_mod.f90

@@ -0,0 +1,81 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+!
+! File: psb_linsolve_mod.f90
+!  Interfaces for linear solvers.
+!
+Module psb_d_linsolve_mod
+
+  use psb_const_mod
+  public 
+
+  interface psb_krylov
+    Subroutine psb_dkrylov_vect(method,a,prec,b,x,eps,desc_a,info,&
+         & itmax,iter,err,itrace,irst,istop,cond,s1,s2)
+      use psb_base_mod, only  : psb_ipk_, psb_desc_type, psb_dspmat_type, &
+           & psb_dpk_, psb_d_vect_type
+      use psb_prec_mod, only : psb_dprec_type
+      character(len=*)                      :: method
+      Type(psb_dspmat_type), Intent(in)     :: a
+      Type(psb_desc_type), Intent(in)       :: desc_a
+      class(psb_dprec_type), intent(inout)  :: prec 
+      type(psb_d_vect_type), Intent(inout)  :: b
+      type(psb_d_vect_type), Intent(inout)  :: x
+      Real(psb_dpk_), Intent(in)             :: eps
+      integer(psb_ipk_), intent(out)            :: info
+      integer(psb_ipk_), Optional, Intent(in)   :: itmax, itrace, irst,istop
+      integer(psb_ipk_), Optional, Intent(out)  :: iter
+      Real(psb_dpk_), Optional, Intent(out) :: err,cond
+      type(psb_d_vect_type), optional  :: s1,s2
+    end Subroutine psb_dkrylov_vect
+  end interface
+
+  interface psb_richardson
+    Subroutine psb_drichardson_vect(a,prec,b,x,eps,desc_a,info,&
+         & itmax,iter,err,itrace,istop)
+      use psb_base_mod, only  : psb_ipk_, psb_desc_type, psb_dspmat_type, &
+           & psb_dpk_, psb_d_vect_type
+      use psb_prec_mod, only : psb_dprec_type
+      Type(psb_dspmat_type), Intent(in)     :: a
+      Type(psb_desc_type), Intent(in)       :: desc_a
+      class(psb_dprec_type), intent(inout)  :: prec 
+      type(psb_d_vect_type), Intent(inout)  :: b
+      type(psb_d_vect_type), Intent(inout)  :: x
+      Real(psb_dpk_), Intent(in)            :: eps
+      integer(psb_ipk_), intent(out)            :: info
+      integer(psb_ipk_), Optional, Intent(in)   :: itmax, itrace, istop
+      integer(psb_ipk_), Optional, Intent(out)  :: iter
+      Real(psb_dpk_), Optional, Intent(out) :: err
+    end Subroutine psb_drichardson_vect
+  end interface
+  
+end module psb_d_linsolve_mod

linsolve/psb_s_linsolve_mod.f90

@@ -0,0 +1,81 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+!
+! File: psb_linsolve_mod.f90
+!  Interfaces for linear solvers.
+!
+Module psb_s_linsolve_mod
+
+  use psb_const_mod
+  public 
+
+  interface psb_krylov
+    Subroutine psb_skrylov_vect(method,a,prec,b,x,eps,desc_a,info,&
+         & itmax,iter,err,itrace,irst,istop,cond,s1,s2)
+      use psb_base_mod, only  : psb_ipk_, psb_desc_type, psb_sspmat_type, &
+           & psb_spk_, psb_s_vect_type
+      use psb_prec_mod, only : psb_sprec_type
+      character(len=*)                      :: method
+      Type(psb_sspmat_type), Intent(in)     :: a
+      Type(psb_desc_type), Intent(in)       :: desc_a
+      class(psb_sprec_type), intent(inout)  :: prec 
+      type(psb_s_vect_type), Intent(inout)  :: b
+      type(psb_s_vect_type), Intent(inout)  :: x
+      Real(psb_spk_), Intent(in)             :: eps
+      integer(psb_ipk_), intent(out)            :: info
+      integer(psb_ipk_), Optional, Intent(in)   :: itmax, itrace, irst,istop
+      integer(psb_ipk_), Optional, Intent(out)  :: iter
+      Real(psb_spk_), Optional, Intent(out) :: err,cond
+      type(psb_s_vect_type), optional  :: s1,s2
+    end Subroutine psb_skrylov_vect
+  end interface
+
+  interface psb_richardson
+    Subroutine psb_srichardson_vect(a,prec,b,x,eps,desc_a,info,&
+         & itmax,iter,err,itrace,istop)
+      use psb_base_mod, only  : psb_ipk_, psb_desc_type, psb_sspmat_type, &
+           & psb_spk_, psb_s_vect_type
+      use psb_prec_mod, only : psb_sprec_type
+      Type(psb_sspmat_type), Intent(in)     :: a
+      Type(psb_desc_type), Intent(in)       :: desc_a
+      class(psb_sprec_type), intent(inout)  :: prec 
+      type(psb_s_vect_type), Intent(inout)  :: b
+      type(psb_s_vect_type), Intent(inout)  :: x
+      Real(psb_spk_), Intent(in)            :: eps
+      integer(psb_ipk_), intent(out)            :: info
+      integer(psb_ipk_), Optional, Intent(in)   :: itmax, itrace, istop
+      integer(psb_ipk_), Optional, Intent(out)  :: iter
+      Real(psb_spk_), Optional, Intent(out) :: err
+    end Subroutine psb_srichardson_vect
+  end interface
+  
+end module psb_s_linsolve_mod

linsolve/psb_z_linsolve_mod.f90

@@ -0,0 +1,81 @@
+!   
+!                Parallel Sparse BLAS  version 3.5
+!      (C) Copyright 2006-2018
+!        Salvatore Filippone    
+!        Alfredo Buttari      
+!   
+!    Redistribution and use in source and binary forms, with or without
+!    modification, are permitted provided that the following conditions
+!    are met:
+!      1. Redistributions of source code must retain the above copyright
+!         notice, this list of conditions and the following disclaimer.
+!      2. Redistributions in binary form must reproduce the above copyright
+!         notice, this list of conditions, and the following disclaimer in the
+!         documentation and/or other materials provided with the distribution.
+!      3. The name of the PSBLAS group or the names of its contributors may
+!         not be used to endorse or promote products derived from this
+!         software without specific prior written permission.
+!   
+!    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+!    ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
+!    TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+!    PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE PSBLAS GROUP OR ITS CONTRIBUTORS
+!    BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+!    CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+!    SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+!    INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+!    CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+!    ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+!    POSSIBILITY OF SUCH DAMAGE.
+!   
+!    
+!
+! File: psb_linsolve_mod.f90
+!  Interfaces for linear solvers.
+!
+Module psb_z_linsolve_mod
+
+  use psb_const_mod
+  public 
+
+  interface psb_krylov
+    Subroutine psb_zkrylov_vect(method,a,prec,b,x,eps,desc_a,info,&
+         & itmax,iter,err,itrace,irst,istop,cond,s1,s2)
+      use psb_base_mod, only  : psb_ipk_, psb_desc_type, psb_zspmat_type, &
+           & psb_dpk_, psb_z_vect_type
+      use psb_prec_mod, only : psb_zprec_type
+      character(len=*)                      :: method
+      Type(psb_zspmat_type), Intent(in)     :: a
+      Type(psb_desc_type), Intent(in)       :: desc_a
+      class(psb_zprec_type), intent(inout)  :: prec 
+      type(psb_z_vect_type), Intent(inout)  :: b
+      type(psb_z_vect_type), Intent(inout)  :: x
+      Real(psb_dpk_), Intent(in)             :: eps
+      integer(psb_ipk_), intent(out)            :: info
+      integer(psb_ipk_), Optional, Intent(in)   :: itmax, itrace, irst,istop
+      integer(psb_ipk_), Optional, Intent(out)  :: iter
+      Real(psb_dpk_), Optional, Intent(out) :: err,cond
+      type(psb_z_vect_type), optional  :: s1,s2
+    end Subroutine psb_zkrylov_vect
+  end interface
+
+  interface psb_richardson
+    Subroutine psb_zrichardson_vect(a,prec,b,x,eps,desc_a,info,&
+         & itmax,iter,err,itrace,istop)
+      use psb_base_mod, only  : psb_ipk_, psb_desc_type, psb_zspmat_type, &
+           & psb_dpk_, psb_z_vect_type
+      use psb_prec_mod, only : psb_zprec_type
+      Type(psb_zspmat_type), Intent(in)     :: a
+      Type(psb_desc_type), Intent(in)       :: desc_a
+      class(psb_zprec_type), intent(inout)  :: prec 
+      type(psb_z_vect_type), Intent(inout)  :: b
+      type(psb_z_vect_type), Intent(inout)  :: x
+      Real(psb_dpk_), Intent(in)            :: eps
+      integer(psb_ipk_), intent(out)            :: info
+      integer(psb_ipk_), Optional, Intent(in)   :: itmax, itrace, istop
+      integer(psb_ipk_), Optional, Intent(out)  :: iter
+      Real(psb_dpk_), Optional, Intent(out) :: err
+    end Subroutine psb_zrichardson_vect
+  end interface
+  
+end module psb_z_linsolve_mod