implementing beamquality example with FFT

program2D/main.cc

@@ -109,6 +109,10 @@ void spectrumPlaneWave(std::complex<double> (*aperture) (double x, double y),
           varFar.inversFFT();
       }
 
+void estimateBeamQuality(VariableFFT& varIn, double lambda)
+{
+
+}
 
 int main(int argc, char** argv) {
 	std::ofstream DATEI;
@@ -145,6 +149,7 @@ int main(int argc, char** argv) {
     
     VariableFFT varE(n,n,geo);
     VariableFFT varIn(varE);
+    VariableFFT varPhase(varE);