Unverified Commit e368ae9f authored by Axel Kohlmeyer's avatar Axel Kohlmeyer
Browse files

fix bug in pair styles lubricate and lubricate/poly 

Closes #1933
parent e2b7054f

src/COLLOID/pair_lubricate.cpp

+1 −1
Original line number Diff line number Diff line
@@ -142,7 +142,7 @@ void PairLubricate::compute(int eflag, int vflag) 
    Ef[2][2] = h_rate[2]/domain->zprd;

    Ef[0][1] = Ef[1][0] = 0.5 * h_rate[5]/domain->yprd;

    Ef[0][2] = Ef[2][0] = 0.5 * h_rate[4]/domain->zprd;

    Ef[1][2] = Ef[2][1] = 0.5 * h_rate[3]/domain->zprd;

    Ef[1][2] = Ef[2][1] = 0.5 * h_rate[3]/domain->xprd;



    // copy updated velocity/omega/angmom to the ghost particles

    // no need to do this if not shearing since comm->ghost_velocity is set

src/COLLOID/pair_lubricate_poly.cpp

+1 −1
Original line number Diff line number Diff line
@@ -124,7 +124,7 @@ void PairLubricatePoly::compute(int eflag, int vflag) 
    Ef[2][2] = h_rate[2]/domain->zprd;

    Ef[0][1] = Ef[1][0] = 0.5 * h_rate[5]/domain->yprd;

    Ef[0][2] = Ef[2][0] = 0.5 * h_rate[4]/domain->zprd;

    Ef[1][2] = Ef[2][1] = 0.5 * h_rate[3]/domain->zprd;

    Ef[1][2] = Ef[2][1] = 0.5 * h_rate[3]/domain->xprd;



    // copy updated omega to the ghost particles

    // no need to do this if not shearing since comm->ghost_velocity is set

CRA Git | Maintained and supported by SUSTech CRA and CCSE