В настоящее время я пытаюсь восстановить некоторые функции функции оценки числа условий R * kappa, которая оценивает число условий матрицы X по:
- Разработка QR-разложения X.
- Вызов
LAPACK dtrcon или LINPACK dtrco (в зависимости от того, каковы основные зависимости в системе пользователей) и вычисление числа условий R, верхней треугольной матрицы, которая должна иметь тот же номер условия, что и X (см. здесь ).
Я пытался понять, что делают алгоритмы LAPACK и LINPACK, поскольку это может быть чрезвычайно полезно для моего собственного кодирования.
Мне удалось найти алгоритм, который LINPACK использует, который был описан здесь , но мне не повезло найти источник алгоритма LAPACK. Комментарии в функции R kappa предполагают, что они используют разные алгоритмы (см. здесь ), но я не уверен ...
Короче говоря, мой вопрос:
Кто-нибудь знает, если LAPACK dtrcon и LINPACK dtrco используют один и тот же алгоритм, и если нет, то какой алгоритм LAPACK dtrcon использует?
Заранее спасибо!