ExtremeOverlapUnitTest.cpp File Reference

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

Go to the source code of this file.

Classes
class	ExtremeOverlapUnitTest
	Makes sure that the behavior is still sensible if the overlap of two particles grows extremely large. More...

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

Function Documentation

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

Definition at line 97 of file ExtremeOverlapUnitTest.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.

{
    ExtremeOverlapUnitTest OverlapProblem;
    
    //set some contact parameters; the contact is elastic (so we can check 
    //energy conservation), and very soft (so the particles repel each other slowly))
    //To check Energy conservation, use gnuplot:
    //  p 'ExtremeOverlapUnitTest.ene' u 1:($2+$3+$4+$5)
    auto species = OverlapProblem.speciesHandler.copyAndAddObject(LinearViscoelasticSlidingFrictionSpecies());
    species->setDensity(2000);
    species->setStiffness(1e1);
    species->setDissipation(0.0);
    species->setSlidingFrictionCoefficient(1.0);
    species->setSlidingStiffness(2.0/7.0*species->getStiffness());
    species->setSlidingDissipation(2.0/7.0*species->getDissipation());
    
    OverlapProblem.setName("ExtremeOverlapUnitTest");
    OverlapProblem.setFileType(FileType::NO_FILE);
    OverlapProblem.dataFile.setFileType(FileType::ONE_FILE);
    OverlapProblem.setSaveCount(10000);
    OverlapProblem.setDimension(2);
    OverlapProblem.setTimeStep(1e-8);
    OverlapProblem.setTimeMax(0.01);
    OverlapProblem.solve(argc,argv);
    
    Vec3D position = OverlapProblem.particleHandler.getObject(0)->getPosition();
    Vec3D positionToCompare = Vec3D(-1.031389999146e-06, 0.009506389407855, 0);
    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.01010314899993, 0.01487096023445, 0);
    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;
    //std::cout.precision(13); std::cout << position << std::endl;
    return 0;
}