Dear deal.ii developers and community,
Recently, our group developed a limited-memory BFGS solver (a type of quasi Newton method) as a monolithic approach to solve the phase-field crack propagation problems via deal.ii. My question is, which is a better way to share the code with the deal.ii community? Should we contribute the code in the deal.ii code-gallery? Or could the code be potentially included as a tutorial step?
Here is more information:
And our source code and input files for all the numerical examples are hosted here:
The unique features of our code comparing with existing tutorial steps:
1. The implementation of the limited-memory BFGS solver (a type of quasi Newton solver) that can be used in various non-convex optimization problems.
2. The limited-memory BFGS solver is combined with an adaptive refinement technique, which involves solution transfer and projection from coarse meshes to refined meshes.
3. A widely adopted phase-field model with various 2D and 3D benchmark problems.
Thanks for your feedback.
Tao Jin, PhD
University of Ottawa