Сижу уже пару часов, пытаюсь нагуглить нормальный вариант решения системы уравнений с комплексными числами: Ax=b, где A — прямоугольная комплексная матрица (ортонормированный базис), x — искомый вектор действительных чисел (коэффициенты разложения b по базису A), b — раскладываемый ряд.
«Просто» разложить не могу, т.к. мало данных (ошибка зверская вылезает), поэтому хочу попробовать метод наименьших квадратов (на том же наборе точек для действительного вектора b он дал превосходные результаты).
Сразу увидел, что применяемый мною метод QR-разложения в GSL не имеет комплексного аналога. Взялся за LU-разложение, но на этапе компиляции увидел «gsl: luc.c:64: ERROR: LU decomposition requires square matrix» — матрица-то должна быть квадратной (а я что-то в мануале на это внимания не обратил).
Вопрос: есть ли в GSL (ну, пусть даже не GSL, а BLAS или даже LAPACK) методики решения таких уравнений?
