Commit 56e84ac5 authored by Tilman Steinweg's avatar Tilman Steinweg
Browse files

fixed copy and paste bug

fixed copy and paste bug in update function of the cpml boundaries.
parent f83a7edb
...@@ -136,11 +136,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n ...@@ -136,11 +136,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n
if((POS[3]==0) && (k<=FW)){ if((POS[3]==0) && (k<=FW)){
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k];} vzz = vzz / K_z[k] + psi_vzz[j][i][k];}
if((POS[3]==NPROCZ-1) && (k>=nz2+1)){ if((POS[3]==NPROCZ-1) && (k>=nz2+1)){
...@@ -273,11 +273,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n ...@@ -273,11 +273,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n
if((POS[3]==0) && (k<=FW)){ if((POS[3]==0) && (k<=FW)){
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k];} vzz = vzz / K_z[k] + psi_vzz[j][i][k];}
if((POS[3]==NPROCZ-1) && (k>=nz2+1)){ if((POS[3]==NPROCZ-1) && (k>=nz2+1)){
...@@ -389,11 +389,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n ...@@ -389,11 +389,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n
if((POS[3]==0) && (k<=FW)){ if((POS[3]==0) && (k<=FW)){
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k];} vzz = vzz / K_z[k] + psi_vzz[j][i][k];}
if((POS[3]==NPROCZ-1) && (k>=nz2+1)){ if((POS[3]==NPROCZ-1) && (k>=nz2+1)){
...@@ -506,11 +506,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n ...@@ -506,11 +506,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n
if((POS[3]==0) && (k<=FW)){ if((POS[3]==0) && (k<=FW)){
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k];} vzz = vzz / K_z[k] + psi_vzz[j][i][k];}
if((POS[3]==NPROCZ-1) && (k>=nz2+1)){ if((POS[3]==NPROCZ-1) && (k>=nz2+1)){
...@@ -613,11 +613,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n ...@@ -613,11 +613,11 @@ double update_s_CPML(int nx1, int nx2, int ny1, int ny2, int nz1, int nz2, int n
vzz = (b1*(vz[j][i][k]-vz[j][i][k-1])+b2*(vz[j][i][k+1]-vz[j][i][k-2]))/DZ; vzz = (b1*(vz[j][i][k]-vz[j][i][k-1])+b2*(vz[j][i][k+1]-vz[j][i][k-2]))/DZ;
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k]; vzz = vzz / K_z[k] + psi_vzz[j][i][k];
......
...@@ -137,11 +137,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz ...@@ -137,11 +137,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz
if((POS[3]==0) && (k<=FW)){ if((POS[3]==0) && (k<=FW)){
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k];} vzz = vzz / K_z[k] + psi_vzz[j][i][k];}
if((POS[3]==NPROCZ-1) && (k>=nz2+1)){ if((POS[3]==NPROCZ-1) && (k>=nz2+1)){
...@@ -231,11 +231,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz ...@@ -231,11 +231,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz
if((POS[3]==0) && (k<=FW)){ if((POS[3]==0) && (k<=FW)){
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k];} vzz = vzz / K_z[k] + psi_vzz[j][i][k];}
if((POS[3]==NPROCZ-1) && (k>=nz2+1)){ if((POS[3]==NPROCZ-1) && (k>=nz2+1)){
...@@ -304,11 +304,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz ...@@ -304,11 +304,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz
if((POS[3]==0) && (k<=FW)){ if((POS[3]==0) && (k<=FW)){
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k];} vzz = vzz / K_z[k] + psi_vzz[j][i][k];}
if((POS[3]==NPROCZ-1) && (k>=nz2+1)){ if((POS[3]==NPROCZ-1) && (k>=nz2+1)){
...@@ -377,11 +377,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz ...@@ -377,11 +377,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz
if((POS[3]==0) && (k<=FW)){ if((POS[3]==0) && (k<=FW)){
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k];} vzz = vzz / K_z[k] + psi_vzz[j][i][k];}
if((POS[3]==NPROCZ-1) && (k>=nz2+1)){ if((POS[3]==NPROCZ-1) && (k>=nz2+1)){
...@@ -440,11 +440,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz ...@@ -440,11 +440,11 @@ double update_s_CPML_elastic(int nx1, int nx2, int ny1, int ny2, int nz1, int nz
vzz = (b1*(vz[j][i][k]-vz[j][i][k-1])+b2*(vz[j][i][k+1]-vz[j][i][k-2]))/DZ; vzz = (b1*(vz[j][i][k]-vz[j][i][k-1])+b2*(vz[j][i][k+1]-vz[j][i][k-2]))/DZ;
psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz; psi_vxz[j][i][k] = b_z_half[k] * psi_vxz[j][i][k] + a_z_half[k] * vxz;
vxz = vxz / K_y_half[k] + psi_vxz[j][i][k]; vxz = vxz / K_z_half[k] + psi_vxz[j][i][k];
psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz; psi_vyz[j][i][k] = b_z_half[k] * psi_vyz[j][i][k] + a_z_half[k] * vyz;
vyz = vyz / K_y_half[k] + psi_vyz[j][i][k]; vyz = vyz / K_z_half[k] + psi_vyz[j][i][k];
psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz; psi_vzz[j][i][k] = b_z[k] * psi_vzz[j][i][k] + a_z[k] * vzz;
vzz = vzz / K_y[k] + psi_vzz[j][i][k]; vzz = vzz / K_z[k] + psi_vzz[j][i][k];
......
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