debugging SSO

This commit is contained in:
James E McClure
2015-06-04 07:26:31 -04:00
parent 965257e9cd
commit 57d3e282cb
3 changed files with 12 additions and 3 deletions

View File

@@ -803,8 +803,12 @@ inline void SSO(DoubleArray &Distance, char *ID, Domain &Dm, int timesteps){
for (q=0; q<Q; q++){
weights[q] = sqrt(1.0*(D3Q27[q][0]*D3Q27[q][0]) + 1.0*(D3Q27[q][1]*D3Q27[q][1]) + 1.0*(D3Q27[q][2]*D3Q27[q][2]));
}
const RankInfoStruct rank_info(Dm.rank,Dm.nprocx,Dm.nprocy,Dm.nprocz);
fillHalo<double> fillData(rank_info,Dm.Nx-2,Dm.Ny-2,Dm.Nz-2,1,1,1,0,1);
int xdim,ydim,zdim;
xdim=Dm.Nx-2;
ydim=Dm.Ny-2;
zdim=Dm.Nz-2;
fillHalo<double> fillData(Dm.rank_info,xdim,ydim,zdim,1,1,1,0,1);
int count = 0;
while (count < timesteps){

View File

@@ -161,7 +161,8 @@ int main(int argc, char **argv)
readRankData( rank, nx+2, ny+2, nz+2, Phase, SignDist );
// Communication the halos
fillHalo<double> fillData(Dm.rank_info,nx,ny,nz,1,1,1,0,1);
const RankInfoStruct rank_info(rank,nprocx,nprocy,nprocz);
fillHalo<double> fillData(rank_info,nx,ny,nz,1,1,1,0,1);
fillData.fill(Phase);
fillData.fill(SignDist);

View File

@@ -45,6 +45,10 @@ int main(int argc, char **argv)
// Get the rank info
Domain Dm(nx,ny,nz,rank,nprocx,nprocy,nprocz,Lx,Ly,Lz,BC);
if (rank==0){
printf()
}
nx+=2; ny+=2; nz+=2;
int count = 0;