Inversion for C2-Benchmark with taup parameterization works: config...

Inversion for C2-Benchmark with taup parameterization works: config FWI/conf/fwiVAcoustic_inversion_C2.conf

FWI/conf/fwiVAcoustic_inversion_C2.conf

+ForwardSolver = VADGTS;
+Distribution = RCB;
+Model= FWIInverse_DGTS;
+Mesh = SeismogramGrid
+#Mesh = square2000x2000_20x20;
+#Mesh = square2000x2000;
+Problem = VAcousticWaveC2Benchmark;
+############################# Receiver Positions - BenchmarkC2 setup
+FWIReceiverCountY = 64;
+FWIReceiverCountX = 1;
+###############################
+FWIReceiverFirstXPosition = 1351;
+FWIReceiverFirstYPosition = 371;
+FWIReceiverLastYPosition  = 1631;
+centerx = 599;
+centery = 369;
+f_source = 25;
+numL =1;
+#source_factor = 100000;
+measurement_space = cellavg;
+measurement_time = point;
+compMeasure = 0;
+sourceSpace = 1; # cellavg
+source_component = 2;
+
+allowPlotting_M =1;
+allowPlotting_WF =0;
+plotfactor = 10;
+#parametrization = native;
+parametrization = taup;
+T=0.8;
+dt = 0.008;
+#############################
+shots = 1;
+
+level  =1;
+plevel =1;
+deg=1;
+
+############################# Time Disc Settings
+TimeVerbose = 1;
+plot_tStep = 100000;
+rkorder = -2;
+Kmax =150;
+Keps = 1e-07;
+Overlap = dG1;
+debugAdjoint = 0;
+debugFrechet = 0;
+
+############################# Solver Settings
+#LinearSolver = CG;
+LinearSolver = GMRES;
+LinearVerbose = -1;
+LinearSteps =100;
+LinearNthIterOutput = 1;
+
+############################# unused Settings (so far)
+FWINonlinearSolver = FWISteepestDescent;
+#FWINonlinearSolver = FWI_CGREGINN;
+#FWIReconstructRho    = 0;
+#FWIReconstructKappa0 = 1;
+#FWIReconstructTaus   = 0;
+FWI_DebugPlottingLevel = 3;
+#weight_rho = 0.005;
+#weight_kappa_0 = 2.42574e+10;
+#weight_kappa_1 = 4.85149e+08;
+#weight_tau_1 = 0.006;
+FWIReconstructRho = 0;
+FWIReconstructV_P = 1;
+FWIReconstructTauP = 0;
+
+source_factor =1;
+FWI_STDCVerbose = 2;
+FWI_STDCSteps = 5;
+FWI_STDCLineSearchSteps = 0;
+FWI_STDCLineSearchFactor = 0.3333333333333333;
+
+FWI_CGREGINNVerbose = 2;
+FWI_CGREGINNSteps = 5;
+FWI_CGREGINN_MaxInnerLoopSteps = 10;
+FWI_CGREGINN_ThetaMax = 0.999
+FWI_CGREGINN_Gamma = 0.9;
+FWIMeasurementKernel = ST_L1;
+
+LinearSteps = 10000;
+#loadconf = linearsteps.conf
+
+StdDamping = 0;
+clearcache = 0;
+# for i in $(seq -w 1 9 96); do echo "LinearSteps = $i;\n" > FWI/conf/linearsteps.conf; ./fwi; mv data/vtk/fwi_FOSLS_forward.vtk data/vtk/fwi_FOSLS_forward_CG$i.vtk; done
+
+SeismogramToPlot = data/FWI/test_seismogram.txt
+SeismogramToPlot = data/FWI/test_seismogram0
+SeismogramToPlot = data/FWI/test_seismogram_linesearch0
+