fix phase field init

This commit is contained in:
James E McClure 2018-09-04 21:40:44 -04:00
parent 11252ecb38
commit 50f26adcca
2 changed files with 2 additions and 2 deletions

View File

@ -2785,7 +2785,7 @@ extern "C" void ScaLBL_PhaseField_Init(int *Map, double *Phi, double *Den, doubl
}
else{
nA=0.5*(phi-1.f);
nB=0.5*(phi+1.f);
nB=0.5*(1.f-phi);
}
Den[idx] = nA;
Den[Np+idx] = nB;

View File

@ -3910,7 +3910,7 @@ __global__ void dvc_ScaLBL_PhaseField_Init(int *Map, double *Phi, double *Den, d
}
else{
nA=0.5*(phi-1.f);
nB=0.5*(phi+1.f);
nB=0.5*(1.f-phi);
}
Den[idx] = nA;
Den[Np+idx] = nB;