Adapt demos

Author: schnellerhase

python/demo/demo_axis.py

@@ -794,6 +794,6 @@ def create_eps_mu(pml, rho, eps_bkg, mu_bkg):
     Es_dg = fem.Function(W)
     Es_expr = fem.Expression(Esh, W.element.interpolation_points)
     Es_dg.interpolate(Es_expr)
-    with VTXWriter(mesh_data.mesh.comm, out_folder / "Es.bp", Es_dg) as f:
+    with VTXWriter(mesh_data.mesh.comm, out_folder / "Es.bp", "w", Es_dg) as f:
         f.write(0.0)
 # -

python/demo/demo_biharmonic.py

@@ -330,7 +330,7 @@
 if has_adios2:
     out_folder = Path("out_biharmonic")
     out_folder.mkdir(parents=True, exist_ok=True)
-    with io.VTXWriter(msh.comm, out_folder / "biharmonic.bp", [uh]) as file:
+    with io.VTXWriter(msh.comm, out_folder / "biharmonic.bp", "w", [uh]) as file:
         file.write(0.0)
 
 # and displayed using [pyvista](https://docs.pyvista.org/).

python/demo/demo_half_loaded_waveguide.py

@@ -469,10 +469,10 @@ def verify_mode(
 
         if has_vtx:
             # Save solutions
-            with VTXWriter(msh.comm, out_folder / f"Et_{i}.bp", Et_dg) as f:
+            with VTXWriter(msh.comm, out_folder / f"Et_{i}.bp", "w", Et_dg) as f:
                 f.write(0.0)
 
-            with VTXWriter(msh.comm, out_folder / f"Ez_{i}.bp", ezh) as f:
+            with VTXWriter(msh.comm, out_folder / f"Ez_{i}.bp", "w", ezh) as f:
                 f.write(0.0)
 
         # Visualize solutions with Pyvista

python/demo/demo_hdg.py

@@ -264,9 +264,9 @@ def u_e(x):
 if dolfinx.has_adios2:
     from dolfinx.io import VTXWriter
 
-    with VTXWriter(msh.comm, "u.bp", u, "bp4") as f:
+    with VTXWriter(msh.comm, "u.bp", "w", u, "bp4") as f:
         f.write(0.0)
-    with VTXWriter(msh.comm, "ubar.bp", ubar, "bp4") as f:
+    with VTXWriter(msh.comm, "ubar.bp", "w", ubar, "bp4") as f:
         f.write(0.0)
 else:
     print("ADIOS2 required for VTX output")

python/demo/demo_interpolation-io.py

@@ -82,7 +82,7 @@
 if has_adios2:
     from dolfinx.io import VTXWriter
 
-    with VTXWriter(msh.comm, out_folder / "output_nedelec.bp", u0) as f:
+    with VTXWriter(msh.comm, out_folder / "output_nedelec.bp", "w", u0) as f:
         f.write(0.0)
 else:
     print("ADIOS2 required for VTX output")

python/demo/demo_mixed-poisson.py

@@ -352,7 +352,7 @@
 if dolfinx.has_adios2:
     from dolfinx.io import VTXWriter
 
-    with VTXWriter(msh.comm, "output_mixed_poisson.bp", u) as f:
+    with VTXWriter(msh.comm, "output_mixed_poisson.bp", "w", u) as f:
         f.write(0.0)
 else:
     print("ADIOS2 required for VTX output.")

python/demo/demo_navier-stokes.py

@@ -367,8 +367,8 @@ def jump(phi, n):
 if has_adios2:
     u_vis = fem.Function(W, name="u_init")
     u_vis.interpolate(u_h)
-    u_file = io.VTXWriter(msh.comm, "u.bp", u_vis)
-    p_file = io.VTXWriter(msh.comm, "p.bp", p_h)
+    u_file = io.VTXWriter(msh.comm, "u.bp", "w", u_vis)
+    p_file = io.VTXWriter(msh.comm, "p.bp", "w", p_h)
     u_file.write(t)
     p_file.write(t)
 else:

python/demo/demo_pml.py

@@ -666,7 +666,7 @@ def create_eps_mu(
 Esh_dg = fem.Function(V_dg)
 Esh_dg.interpolate(Esh)
 
-with VTXWriter(mesh_data.mesh.comm, out_folder / "Esh.bp", Esh_dg) as vtx:
+with VTXWriter(mesh_data.mesh.comm, out_folder / "Esh.bp", "w", Esh_dg) as vtx:
     vtx.write(0.0)
 # -
 
@@ -702,7 +702,7 @@ def create_eps_mu(
 E_dg = fem.Function(V_dg)
 E_dg.interpolate(E)
 
-with VTXWriter(mesh_data.mesh.comm, out_folder / "E.bp", E_dg) as vtx:
+with VTXWriter(mesh_data.mesh.comm, out_folder / "E.bp", "w", E_dg) as vtx:
     vtx.write(0.0)
 # -
 

python/demo/demo_scattering_boundary_conditions.py

@@ -655,7 +655,7 @@ def curl_2d(f: fem.Function):
 assert isinstance(Esh, fem.Function)
 Esh_dg.interpolate(Esh)
 
-with io.VTXWriter(mesh_data.mesh.comm, out_folder / "Esh.bp", Esh_dg) as vtx:
+with io.VTXWriter(mesh_data.mesh.comm, out_folder / "Esh.bp", "w", Esh_dg) as vtx:
     vtx.write(0.0)
 # -
 
@@ -694,7 +694,7 @@ def curl_2d(f: fem.Function):
 E.x.array[:] = Eb.x.array[:] + Esh.x.array[:]
 E_dg = fem.Function(V_dg)
 E_dg.interpolate(E)
-with io.VTXWriter(mesh_data.mesh.comm, "E.bp", E_dg) as vtx:
+with io.VTXWriter(mesh_data.mesh.comm, "E.bp", "w", E_dg) as vtx:
     vtx.write(0.0)
 # -