A massively parallel program to solve the phase field formulation for crack propagation