Can I use a stabilized method such as SUPG or Galerkin Least Squares in Sundance?
Yes; stabilization terms can be represented by the expression language just like any other term. There are two ways in which stabilization terms differ from the usual terms in weak forms: they include mesh-specific information (usually a cell diameter), and they drop surface terms that arise during integrations by parts. The cell diameter can be represented with a CellDiameterExpr object, and the surface terms are left out implicitly. The tutorial contains an example of pressure stabilization of the Navier-Stokes equations.