Оценка числа обусловленности симметричной матрицы

Алгоритм оценки числа обусловленности симметричной матрицы аналогичен алгоритму, оценивающему число обусловленности матрицы общего вида, за исключением того, что в ходе оценки используется не LU-разложение матрицы, а LDLT-разложение. В связи с этим я не буду повторять здесь описание принципов работы алгоритма, поскольку желающие могут познакомиться с ним по ссылке, указанной выше.

Модуль включает в себя две подпрограммы: SMatrixLDLTRCond и SMatrixRCond. Первая подпрограмма оценивает число обусловленности матрицы, заданной LDLT-разложением, а вторая - симметричной матрицы, чье LDLT-разложение ещё не построено. В связи с тем, что у симметричных матриц 1-норма и ∞-норма совпадают, различные варианты подпрограмм для каждой из этих норм не приводятся.

Этот алгоритм перенесен из библиотеки LAPACK.