test dphi_para * (2pi/lambda)^2
changing cMat to D+