ConsistencyInitialParticleWallOverlapUnitTest.cpp File Reference

#include <cmath>
#include "Mercury3D.h"
#include "Walls/InfiniteWall.h"
#include "Species/LinearViscoelasticFrictionSpecies.h"

Classes
class	ConsistencyInitialParticleWallOverlapUnitTest

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

Function Documentation

main()

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

{
    // Problem setup
    ConsistencyInitialParticleWallOverlapUnitTest problem;
    problem.setXMax(0.4);
    problem.setYMax(0.6);
    problem.setZMax(0.1);
    problem.setXMin(0.0);
    problem.setYMin(0.2);
    problem.setZMin(0.0);
 
    problem.setGravity(Vec3D{0.,0.,-9.81});
 
    problem.setRadius(1.e-2);
    problem.setInitOverlap(1.e-3);
    problem.setInitPos0({0., 0.4, problem.getRadius() - problem.getInitOverlap()});
    problem.setInitPos1({0.2, 0.4, problem.getRadius() - problem.getInitOverlap()});
 
    problem.setTimeMax(1.e-1);
    problem.setTimeStep(0.005/50.);
 
    problem.setName("ConsistencyInitialParticleWallOverlapUnitTest");
    // Uncomment to check vtu files
    //problem.setParticlesWriteVTK(true);
    //problem.wallHandler.setWriteVTK(true);
    problem.dataFile.setFileType(FileType::NO_FILE);
    problem.restartFile.setFileType(FileType::NO_FILE);
    problem.eneFile.setFileType(FileType::NO_FILE);
    problem.fStatFile.setFileType(FileType::NO_FILE);
 
    problem.solve();
 
    Vec3D initDiff = problem.getInitPos1() - problem.getInitPos0();
    Vec3D currDiff = problem.particleHandler.getObject(1)->getPosition() - problem.particleHandler.getObject(0)->getPosition();
 
    // Check if particle 0 and 1 only moved in z-dir
    helpers::check(problem.getInitPos0().X, problem.particleHandler.getObject(0)->getPosition().X, 1.e-15, "Particle 0 x-movement");
    helpers::check(problem.getInitPos0().Y, problem.particleHandler.getObject(0)->getPosition().Y, 1.e-15, "Particle 0 y-movement");
 
    helpers::check(problem.getInitPos1().X, problem.particleHandler.getObject(1)->getPosition().X, 1.e-15, "Particle 1 x-movement");
    helpers::check(problem.getInitPos1().Y, problem.particleHandler.getObject(1)->getPosition().Y, 1.e-15, "Particle 1 y-movement");
 
    helpers::check((currDiff.X-initDiff.X),0.,1.e-15,"Difference in x-movement");
    helpers::check((currDiff.Y-initDiff.Y),0.,1.e-15,"Difference in y-movement");
    helpers::check((currDiff.Z-initDiff.Z),0.,1.e-15,"Difference in z-movement");
 
    return 0;
}

References helpers::check(), NO_FILE, problem, Vec3D::X, Vec3D::Y, and Vec3D::Z.