High-order mixed finite elements for an energy-based model of the polarization process in ferroelectric materials. (February 2021)