From 0f2594dcc5433880a96396128f6d98110f1eeb6b Mon Sep 17 00:00:00 2001
From: Alfredo Buttari <alfredo.buttari@enseeiht.fr>
Date: Tue, 29 Nov 2005 15:37:47 +0000
Subject: [PATCH] *** empty log message ***

---
 src/modules/psb_spsb_mod.f90 | 143 +++++++++++++++++++++++++++++++++++
 1 file changed, 143 insertions(+)
 create mode 100644 src/modules/psb_spsb_mod.f90

diff --git a/src/modules/psb_spsb_mod.f90 b/src/modules/psb_spsb_mod.f90
new file mode 100644
index 00000000..b1893d87
--- /dev/null
+++ b/src/modules/psb_spsb_mod.f90
@@ -0,0 +1,143 @@
+module psb_spsb_mod
+
+  interface
+      real(kind(1.d0)) function ddot(n,dx,incx,dy,incy)
+      real(kind(1.d0)) :: dx(*),dy(*),dtemp
+      integer          :: i,incx,incy,ix,iy,m,mp1,n
+    end function ddot
+ end interface
+
+  interface
+      subroutine  daxpby(m, n, alpha, X, lldx, beta, Y, lldy, info)
+      real(kind(1.d0)) :: X(lldx,*), Y(lldy,*)
+      real(kind(1.d0)) :: alpha, beta
+      integer          :: n, m, lldx, lldy, info
+    end subroutine daxpby
+ end interface
+
+ interface
+      subroutine dcsmm(trans,m,n,k,alpha,pl,fida,descra,a,ia1,ia2,&
+     &   infoa,pr,b,ldb,beta,c,ldc,work,lwork,ierror)
+      integer          ::  m,n,k,ldb,ldc,lwork, ierror
+      character        ::  trans
+      real(kind(1.d0)) ::  alpha,beta
+      integer          ::  ia1(*),ia2(*),infoa(*),pl(*),pr(*)
+      character        ::  descra*11, fida*5
+      real(kind(1.d0)) ::  a(*),b(ldb,*),c(ldc,*),work(*)
+    end subroutine dcsmm
+ end interface
+
+ interface
+    real(kind(1.d0)) function dcsnmi(trans,m,n,fida,descra,a,ia1,ia2,&
+         &                 infoa,ierror)
+      integer          ::  m,n, ierror
+      character        ::  trans
+      integer          ::  ia1(*),ia2(*),infoa(*)
+      character        ::  descra*11, fida*5
+      real(kind(1.d0)) ::  a(*)
+    end function dcsnmi
+ end interface
+
+ interface
+      subroutine dcsrp(trans,m,n,fida,descra,ia1,ia2,infoa,&
+     &  p,work,lwork,ierror)
+      implicit none                                                      
+      integer          :: lwork, m, n, ierror
+      character        :: trans
+      real(kind(1.d0)) :: work(lwork)
+      integer          :: ia1(*), ia2(*), infoa(*), p(*), int_val(5)
+      character        :: descra*11, fida*5
+    end subroutine dcsrp
+ end interface
+
+ interface
+      subroutine dcssm(trans,m,n,alpha,unitd,d,&
+     &   pl,fidt,descrt,t,it1,it2,infot,pr,&
+     &   b,ldb,beta,c,ldc,work,lwork,ierror)
+      real(kind(1.d0)) :: alpha, beta
+      integer          :: n, ldb, ldc, m, lwork, ierror
+      character        :: unitd, trans
+      real(kind(1.d0)) :: t(*), b(ldb,*), c(ldc,*), d(*), work(*)
+      integer          :: it1(*), it2(*), infot(*), pl(*), pr(*)
+      character        :: descrt*11, fidt*5
+    end subroutine dcssm
+ end interface
+
+ interface
+    subroutine dcsupd(m, n, fida, descra, a, ia1,&
+         &  ia2, infoa, ia, ja, fidh, descrh, h, ih1, ih2,&
+         &  infoh, ih, jh, flag, glob_to_loc,&
+         &  iwork, liwork, ierror)
+      integer          :: ia, ja, ih, jh, m, n,&
+           &  ierror, flag, liwork 
+      integer          :: ia1(*),ia2(*),ih1(*),ih2(*),&
+           &  infoa(*),infoh(*),iwork(*),&
+           &  glob_to_loc(*)
+      character        :: descra*11,descrh*11, fida*5, fidh*5
+      real(kind(1.d0)) :: a(*),h(*)    
+    end subroutine dcsupd
+ end interface
+ 
+ interface
+    subroutine dgelp(trans,m,n,p,b,ldb,work,lwork,ierror)
+      integer          :: ldb, m, n, lwork, ierror
+      character        :: trans
+      real(kind(1.d0)) :: b(ldb,*), work(*)
+      integer          :: p(*)
+    end subroutine dgelp
+ end interface
+ 
+ interface
+    subroutine dlpupd(m,n,perm,b,ldb,beta,c,ldc)
+      integer          :: m, n, ldb, ldc
+      real(kind(1.d0)) :: beta
+      integer          :: perm(*)
+      real(kind(1.d0)) :: b(ldb,*), c(ldc,*)
+    end subroutine dlpupd
+ end interface
+ 
+ interface
+      subroutine dswmm(trans,m,n,k,alpha,fida,descra,a,ia1,ia2,&
+           & infoa,b,ldb,beta,c,ldc,work,lwork,ierror)
+      integer          :: m,n,k,ldb,ldc,lwork,ierror
+      character        :: trans
+      real(kind(1.d0)) :: alpha,beta
+      integer          :: ia1(*),ia2(*),infoa(*), int_val(5)
+      character        :: descra*11, fida*5
+      real(kind(1.d0)) :: a(*),b(ldb,*),c(ldc,*),work(*)
+    end subroutine dswmm
+ end interface
+
+
+ interface
+    subroutine dswprt(m,n,fida,descra,a,ia1,ia2,infoa,title,&
+         & iout,ierror)
+      integer       :: m,n,iout,ierror
+      integer       :: ia1(*),ia2(*),infoa(*)
+      character     :: descra*11, fida*5, title*(*)
+    end subroutine dswprt
+ end interface
+ 
+ interface
+      subroutine dswsm(trans,m,n,alpha,unitd,d,fidt,descrt,t,it1,it2,& 
+           & infot,b,ldb,beta,c,ldc,work,lwork,ierror)
+      integer          :: m, n, ldb, ldc, lwork, ierror
+      character        :: unitd, trans
+      real(kind(1.d0)) :: alpha, beta
+      integer          :: it1(*), it2(*), infot(*)
+      character        :: descrt*11, fidt*5
+      real(kind(1.d0)) :: t(*), b(ldb,*), c(ldc,*), d(*), work(*)
+    end subroutine dswsm
+ end interface
+
+ interface
+      subroutine symbmm (n, m, l, ia, ja, diaga, ib,&
+           & jb, diagb, ic, jc, diagc, index)
+      integer          :: ia(*), ja(*), diaga, ib(*),&
+           & jb(*), diagb, diagc, index(*)
+      integer, pointer :: ic(:),jc(:)
+      integer          :: nze, info
+    end subroutine symbmm
+ end interface
+
+end module psb_spsb_mod