Currently, the initialization of the crack field used for anisotropic models only initializes the field on the crack boundaries. This could be enhanced.
One must be cautious not to set the lower bounds with the predicted initial field when initializing. Indeed, there is an issue if the initial field is overestimated and can not be reduced during the first solver step.