Homec4science
Diffusion GooseFEM (main) Browse

rGOOSEFEM/include/GooseFEM/ParaView.h
6adf6418a717main

Last Change
Raw File

ParaView.h
View Options

/* =================================================================================================
(c - GPLv3) T.W.J. de Geus (Tom) | tom@geus.me | www.geus.me | github.com/tdegeus/GooseFEM
================================================================================================= */
#ifndef GOOSEFEM_PARAVIEW_H
#define GOOSEFEM_PARAVIEW_H
// -------------------------------------------------------------------------------------------------
#include "config.h"
#include <fstream>
#ifndef GOOSEFEM_NO_HIGHFIVE
#include <highfive/H5Easy.hpp>
#endif
// =================================================================================================
namespace GooseFEM {
namespace ParaView {
namespace HDF5 {
// -------------------------------------------------------------------------------------------------
inline std::string join(const std::vector<std::string>& lines, const std::string& sep="\n");
// -------------------------------------------------------------------------------------------------
inline std::string indent(size_t n);
// -------------------------------------------------------------------------------------------------
xt::xtensor<double,2> as3d(const xt::xtensor<double,2> &data);
// -------------------------------------------------------------------------------------------------
// http://xdmf.org/index.php/XDMF_Model_and_Format
enum class ElementType {
Triangle,
Quadrilateral,
Hexahedron
};
// -------------------------------------------------------------------------------------------------
enum class AttributeType {
Cell,
Node
};
// -------------------------------------------------------------------------------------------------
inline ElementType convert(GooseFEM::Mesh::ElementType type);
// -------------------------------------------------------------------------------------------------
inline std::string to_string(ElementType type);
inline std::string to_string(AttributeType type);
// -------------------------------------------------------------------------------------------------
class Connectivity
{
public:
// constructors
Connectivity() = default;
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
Connectivity(
const H5Easy::File& data,
const std::string& dataset,
GooseFEM::Mesh::ElementType type);
#endif
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires opening & reading the HDF5-file)
Connectivity(
const std::string& fname,
const std::string& dataset,
GooseFEM::Mesh::ElementType type);
#endif
Connectivity(
const std::string& fname,
const std::string& dataset,
GooseFEM::Mesh::ElementType type,
const std::vector<size_t>& shape);
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
Connectivity(
const H5Easy::File& data,
const std::string& dataset,
ElementType type);
#endif
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires opening & reading the HDF5-file)
Connectivity(
const std::string& fname,
const std::string& dataset,
ElementType type);
#endif
Connectivity(
const std::string& fname,
const std::string& dataset,
ElementType type,
const std::vector<size_t>& shape);
// info
size_t nelem() const;
size_t nne() const;
std::vector<size_t> shape() const;
std::string fname() const;
// check
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires opening & reading the HDF5-file)
void checkShape();
#endif
// return lines of XDMF file (indented, indentation starts at 0)
std::vector<std::string> xdmf(size_t indent=4) const;
private:
// get shape from opened HDF5-file
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
void readShape(const H5Easy::File& data);
#endif
// back-end for constructors (allows single implementation for several overloads)
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
void init(
const H5Easy::File& data,
const std::string& dataset,
ElementType type);
#endif
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires opening & reading the HDF5-file)
void init(
const std::string& fname,
const std::string& dataset,
ElementType type);
#endif
void init(
const std::string& fname,
const std::string& dataset,
ElementType type,
const std::vector<size_t>& shape);
// internal data
ElementType m_type;
std::string m_fname;
std::string m_dataset;
std::vector<size_t> m_shape;
#ifndef GOOSEFEM_NO_HIGHFIVE
bool m_verified=false; // if true: shape read from file, not need to check
#endif
};
// -------------------------------------------------------------------------------------------------
class Coordinates
{
public:
// constructors
Coordinates() = default;
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
Coordinates(
const H5Easy::File& data,
const std::string& dataset);
#endif
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
Coordinates(
const std::string& fname,
const std::string& dataset); // (requires reading the HDF5)
#endif
Coordinates(
const std::string& fname,
const std::string& dataset,
const std::vector<size_t>& shape);
// info
size_t nnode() const;
size_t ndim() const;
std::vector<size_t> shape() const;
std::string fname() const;
// check
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires opening & reading the HDF5-file)
void checkShape();
#endif
// return lines of XDMF file (indented, indentation starts at 0)
std::vector<std::string> xdmf(size_t indent=4) const;
private:
// get shape from opened HDF5-file
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
void readShape(const H5Easy::File& data);
#endif
// internal data
std::string m_fname;
std::string m_dataset;
std::vector<size_t> m_shape;
#ifndef GOOSEFEM_NO_HIGHFIVE
bool m_verified=false; // if true: shape read from file, not need to check
#endif
};
// -------------------------------------------------------------------------------------------------
class Attribute
{
public:
// constructors
Attribute() = default;
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
Attribute(
const H5Easy::File& data,
const std::string& dataset,
const std::string& name,
AttributeType type);
#endif
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires opening & reading the HDF5-file)
Attribute(
const std::string& fname,
const std::string& dataset,
const std::string& name,
AttributeType type);
#endif
Attribute(
const std::string& fname,
const std::string& dataset,
const std::string& name,
AttributeType type,
const std::vector<size_t>& shape);
// info
std::vector<size_t> shape() const;
std::string fname() const;
// check
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires opening & reading the HDF5-file)
void checkShape();
#endif
std::vector<std::string> xdmf(size_t indent=4) const;
private:
// get shape from opened HDF5-file
#ifndef GOOSEFEM_NO_HIGHFIVE // (requires reading the HDF5-file)
void readShape(const H5Easy::File& data);
#endif
// internal data
AttributeType m_type;
std::string m_fname;
std::string m_dataset;
std::string m_name;
std::vector<size_t> m_shape;
#ifndef GOOSEFEM_NO_HIGHFIVE
bool m_verified=false; // if true: shape read from file, not need to check
#endif
};
// -------------------------------------------------------------------------------------------------
class Mesh
{
public:
// constructors
Mesh() = default;
Mesh(
const Connectivity& conn,
const Coordinates& coor);
// return lines of XDMF file (indented, indentation starts at 0)
std::vector<std::string> xdmf(size_t indent=4) const;
// write readable XDMF file
void write(const std::string& fname, size_t n_indent=4) const;
private:
// internal data
Connectivity m_conn;
Coordinates m_coor;
};
// -------------------------------------------------------------------------------------------------
class Increment
{
public:
// constructors
Increment() = default;
Increment(
const Connectivity& conn,
const Coordinates& coor);
Increment(
const Connectivity& conn,
const Coordinates& coor,
const std::vector<Attribute>& attr);
// add attribute to this increment
void push_back(const Connectivity& data);
void push_back(const Coordinates& data);
void push_back(const Attribute& data);
// return lines of XDMF file (indented, indentation starts at 0)
std::vector<std::string> xdmf(size_t indent=4) const;
private:
// internal data
std::vector<Connectivity> m_conn;
std::vector<Coordinates> m_coor;
std::vector<Attribute> m_attr;
};
// -------------------------------------------------------------------------------------------------
class TimeSeries
{
public:
// constructors
TimeSeries() = default;
TimeSeries(
const std::vector<Increment>& data);
// add increment
void push_back(const Increment& data);
// return lines of XDMF file (indented, indentation starts at 0)
std::vector<std::string> xdmf(size_t indent=4) const;
// write readable XDMF file
void write(const std::string& fname, size_t n_indent=4) const;
private:
// internal data
std::vector<Increment> m_data;
};
// -------------------------------------------------------------------------------------------------
}}} // namespace ...
// =================================================================================================
#include "ParaView.hpp"
// =================================================================================================
#endif
c4science · Help