Function complex_gaussian_elimination 

pub fn complex_gaussian_elimination(
    a: &[Vec<Complex>],
    b: &[Complex],
) -> Option<Vec<Complex>>

Solve the complex linear system A·x = b using Gaussian elimination with partial pivoting.

a is an N×N matrix in row-major format (Vec of rows). Returns the solution vector x, or None if the system is singular.