[Getfem-commits] (no subject)

[Tetsuo Koyama]

branch: devel-tetsuo-xml
commit f2fedea1eafb8b5dca5dc12bbde7bc5216024560
Author: Tetsuo Koyama <address@hidden>
AuthorDate: Sat May 2 06:05:17 2020 +0000

    :test_tube: export_to_vtu mesh_fem
---
 tests/thermo_elasticity_electrical_coupling.cc | 36 +++++++++++++++++++-------
 1 file changed, 26 insertions(+), 10 deletions(-)

diff --git a/tests/thermo_elasticity_electrical_coupling.cc 
b/tests/thermo_elasticity_electrical_coupling.cc
index 18a8ebe..d545937 100644
--- a/tests/thermo_elasticity_electrical_coupling.cc
+++ b/tests/thermo_elasticity_electrical_coupling.cc
@@ -309,21 +309,37 @@ int main(int argc, char *argv[]) {
   plain_vector CO(mfvm.nb_dof() * 2);
   getfem::ga_interpolation_Lagrange_fem(md, "-"+sigmaeps+"*Grad_V",  mfvm, CO);
   
-  getfem::vtk_export exp("displacement_with_von_mises.vtk", false);
-  exp.exporting(mfu);
-  exp.write_point_data(mfu, U, "elastostatic displacement");
-  exp.write_point_data(mfvm, VM, "Von Mises stress");
+  getfem::vtk_export vtk_exp("displacement_with_von_mises.vtk", false);
+  vtk_exp.exporting(mfu);
+  vtk_exp.write_point_data(mfu, U, "elastostatic displacement");
+  vtk_exp.write_point_data(mfvm, VM, "Von Mises stress");
+
+  getfem::vtu_export vtu_exp("displacement_with_von_mises.vtu", true);
+  vtu_exp.exporting(mfu);
+  vtu_exp.write_point_data(mfu, U, "elastostatic displacement");
+  vtu_exp.write_point_data(mfvm, VM, "Von Mises stress");
+
   cout << "\nYou can view solutions with for instance:\n\nmayavi2 "
     "-d displacement_with_von_mises.vtk -f WarpVector -m Surface\n" << endl;
   
-  getfem::vtk_export exp2("temperature.vtk", false);
-  exp2.exporting(mft);
-  exp2.write_point_data(mft, THETA, "Temperature");
+  getfem::vtk_export vtk_exp2("temperature.vtk", false);
+  vtk_exp2.exporting(mft);
+  vtk_exp2.write_point_data(mft, THETA, "Temperature");
+
+  getfem::vtu_export vtu_exp2("temperature.vtu", true);
+  vtu_exp2.exporting(mft);
+  vtu_exp2.write_point_data(mft, THETA, "Temperature");
+
   cout << "mayavi2 -d temperature.vtk -f WarpScalar -m Surface\n" << endl;
 
-  getfem::vtk_export exp3("electric_potential.vtk", false);
-  exp3.exporting(mft);
-  exp3.write_point_data(mft, V, "Electric potential");
+  getfem::vtk_export vtk_exp3("electric_potential.vtk", false);
+  vtk_exp3.exporting(mft);
+  vtk_exp3.write_point_data(mft, V, "Electric potential");
+
+  getfem::vtu_export vtu_exp3("electric_potential.vtu", true);
+  vtu_exp3.exporting(mft);
+  vtu_exp3.write_point_data(mft, V, "Electric potential");
+
   cout << "mayavi2 -d electric_potential.vtk -f WarpScalar -m Surface\n"
        << endl;