unkind05
Newbie | Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору Здравствуйте! Надеюсь на помощь) Только только осваиваю фортран (по учёбе приходится), и начало освоения началось с копирования нужной мне подпрограммы библиотеки IMSL. Данная находится на с.216 IMSL Fortran 90 Library Help ,решает,соответственно,разреженные матрицы комплексные...
При запуске текста программы описанной на этой странице винда выдаёт критическую ошибку и программа не работает. Сам найти причину не в состоянии,надеюсь на вашу помощь)
Текст программы:
INTEGER N, NFAC, NZ PARAMETER (N=6, NZ=15, NFAC=3*NZ) INTEGER IPARAM(6), IPVT(N), IRFAC(NFAC), IROW(NZ), JCFAC(NFAC), & & JCOL(NZ), JPVT(N), NL REAL RPARAM(5) COMPLEX A(NZ), FAC(NFAC) DATA A/(3.0,7.0), (3.0,2.0), (-3.0,0.0), (-1.0,3.0), (4.0,2.0), & & (10.0,7.0), (-5.0,4.0), (1.0,6.0), (-1.0,12.0), (-5.0,0.0), & & (12.0,2.0), (-2.0,8.0), (-2.0,-4.0), (-1.0,2.0), (-7.0,7.0)/ DATA IROW/6, 2, 2, 4, 3, 1, 5, 4, 6, 5, 5, 6, 4, 2, 5/ DATA JCOL/6, 2, 3, 5, 3, 1, 1, 4, 1, 4, 5, 2, 1, 4, 6/ IPARAM(1) = 0 CALL LFTZG (N, NZ, A, IROW, JCOL, IPARAM, RPARAM, NFAC, NL, FAC, & & IRFAC, JCFAC, IPVT, JPVT) CALL WRCRN ('fac',NFAC,1,FAC,NFAC,0) CALL WRIRN (' irfac ', 1, NFAC, IRFAC, 1, 0) CALL WRIRN (' jcfac ', 1, NFAC, JCFAC, 1, 0) CALL WRIRN (' p ', 1, N, IPVT, 1, 0) CALL WRIRN (' q ', 1, N, JPVT, 1, 0) END |