GeometryConcrete.h

Go to the documentation of this file.

/*
 * GeometryConcrete.h
 *
 *  Created on: 4 Sep 2015
 *      Author: ttitsche
 */

#pragma once

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

namespace NuTo
{

class Specimen;
class ParticleHandler;

class GeometryConcrete
{
public:
    enum eGradingCurve
    {
        A16,
        B16,
        C16
    };

    GeometryConcrete();

    virtual ~GeometryConcrete();


    void MaximizeParticleDistance(double rParticleDistance);


    void MaximizeParticleVolumeFraction(double rShrinkage);


    void ExportGmshGeo2D(std::string rGmshFile, double rMeshSize, double rZSlice, double rMinRadius);


    void ExportGmshGeo3D(std::string rGmshFile, double rMeshSize);


    void SetSpecimenBox(double rXs, double rXe, double rYs, double rYe, double rZs, double rZe);
    void SetSpecimenCylinder(double rXs, double rXe, double rYs, double rYe, double rZs, double rZe);
    void SetSpecimenCylinder(double radius, double height);

    void SetGradingCurve(eGradingCurve rGradingCurveEnum, int rNumClasses);

    void SetGradingCurve(const Eigen::MatrixXd& rGradingCurve);


    Eigen::MatrixXd GetParticles(bool rBeforeEDMD);

    void SetParticles(Eigen::MatrixXd);

    void SetAbsoluteGrowthRate(double rAbsoluteGrowthRate);

    void SetRelativeGrowthRate(double rRelativeGrowthRate);

    void SetGmshBinary(const std::string& rGmshBinary);

    void SetInitialTimeBarrier(double rInitialTimeBarrier);

    void SetNumEventsMax(double rNumEventsMax);

    void SetParticleVolumeFraction(double rParticleVolumeFraction);

    void SetRandomVelocityRange(double rRandomVelocityRange);

    void SetSecondsPrint(double rSecondsPrint);

    void SetSecondsWallTimeMax(double rSecondsWallTimeMax);

    void SetSeed(double rSeed);

    bool ContinueOnException() const;

    void SetContinueOnException(bool rContinueOnException);

private:
    void CheckParameters();

    Specimen* mSpecimen = nullptr;
    ParticleHandler* mParticleHandler = nullptr;
    Eigen::MatrixXd mGradingCurve;

    double mSeed = 42;
    double mParticleVolumeFraction = 0.5;

    double mNumEventsMax = 2.e6;
    double mSecondsWallTimeMax = 60.;
    double mSecondsPrint = 5.;
    double mInitialTimeBarrier = 10;

    double mRandomVelocityRange = 1.;
    double mRelativeGrowthRate = 0.0;
    double mAbsoluteGrowthRate = 0.0;

    bool mContinueOnException = false;
};

} /* namespace NuTo */