Added elastic VTI PSV RTM module

include/fd.h

@@ -258,10 +258,17 @@ float ** psi_vyy, float ** psi_vxy, float ** psi_vxxs, float ***pr, float ***pp,
 
 void alloc_matVTI(struct matVTI *matVTI);
 
+void ass_gradVTI(struct fwiPSV *fwiPSV, struct matVTI *matVTI, int iter);
+
 void checkfd_ssg_VTI(FILE *fp, float ** prho, float ** c11, float ** c13, float ** c33, float ** c44, float *hc);
 
 void FD_VTI();
 
+float grad_obj_VTI(struct wavePSV *wavePSV, struct wavePSV_PML *wavePSV_PML, struct matVTI *matVTI, struct fwiPSV *fwiPSV, struct mpiPSV *mpiPSV, 
+         struct seisPSV *seisPSV, struct seisPSVfwi *seisPSVfwi, struct acq *acq, float *hc, int iter, int nsrc, int ns, int ntr, int ntr_glob, int nsrc_glob, 
+         int nsrc_loc, int ntr_loc, int nstage, float **We, float **Ws, float **Wr, float ** taper_coeff, int hin, int *DTINV_help, 
+         MPI_Request * req_send, MPI_Request * req_rec);
+
 void matcopy_elastic_VTI(float ** rho, float ** pi, float ** u);
 
 void model_elastic_VTI(float  **  rho, float **  c11, float **  c13, float **  c33, float **  c44);
@@ -270,6 +277,8 @@ void physics_VTI();
 
 void readmod_elastic_VTI(float  **  rho, float **  c11, float **  c13, float **  c33, float **  c44);
 
+void RTM_VTI();
+
 void seismo_ssg_VTI(int lsamp, int ntr, int **recpos, float **sectionvx, 
 float **sectionvy, float **sectionp, float **sectioncurl, float **sectiondiv,
 float **vx, float **vy, float **sxx, float **syy, float *hc);

src/Makefile

@@ -115,12 +115,15 @@ DENISE += alloc_fwiPSV.c \
 
 # add VTI source files
 DENISE += alloc_matVTI.c \
+	  ass_gradVTI.c \
 	  checkfd_ssg_VTI.c \
 	  FD_VTI.c \
+	  grad_obj_VTI.c \
 	  matcopy_elastic_VTI.c \
 	  $(MODEL_VTI) \
 	  physics_VTI.c \
 	  readmod_elastic_VTI.c \
+	  RTM_VTI.c \
 	  seismo_ssg_VTI.c \
 	  snap_ssg_VTI.c \
 	  update_s_elastic_PML_VTI.c \