Contact force

prims.get_contact_force_matrix() can only provide the net force between two rigid bodies, but how can I get a more dense contact force distribution?