OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

membrane test

Browse Source
This commit is contained in:
James McClure 2021-12-28 17:21:11 -05:00
parent 74bf7dbec4
commit ec76dc9a3b
1 changed files with 11 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
tests/TestMembrane.cpp
View File

@ -59,11 +59,19 @@ int main(int argc, char **argv)
int N = Nx*Ny*Nz; int N = Nx*Ny*Nz;
//....................................................................... //.......................................................................
int Np = 0; int Np = 0;
for (k=1;k<Nz-1;k++){ double distance,radius;
for (j=1;j<Ny-1;j++){ DoubleArray Distance(Nx,Ny,Nz);
for (i=1;i<Nx-1;i++){ for (k=0;k<Nz;k++){
for (j=0;j<Ny;j++){
for (i=0;i<Nx;i++){
n = k*Nx*Ny+j*Nx+i; n = k*Nx*Ny+j*Nx+i;
Dm->id[n] = 1; Dm->id[n] = 1;
radius = double(Nx)/6;
distance = sqrt(double((i-0.5*Nx)*(i-0.5*Nx)+ (j-0.5*Ny)*(j-0.5*Ny)+ (k-0.5*Nz)*(k-0.5*Nz)))-radius;
if (distance < 0.0 ){
Dm->id[n] = 1;
}
Distance(i,j,k) = distance;
Np++; Np++;
} }
} }