Skip to content

Re: multicomponent with different densities

#2737
mathias
Keymaster

Dear John,

basically you set the pressure! You may leave it to the values 1 and chance the charRho in..

std::vector<T> rho0;
rho0.push_back( 1 ); //first lattice
rho0.push_back( 1 ); //second lattice -> e.g. change it to 2
PsiEqualsRho<T,T> interactionPotential;
ShanChenForcedGenerator3D<T,DESCRIPTOR> coupling( G,rho0,interactionPotential );

Best
Mathias