Commit 96f22a36 authored by laura.gassner's avatar laura.gassner

rearranged all input files to match the order of the parameter description in the manual

parent 2f167825
......@@ -550,9 +550,34 @@ A step length search which is based on the wolfe condition can be used which the
However, it is possible to chose a second step length search which the L-BFGS method. This one is available by the switch LBFGS\_STEP\_LENGTH=0. If you use LBFGS\_STEP\_LENGTH the wolfe condition based step length search will be deactivated automatically. Whit this search the step length 0.1, 0.5 and 1.0 will be tried and with a parabolic fit the best step length will be estimated. This search algorithm is implemented to save computation time, however the wolfe condition will not be checked, so the L-BFGS stability will not be satisfied. It is highly recommended to use the step length search based on the wolfe condition. It is likely that this option will be removed in a future releas.
\section{Step length estimation}
{\color{blue}{\begin{verbatim}
"Step length estimation" : "comment",
"EPS_SCALE" : "0.01",
"STEPMAX" : "4",
"SCALEFAC" : "4.0",
"TESTSHOT_START , TESTSHOT_END , TESTSHOT_INCR" : "1 , 2 , 1",
\end{verbatim}}}
For the step length estimation a parabolic line search method proposed by \cite{sourbier:09,sourbier:09b}, \cite{brossier:2009} and \cite{nocedal:1999} is implemented. For this step length estimation only two further test forward modelings are needed. The vector L2t contains the misfit values and the vector epst contains the corresponding step length. During the forward modeling of an iteration step the misfit norm of the data residuals is calculated for the shots defined by TESTSHOT\_START, TESTSHOT\_END and TESTSHOT\_INC. The value L2t(1) then contains the misfit from the forward modeling and the corresponding epst(1) value is 0.0.\\
The step lengths for the different parameters are defined as:\\
EPSILON = EPS\_SCALE * m\_max/grad\_max
EPSILON = epst[i] * m\_max/grad\_max\\
where m\_max is the maximum value of the corresponding model parameter in the whole model and grad\_max is the maximum absolute value of the gradient.\\
For a better definition of the parabola the improved line search is now trying to estimate a steplength epst(2) with L2t(2)<L2t(1). If the code is not able to find an appropiate steplength using the user-defined value EPS\_SCALE (f.e. EPS\_SCALE = 0.01 = 1\% change in terms of m\_max/grad\_max), the code divides this steplength by the variable SCALEFAC and calculates the misfit norm again. If this search fails after STEPMAX attempts DENISE exits with an error message. If the algorithm has found an appropriate value for epst(2), it is trying to estimate a steplength epst(3) with L2t(3)> L2t(2), by increasing the steplength\\
EPS\_SCALE += EPS\_SCALE/SCALEFAC.\\
If a corresponding value epst(3) can be found after STEPMAX forward modellings, DENISE can fit a parabola through the 3 points (L2t(i),epst(i)) and estimates an optimum step length at the minimum of the parabola. If the L2-value L2t(3) after STEPMAX forward models is still smaller than L2t(2) the optimum steplength estimated by parabolic fitting will be not larger than epst(3).\\
Please note: This step length search is only available wenn PCG method (GRAD\_METHOD==1) is used and will be automatically selected.
\section{Misfit definition}
{\color{blue}{\begin{verbatim}
"Gradient calculation" : "comment",
"Misfit Definition" : "comment",
"LNORM" : "2",
"NORMALIZE" : "0",
"DTINV" : "2",
......@@ -582,31 +607,6 @@ If NORMALIZE is set to 1, the synthetic data and the measured data will be norma
To reduce the memory requirements during an inversion one can define that only every DTINV time sample is used for the calculation of the gradients. To set this parameter appropriately one has to keep in mind the Nyquist criterion to avoid aliasing effects.
\newpage
\section{Step length estimation}
{\color{blue}{\begin{verbatim}
"Step length estimation" : "comment",
"EPS_SCALE" : "0.01",
"STEPMAX" : "4",
"SCALEFAC" : "4.0",
"TESTSHOT_START , TESTSHOT_END , TESTSHOT_INCR" : "1 , 2 , 1",
\end{verbatim}}}
For the step length estimation a parabolic line search method proposed by \cite{sourbier:09,sourbier:09b}, \cite{brossier:2009} and \cite{nocedal:1999} is implemented. For this step length estimation only two further test forward modelings are needed. The vector L2t contains the misfit values and the vector epst contains the corresponding step length. During the forward modeling of an iteration step the misfit norm of the data residuals is calculated for the shots defined by TESTSHOT\_START, TESTSHOT\_END and TESTSHOT\_INC. The value L2t(1) then contains the misfit from the forward modeling and the corresponding epst(1) value is 0.0.\\
The step lengths for the different parameters are defined as:\\
EPSILON = EPS\_SCALE * m\_max/grad\_max
EPSILON = epst[i] * m\_max/grad\_max\\
where m\_max is the maximum value of the corresponding model parameter in the whole model and grad\_max is the maximum absolute value of the gradient.\\
For a better definition of the parabola the improved line search is now trying to estimate a steplength epst(2) with L2t(2)<L2t(1). If the code is not able to find an appropiate steplength using the user-defined value EPS\_SCALE (f.e. EPS\_SCALE = 0.01 = 1\% change in terms of m\_max/grad\_max), the code divides this steplength by the variable SCALEFAC and calculates the misfit norm again. If this search fails after STEPMAX attempts DENISE exits with an error message. If the algorithm has found an appropriate value for epst(2), it is trying to estimate a steplength epst(3) with L2t(3)> L2t(2), by increasing the steplength\\
EPS\_SCALE += EPS\_SCALE/SCALEFAC.\\
If a corresponding value epst(3) can be found after STEPMAX forward modellings, DENISE can fit a parabola through the 3 points (L2t(i),epst(i)) and estimates an optimum step length at the minimum of the parabola. If the L2-value L2t(3) after STEPMAX forward models is still smaller than L2t(2) the optimum steplength estimated by parabolic fitting will be not larger than epst(3).\\
Please note: This step length search is only available wenn PCG method (GRAD\_METHOD==1) is used and will be automatically selected.
\section{Abort criterion}
\label{json:abort_criterion}
......@@ -715,7 +715,7 @@ If you are using frequeny filtering (TIME\_FILT==1) during the inversion, you ca
\section{Data manipulation}
\subsection{Time windowing}
{\color{blue}{\begin{verbatim}
"Time windowing and damping" : "comment",
"Time windowing" : "comment",
"TIMEWIN" : "0",
"TW_IND" : "0",
"PICKS_FILE" : "./picked_times/picks"
......
......@@ -68,7 +68,6 @@
"SEIS_FILE_P" : "su/DENISE_p.su",
"General inversion parameters" : "comment",
"INVMAT1" : "1",
"INVMAT" : "10",
......
......@@ -6,7 +6,6 @@
#
{
"Domain Decomposition" : "comment",
"NPROCX" : "4",
"NPROCY" : "2",
......@@ -39,6 +38,9 @@
"Acoustic Computation" : "comment",
"ACOUSTIC" : "0",
"Wavetype Computation" : "comment",
"WAVETYPE" : "1",
"Model" : "comment",
"READMOD" : "0",
"MFILE" : "model/model_Test",
......
......@@ -68,7 +68,6 @@
"SEIS_FILE_P" : "su/DENISE_p.su",
"General inversion parameters" : "comment",
"ITERMAX" : "10",
"DATA_DIR" : "su/measured_data/DENISE_real",
......
......@@ -38,6 +38,9 @@
"Acoustic Computation" : "comment",
"ACOUSTIC" : "0",
"Wavetype Computation" : "comment",
"WAVETYPE" : "1",
"Model" : "comment",
"READMOD" : "0",
"MFILE" : "model/model_Test",
......@@ -102,48 +105,11 @@
"MISFIT_LOG_FILE" : "L2_LOG.dat",
"VELOCITY" : "0",
"Gradient-Method" : "comment",
"GRAD_METHOD" : "1",
"LBFGS_STEP_LENGTH" : "1",
"N_LBFGS" : "5",
"Wolfe Condition" : "comment",
"WOLFE_CONDITION" : "1",
"WOLFE_NUM_TEST" : "5",
"WOLFE_TRY_OLD_STEPLENGTH" : "0",
"WOLFE_C1_SL" : "1e-4",
"WOLFE_C2_SL" : "0.9",
"Approx. Hessian" : "comment",
"EPRECOND" : "3",
"EPSILON_WE" : "0.005",
"EPRECOND_ITER" : "0",
"EPRECOND_PER_SHOT" : "1",
"Workflow" : "comment",
"USE_WORKFLOW" : "0",
"FILE_WORKFLOW" : "workflow.txt",
"Inversion for density" : "comment",
"Inversion for parameter starting from iteration" : "comment",
"INV_RHO_ITER" : "0",
"Inversion for Vp" : "comment",
"INV_VP_ITER" : "0",
"Inversion for Vs" : "comment",
"INV_VS_ITER" : "0",
"Minimum Vp/Vs-ratio" : "comment",
"VP_VS_RATIO" : "0.0",
"Limited update of model parameters in reference to the starting model" : "comment",
"S" : "0",
"S_VS" : "0.0",
"S_VP" : "0.0",
"S_RHO" : "0.0",
"Output of inverted models" : "comment",
"INV_MODELFILE" : "model/modelTest",
"nfstart" : "1",
......@@ -154,12 +120,27 @@
"nfstart_jac" : "1",
"nf_jac" : "1",
"Gradient calculation" : "comment",
"LNORM" : "2",
"LNORM values: L1 norm=1, L2 norm=2, Cauchy=3, SECH=4, Global correlation=5, normalized L2 norm=7, Envelope-based norm=8" : "comment",
"NORMALIZE" : "0",
"DTINV" : "2",
"WATERLEVEL_LNORM8" : "0.0",
"Workflow" : "comment",
"USE_WORKFLOW" : "0",
"FILE_WORKFLOW" : "workflow.txt",
"Approx. Hessian" : "comment",
"EPRECOND" : "3",
"EPSILON_WE" : "0.005",
"EPRECOND_ITER" : "0",
"EPRECOND_PER_SHOT" : "1",
"Gradient-Method" : "comment",
"GRAD_METHOD" : "1",
"LBFGS_STEP_LENGTH" : "1",
"N_LBFGS" : "5",
"Wolfe Condition" : "comment",
"WOLFE_CONDITION" : "1",
"WOLFE_NUM_TEST" : "5",
"WOLFE_TRY_OLD_STEPLENGTH" : "0",
"WOLFE_C1_SL" : "1e-4",
"WOLFE_C2_SL" : "0.9",
"Step length estimation" : "comment",
"EPS_SCALE" : "0.01",
......@@ -167,13 +148,52 @@
"SCALEFAC" : "4.0",
"TESTSHOT_START , TESTSHOT_END , TESTSHOT_INCR" : "1 , 2 , 1",
"Misfit Definition" : "comment",
"LNORM" : "2",
"LNORM values: L1 norm=1, L2 norm=2, Cauchy=3, SECH=4, Global correlation=5, normalized L2 norm=7, Envelope-based norm=8" : "comment",
"NORMALIZE" : "0",
"DTINV" : "2",
"WATERLEVEL_LNORM8" : "0.0",
"Termination of the programmme" : "comment",
"PRO" : "0.01",
"Definition of inversion for source time function" : "comment",
"INV_STF" : "0",
"PARA" : "fdlsq:tshift=0.0",
"N_STF" : "10",
"N_STF_START" : "1",
"TAPER_STF" : "0",
"TRKILL_STF" : "0",
"TRKILL_FILE_STF" : "./trace_kill/trace_kill.dat",
"Frequency filtering during inversion" : "comment",
"TIME_FILT" : "0",
"F_HP" : "1",
"FC_START" : "10.0",
"FC_END" : "75.0",
"FC_INCR" : "10.0",
"ORDER" : "2",
"ZERO_PHASE" : "0",
"FREQ_FILE" : "frequencies.dat",
"Minimum number of iteration per frequency" : "comment",
"MIN_ITER" : "0",
"Definition of gradient taper geometry" : "comment",
"Time windowing" : "comment",
"TIMEWIN" : "0",
"TW_IND" : "0",
"PICKS_FILE" : "./picked_times/picks",
"TWLENGTH_PLUS" : "0.01",
"TWLENGTH_MINUS" : "0.01",
"GAMMA" : "100000",
"Trace killing" : "comment",
"TRKILL" : "0",
"TRKILL_FILE" : "./trace_kill/trace_kill.dat",
"Definition of a gradient taper" : "comment",
"SWS_TAPER_GRAD_VERT" : "0",
"SWS_TAPER_GRAD_HOR" : "0",
"GRADT1 , GRADT2 , GRADT3 , GRADT4" : "5 , 15 , 490 , 500",
......@@ -188,7 +208,7 @@
"TAPER_FILE_NAME_U" : "taper_u.bin",
"TAPER_FILE_NAME_RHO" : "taper_rho.bin",
"Definition of smoothing (spatial filtering) of the gradients" : "comment",
"Definition of spatial filtering of the gradients" : "comment",
"SPATFILTER" : "0",
"SPAT_FILT_SIZE" : "40",
"SPAT_FILT_1" : "1",
......@@ -208,41 +228,19 @@
"RHOUPPERLIM" : "5000",
"RHOLOWERLIM" : "0",
"Limited update of model parameters in reference to the starting model" : "comment",
"S" : "0",
"S_VS" : "0.0",
"S_VP" : "0.0",
"S_RHO" : "0.0",
"Minimum Vp/Vs-ratio" : "comment",
"VP_VS_RATIO" : "0.0",
"Definition of smoothing the models vp and vs" : "comment",
"MODEL_FILTER" : "0",
"FILT_SIZE" : "5",
"Definition of inversion for source time function" : "comment",
"INV_STF" : "0",
"PARA" : "fdlsq:tshift=0.0",
"N_STF" : "10",
"N_STF_START" : "1",
"TAPER_STF" : "0",
"TRKILL_STF" : "0",
"TRKILL_FILE_STF" : "./trace_kill/trace_kill.dat",
"Frequency filtering during inversion" : "comment",
"TIME_FILT" : "0",
"F_HP" : "1",
"FC_START" : "10.0",
"FC_END" : "75.0",
"FC_INCR" : "10.0",
"ORDER" : "2",
"ZERO_PHASE" : "0",
"FREQ_FILE" : "frequencies.dat",
"Trace killing" : "comment",
"TRKILL" : "0",
"TRKILL_FILE" : "./trace_kill/trace_kill.dat",
"Time windowing and damping" : "comment",
"TIMEWIN" : "0",
"TW_IND" : "0",
"PICKS_FILE" : "./picked_times/picks",
"TWLENGTH_PLUS" : "0.01",
"TWLENGTH_MINUS" : "0.01",
"GAMMA" : "100000",
"Verbose mode" : "comment",
"VERBOSE" : "0",
}
......@@ -40,6 +40,9 @@
"Acoustic Computation" : "comment",
"ACOUSTIC" : "0",
"Wavetype Computation" : "comment",
"WAVETYPE" : "1",
"Model" : "comment",
"READMOD" : "0",
"MFILE" : "model/mod_toy_example_true",
......@@ -79,7 +82,7 @@
"ITERMAX" : "1",
"INVMAT1" : "1",
"INVMAT" : "10",
"Verbose mode" : "comment",
"VERBOSE" : "0",
}
......@@ -40,26 +40,12 @@
"SOURCE_FILE" : "./source/src_toy_example.dat",
"RUN_MULTIPLE_SHOTS" : "1",
"Snapshots" : "comment",
"SNAP" : "0",
"TSNAP1" : "5e-5",
"TSNAP2" : "0.6",
"TSNAPINC" : "0.006",
"IDX" : "1",
"IDY" : "1",
"SNAP_FORMAT" : "3",
"SNAPSHOT_START , SNAPSHOT_END , SNAPSHOT_INCR" : "2, 3 , 2",
"SNAP_FILE" : "./snap/waveform_forward",
"Acoustic Computation" : "comment",
"ACOUSTIC" : "0",
"Wavetype Computation" : "comment",
"WAVETYPE" : "3",
"Verbose" : "comment",
"VERBOSE" : "0",
"Model" : "comment",
"READMOD" : "0",
"MFILE" : "model/mod_toy_example_true",
......@@ -96,10 +82,23 @@
"FL3" : "72.6774",
"TAU" : "0.0966",
"Snapshots" : "comment",
"SNAP" : "0",
"TSNAP1" : "5e-5",
"TSNAP2" : "0.6",
"TSNAPINC" : "0.006",
"IDX" : "1",
"IDY" : "1",
"SNAP_FORMAT" : "3",
"SNAPSHOT_START , SNAPSHOT_END , SNAPSHOT_INCR" : "2, 3 , 2",
"SNAP_FILE" : "./snap/waveform_forward",
"General inversion parameters" : "comment",
"ITERMAX" : "1",
"INVMAT1" : "1",
"INVMAT" : "10",
"Verbose" : "comment",
"VERBOSE" : "0",
}
......@@ -83,14 +83,11 @@
"QUELLTYPB" : "1",
"MISFIT_LOG_FILE" : "LOG_toy_example.dat",
"Workflow" : "comment",
"USE_WORKFLOW" : "0",
"FILE_WORKFLOW" : "workflow.txt",
"Inversion for density" : "comment",
"Inversion for parameter starting from iteration" : "comment",
"INV_RHO_ITER" : "400",
"INV_VP_ITER" : "400",
"INV_VS_ITER" : "0",
"Output of inverted models" : "comment",
"INV_MODELFILE" : "model/toy_example/mod_toy_example",
"nfstart" : "1",
......@@ -101,6 +98,15 @@
"nfstart_jac" : "1",
"nf_jac" : "1",
"Workflow" : "comment",
"USE_WORKFLOW" : "0",
"FILE_WORKFLOW" : "workflow.txt",
"Approx. Hessian" : "comment",
"EPRECOND" : "0",
"EPSILON_WE" : "0.005",
"EPRECOND_ITER" : "0",
"Gradient-Method" : "comment",
"GRAD_METHOD" : "1",
"LBFGS_STEP_LENGTH" : "1",
......@@ -111,12 +117,7 @@
"WOLFE_NUM_TEST" : "5",
"WOLFE_TRY_OLD_STEPLENGTH" : "1",
"Approx. Hessian" : "comment",
"EPRECOND" : "0",
"EPSILON_WE" : "0.005",
"EPRECOND_ITER" : "0",
"Gradient calculation" : "comment",
"Misfit Definition" : "comment",
"LNORM" : "7",
"LNORM values: L1 norm=1, L2 norm=2, Cauchy=3, SECH=4, Global correlation=5, L2 norm with time windowing=6" : "comment",
"NORMALIZE" : "0",
......@@ -131,6 +132,13 @@
"Termination of the programmme" : "comment",
"PRO" : "0.01",
"Frequency filtering during inversion" : "comment",
"TIME_FILT" : "1",
"FC_START" : "10.0",
"FC_END" : "70.0",
"FC_INCR" : "10.0",
"ORDER" : "4",
"Minimum number of iteration per frequency" : "comment",
"MIN_ITER" : "10",
......@@ -141,7 +149,6 @@
"FILTSIZE" : "1",
"SWS_TAPER_FILE" : "0",
"Upper and lower limits for model parameters" : "comment",
"VPUPPERLIM" : "2420",
"VPLOWERLIM" : "0",
......@@ -154,13 +161,6 @@
"MODEL_FILTER" : "1",
"FILT_SIZE" : "3",
"Frequency filtering during inversion" : "comment",
"TIME_FILT" : "1",
"FC_START" : "10.0",
"FC_END" : "70.0",
"FC_INCR" : "10.0",
"ORDER" : "4",
"Verbose mode" : "comment",
"VERBOSE" : "0",
}
......@@ -46,9 +46,6 @@
"Wavetype Computation" : "comment",
"WAVETYPE" : "2",
"Verbose" : "comment",
"VERBOSE" : "0",
"Model" : "comment",
"READMOD" : "0",
"MFILE" : "model_true/mod_toy_example",
......@@ -84,7 +81,7 @@
"FL2" : "7.6660",
"FL3" : "72.6774",
"TAU" : "0.0966",
"General inversion parameters" : "comment",
"ITERMAX" : "5",
......@@ -95,30 +92,11 @@
"QUELLTYPB values: 1=x y; 2=x; 3=y; 4=z; 5=x y z" : "comment",
"MISFIT_LOG_FILE" : "LOG_toy_example.dat",
"Gradient-Method" : "comment",
"GRAD_METHOD" : "1",
"LBFGS_STEP_LENGTH" : "1",
"N_LBFGS" : "5",
"Wolfe Condition" : "comment",
"WOLFE_CONDITION" : "1",
"WOLFE_NUM_TEST" : "5",
"WOLFE_TRY_OLD_STEPLENGTH" : "1",
"Approx. Hessian" : "comment",
"EPRECOND" : "3",
"EPSILON_WE" : "0.005",
"EPRECOND_ITER" : "0",
"Workflow" : "comment",
"USE_WORKFLOW" : "0",
"FILE_WORKFLOW" : "workflow.txt",
"Inversion for density" : "comment",
"Inversion for parameter starting from iteration" : "comment",
"INV_RHO_ITER" : "400",
"INV_VP_ITER" : "400",
"INV_VS_ITER" : "0",
"Output of inverted models" : "comment",
"INV_MODELFILE" : "model/toy_example/mod_toy_example",
"nfstart" : "1",
......@@ -129,7 +107,26 @@
"nfstart_jac" : "1",
"nf_jac" : "1",
"Gradient calculation" : "comment",
"Workflow" : "comment",
"USE_WORKFLOW" : "0",
"FILE_WORKFLOW" : "workflow.txt",
"Approx. Hessian" : "comment",
"EPRECOND" : "3",
"EPSILON_WE" : "0.005",
"EPRECOND_ITER" : "0",
"Gradient-Method" : "comment",
"GRAD_METHOD" : "1",
"LBFGS_STEP_LENGTH" : "1",
"N_LBFGS" : "5",
"Wolfe Condition" : "comment",
"WOLFE_CONDITION" : "1",
"WOLFE_NUM_TEST" : "5",
"WOLFE_TRY_OLD_STEPLENGTH" : "1",
"Misfit Definition" : "comment",
"LNORM" : "2",
"LNORM values: L1 norm=1, L2 norm=2, Cauchy=3, SECH=4, Global correlation=5, L2 norm with time windowing=6" : "comment",
"NORMALIZE" : "0",
......@@ -144,6 +141,22 @@
"Termination of the programmme" : "comment",
"PRO" : "0.01",
"Definition of inversion for source time function" : "comment",
"INV_STF" : "0",
"PARA" : "fdlsq:tshift=0.0",
"N_STF" : "10",
"N_STF_START" : "1",
"Frequency filtering during inversion" : "comment",
"TIME_FILT" : "0",
"FC_START" : "5.0",
"FC_END" : "5.0",
"FC_INCR" : "10.0",
"ORDER" : "4",
"Minimum number of iteration per frequency" : "comment",
"MIN_ITER" : "10",
"Definition of gradient taper geometry" : "comment",
"SWS_TAPER_CIRCULAR_PER_SHOT" : "1",
"SRTSHAPE" : "1",
......@@ -152,8 +165,6 @@
"SWS_TAPER_FILE" : "0",
"TAPER_FILE" : "taper/toy_example/taper_toy_example",
"Definition of inversion for source time function" : "comment", "INV_STF" : "0", "PARA" : "fdlsq:tshift=0.0", "N_STF" : "10", "N_STF_START" : "1",
"Upper and lower limits for model parameters" : "comment",
"VPUPPERLIM" : "2420",
"VPLOWERLIM" : "0",
......@@ -166,12 +177,6 @@
"MODEL_FILTER" : "1",
"FILT_SIZE" : "3",
"Minimum number of iteration per frequency" : "comment",
"MIN_ITER" : "10",
"Frequency filtering during inversion" : "comment",
"TIME_FILT" : "0",
"FC_START" : "5.0",
"FC_END" : "5.0",
"FC_INCR" : "10.0",
"ORDER" : "4",
"Verbose" : "comment",
"VERBOSE" : "0",
}
......@@ -91,5 +91,5 @@
"INVMAT" : "10",
"Verbose mode" : "comment",
"VERBOSE" : "0",
"VERBOSE" : "0",
}
......@@ -75,14 +75,39 @@
"FL3" : "72.6774",
"TAU" : "0.0966",
"General inversion parameters" : "comment",
"ITERMAX" : "100",
"ITERMAX" : "400",
"DATA_DIR" : "su/measured_data/toy_example_ac",
"INVMAT1" : "1",
"INVMAT" : "0",
"QUELLTYPB" : "4",
"MISFIT_LOG_FILE" : "LOG_toy_example_ac.dat",
"Inversion for parameter starting from iteration" : "comment",
"INV_RHO_ITER" : "400",
"INV_VP_ITER" : "1",
"Output of inverted models" : "comment",
"INV_MODELFILE" : "model/toy_example/mod_toy_example_ac",
"nfstart" : "1",
"nf" : "1",
"Output of gradients" : "comment",
"JACOBIAN" : "jacobian/toy_example/jac_toy_example_ac",
"nfstart_jac" : "1",
"nf_jac" : "1",
"Workflow" : "comment",
"USE_WORKFLOW" : "0",
"FILE_WORKFLOW" : "workflow.txt",
"Approx. Hessian" : "comment",
"EPRECOND" : "3",
"EPSILON_WE" : "0.005",
"EPRECOND_ITER" : "0",
"EPRECOND_PER_SHOT" : "0",
"Gradient-Method" : "comment",
"GRAD_METHOD" : "2",
"LBFGS_STEP_LENGTH" : "1",
......@@ -94,38 +119,13 @@
"WOLFE_TRY_OLD_STEPLENGTH" : "1",
"WOLFE_C1_SL" : "0.0",
"Approx. Hessian" : "comment",
"EPRECOND" : "3",
"EPSILON_WE" : "0.005",
"EPRECOND_ITER" : "0",
"EPRECOND_PER_SHOT" : "0",
"Workflow" : "comment",
"USE_WORKFLOW" : "0",
"FILE_WORKFLOW" : "workflow.txt",
"Inversion for density" : "comment",
"INV_RHO_ITER" : "100",
"INV_VP_ITER" : "1",
"Step length estimation" : "comment",
"EPS_SCALE" : "0.01",
"STEPMAX" : "4",
"SCALEFAC" : "5.0",
"TESTSHOT_START , TESTSHOT_END , TESTSHOT_INCR" : "1 , 5 , 2",
"Output of inverted models" : "comment",
"INV_MODELFILE" : "model/toy_example/mod_toy_example_ac",
"nfstart" : "1",
"nf" : "1",
"Output of gradients" : "comment",
"JACOBIAN" : "jacobian/toy_example/jac_toy_example_ac",
"nfstart_jac" : "1",
"nf_jac" : "1",
"Gradient calculation" : "comment",
"Misfit Definition" : "comment",
"LNORM" : "2",
"LNORM values: L1 norm=1, L2 norm=2, Cauchy=3, SECH=4, Global correlation=5, L2 norm with time windowing=6" : "comment",
"NORMALIZE" : "0",
......@@ -134,19 +134,6 @@
"Termination of the programmme" : "comment",
"PRO" : "0.01",
"Minimum number of iteration per frequency" : "comment",
"MIN_ITER" : "0",
"Definition of gradient taper geometry" : "comment",
"SWS_TAPER_CIRCULAR_PER_SHOT" : "0",
"SRTSHAPE" : "1",
"SRTRADIUS" : "3.0",
"FILTSIZE" : "1",
"SWS_TAPER_FILE" : "1",
"TAPER_FILE_NAME" : "taper_te_ac.bin",
"TAPER_FILE_NAME_U" : "taper_te_ac_u.bin",
"TAPER_FILE_NAME_RHO" : "taper_te_ac_rho.bin",
"Definition of inversion for source time function" : "comment",
"INV_STF" : "0",
"PARA" : "fdlsq:tshift=0.0",
......@@ -155,18 +142,6 @@
"TAPER_STF" : "