wrote several complex mult routines, changed B matrices to eps vectors, confirmed residual = 1e-10 using complex pointwise mult of eps and psi