Find here a solver for generic systems of PDEs, deterministic or stochastic, written in C++ and CUDA.
See tutorials on how to write solvers here, and examples here.