“…The results require some regularity assumption and that the mesh size of the discretization is small enough. On the other hand, once the coerciveness of the least-squares method is proved, see [18,6,9,11,45,59], the LSFEM with H 1 -H(div) conforming finite element approximations is automatic coercive, without the restriction on the mesh size.…”