SquarePackingSelfTest.cpp File Reference

#include "Mercury3D.h"
#include "StatisticsVector.h"
#include "Particles/BaseParticle.h"
#include "Walls/InfiniteWall.h"
#include <cmath>
#include <iostream>
#include <iomanip>
#include <Species/LinearViscoelasticSpecies.h>

Go to the source code of this file.

Classes
class	SquarePacking

Functions
int	main (int argc UNUSED, char *argv[] UNUSED)

Function Documentation

int main	(	int argc	UNUSED,
		char *argv[]	UNUSED
	)

Definition at line 82 of file SquarePackingSelfTest.cpp.

References BaseHandler< T >::copyAndAddObject(), DPMBase::getTimeMax(), SquarePacking::N, constants::pi, StatisticsVector< T >::setCGTimeMin(), StatisticsVector< T >::setCGWidth(), DPMBase::setGravity(), StatisticsVector< T >::setN(), File::setName(), Files::setName(), DPMBase::setParticleDimensions(), Files::setSaveCount(), DPMBase::setSystemDimensions(), DPMBase::setTimeMax(), StatisticsVector< T >::setTimeMaxStat(), DPMBase::setTimeStep(), DPMBase::solve(), DPMBase::speciesHandler, Files::statFile, and StatisticsVector< T >::statistics_from_fstat_and_data().

{
    SquarePacking problem;
    problem.setName("SquarePackingSelfTest");
    auto species=problem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
    //set the number of particles
    problem.N=5;
    problem.setSystemDimensions(2);
    problem.setParticleDimensions(2);
    species->setDensity(6.0/constants::pi*2./3.);
    //species->setDensity(1.9098593*2/3);
    problem.setGravity(Vec3D(0.,-1,0.));
    species->setCollisionTimeAndRestitutionCoefficient(.01,.1,1.);
    problem.setTimeStep(.0002);
    problem.setTimeMax(1.0);
    problem.setSaveCount(1000);

    problem.setCGTimeMin(problem.getTimeMax()*.98);
    problem.setTimeMaxStat(1e20);

    std::cout << "Relax the packing" << std::endl;
    problem.solve();

    std::cout << "Get statistics" << std::endl;
    StatisticsVector<Y> stats("SquarePackingSelfTest");
    stats.statFile.setName("SquarePackingSelfTest.Y.stat");
    double n = 500;
    stats.setN(n);
    stats.setCGWidth(.1);
    stats.setCGTimeMin(problem.getTimeMax()*.98);
    stats.setTimeMaxStat(1e20);
    //stats.verbose();
    //stats.set_boundedDomain(true);
    stats.statistics_from_fstat_and_data();

    std::cout << "Get fully averaged statistics" << std::endl;
    StatisticsVector<O> statsO("SquarePackingSelfTest");
    statsO.statFile.setName("SquarePackingSelfTest.O.stat");
    statsO.setCGTimeMin(problem.getTimeMax()*.98);
    statsO.setTimeMaxStat(1e20);
    statsO.statistics_from_fstat_and_data();
    // should give you Density 1
}