CGHandlerRadialSelfTest.cpp File Reference

#include "Mercury3D.h"
#include "StatisticsVector.h"
#include "Walls/InfiniteWall.h"
#include "CG/CG.h"
#include "CG/TimeSmoothedCG.h"
#include "CG/TimeAveragedCG.h"
#include "CG/Functions/Lucy.h"
#include "CG/Functions/Linear.h"
#include "CG/Functions/Heaviside.h"
#include "CG/Functions/Gauss.h"
#include <cmath>
#include <iostream>
#include <iomanip>
#include <Species/LinearViscoelasticSpecies.h>

Classes
class	CGHandlerSelfTest
	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, char *argv[])

Function Documentation

main()

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

{
    //declare the DPM problem and set the name
    CGHandlerSelfTest problem;
    problem.setTimeMax(1.0);
    problem.setSaveCount(999999);
 
    problem.cgHandler.copyAndAddObject(CG<CGCoordinates::X,CGFunctions::Gauss>());
    problem.cgHandler.copyAndAddObject(CG<CGCoordinates::R,CGFunctions::Gauss>());
//    problem.cgHandler.copyAndAddObject(CG<CGCoordinates::XZ,CGFunctions::Gauss>());
//    problem.cgHandler.copyAndAddObject(CG<CGCoordinates::RZ,CGFunctions::Gauss>());
 
    for (auto cg : problem.cgHandler)
    {
        cg->setN(100);
        cg->setWidth(0.5/3.01);
        cg->setTimeMin(problem.getTimeMax());
    }
 
    //run the simulation
    problem.solve(argc, argv);
 
    helpers::writeToFile("CGHandlerRadialSelfTest.0.gnu","p 'CGHandlerRadialSelfTest.0.stat' u 2:($3*16), 'CGHandlerRadialSelfTest.1.stat' u 2:($3*2*pi*$2*4)");
    helpers::writeToFile("CGHandlerRadialSelfTest.1.gnu","p 'CGHandlerRadialSelfTest.0.stat' u 2:($14*16), 'CGHandlerRadialSelfTest.1.stat' u 2:($14*2*pi*$2*4)");
 
    //helpers::gnuplot("p 'CGHandlerRadialSelfTest.0.stat' u 2:($3*3), 'CGHandlerRadialSelfTest.1.stat' u 2:($3*2*pi*$2)\n");
}

References problem, and helpers::writeToFile().