CGForceDensityUnitTest.cpp File Reference

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

Classes
class	CGForceDensityUnitTest
	Tests if the force density between two species has the correct sign assigned. 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
    CGForceDensityUnitTest problem;
    problem.setTimeMax(0.1); //a total of 5000 time steps
    problem.setTimeStep(0.1);
    problem.setFileType(FileType::NO_FILE);
 
    //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.1);
    cg0->selectSpecies(0);
    auto cg1 = problem.cgHandler.copyAndAddObject(CG<CGCoordinates::O>());
    cg1->setTimeMin(0.1);
    cg1->selectSpecies(1);
    
    //run the simulation
    problem.solve(argc, argv);
 
    Vec3D cg0ForceDensityResult = static_cast<CG<CGCoordinates::O,CGFunctions::Lucy,CGFields::StandardFields>*>(cg0)->evaluateAverage().getInteractionForceDensity();
    Vec3D cg1ForceDensityResult = static_cast<CG<CGCoordinates::O,CGFunctions::Lucy,CGFields::StandardFields>*>(cg1)->evaluateAverage().getInteractionForceDensity();
 
    logger(INFO, "ForceDensity species 0 %, forceDensity species 1 %, sum %", cg0ForceDensityResult, cg1ForceDensityResult, (cg0ForceDensityResult + cg1ForceDensityResult));
    logger.assert_always(cg0ForceDensityResult.getLength() > 10 && cg1ForceDensityResult.getLength() > 10, "evaluated force density is too small");
    logger.assert_always(cg0ForceDensityResult.X < -1 && cg1ForceDensityResult.X > 1, "force density does not point in the right direction");
    logger.assert_always( (cg0ForceDensityResult + cg1ForceDensityResult).getLength() < 1e-10, "evaluated force densities do not sum to 0");
}

References e(), Vec3D::getLength(), INFO, logger, NO_FILE, problem, and Vec3D::X.