ZykovAlex
Newbie | Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
program STBSV_MAIN
character*1 uplo, trans, diag
integer n, k, lda, incx
integer mmax, nmax, xmax
parameter (mmax=5, nmax=7, xmax=5)
parameter (lda=mmax)
real(8) a(mmax, nmax), x(xmax)
integer i, j, k1
! Intrinsic Functions
intrinsic abs
! External Subroutines
external STBSV
! Executable Statements
!
print*
print*, ' S T B S V EXAMPLE PROGRAM'
!
! Read input data from input file
open(1, file = 'stbsvx.d', status = 'old')
read(1, *)
read(1, *) n, k
read(1, *) incx
read(1, 100) uplo, trans, diag
if ( (1+(n-1)*abs(incx))>xmax ) then
print*, ' Insufficient memory for arrays'
goto 99
end if
read(1, *) (x(i),i=1, 1+(n-1)*abs(incx))
if ( (k+1)>mmax .or. n>nmax ) then
print*, ' Insufficient memory for arrays'
goto 99
end if
if (k >= n) then
k1 = n-1
else
k1 = k
end if
if ((uplo == 'U').or.(uplo == 'u')) then
read(1, *) ((a(k-k1+i,j),j=k1+2-i,n),i=1,k1+1)
else
do i = 1, k1+1
read(1, *) (a(i,j),j=1,n+1-i)
end do
end if
!
! Print input data
print*
print*, ' INPUT DATA'
print 101, n, k
print 102, uplo, trans, diag
!call PrintVectorS(0,n,x,incx,'X ')
if ((uplo == 'U').or.(uplo == 'u')) then
!call PrintBandArrayS(1,0,k,n,n,a,lda,'A')
else
!call PrintBandArrayS(1,k,0,n,n,a,lda,'A')
end if
!
! Call STBSV subroutine
call STBSV( uplo, trans, diag, n, k, a, lda, x, incx )
!
print*
print*, ' OUTPUT DATA'
!call PrintVectorS(1,n,x,incx,'X ')
99 continue
print*, 'ЖОПА'
100 format(3(a1,1x))
101 format(7x,'N=',i1,' K=',i1)
102 format(7x,'UPLO=',a1,' TRANS=',a1,' DIAG=',a1)
!stop
end
не работает(
Что не так? |
Сомнительность этого утверждения элементарно проверяется уменьшением размера массива: при 32x32 - работает, при 2048x2048 - не работает. 
