Implement also finite-element solver for AE?

check [featom](https://github.com/atomic-solvers/featom). The fe solver can be very stable and numerically robust according to the featom paper, which is useful for meta-gga, or even range-seperated functionals.