Сообщения об ошибках 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.

0 ответов

Другие вопросы по тегам