ExtremeOverlapWithWallsUnitTest.cpp File Reference

#include "DPMBase.h"
#include "Particles/BaseParticle.h"
#include "Walls/InfiniteWall.h"
#include <iostream>
#include <Species/LinearViscoelasticSpecies.h>
#include <Logger.h>

Go to the source code of this file.

Classes
class	ExtremeOverlapWithWallsUnitTest
	Compresses 2 particles (vertically) until they have an extreme overlap. More...

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

Function Documentation

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

Start off my solving the default problem

Definition at line 105 of file ExtremeOverlapWithWallsUnitTest.cpp.

References BaseHandler< T >::copyAndAddObject(), Files::dataFile, FATAL, BaseHandler< T >::getObject(), BaseInteractable::getPosition(), Vec3D::isEqualTo(), logger, NO_FILE, ONE_FILE, DPMBase::particleHandler, DPMBase::setDimension(), File::setFileType(), Files::setFileType(), Files::setName(), Files::setSaveCount(), DPMBase::setTimeMax(), DPMBase::setTimeStep(), DPMBase::solve(), and DPMBase::speciesHandler.

{
    ExtremeOverlapWithWallsUnitTest OverlapProblem;
    auto species = OverlapProblem.speciesHandler.copyAndAddObject(LinearViscoelasticSpecies());
    species->setDensity(2000);
    species->setStiffness(1e5);
    species->setDissipation(1e-2);
    OverlapProblem.setName("ExtremeOverlapWithWallsUnitTest");
    OverlapProblem.setFileType(FileType::NO_FILE);
    OverlapProblem.dataFile.setFileType(FileType::ONE_FILE);
    OverlapProblem.setSaveCount(1000);
    OverlapProblem.setDimension(3);
    OverlapProblem.setTimeStep(1e-7);
    OverlapProblem.setTimeMax(0.01);
    // OverlapProblem.dataFile.getFstream().precision(10);
    OverlapProblem.solve(argc,argv);
    
    Vec3D position = OverlapProblem.particleHandler.getObject(0)->getPosition();
    Vec3D positionToCompare = Vec3D(0, -0.0004900098000491, 0.004008957617341);
    if (!position.isEqualTo(positionToCompare, 1e-10))
        logger(FATAL,"Large particle is in the wrong position. It is at % and should be %",position,positionToCompare);
    else
        std::cout << "Test passed" << std::endl;
    position = OverlapProblem.particleHandler.getObject(1)->getPosition();
    positionToCompare = Vec3D(0, -0.0004900098000491, 0.009004578797915);
    //std::cout.precision(13); std::cout << position << std::endl;
    if (!position.isEqualTo(positionToCompare, 1e-10))
        logger(FATAL,"Large particle is in the wrong position. It is at % and should be %",position,positionToCompare);
    else
        std::cout << "Test passed" << std::endl;

    return 0;
}