fix phase field init

This commit is contained in:
James E McClure 2018-09-04 22:00:33 -04:00
parent b04120e8f0
commit b4936bfa28

View File

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