Numerical modeling of anisotropic ferroelectric materials with hybridizable discontinuous Galerkin methods