SquarePackingSelfTest.cpp File Reference

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

Classes
class	SquarePacking

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

Function Documentation

main()

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

{
    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);
    
    logger(INFO, "Relax the packing");
    problem.solve();
    
    logger(INFO, "Get statistics");
    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();
    
    logger(INFO, "Get fully averaged statistics");
    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
}

References INFO, logger, n, constants::pi, problem, StatisticsVector< T >::setCGTimeMin(), StatisticsVector< T >::setCGWidth(), StatisticsVector< T >::setN(), File::setName(), StatisticsVector< T >::setTimeMaxStat(), DPMBase::statFile, and StatisticsVector< T >::statistics_from_fstat_and_data().