CubicPackingSelfTest.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	NewtonsCradleSelftest
	In this file a cubic packing of 5^3 particles in a tri-axial box is created and allowed to settle under small gravity. After that Z statistics are calculated. More...

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

Function Documentation

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

should give you Density 1

Definition at line 84 of file CubicPackingSelfTest.cpp.

References BaseHandler< T >::copyAndAddObject(), BaseHandler< T >::getObject(), DPMBase::getTimeMax(), NewtonsCradleSelftest::N, StatisticsVector< T >::setCGShape(), StatisticsVector< T >::setCGTimeMin(), StatisticsVector< T >::setCGWidth(), ParticleSpecies::setDensity(), DPMBase::setGravity(), StatisticsVector< T >::setN(), Files::setName(), Files::setSaveCount(), StatisticsVector< T >::setSuperExact(), DPMBase::setSystemDimensions(), DPMBase::setTimeMax(), StatisticsVector< T >::setTimeMaxStat(), DPMBase::setTimeStep(), StatisticsVector< T >::setVerbosityLevel(), DPMBase::solve(), DPMBase::speciesHandler, and StatisticsVector< T >::statistics_from_fstat_and_data().

{
    NewtonsCradleSelftest problem;
    problem.setName("CubicPackingSelfTest");
    auto species=problem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());

    //set the number of particles
    problem.N = 5;
    problem.setSystemDimensions(3);
    problem.speciesHandler.getObject(0)->setDensity(1.9098593);
    problem.setGravity(Vec3D(0., 0., -1.));
    species->setCollisionTimeAndRestitutionCoefficient(.01, .1, 1.);
    problem.setTimeStep(.0002);
    problem.setTimeMax(0.2);
    problem.setSaveCount(1000);
    problem.solve();

    StatisticsVector<Z> stats("CubicPackingSelfTest");
    double n = 500;
    stats.setN(n);
    stats.setCGWidth(.1);
    stats.setSuperExact(false);
    stats.setCGShape("Gaussian");
    stats.setCGTimeMin(problem.getTimeMax() * .999999);
    stats.setTimeMaxStat(1e20);
    stats.setVerbosityLevel(0);
    stats.statistics_from_fstat_and_data();

}