Всем привет!
У меня есть задача решить систему линейных уравнений в комплексных числах. С помощью специальных библиотек я это делаю, все считается. Самописное решение для double матриц тоже находит решение, раскладывая [B +Ci] в [B -C; C B].
Но вот вопрос - как хорошие быстрые библиотеки типа mkl решают комплексные матрицы? Тоже переводят их в вещественную форму или решают их прямо в комплексной форме? И тут я уже совершенно плаваю - можно ли такие матрицы раскладывать/решать просто заменяя арифметические операции над вещественными числами - соответствующими арифметическими операциями над комплексными числами?
Попробовал решалку из octave - там chol-разложение для матриц в комплексной форме (и с перестановками) работает в несколько раз быстрее, чем для «удвоенной» вещественной
P.S. Прошу прощения за «тупняк», исходники SuiteSparse уже скачал, пытаюсь разобраться (нужен еще и разреженный вариант)
