Skip to content

SOFI2D
SOFI2D
Commit c1812267 authored by Tilman Steinweg's avatar Tilman Steinweg
Browse files
Options

Cleanup

parent 3b01f88e
Hide whitespace changes
Inline Side-by-side
Showing with 12 additions and 11 deletions
src/exchange_s.c
View file @ c1812267
......@@ -31,8 +31,9 @@ void exchange_s(int nd, float ** sxx, float ** syy,
MPI_Request *req_send, MPI_Request *req_rec){
extern int NX, NY, POS[3], NPROCX, NPROCY, BOUNDARY, FDORDER;
extern int INDEX[5],MYID;
extern int NX, NY, POS[3], NPROCX, NPROCY, BOUNDARY;
// extern int FDORDER;
extern int INDEX[5];
extern const int TAG1,TAG2,TAG5,TAG6;
MPI_Status status;
int i, j, fdo, fdo2, n, l;
......
src/exchange_v.c
View file @ c1812267
......@@ -31,7 +31,8 @@ void exchange_v(int nd, float ** vx, float ** vy,
MPI_Request * req_send, MPI_Request * req_rec){
extern int NX, NY, POS[3], NPROCX, NPROCY, BOUNDARY, FDORDER;
extern int NX, NY, POS[3], NPROCX, NPROCY, BOUNDARY;
// extern int FDORDER;
extern int INDEX[5];
extern const int TAG1,TAG2,TAG5,TAG6;
MPI_Status status;
......
src/psource.c
View file @ c1812267
......@@ -26,8 +26,9 @@
void psource(int nt, float ** sxx, float ** syy,
float ** srcpos_loc, float ** signals, int nsrc){
extern float DH, DT;
extern int RSG, NT;
extern float DH;
// extern int RSG;
extern int NT;
int i, j, l;
float amp=0;
......
src/surface.c
View file @ c1812267
......@@ -29,14 +29,14 @@ float ** pi, float ** u, float ** taup,
float ** taus, float * etajm, float * peta, float *hc, float * K_x, float * a_x, float * b_x, float ** psi_vxx){
int i,j,l,m,h,h1,fdoh;
int i,j,l,m,h1,fdoh;
float bjm, djm, e, fjm, g, dhi;
float vxx, vyy, sump=0.0;
float dthalbe;
extern float DT, DH;
extern int NX, L, FDORDER;
extern int FW, BOUNDARY;
extern int NPROCX, NPROCY, POS[3], MYID;
extern int NPROCX,POS[3];
extern int ABS_TYPE;
fdoh = FDORDER/2;
dhi = 1.0/DH;
......@@ -98,7 +98,6 @@ float ** taus, float * etajm, float * peta, float *hc, float * K_x, float * a_x,
if((!BOUNDARY) && (POS[1]==NPROCX-1) && (i>=NX-FW+1)){
h1 = (i-NX+2*FW);
h = i;
psi_vxx[j][h1] = b_x[h1] * psi_vxx[j][h1] + a_x[h1] * vxx;
vxx = vxx / K_x[h1] + psi_vxx[j][h1];
......
src/surface_elastic.c
View file @ c1812267
......@@ -27,13 +27,13 @@ void surface_elastic(int ndepth, int * gx, float ** vx, float ** vy, float ** sx
float **sxy, float **pi, float **u, float *hc, float * K_x, float * a_x, float * b_x, float ** psi_vxx){
int i,j,m,fdoh,h,h1;
int i,j,m,fdoh,h1;
float fjm, g, dhi;
float vxx, vyy;
extern float DT, DH;
extern int NX, FDORDER;
extern int FW, BOUNDARY;
extern int NPROCX, NPROCY, POS[3], MYID;
extern int NPROCX, POS[3];
extern int ABS_TYPE;
fdoh = FDORDER/2;
dhi = 1.0/DH;
......@@ -77,7 +77,6 @@ float **sxy, float **pi, float **u, float *hc, float * K_x, float * a_x, float
if((!BOUNDARY) && (POS[1]==NPROCX-1) && (i>=NX-FW+1)){
h1 = (i-NX+2*FW);
h = i;
psi_vxx[j][h1] = b_x[h1] * psi_vxx[j][h1] + a_x[h1] * vxx;
vxx = vxx / K_x[h1] + psi_vxx[j][h1];
......
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