NRPy+ ("Python-based code generation for numerical relativity... and beyond!") is a Python package aimed at accelerating the process of encoding partial differential equations (PDEs) into high-performance software. Users may input complex expressions in clean tensorial form (Einstein notation), which acts to minimize human error. NRPy+ processes this input and outputs highly optimized C code with CSE and/or SIMD compiler intrinsics if desired.
Link to Source Code