try cg

FWI/src/ForwardModels/ForwardModelAAO.cpp

@@ -336,7 +336,7 @@ void ForwardModelAAO::ApplyPrecon(const FWIDataContainerDomain &m,
   const Vector &NonlinearPointMat = m.getMat();
 
   Vectors &UpdateWave = m_update.getVectorsWrite();
-
+  mout << "First AWE: " << sqrt(L2TScalarAAO(UpdateWave, UpdateWave)) << endl;
   bool usePenaltyForPrecon = true;
   Config::Get("usePenaltyForPrecon", usePenaltyForPrecon);
   int CGImprovement = 0;
@@ -367,7 +367,7 @@ void ForwardModelAAO::ApplyPrecon(const FWIDataContainerDomain &m,
     Vectors WaveTmp2(UpdateWave);
     WaveTmp2.Clear();
     AdjointMeasurement(meas, WaveTmp2);
-    mout << " residual to AWE: " << sqrt(L2TScalarAAO(WaveTmp2, WaveTmp2)) << endl;
+    mout << " residual to  AWE second: " << sqrt(L2TScalarAAO(WaveTmp2, WaveTmp2)) << endl;
     WaveTmp2.Clear();
     DGSeismogram result =
         ApplyDataSpaceCG(WavesolutionPrecon, meas, WaveTmp, WaveTmp2, CGImprovement);
@@ -386,7 +386,7 @@ void ForwardModelAAO::ApplyPrecon(const FWIDataContainerDomain &m,
     Obs->MeasureFullWave(UpdateWave, meas);
     WaveTmp.Clear();
     AdjointMeasurement(meas, WaveTmp);
-    mout << " residual to AWE: " << sqrt(L2TScalarAAO(WaveTmp, WaveTmp)) << endl;
+    mout << " residual to AWE third: " << sqrt(L2TScalarAAO(WaveTmp, WaveTmp)) << endl;
   }
 }