Function complex_LU_solve 

pub fn complex_LU_solve(
    lu: &MatrixComplexF64,
    p: &Permutation,
    b: &VectorComplexF64,
    x: &mut VectorComplexF64,
) -> Result<(), Value>

This function solves the square system A x = b using the LU decomposition of A into (LU, p) given by LU_decomp or LU_decomp as input.