I am interested in study critical phenomena for the Eisntein-Vlasov system. I am working in maximal slicing condition with areal coordinates. I use two methods to attack this problem: particle-mesh (PM) method and solve the Boltzmman equation directly (using RNPL) on the tangent bundle of the spacetime. In the next PS file you can have a look at the equations in both cases.

Particle-Mesh Code:

Solving Boltzmman Equation: