Write VTU output with higher order elements

diff --git a/examples/step-44/step-44.cc b/examples/step-44/step-44.cc
index 2cb1b0c4773db08ee5d826ea5c277300bc06e039..f57067007ea35fc5f9ecbcc0776bd484d4a14aa5 100644 (file)
--- a/examples/step-44/step-44.cc
+++ b/examples/step-44/step-44.cc
@@ -3373,6 +3373,10 @@ namespace Step44
     solution_name.emplace_back("pressure");
     solution_name.emplace_back("dilatation");
 
+    DataOutBase::VtkFlags output_flags;
+    output_flags.write_higher_order_cells = true;
+    data_out.set_flags(output_flags);
+
     data_out.attach_dof_handler(dof_handler);
     data_out.add_data_vector(solution_n,
                              solution_name,
@@ -3394,8 +3398,8 @@ namespace Step44
     data_out.build_patches(q_mapping, degree);
 
     std::ofstream output("solution-" + std::to_string(dim) + "d-" +
-                         std::to_string(time.get_timestep()) + ".vtk");
-    data_out.write_vtk(output);
+                         std::to_string(time.get_timestep()) + ".vtu");
+    data_out.write_vtu(output);
   }
 
 } // namespace Step44