Блог

Я хочу решить Y = X * A для A (где Y-NxO, X-NXi, а A-IxO с N>>>O,I)

Я хочу использовать инверсию Мура-Пенроуза A = inv(транспонировать(X) * X) * (транспонировать(X) * Y)

квадратная матрица: (транспонирование(X) * X) является IxI, симметричной и полуположительной определенной.

Я думал об использовании DSYSV, который решает A*X=B для X, где A симметрично и положительно определено. Я мог бы заменить транспонирование(X)*X на A, а транспонирование(X)*Y на B, тогда решение X будет матрицей коэффициентов, которую я назвал выше.

Мне интересно, существует ли более эффективная комбинация процедур драйвера LAPACK для выполнения этой операции?