Skip to content

GitLab

Commit 68b7d14c authored by laura.gassner's avatar laura.gassner
Browse files

Merge branch 'develop' into velocity_acoustic

parents a1bfa1e9 7b3cc0ff
Hide whitespace changes
Inline Side-by-side
doc/latex/6_Parameter_Definition.tex
View file @ 68b7d14c
...@@ -487,12 +487,12 @@ Default values are: ...@@ -487,12 +487,12 @@ Default values are:
With the use of a workflow file, you can define different FWI stages. For instance, one FWI stage could refer to one corner frequency of a low-pass filter or to a specific time window. Every line in the workflow file corresponds to one FWI stage. To use a workflow file the switch USE\_WORKFLOW have to be set to 1. The algorithm will automatically change to the next line of the workflow file if the abort criterium of the current line is reached or if no step length could be found which reduces the misfit. The structure of the variables inside the workflow file is as follow: \\ With the use of a workflow file, you can define different FWI stages. For instance, one FWI stage could refer to one corner frequency of a low-pass filter or to a specific time window. Every line in the workflow file corresponds to one FWI stage. To use a workflow file the switch USE\_WORKFLOW have to be set to 1. The algorithm will automatically change to the next line of the workflow file if the abort criterium of the current line is reached or if no step length could be found which reduces the misfit. The structure of the variables inside the workflow file is as follow: \\
\noindent \noindent
\begin{tabular}{llllllllllll} \begin{tabular}{lllllllllllll}
\# & INV\_VS & INV\_VP & INV\_RHO & PRO & TIME\_FILT & FC & WAVETYPE & 0 & 0 & EPRECOND & EPSILON\_WE\\ \# & INV\_VS & INV\_VP & INV\_RHO & PRO & TIME\_FILT & F\_HIGH\_PASS & F\_LOW\_PASS & WAVETYPE & 0 & 0 & EPRECOND & EPSILON\_WE\\
\end{tabular} \end{tabular}
\ \\ \ \\
The first column is the number of the line. With INV\_* etc. you can activate the inversion for VS, VP or RHO, respectively. The abort criterium in percent for this FWI stage will be the declared in the variable PRO. With TIME\_FILT you can activate the frequency filtering with the corner frequency FC. WAVETYPE can be used to switch between PSV and SH modeling, however it is recommended to use PSV modeling (WAVETYPE==1) due to the current development on SH modeling. The following to zeros are placeholders for an upcoming update. With EPRECOND and EPSILON\_WE you can control the approx. Hessian. Please note, that all features which are used eg. TIME\_FILT (see section \ref{sec:filtering}) within the workflow have to be activated in the .JSON file. The first column is the number of the line. With INV\_* etc. you can activate the inversion for VS, VP or RHO, respectively. The abort criterium in percent for this FWI stage will be the declared in the variable PRO. With TIME\_FILT you can activate the frequency filtering with the corner frequency F\_HIGH\_PASS of the high-pass and F\_LOW\_PASS of the low-pass filter. WAVETYPE can be used to switch between PSV and SH modeling, however it is recommended to use PSV modeling (WAVETYPE==1) due to the current development on SH modeling. The following to zeros are placeholders for an upcoming update. With EPRECOND and EPSILON\_WE you can control the approx. Hessian. Please note, that all features which are used eg. TIME\_FILT (see section \ref{sec:filtering}) within the workflow have to be activated in the .JSON file.
For an example of a workflow file, have a look in the par/ folder. For an example of a workflow file, have a look in the par/ folder.
...@@ -840,7 +840,7 @@ For GRAD\_FILT\_WAVELENGTH = 1 (and TIME\_FILT=1) a new wavelength dependent fil ...@@ -840,7 +840,7 @@ For GRAD\_FILT\_WAVELENGTH = 1 (and TIME\_FILT=1) a new wavelength dependent fil
\begin{equation} \begin{equation}
\mbox{FILT\_SIZE\_GRAD} = \frac{V_{s,\text{average}}\cdot \text{A}}{\text{FC}} \mbox{FILT\_SIZE\_GRAD} = \frac{V_{s,\text{average}}\cdot \text{A}}{\text{FC}}
\end{equation} \end{equation}
where FC is the corner frequency of TIME\_FILT and A is a weighting factor. where F\_LOW\_PASS is the corner frequency of TIME\_FILT and A is a weighting factor.
\section{Model manipulation} \section{Model manipulation}
......
par/workflow.txt
View file @ 68b7d14c
# VS VP RHO PRO FIL FC WT J J EPRE EPSI # VS VP RHO PRO F_FILT F_HIGH F_LOW WT J J EPRE EPSI
1 1 0 0 0.01 0 0 2 0 0 0 0.005 1 1 0 0 0.01 0 0 0 2 0 0 0 0.005
2 1 0 0 0.01 0 0 2 0 0 0 0.005 2 1 0 0 0.01 0 0 0 2 0 0 0 0.005
3 1 0 0 0.01 0 0 2 0 0 0 0.005 3 1 0 0 0.01 0 0 0 2 0 0 0 0.005
4 1 0 0 0.01 0 0 2 0 0 0 0.005 4 1 0 0 0.01 0 0 0 2 0 0 0 0.005
src/IFOS2D.c
View file @ 68b7d14c
...@@ -952,6 +952,8 @@ int main(int argc, char **argv){ ...@@ -952,6 +952,8 @@ int main(int argc, char **argv){
gradient_optimization=1; gradient_optimization=1;
} }
/* Reset fail status of parabolic step length search */
step3=0;
} }
if (MYID==0){ if (MYID==0){
...@@ -2424,11 +2426,13 @@ int main(int argc, char **argv){ ...@@ -2424,11 +2426,13 @@ int main(int argc, char **argv){
muss = 0; muss = 0;
lamss = prho[j][i] * ppi[j][i] * ppi[j][i] - 2.0 * muss; lamss = prho[j][i] * ppi[j][i] * ppi[j][i] - 2.0 * muss;
waveconv_lam[j][i] = (1.0/(4.0 * (lamss+muss) * (lamss+muss))) * waveconv_lam[j][i]; if(!ACOUSTIC)
waveconv_lam[j][i] = (1.0/(4.0 * (lamss+muss) * (lamss+muss))) * waveconv_lam[j][i];
else
waveconv_lam[j][i] = (1.0/((lamss+muss) * (lamss+muss))) * waveconv_lam[j][i];
/* calculate Vp gradient */ /* calculate Vp gradient */
waveconv_shot[j][i] = 2.0 * ppi[j][i] * prho[j][i] * waveconv_lam[j][i]; waveconv_shot[j][i] = 2.0 * ppi[j][i] * prho[j][i] * waveconv_lam[j][i];
// waveconv_shot[j][i] = waveconv_lam[j][i];
} }
if(PARAMETERIZATION==2){ if(PARAMETERIZATION==2){
...@@ -3885,7 +3889,7 @@ int main(int argc, char **argv){ ...@@ -3885,7 +3889,7 @@ int main(int argc, char **argv){
if(MYID==0){ if(MYID==0){
printf("\n Did not find a step length which decreases the misfit.\n"); printf("\n Did not find a step length which decreases the misfit.\n");
} }
step3=0;
break; break;
} }
...@@ -3924,6 +3928,8 @@ int main(int argc, char **argv){ ...@@ -3924,6 +3928,8 @@ int main(int argc, char **argv){
PCG_iter_start=iter+1; PCG_iter_start=iter+1;
wolfe_SLS_failed=0; wolfe_SLS_failed=0;
step3=0;
} }
/* ------------------------------------------------- */ /* ------------------------------------------------- */
...@@ -3959,6 +3965,8 @@ int main(int argc, char **argv){ ...@@ -3959,6 +3965,8 @@ int main(int argc, char **argv){
wolfe_SLS_failed=0; wolfe_SLS_failed=0;
alpha_SL_old=1; alpha_SL_old=1;
step3=0;
} }
/* ------------------------------------------------- */ /* ------------------------------------------------- */
...@@ -3994,6 +4002,8 @@ int main(int argc, char **argv){ ...@@ -3994,6 +4002,8 @@ int main(int argc, char **argv){
wolfe_SLS_failed=0; wolfe_SLS_failed=0;
alpha_SL_old=1; alpha_SL_old=1;
step3=0;
} }
} }
......
src/apply_workflow.c
View file @ 68b7d14c
...@@ -31,6 +31,7 @@ void apply_workflow(float ** workflow,int workflow_lines,char workflow_header[ST ...@@ -31,6 +31,7 @@ void apply_workflow(float ** workflow,int workflow_lines,char workflow_header[ST
/* extern variables */ /* extern variables */
extern int INV_RHO_ITER,INV_VS_ITER,INV_VP_ITER; extern int INV_RHO_ITER,INV_VS_ITER,INV_VP_ITER;
extern int TIME_FILT,MYID; extern int TIME_FILT,MYID;
extern float F_HIGH_PASS;
extern float PRO; extern float PRO;
extern int WAVETYPE; extern int WAVETYPE;
extern float JOINT_INVERSION_PSV_SH_ALPHA_VS; extern float JOINT_INVERSION_PSV_SH_ALPHA_VS;
...@@ -88,42 +89,45 @@ void apply_workflow(float ** workflow,int workflow_lines,char workflow_header[ST ...@@ -88,42 +89,45 @@ void apply_workflow(float ** workflow,int workflow_lines,char workflow_header[ST
PRO=workflow[WORKFLOW_STAGE][5]; PRO=workflow[WORKFLOW_STAGE][5];
/* Frequency filtering */ /* Frequency filtering */
if(TIME_FILT==1) { if( TIME_FILT == 1 ) {
TIME_FILT=workflow[WORKFLOW_STAGE][6]; TIME_FILT=workflow[WORKFLOW_STAGE][6];
if(*F_LOW_PASS>workflow[WORKFLOW_STAGE][7]&&(workflow[WORKFLOW_STAGE][6]>0)) {
if(MYID==0)printf("\n Due to the abort criteriom F_LOW_PASS is already higher than specified in workflow\n"); if( TIME_FILT > 0 ) {
if(MYID==0)printf(" therefore instead of %.2f HZ F_LOW_PASS=%.2f HZ is used\n",workflow[WORKFLOW_STAGE][7],*F_LOW_PASS); if( F_HIGH_PASS != workflow[WORKFLOW_STAGE][7] ) *LBFGS_iter_start=*iter;
} else { F_HIGH_PASS=workflow[WORKFLOW_STAGE][7];
if(*F_LOW_PASS!=workflow[WORKFLOW_STAGE][7]) *LBFGS_iter_start=*iter;
*F_LOW_PASS=workflow[WORKFLOW_STAGE][7]; if( *F_LOW_PASS != workflow[WORKFLOW_STAGE][8] ) *LBFGS_iter_start=*iter;
*F_LOW_PASS=workflow[WORKFLOW_STAGE][8];
} }
} else { } else {
if(MYID==0&&(workflow[WORKFLOW_STAGE][6]>0))printf("\n TIME_FILT cannot be activated due to it is not activated in the JSON File \n"); if(MYID==0&&(workflow[WORKFLOW_STAGE][6]>0))printf("\n TIME_FILT cannot be activated due to it is not activated in the JSON File \n");
} }
/* Change of wavetype */ /* Change of wavetype */
if(wavetype_start!=3&&(WAVETYPE!=workflow[WORKFLOW_STAGE][8])){ if(wavetype_start!=3&&(WAVETYPE!=workflow[WORKFLOW_STAGE][9])){
if(MYID==0)printf("\n Sorry, change of WAVETYPE with workflow only possible if WAVETYPE==3 in *.json"); if(MYID==0)printf("\n Sorry, change of WAVETYPE with workflow only possible if WAVETYPE==3 in *.json");
if(MYID==0)printf("\n WAVETYPE will remain unchanged %i",WAVETYPE); if(MYID==0)printf("\n WAVETYPE will remain unchanged %i",WAVETYPE);
} else { } else {
/* detect change and reset some things */ /* detect change and reset some things */
if(WAVETYPE!=workflow[WORKFLOW_STAGE][8]) { if(WAVETYPE!=workflow[WORKFLOW_STAGE][9]) {
*change_wavetype_iter=*iter; *change_wavetype_iter=*iter;
*LBFGS_iter_start=*iter; *LBFGS_iter_start=*iter;
} }
WAVETYPE=workflow[WORKFLOW_STAGE][8]; WAVETYPE=workflow[WORKFLOW_STAGE][9];
} }
/* Joint inversion PSV and SH */ /* Joint inversion PSV and SH */
JOINT_INVERSION_PSV_SH_ALPHA_VS=workflow[WORKFLOW_STAGE][9]; JOINT_INVERSION_PSV_SH_ALPHA_VS=workflow[WORKFLOW_STAGE][10];
JOINT_INVERSION_PSV_SH_ALPHA_RHO=workflow[WORKFLOW_STAGE][10]; JOINT_INVERSION_PSV_SH_ALPHA_RHO=workflow[WORKFLOW_STAGE][11];
/* Approx. Hessian */ /* Approx. Hessian */
if(EPRECOND==0 && workflow[WORKFLOW_STAGE][11]!=0){ if(EPRECOND==0 && workflow[WORKFLOW_STAGE][12]!=0){
if(MYID==0) printf(" WARNING: EPRECOND have to be set >0 in JSON (if so, ignore this message)"); if(MYID==0) printf(" WARNING: EPRECOND have to be set >0 in JSON (if so, ignore this message)");
} }
EPRECOND=workflow[WORKFLOW_STAGE][11]; EPRECOND=workflow[WORKFLOW_STAGE][12];
EPSILON_WE=workflow[WORKFLOW_STAGE][12]; EPSILON_WE=workflow[WORKFLOW_STAGE][13];
if(*LBFGS_iter_start==*iter && GRAD_METHOD==2){ if(*LBFGS_iter_start==*iter && GRAD_METHOD==2){
if(MYID==0)printf("\n L-BFGS will be used from iteration %d on.",*LBFGS_iter_start+1); if(MYID==0)printf("\n L-BFGS will be used from iteration %d on.",*LBFGS_iter_start+1);
......
src/fd.h
View file @ 68b7d14c
...@@ -22,7 +22,7 @@ ...@@ -22,7 +22,7 @@
#define STRING_SIZE 74 #define STRING_SIZE 74
#define STRING_SIZE2 256 #define STRING_SIZE2 256
#define REQUEST_COUNT 4 #define REQUEST_COUNT 4
#define WORKFLOW_MAX_VAR 12 #define WORKFLOW_MAX_VAR 13
/* declaration of functions */ /* declaration of functions */
......
src/taper_grad.c
View file @ 68b7d14c
...@@ -504,11 +504,13 @@ void taper_grad(float ** waveconv,float ** taper_coeff, float **srcpos, int nsho ...@@ -504,11 +504,13 @@ void taper_grad(float ** waveconv,float ** taper_coeff, float **srcpos, int nsho
} }
fclose(fp_taper); fclose(fp_taper);
sprintf(modfile,"taper_coeff_file.bin"); if(VERBOSE == 1) {
writemod(modfile,taper_coeff,3); sprintf(modfile,"taper_coeff_file.bin");
writemod(modfile,taper_coeff,3);
MPI_Barrier(MPI_COMM_WORLD); MPI_Barrier(MPI_COMM_WORLD);
if (MYID==0) mergemod(modfile,3); if (MYID==0) mergemod(modfile,3);
}
} }
} }
\ No newline at end of file
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Imprint