GmshWriter.h

Go to the documentation of this file.

#pragma once

#include <iosfwd>
#include <string>
#include <Eigen/Core>

namespace NuTo
{

class GmshWriter
{
public:
    struct Box2D
    {
        Box2D(double lx, double ly)
            : mStart(Eigen::Vector2d::Zero())
            , mEnd(Eigen::Vector2d(lx, ly))
        {
        }
        Box2D(Eigen::Vector2d start, Eigen::Vector2d end)
            : mStart(start)
            , mEnd(end)
        {
        }

        Eigen::Vector2d mStart;
        Eigen::Vector2d mEnd;
    };

    struct Box3D
    {
        Box3D(double lx, double ly, double lz)
            : mStart(Eigen::Vector3d::Zero())
            , mEnd(Eigen::Vector3d(lx, ly, lz))
        {
        }
        Box3D(Eigen::Vector3d start, Eigen::Vector3d end)
            : mStart(start)
            , mEnd(end)
        {
        }

        Eigen::Vector3d mStart;
        Eigen::Vector3d mEnd;
    };

    struct Cylinder
    {
        double mRadius;
        double mHeight;
    };

    struct Options
    {
        Options(double meshSizeAggregates, double meshSizeMatrix, double interfaceThickness = 0.)
            : interfaceThickness(interfaceThickness)
            , meshSizeAggregates(meshSizeAggregates)
            , meshSizeMatrix(meshSizeMatrix)
        {
        }

        double interfaceThickness = 0;
        double meshSizeAggregates = 1;
        double meshSizeMatrix = 1;

        int meshAlg = 6;
        int meshAlg3D = 4;
        int meshRecombinationAlgorithm = 0;
        int meshOptimize = 2;
        int meshSmoothing = 2;

        std::string physicalGroupMatrix = "Matrix";
        std::string physicalGroupAggregates = "Aggregates";
        std::string physicalGroupInterfaces = "Interfaces";
        std::string additionalOptions = "";
    };

    // Write to output stream

    static void Write(std::ostream& out, Box2D box, const Eigen::MatrixX3d& aggregates, Options opt);
    static void Write(std::ostream& out, Box3D box, const Eigen::MatrixX4d& aggregates, Options opt);
    static void Write(std::ostream& out, Cylinder cylinder, const Eigen::MatrixX4d& aggregates, Options opt);

    // Write to file

    static void Write(std::string filename, Box2D box, const Eigen::MatrixX3d& aggregates, Options opt);
    static void Write(std::string filename, Box3D box, const Eigen::MatrixX4d& aggregates, Options opt);
    static void Write(std::string filename, Cylinder cylinder, const Eigen::MatrixX4d& aggregates, Options opt);
};
} /* NuTo */