Conventional crystal plasticity (CP) solvers are based on a Newton-Raphson (NR) approach which use an initial guess for the free variables (often stress) to be solved. These solvers are limited by a finite interval of convergence and often fail when the free variable falls outside this interval. Solution failure results in the reduction of the t…