fixes for vtu

Author: gtheler

src/mesh/vtu.c

@@ -24,7 +24,7 @@
 
 int feenox_mesh_write_header_vtu(mesh_t *mesh, FILE *fp) {
   fprintf(fp, "<?xml version=\"1.0\"?>\n");
-  fprintf(fp, "<VTKFile type=\"UnstructuredGrid\" version=\"1.0\" byte_order=\"LittleEndian\">\n");
+  fprintf(fp, "<VTKFile type=\"UnstructuredGrid\" version=\"1.0\">\n");
   fprintf(fp, "  <UnstructuredGrid>\n");
   if (mesh->n_cells == 0) {
     feenox_call(feenox_mesh_element2cell(mesh));

src/mesh/write.c

@@ -28,6 +28,10 @@ int feenox_instruction_mesh_write(void *arg) {
   }
 
   mesh_write_t *mesh_write = (mesh_write_t *)arg;
+  
+  // re-initialize flags
+  mesh_write->point_init = 0;
+  mesh_write->cell_init = 0;
 
   // TODO: in parallel runs only print from first processor
   if (feenox.mpi_rank != 0) {

src/parser/parser.c

@@ -2654,8 +2654,8 @@ int feenox_parse_write_results(void) {
       feenox_check_alloc(strcat(file_path, feenox.argv[i]));
     }
     feenox_check_alloc(strcat(file_path, "."));
-    feenox_check_alloc(strcat(file_path, (mesh_write->post_format == post_format_vtu) ? "vtu" :
-                                           ((mesh_write->post_format == post_format_vtk) ? "vtk" : "msh")));
+    feenox_check_alloc(strcat(file_path, (mesh_write->post_format == post_format_gmsh) ? "msh" :
+                                           ((mesh_write->post_format == post_format_vtk) ? "vtk" : "vtu")));
 
     feenox_call(feenox_define_file(file_path, file_path, 0, "w"));
     if ((mesh_write->file = feenox_get_file_ptr(file_path)) == NULL) {