*** empty log message ***
Alfredo Buttari
19 years ago
parent
51046f0c31
commit
0f2594dcc5
@ -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
|
||||
Loading…
Reference in New Issue