CGBasicSelfTest.cpp File Reference

#include "Mercury3D.h"
#include "Species/LinearViscoelasticSpecies.h"
#include "CG/CG.h"
#include "CG/TimeSmoothedCG.h"
#include "CG/TimeAveragedCG.h"

Classes
class	CGBasicSelfTest
	Tests if the different CG templates work correctly. 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
    CGBasicSelfTest problem;
    problem.setTimeMax(0.21); //a total of 5000 time steps
 
    //int n1 = 1000, n2 = 100, n3=20;
    size_t n1 = 1000, n2 = 33, n3=10;
    //define different coarse-graining objects (which should all result in the same mean values)
    auto cg0 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::O>());
    cg0->setTimeMin(0.21);
    auto cg1 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::Z,CGFunctions::Gauss>());
    cg1->setN(n1);
    cg1->setWidth(0.5/3.01);
    cg1->setTimeMin(0.21);
    auto cg2 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::YZ,CGFunctions::Gauss>());
    cg2->setN(n2);
    cg2->setWidth(0.5/3.01);
    cg2->setTimeMin(0.21);
    auto cg3 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::XYZ,CGFunctions::Gauss>());
    cg3->setN(n3);
    cg3->setWidth(0.5/3.01);
    cg3->setTimeMin(0.21);
    auto cg4 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::Z,CGFunctions::Heaviside>());
    cg4->setNZ(n1);
    cg4->setWidth(0.5/2.01);
    cg4->setTimeMin(0.21);
    auto cg5 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::YZ,CGFunctions::Heaviside>());
    cg5->setN(n2);
    cg5->setWidth(0.5/2.01);
    cg5->setTimeMin(0.21);
    auto cg6 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::XYZ,CGFunctions::Heaviside>());
    cg6->setN(n3);
    cg6->setWidth(0.5/2.01);
    cg6->setTimeMin(0.21);
    auto cg7 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::Z,CGFunctions::Lucy>());
    cg7->setNZ(n1);
    cg7->setWidth(0.5/2.01);
    cg7->setTimeMin(0.21);
    auto cg8 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::YZ,CGFunctions::Lucy>());
    cg8->setN(n2);
    cg8->setWidth(0.5/2.01);
    cg8->setTimeMin(0.21);
    auto cg9 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::XYZ,CGFunctions::Lucy>());
    cg9->setN(n3);
    cg9->setWidth(0.5/2.01);
    cg9->setTimeMin(0.21);
    auto cg10 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::X,CGFunctions::Lucy>());
    cg10->setN(n1);
    cg10->setWidth(0.5/2.01);
    cg10->setTimeMin(0.21);
    auto cg11 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::X,CGFunctions::Gauss>());
    cg11->setN(n1);
    cg11->setWidth(0.5/3.01);
    cg11->setTimeMin(0.21);
    auto cg12 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::Y,CGFunctions::Lucy>());
    cg12->setN(n1);
    cg12->setWidth(0.5/2.01);
    cg12->setTimeMin(0.21);
    auto cg13 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::Y,CGFunctions::Gauss>());
    cg13->setN(n1);
    cg13->setWidth(0.5/3.01);
    cg13->setTimeMin(0.21);
 
    //run the simulation
    problem.solve(argc, argv);
 
    helpers::writeToFile("CGBasicSelfTest.1D.gnu","p 'CGBasicSelfTest.1.stat' u 2:22 w d, 'CGBasicSelfTest.4.stat' u 2:22 w d, 'CGBasicSelfTest.7.stat' u 2:22 w d");
    helpers::writeToFile("CGBasicSelfTest.2D.gnu","sp 'CGBasicSelfTest.8.stat' u 2:3:23 w d");
}

References problem, and helpers::writeToFile().