text: from firedrake import *
  
  # Stabilizing parameters
  delta_0 = Constant(1)
  delta_1 = Constant(-0.5)
  delta_2 = Constant(0.5)
  delta_3 = Constant(0.5)
  eta_u = Constant(50.0)
  eta_p = 10 * eta_u
  beta_0 = Constant(1.0e-15)
  mesh_parameter = True
  
  solvers_args = {
      "cgls_full": {
          "delta_0": Constant(1),
          "delta_1": Constant(-0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
      },
      "cgls_div": {
          "delta_0": Constant(1),
          "delta_1": Constant(-0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
      },
      "mgls_full": {
          "delta_0": Constant(1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
      },
      "mgls": {
          "delta_0": Constant(1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
      },
      "mvh_full": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
      },
      "mvh_div": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
      },
      "mvh": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.0),
          "delta_3": Constant(0.0),
      },
      ###############################################
      "dgls_full": {
          "delta_0": Constant(1),
          "delta_1": Constant(-0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
          "eta_u": eta_u,
          "eta_p": eta_p,
      },
      "dgls_div": {
          "delta_0": Constant(1),
          "delta_1": Constant(-0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
          "eta_u": eta_u,
          "eta_p": eta_p,
      },
      "dmgls_full": {
          "delta_0": Constant(1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
          "eta_u": eta_u,
          "eta_p": eta_p,
      },
      "dmgls": {
          "delta_0": Constant(1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
          "eta_u": eta_u,
          "eta_p": eta_p,
      },
      "dmvh_full": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
          "eta_u": eta_u,
          "eta_p": eta_p,
      },
      "dmvh_div": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
          "eta_u": eta_u,
          "eta_p": eta_p,
      },
      "dmvh": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.0),
          "delta_3": Constant(0.0),
          "eta_u": eta_u,
          "eta_p": eta_p,
      },
      ###############################################
      "sdhm_full": {
          "delta_0": Constant(1),
          "delta_1": Constant(-0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
          "beta_0": beta_0,
      },
      "sdhm_div": {
          "delta_0": Constant(1),
          "delta_1": Constant(-0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
          "beta_0": beta_0,
      },
      "hmgls_full": {
          "delta_0": Constant(1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
          "beta_0": beta_0,
      },
      "hmgls": {
          "delta_0": Constant(1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
          "beta_0": beta_0,
      },
      "hmvh_full": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.5),
          "beta_0": beta_0,
      },
      "hmvh_div": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.5),
          "delta_3": Constant(0.0),
          "beta_0": beta_0,
      },
      "hmvh": {
          "delta_0": Constant(-1),
          "delta_1": Constant(0.5),
          "delta_2": Constant(0.0),
          "delta_3": Constant(0.0),
          "beta_0": beta_0,
      },
  }
  