Working on implementing ADM and BSSN 3+1 formalism. I coded this up after reading through a few books and papers on Numerical Relativity. These books and papers all tend to jump right into the math required to write the variables in terms of their partial time derivatives, but don't seem to get to the hyperbolicity requirements for numerically modeling the equations until later. Key note.

I implemented the equations in FTCS and wrote this simulation. I then learned the importance of hyperbolic formalisms, and am currently educating myself in them. To keep tabs on that, check out my hydrodynamics projects: in Javascript, in C++, and in OpenCL.

Browse the FTCS BSSN source code here.