Сообщения об ошибках Slatec
Кто отвечает за библиотеку Slatec? Я нашел ошибку и хотел бы сообщить об этом.
Ошибка заключается в том, что в подпрограмме SSMMI2 переменная U должна иметь размер, а не размер:
REAL B(N), DINV(N), L(*), U(n), X(N)
должно быть
REAL B(N), DINV(N), L(*), U(*), X(N)
Если размер определен явно, позже в подпрограмме в этом цикле:
DO 50 J = JBGN, JEND
X(IU(J)) = X(IU(J)) - U(J)*X(ICOL)
50 CONTINUE
J может быть больше, чем N, например, в тестовых программах Джона Беркардта для slatec.