vtu_writer.hpp

Go to the documentation of this file.

// Creates a .vtu file of the solution, to be visualised by paraview
//
// Author: Jerome Droniou (jerome.droniou@monash.edu)
//
 
/*
*
*  This library was developed around HHO methods, although some parts of it have a more
* general purpose. If you use this code or part of it in a scientific publication, 
* please mention the following book as a reference for the underlying principles
* of HHO schemes:
*
* The Hybrid High-Order Method for Polytopal Meshes: Design, Analysis, and Applications. 
*  D. A. Di Pietro and J. Droniou. Modeling, Simulation and Applications, vol. 19. 
*  Springer International Publishing, 2020, xxxi + 525p. doi: 10.1007/978-3-030-37203-3. 
*  url: https://hal.archives-ouvertes.fr/hal-02151813.
*
*/
 
 
 
#ifndef VTU_WRITER_HPP
#define VTU_WRITER_HPP
#include <mesh.hpp>
#include <cstring>
#include <vector>
#include <Eigen/Dense>
using Eigen::Vector2d;
 
namespace HArDCore2D {
 
class VtuWriter {
public:
    VtuWriter(const Mesh* mesh); 
 
    bool write_to_vtu(
          const std::string & filename, 
          const std::vector<Eigen::VectorXd> & sol_vrtx,   
                const std::vector<std::string> & sol_names = {}     
          );
 
    bool write_to_vtu(std::string file_name); 
 
    inline bool write_to_vtu(
                  const std::string & filename, 
                  const Eigen::VectorXd & sol_vrtx,  
                  const std::string & sol_name = "solution"  
                  )
          {
            return write_to_vtu(filename, std::vector<Eigen::VectorXd> {sol_vrtx}, std::vector<std::string> {sol_name});
          }
 
private:
    void write_vertices(FILE* pFile);
    void write_vertices(FILE* pFile, const Eigen::VectorXd & data);
    void write_header(FILE* pFile); 
    void write_cells(FILE* pFile); 
    void write_solution(FILE* pFile, const Eigen::VectorXd & sol, const std::string & name); 
    void write_footer(FILE* pFile); 
 
    const Mesh* _mesh;
 
};
 
}  // namespace HArDCore2D
 
#endif /* VTU_WRITER_HPP */