public void AddPressureForces(float timeStep)
{
if (pressure == 0.0f || volume == 0.0f) return;
float invVolume = 1.0f / volume;
foreach (Triangle t in triangles)
{
JVector v1 = points[t.indices.I0].position;
JVector v2 = points[t.indices.I1].position;
JVector v3 = points[t.indices.I2].position;
JVector cross = (v3 - v1) % (v2 - v1);
JVector center = (v1 + v2 + v3) * (1.0f / 3.0f);
points[t.indices.I0].AddForce(invVolume * cross * pressure);
points[t.indices.I1].AddForce(invVolume * cross * pressure);
points[t.indices.I2].AddForce(invVolume * cross * pressure);
}
}
#endregion