Proposed changes

This variant of the elliptic DG formulation moves the Jacobian into the divergence like this:

$$ \partial_i F^i = \frac{1}{J} \partial_\hat{\imath} J {J^{-1}}^{\hat{\imath}}_j F^j $$

This formulation makes the strong scheme work as well as the weak scheme for extremely stretched grids (initial data with large outer boundaries).

Upgrade instructions

Code review checklist

• The code is documented and the documentation renders correctly. Run
  make doc to generate the documentation locally into BUILD_DIR/docs/html.
  Then open index.html.
• The code follows the stylistic and code quality guidelines listed in the
  code review guide.
• The PR lists upgrade instructions and is labeled bugfix or
  new feature if appropriate.

Further comments