Formation of grain boundaries in ductile single crystals under plane-strain simple shear: MATLAB implementation

MATLAB implementation and the scripts that reproduce every figure and table of the paper "Formation of grain boundaries in ductile single crystals under plane-strain simple shear: a block-coordinate finite element method" by K. C. Le and T. D. Nguyen.

The code minimises the continuum-dislocation energy over the deformation and the plastic slip by a block-coordinate scheme: an ADMM solve for the non-smooth slip block and a Levenberg-regularised Newton solve for the deformation block, with a polyconvex Ciarlet-Geymonat elastic energy.

Running reproduce_all.m regenerates every figure and prints every table in a few seconds from the saved results. Twelve further scripts verify the propositions and the algebra of the paper, each printing PASS or FAIL.
Requires MATLAB R2025b; no toolboxes.