Major Revision of DENISE-BE SH Module

daniel-koehn · Jan 29, 2023 · b4246e2 · b4246e2
1 parent 540e8c3
commit b4246e2
Show file tree

Hide file tree

Showing 31 changed files with 4,659 additions and 534 deletions.
diff --git a/include/fd.h b/include/fd.h
@@ -607,6 +607,11 @@ void outseis_SHres(struct seisSH *seisSH, struct seisSHfwi *seisSHfwi, int *recs
 
 void physics_SH();
 
+void PML_pro_SH(float * d_x, float * K_x, float * alpha_prime_x, float * a_x, float * b_x, 
+            float * d_x_half, float * K_x_half, float * alpha_prime_x_half, float * a_x_half, float * b_x_half,
+            float * d_y, float * K_y, float * alpha_prime_y, float * a_y, float * b_y, 
+            float * d_y_half, float * K_y_half, float * alpha_prime_y_half, float * a_y_half, float * b_y_half);
+
 void precond_SH(struct fwiSH *fwiSH, struct acq *acq, int nsrc, int ntr_glob, float ** taper_coeff, FILE *FP_GRAV);
 
 void prepare_update_s_visc_SH(float *etajm, float *etaip, float *peta, float **fipjp, float **pujp, 
@@ -623,7 +628,7 @@ void saveseis_glob_SH(FILE *fp, float **sectionvz, int  **recpos, int  **recpos_
 
 void sh(struct waveSH *waveSH, struct waveSH_PML *waveSH_PML, struct matSH *matSH, struct fwiSH *fwiSH, struct mpiPSV *mpiPSV, 
          struct seisSH *seisSH, struct seisSHfwi *seisSHfwi, struct acq *acq, float *hc, int ishot, int nshots, int nsrc_loc, 
-         int ns, int ntr, float **Ws, float **Wr, int hin, int *DTINV_help, int mode, MPI_Request * req_send, MPI_Request * req_rec);
+         int ns, int ntr, float **Ws, float **Wr, int hin, int *DTINV_help, int mode, MPI_Request * req_send, MPI_Request * req_rec);	  
 
 float step_length_est_sh(struct waveSH *waveSH, struct waveSH_PML *waveSH_PML, struct matSH *matSH, struct fwiSH *fwiSH, struct mpiPSV *mpiPSV, 
          struct seisSH *seisSH, struct seisSHfwi *seisSHfwi, struct acq *acq, float *hc, int iter, int nsrc, int ns, int ntr, int ntr_glob, float * epst1, 
@@ -775,6 +780,8 @@ double LU_decomp(double  **A, double *x, double *b,int n);
 
 float maximum_m(float **mat, int nx, int ny);
 
+void median_model(float ** waveconv, int filt_size);
+
 void median_src(float ** waveconv,float ** taper_coeff, float **srcpos, int nshots, int **recpos, int ntr, int iter, int sws);
 
 float minimum_m(float **mat, int nx, int ny);

diff --git a/src/SH/FD_SH.c b/src/SH/FD_SH.c
@@ -141,7 +141,7 @@ mem_SH(nseismograms,ntr,ns,fdo3,nd,buffsize);
 alloc_SH(&waveSH,&waveSH_PML);
 
 /* calculate damping coefficients for CPMLs (SH problem)*/
-if(FW>0){PML_pro(waveSH_PML.d_x, waveSH_PML.K_x, waveSH_PML.alpha_prime_x, waveSH_PML.a_x, waveSH_PML.b_x, waveSH_PML.d_x_half, waveSH_PML.K_x_half, waveSH_PML.alpha_prime_x_half, waveSH_PML.a_x_half, 
+if(FW>0){PML_pro_SH(waveSH_PML.d_x, waveSH_PML.K_x, waveSH_PML.alpha_prime_x, waveSH_PML.a_x, waveSH_PML.b_x, waveSH_PML.d_x_half, waveSH_PML.K_x_half, waveSH_PML.alpha_prime_x_half, waveSH_PML.a_x_half, 
                  waveSH_PML.b_x_half, waveSH_PML.d_y, waveSH_PML.K_y, waveSH_PML.alpha_prime_y, waveSH_PML.a_y, waveSH_PML.b_y, waveSH_PML.d_y_half, waveSH_PML.K_y_half, waveSH_PML.alpha_prime_y_half, 
                  waveSH_PML.a_y_half, waveSH_PML.b_y_half);
 }

diff --git a/src/SH/FD_grad_SH.c b/src/SH/FD_grad_SH.c
@@ -233,7 +233,7 @@ FD_GRAD_MAT = 2;
 alloc_SH(&waveSH,&waveSH_PML);
 
 /* calculate damping coefficients for CPMLs (SH problem)*/
-if(FW>0){PML_pro(waveSH_PML.d_x, waveSH_PML.K_x, waveSH_PML.alpha_prime_x, waveSH_PML.a_x, waveSH_PML.b_x, waveSH_PML.d_x_half, waveSH_PML.K_x_half, waveSH_PML.alpha_prime_x_half, waveSH_PML.a_x_half, 
+if(FW>0){PML_pro_SH(waveSH_PML.d_x, waveSH_PML.K_x, waveSH_PML.alpha_prime_x, waveSH_PML.a_x, waveSH_PML.b_x, waveSH_PML.d_x_half, waveSH_PML.K_x_half, waveSH_PML.alpha_prime_x_half, waveSH_PML.a_x_half, 
                  waveSH_PML.b_x_half, waveSH_PML.d_y, waveSH_PML.K_y, waveSH_PML.alpha_prime_y, waveSH_PML.a_y, waveSH_PML.b_y, waveSH_PML.d_y_half, waveSH_PML.K_y_half, waveSH_PML.alpha_prime_y_half, 
                  waveSH_PML.a_y_half, waveSH_PML.b_y_half);
 }