Текст подпрограммы и версий ( Фортран ) afh6r.zip |
Тексты тестовых примеров ( Фортран ) tafh6r.zip |
Текст подпрограммы и версий ( Си ) afh6r_c.zip |
Тексты тестовых примеров ( Си ) tafh6r_c.zip |
Текст подпрограммы и версий ( Паскаль ) afh6r_p.zip |
Тексты тестовых примеров ( Паскаль ) tafh6r_p.zip |
Символическое треугольное разложение симметричной положительно определенной разреженной матрицы, заданной в формате RR (U) U.
Описание формата RR (U) U приведено в описании подпрограммы AM21R .
Пусть симметричная положительно определенная разреженная матрица А порядка N задана в формате RR (U) U. Подпрограмма AFH6R выполняет символическое треугольное разложение матрицы А, т.е. формирует портрет матрицы U в разложении матрицы А в произведение трех матриц: A = UTDU, где U - верхняя треугольная матрица с единичной диагональю, а D - диагональная матрица. Портрет результирующей матрицы U формируется также в формате RR (U) U.
С.Писсанецки. Технология разреженных матриц. - М.: Мир, 1988.
SUBROUTINE AFH6R ( IA, JA, N, IU, JU, IP)
Параметры
IA, JA - | заданный портрет матрицы А; |
N - | заданный порядок матрицы А (тип: целый); |
IU, JU - | вычисленный портрет матрицы U; |
IP - | целый одномерный массив длины N, используемый в подпрограмме в качестве рабочего |
Версии: нет
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION IA(6), JA(4), IU(6), JU(4), IP(5) DATA IA /1, 2, 3, 4, 5, 5/ DATA JA /5, 5, 5, 5/ N = 5 CALL AFH6R (IA, JA, N, IU, JU, IP) Результаты: IU = ( 1, 2, 3, 4, 5, 5 ) JU = ( 5, 5, 5, 5 )