ParticleCreator.h

Go to the documentation of this file.

/*
 * ParticleCreator.h
 *
 *  Created on: 27 Feb 2014
 *      Author: ttitsche
 */

#pragma once


#include <vector>
#include "nuto/geometryConcrete/Specimen.h"


namespace NuTo
{
class InputReader;

class ParticleCreator
{
public:
    ParticleCreator(NuTo::Specimen rSpecimen, const double rShrinkage, const long rNumMaxTries = 10000000);

    Eigen::MatrixXd CreateSpheresInSpecimen(const double rRelParticleVolume, const Eigen::MatrixXd& rGradingCurve,
                                            const double rRelativeDistance, const double rAbsoluteDistance,
                                            const int rSeed, const Eigen::MatrixXd& rSpheresBoundary) const;

    Eigen::MatrixXd PerformTakePhase(const Eigen::MatrixXd& rGradingCurve, const Eigen::MatrixXd& rSpheresBoundary,
                                     const double rRelParticleVolume) const;

    void PerformPlacePhase(Eigen::MatrixXd& rParticles, const double rRelativeDistance,
                           const double rAbsoluteDistance) const;

private:
    void InsertParticleIntoBox(const Eigen::MatrixXd& rParticles, const int rTheParticle,
                               std::vector<std::vector<int>>& rSubBox, const Eigen::Vector3i& rNSubBox,
                               const Eigen::Vector3d& rLSubBox) const;

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

    bool CollidesWithBoundary(const Eigen::Vector4d& rParticle, const double rRelativeDistance,
                              const double rAbsoluteDistance) const;

    const std::vector<double> GetSizeClasses(const Eigen::MatrixXd& rParticles) const;

    const std::vector<double> GetNumParticlesPerSizeClass(const Eigen::MatrixXd& rParticles,
                                                          const std::vector<double>& rSizes) const;
    double GetVolume(double radius) const;

    NuTo::Specimen mSpecimen;
    const double mShrinkage;
    const long mNumMaxTries;

    double mVolume;
};

} /* namespace NuTo */