Python least squares solve with complex matrices
#!/usr/bin/env python import os, sys import numpy as np from scipy import linalg A_real = np.random.random((20,5)) A_imag = np.random.random((20,5)) A = A_real + 1j*A_imag; b_real = np.random.random((20,1)) b_imag = np.random.random((20,1)) b = b_real + 1j*b_imag; c,resid,rank,sigma = np.linalg.lstsq(A, b) print c










