Recursive multilevel trust region method with application to fully monolithic phase-field models of brittle fracture.