updated TestBlobAnalyze

This commit is contained in:
James E McClure
2015-06-05 12:56:29 -04:00
parent fb67c949a9
commit bf3d5616a9
2 changed files with 8 additions and 1 deletions

View File

@@ -549,7 +549,7 @@ void TwoPhase::ComputeLocal(){
void TwoPhase::ComputeLocalBlob(){ void TwoPhase::ComputeLocalBlob(){
int i,j,k,n,label; int i,j,k,n,label;
double vF,vS; double vF,vS;
vF = 0.0; vS= -1.0; vF = 0.0; vS=0.0;
// const RankInfoStruct rank_info(Dm.rank,Dm.nprocx,Dm.nprocy,Dm.nprocz); // const RankInfoStruct rank_info(Dm.rank,Dm.nprocx,Dm.nprocy,Dm.nprocz);
int cube[8][3] = {{0,0,0},{1,0,0},{0,1,0},{1,1,0},{0,0,1},{1,0,1},{0,1,1},{1,1,1}}; int cube[8][3] = {{0,0,0},{1,0,0},{0,1,0},{1,1,0},{0,0,1},{1,0,1},{0,1,1},{1,1,1}};

View File

@@ -207,6 +207,13 @@ int main(int argc, char **argv)
else{ else{
Dm.id[n] = 1; Dm.id[n] = 1;
} }
Averages.SDn(i,j,k) = Averages.Phase(i,j,k);
Averages.Phase_tplus(i,j,k) = Averages.Phase(i,j,k);
Averages.Phase_tminus(i,j,k) = Averages.Phase(i,j,k);
Averages.Press(i,j,k) = 0.0;
Averages.Vel_x(i,j,k) = 0.0;
Averages.Vel_y(i,j,k) = 0.0;
Averages.Vel_z(i,j,k) = 0.0;
} }
} }
} }