nag_complex_lu_solve_mult_rhs (f04akc) calculates the approximate solution of a set of complex linear equations with multiple right-hand sides
has been factorized by nag_complex_lu (f03ahc)
To solve a set of complex linear equations
, the function must be preceded by a call to nag_complex_lu (f03ahc)
which computes an
with partial pivoting,
is a permutation matrix,
is lower triangular and
is unit upper triangular. The columns
of the solution
are found by forward and backward substitution in
is a column of the right-hand side.
The accuracy of the computed solution depends on the conditioning of the original matrix. For a detailed error analysis see page 106 of Wilkinson and Reinsch (1971)
The time taken by nag_complex_lu_solve_mult_rhs (f04akc) is approximately proportional to .
To solve the set of linear equations