Update use of bounds constraints

In newer versions of Irksome, bounds constraints are not passed to the
.advance method but rather in the initialization of the time stepper.

Author: danshapero

notebooks/rainier-monolithic.ipynb

@@ -610,7 +610,7 @@
     "    },\n",
     "    \"stage_type\": \"value\",\n",
     "    \"basis_type\": \"Bernstein\",\n",
-    "    #\"bounds\": bounds,\n",
+    "    \"bounds\": bounds,\n",
     "}\n",
     "\n",
     "solver = irksome.TimeStepper(F, tableau, t, dt, w, **bparams)"
@@ -637,7 +637,7 @@
     "final_time = 500.0\n",
     "num_steps = int(final_time / float(dt))\n",
     "for step in trange(num_steps):\n",
-    "    solver.advance(bounds=bounds)\n",
+    "    solver.advance()\n",
     "    h = w.subfunctions[3]\n",
     "    a.interpolate(smb(b + h))\n",
     "\n",