public unsafe override void Initialize(ContentArchive content, Camera camera)
{
camera.Position = new Vector3(-13f, 6, -13f);
camera.Yaw = MathHelper.Pi * 3f / 4;
camera.Pitch = MathHelper.Pi * 0.05f;
Simulation = Simulation.Create(BufferPool, new DemoNarrowPhaseCallbacks(), new DemoPoseIntegratorCallbacks(new Vector3(0, -10, 0)));
using (var compoundBuilder = new CompoundBuilder(BufferPool, Simulation.Shapes, 8))
{
{
//Note that, in bepuphysics v2, there is no 'recentering' when constructing a shape. The pose you pass in for a child is exactly the pose that the compound will use,
//even if the 'true' center of mass isn't at the local origin.
//Instead, if recentering is desired, it should performed ahead of time. The CompoundBuilder can help with this.
//We'll construct this compound using shapes far from the origin, and then use the CompoundBuilder overload that recenters the children and outputs the computed center.
var capsuleChildShape = new Capsule(0.5f, 0.5f);
var capsuleLocalPose = new RigidPose {
Position = new Vector3(-0.5f, 4, 4), Orientation = Quaternion.Identity
};
var boxChildShape = new Box(0.5f, 1f, 1.5f);
var boxLocalPose = new RigidPose {
Position = new Vector3(0.5f, 4, 4), Orientation = Quaternion.Identity
};
//All allocations from the buffer pool used for the final compound shape will be disposed when the demo is disposed. Don't have to worry about leaks in these demos.
compoundBuilder.Add(capsuleChildShape, capsuleLocalPose, 1);
compoundBuilder.Add(boxChildShape, boxLocalPose, 1);
compoundBuilder.BuildDynamicCompound(out var compoundChildren, out var compoundInertia, out var compoundCenter);
compoundBuilder.Reset();
var compound = new Compound(compoundChildren);
var compoundDescription = new BodyDescription
{
Activity = new BodyActivityDescription {
SleepThreshold = 0.01f, MinimumTimestepCountUnderThreshold = 32
},
Collidable = new CollidableDescription
{
Shape = Simulation.Shapes.Add(compound),
SpeculativeMargin = 0.1f,
},
LocalInertia = compoundInertia,
Pose = new RigidPose {
Position = compoundCenter, Orientation = Quaternion.Identity
},
};
Simulation.Bodies.Add(compoundDescription);
}
//Build a stack of sphere grids to stress manifold reduction heuristics in a convex-ish situation.
{
var gridShape = new Sphere(0.5f);
const float gridSpacing = 1.5f;
const int gridWidth = 3;
var gridShapeIndex = Simulation.Shapes.Add(gridShape);
gridShape.ComputeInertia(1, out var gridBoxInertia);
float localPoseOffset = -0.5f * gridSpacing * (gridWidth - 1);
for (int i = 0; i < gridWidth; ++i)
{
for (int j = 0; j < gridWidth; ++j)
{
var localPose = new RigidPose
{
Orientation = Quaternion.Identity,
Position = new Vector3(localPoseOffset, 0, localPoseOffset) + new Vector3(gridSpacing) * new Vector3(i, 0, j)
};
compoundBuilder.Add(gridShapeIndex, localPose, gridBoxInertia.InverseInertiaTensor, 1);
}
}
compoundBuilder.BuildDynamicCompound(out var gridChildren, out var gridInertia, out var center);
compoundBuilder.Reset();
var gridCompound = new Compound(gridChildren);
var bodyDescription = new BodyDescription
{
Activity = new BodyActivityDescription {
SleepThreshold = 0.01f, MinimumTimestepCountUnderThreshold = 32
},
Collidable = new CollidableDescription
{
Shape = Simulation.Shapes.Add(gridCompound),
SpeculativeMargin = 0.1f,
},
LocalInertia = gridInertia,
Pose = new RigidPose {
Orientation = Quaternion.Identity
}
};
for (int i = 0; i < 4; ++i)
{
bodyDescription.Pose.Position = new Vector3(0, 2 + i * 3, 0);
//if (i == 0)
// gridDescription.LocalInertia = new BodyInertia();
//else
// gridDescription.LocalInertia = gridInertia;
Simulation.Bodies.Add(bodyDescription);
}
}
//Build a table and use it for a couple of different tests.
{
var legShape = new Box(0.2f, 1, 0.2f);
legShape.ComputeInertia(1f, out var legInverseInertia);
var legShapeIndex = Simulation.Shapes.Add(legShape);
var legPose0 = new RigidPose {
Position = new Vector3(-1.5f, 0, -1.5f), Orientation = Quaternion.Identity
};
var legPose1 = new RigidPose {
Position = new Vector3(-1.5f, 0, 1.5f), Orientation = Quaternion.Identity
};
var legPose2 = new RigidPose {
Position = new Vector3(1.5f, 0, -1.5f), Orientation = Quaternion.Identity
};
var legPose3 = new RigidPose {
Position = new Vector3(1.5f, 0, 1.5f), Orientation = Quaternion.Identity
};
compoundBuilder.Add(legShapeIndex, legPose0, legInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.Add(legShapeIndex, legPose1, legInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.Add(legShapeIndex, legPose2, legInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.Add(legShapeIndex, legPose3, legInverseInertia.InverseInertiaTensor, 1);
var tableTopPose = new RigidPose {
Position = new Vector3(0, 0.6f, 0), Orientation = Quaternion.Identity
};
var tableTopShape = new Box(3.2f, 0.2f, 3.2f);
compoundBuilder.Add(tableTopShape, tableTopPose, 3);
compoundBuilder.BuildDynamicCompound(out var tableChildren, out var tableInertia, out var tableCenter);
compoundBuilder.Reset();
var table = new Compound(tableChildren);
var tableDescription = new BodyDescription
{
Activity = new BodyActivityDescription {
SleepThreshold = 0.01f, MinimumTimestepCountUnderThreshold = 32
},
Collidable = new CollidableDescription
{
Shape = Simulation.Shapes.Add(table),
SpeculativeMargin = 0.1f,
},
LocalInertia = tableInertia,
Pose = new RigidPose {
Orientation = Quaternion.Identity
}
};
//Stack some tables.
{
for (int i = 0; i < 10; ++i)
{
tableDescription.Pose.Position = new Vector3(10, 3 + i * 1.4f, 10);
Simulation.Bodies.Add(tableDescription);
}
}
{
for (int k = 0; k < 5; ++k)
{
tableDescription.Pose.Position = new Vector3(64 + k * 3, 6 + k * 1.4f, 32);
Simulation.Bodies.Add(tableDescription);
}
//for (int i = 0; i < 10; ++i)
//{
// for (int j = 0; j < 20; ++j)
// {
// for (int k = 0; k < 10; ++k)
// {
// tableDescription.Pose.Position = new Vector3(32 + i * 6, 6 + j * 1.4f, 16 + k * 6);
// Simulation.Bodies.Add(tableDescription);
// }
// }
//}
}
//Put a table on top of a sphere to stress out nonconvex reduction for divergent normals.
{
tableDescription.Pose.Position = new Vector3(10, 6, 0);
Simulation.Bodies.Add(tableDescription);
var sphereShape = new Sphere(3);
var sphereIndex = Simulation.Shapes.Add(sphereShape);
var sphereDescription = new StaticDescription
{
Collidable = new CollidableDescription
{
Shape = sphereIndex,
SpeculativeMargin = 0.1f,
},
Pose = new RigidPose {
Position = new Vector3(10, 2, 0), Orientation = Quaternion.Identity
}
};
Simulation.Statics.Add(sphereDescription);
}
//Put another table on the ground, but with a clamp-ish thing on it that generates opposing normals.
{
tableDescription.Pose.Position = new Vector3(10, 3, -10);
Simulation.Bodies.Add(tableDescription);
var clampPieceShape = new Box(2f, 0.1f, 0.3f);
clampPieceShape.ComputeInertia(1f, out var clampPieceInverseInertia);
var clampPieceShapeIndex = Simulation.Shapes.Add(clampPieceShape);
var clamp0 = new RigidPose {
Position = new Vector3(0, -0.2f, -1.1f), Orientation = Quaternion.Identity
};
var clamp1 = new RigidPose {
Position = new Vector3(0, 0.2f, -1.1f), Orientation = Quaternion.Identity
};
var clamp2 = new RigidPose {
Position = new Vector3(0, -0.2f, 0), Orientation = Quaternion.Identity
};
var clamp3 = new RigidPose {
Position = new Vector3(0, 0.2f, 0), Orientation = Quaternion.Identity
};
var clamp4 = new RigidPose {
Position = new Vector3(0, -0.2f, 1.1f), Orientation = Quaternion.Identity
};
var clamp5 = new RigidPose {
Position = new Vector3(0, 0.2f, 1.1f), Orientation = Quaternion.Identity
};
compoundBuilder.Add(clampPieceShapeIndex, clamp0, clampPieceInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.Add(clampPieceShapeIndex, clamp1, clampPieceInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.Add(clampPieceShapeIndex, clamp2, clampPieceInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.Add(clampPieceShapeIndex, clamp3, clampPieceInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.Add(clampPieceShapeIndex, clamp4, clampPieceInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.Add(clampPieceShapeIndex, clamp5, clampPieceInverseInertia.InverseInertiaTensor, 1);
compoundBuilder.BuildDynamicCompound(out var clampChildren, out var clampInertia, out var clampCenter);
compoundBuilder.Reset();
var clamp = new Compound(clampChildren);
var clampDescription = new BodyDescription
{
Activity = new BodyActivityDescription {
SleepThreshold = 0.01f, MinimumTimestepCountUnderThreshold = 32
},
Collidable = new CollidableDescription
{
Shape = Simulation.Shapes.Add(clamp),
SpeculativeMargin = 0.1f,
},
LocalInertia = clampInertia,
Pose = new RigidPose {
Position = tableDescription.Pose.Position + new Vector3(2f, 0.3f, 0), Orientation = Quaternion.Identity
}
};
Simulation.Bodies.Add(clampDescription);
}
}
//Create a tree-accelerated big compound.
{
var random = new Random(5);
var treeCompoundBoxShape = new Box(0.5f, 1.5f, 1f);
var treeCompoundBoxShapeIndex = Simulation.Shapes.Add(treeCompoundBoxShape);
treeCompoundBoxShape.ComputeInertia(1, out var childInertia);
for (int i = 0; i < 128; ++i)
{
RigidPose localPose;
localPose.Position = new Vector3(12, 6, 12) * (0.5f * new Vector3((float)random.NextDouble(), (float)random.NextDouble(), (float)random.NextDouble()) - Vector3.One);
float orientationLengthSquared;
do
{
localPose.Orientation = new Quaternion((float)random.NextDouble(), (float)random.NextDouble(), (float)random.NextDouble(), (float)random.NextDouble());
}while ((orientationLengthSquared = localPose.Orientation.LengthSquared()) < 1e-9f);
Quaternion.Scale(localPose.Orientation, 1f / MathF.Sqrt(orientationLengthSquared), out localPose.Orientation);
//Quaternion.CreateFromAxisAngle(new Vector3(1, 0, 0), MathF.PI, out localPose.Orientation);
compoundBuilder.Add(treeCompoundBoxShapeIndex, localPose, childInertia.InverseInertiaTensor, 1);
}
compoundBuilder.BuildDynamicCompound(out var children, out var inertia, out var center);
compoundBuilder.Reset();
var compound = new BigCompound(children, Simulation.Shapes, BufferPool);
//var compound = new Compound(children);
var compoundIndex = Simulation.Shapes.Add(compound);
for (int i = 0; i < 8; ++i)
{
Simulation.Bodies.Add(BodyDescription.CreateDynamic(new Vector3(0, 4 + 5 * i, 32), inertia, new CollidableDescription(compoundIndex, 0.1f), new BodyActivityDescription(0.01f)));
}
}
}
//Prevent stuff from falling into the infinite void.
{
var boxShape = new Box(256, 1, 256);
var groundShapeIndex = Simulation.Shapes.Add(boxShape);
var groundDescription = new StaticDescription
{
Collidable = new CollidableDescription
{
Shape = groundShapeIndex,
SpeculativeMargin = 0.1f,
},
Pose = new RigidPose {
Position = new Vector3(0, 0, 0), Orientation = Quaternion.Identity
}
};
Simulation.Statics.Add(groundDescription);
}
const int planeWidth = 48;
const int planeHeight = 48;
MeshDemo.CreateDeformedPlane(planeWidth, planeHeight,
(int x, int y) =>
{
Vector2 offsetFromCenter = new Vector2(x - planeWidth / 2, y - planeHeight / 2);
return(new Vector3(offsetFromCenter.X, MathF.Cos(x / 4f) * MathF.Sin(y / 4f) - 0.01f * offsetFromCenter.LengthSquared(), offsetFromCenter.Y));
}, new Vector3(2, 1, 2), BufferPool, out var planeMesh);
Simulation.Statics.Add(new StaticDescription(new Vector3(64, 4, 32), Quaternion.CreateFromAxisAngle(new Vector3(0, 1, 0), MathF.PI / 2),
new CollidableDescription(Simulation.Shapes.Add(planeMesh), 0.1f)));
}
public unsafe override void Initialize(ContentArchive content, Camera camera)
{
camera.Position = new Vector3(0, 10, 40);
camera.Yaw = 0;
camera.Pitch = 0;
Simulation = Simulation.Create(BufferPool, new TestCallbacks());
Simulation.PoseIntegrator.Gravity = new Vector3(0, -10, 0);
var box = new Box(2f, 2f, 2f);
var capsule = new Capsule(1f, 1f);
var sphere = new Sphere(1.5f);
box.ComputeInertia(1, out var boxInertia);
capsule.ComputeInertia(1, out var capsuleInertia);
sphere.ComputeInertia(1, out var sphereInertia);
var boxIndex = Simulation.Shapes.Add(box);
var capsuleIndex = Simulation.Shapes.Add(capsule);
var sphereIndex = Simulation.Shapes.Add(sphere);
const int width = 12;
const int height = 3;
const int length = 12;
for (int i = 0; i < width; ++i)
{
for (int j = 0; j < height; ++j)
{
for (int k = 0; k < length; ++k)
{
var location = new Vector3(5, 5, 5) * new Vector3(i, j, k) + new Vector3(-width * 2.5f, 2.5f, -length * 2.5f);
var bodyDescription = new BodyDescription
{
Activity = new BodyActivityDescription {
MinimumTimestepCountUnderThreshold = 32, SleepThreshold = 0.1f
},
Pose = new RigidPose
{
Orientation = Quaternion.Identity,
Position = location
},
Collidable = new CollidableDescription
{
Continuity = new ContinuousDetectionSettings {
Mode = ContinuousDetectionMode.Discrete
},
SpeculativeMargin = 0.1f
}
};
switch (j % 3)
{
case 0:
bodyDescription.Collidable.Shape = boxIndex;
bodyDescription.LocalInertia = boxInertia;
break;
case 1:
bodyDescription.Collidable.Shape = capsuleIndex;
bodyDescription.LocalInertia = capsuleInertia;
break;
case 2:
bodyDescription.Collidable.Shape = sphereIndex;
bodyDescription.LocalInertia = sphereInertia;
break;
}
Simulation.Bodies.Add(bodyDescription);
}
}
}
//var staticShapeIndex = Simulation.Shapes.Add(new Box(100, 1, 100));
//var staticDescription = new StaticDescription
//{
// Collidable = new CollidableDescription
// {
// Continuity = new ContinuousDetectionSettings { Mode = ContinuousDetectionMode.Discrete },
// Shape = Simulation.Shapes.Add(new Box(100, 1, 100)),
// SpeculativeMargin = 0.1f
// },
// Pose = new RigidPose { Position = new Vector3(0, -1, 0), Orientation = Quaternion.Identity }
//};
//Simulation.Statics.Add(staticDescription);
const int planeWidth = 64;
const int planeHeight = 64;
MeshDemo.CreateDeformedPlane(planeWidth, planeHeight,
(int x, int y) =>
{
return(new Vector3(x, 1 * MathF.Cos(x / 4f) * MathF.Sin(y / 4f), y));
}, new Vector3(2, 3, 2), BufferPool, out var planeMesh);
Simulation.Statics.Add(new StaticDescription(new Vector3(-64, -10, -64), new CollidableDescription(Simulation.Shapes.Add(planeMesh), 0.1f)));
}
public unsafe override void Initialize(ContentArchive content, Camera camera)
{
camera.Position = new Vector3(-15f, 20, -15f);
camera.Yaw = MathHelper.Pi * 3f / 4;
camera.Pitch = MathHelper.Pi * 0.1f;
//Using minimum sized allocations forces as many resizes as possible.
Simulation = Simulation.Create(BufferPool, new DemoNarrowPhaseCallbacks(), new DemoPoseIntegratorCallbacks(new Vector3(0, -10, 0)), initialAllocationSizes:
new SimulationAllocationSizes
{
Bodies = 1,
ConstraintCountPerBodyEstimate = 1,
Constraints = 1,
ConstraintsPerTypeBatch = 1,
Islands = 1,
ShapesPerType = 1,
Statics = 1
});
Simulation.Deterministic = false;
const int planeWidth = 8;
const int planeHeight = 8;
MeshDemo.CreateDeformedPlane(planeWidth, planeHeight,
(int x, int y) =>
{
Vector2 offsetFromCenter = new Vector2(x - planeWidth / 2, y - planeHeight / 2);
return(new Vector3(offsetFromCenter.X, MathF.Cos(x / 4f) * MathF.Sin(y / 4f) - 0.2f * offsetFromCenter.LengthSquared(), offsetFromCenter.Y));
}, new Vector3(2, 1, 2), BufferPool, out var staticShape);
var staticShapeIndex = Simulation.Shapes.Add(staticShape);
const int staticGridWidthInSpheres = 128;
const float staticSpacing = 8;
for (int i = 0; i < staticGridWidthInSpheres; ++i)
{
for (int j = 0; j < staticGridWidthInSpheres; ++j)
{
var staticDescription = new StaticDescription
{
Collidable = new CollidableDescription
{
Continuity = new ContinuousDetectionSettings {
Mode = ContinuousDetectionMode.Discrete
},
Shape = staticShapeIndex,
SpeculativeMargin = 0.1f
},
Pose = new RigidPose
{
Position = new Vector3(
-staticGridWidthInSpheres * staticSpacing * 0.5f + i * staticSpacing,
-4 + 4 * (float)Math.Cos(i * 0.3) + 4 * (float)Math.Cos(j * 0.3),
-staticGridWidthInSpheres * staticSpacing * 0.5f + j * staticSpacing),
Orientation = BepuUtilities.Quaternion.Identity
}
};
Simulation.Statics.Add(staticDescription);
}
}
//A bunch of kinematic balls do acrobatics as an extra stressor.
var kinematicShape = new Sphere(8);
var kinematicShapeIndex = Simulation.Shapes.Add(kinematicShape);
var kinematicCount = 64;
var anglePerKinematic = MathHelper.TwoPi / kinematicCount;
var startingRadius = 256;
kinematicHandles = new int[kinematicCount];
for (int i = 0; i < kinematicCount; ++i)
{
var angle = anglePerKinematic * i;
var description = new BodyDescription
{
Collidable = new CollidableDescription
{
Continuity = new ContinuousDetectionSettings {
Mode = ContinuousDetectionMode.Discrete
},
Shape = kinematicShapeIndex,
SpeculativeMargin = 0.1f
},
Pose = new RigidPose
{
Position = new Vector3(
startingRadius * (float)Math.Cos(angle),
0,
startingRadius * (float)Math.Sin(angle)),
Orientation = BepuUtilities.Quaternion.Identity
},
Activity = new BodyActivityDescription {
SleepThreshold = 0, MinimumTimestepCountUnderThreshold = 4
},
};
kinematicHandles[i] = Simulation.Bodies.Add(description);
}
dynamicHandles = new QuickQueue <int>(65536, BufferPool);
removedStatics = new QuickQueue <StaticDescription>(512, BufferPool);
random = new Random(5);
}