PeriodicWallsWithSlidingFrictionUnitTest.cpp File Reference

#include "Species/LinearViscoelasticSlidingFrictionSpecies.h"
#include "Species/LinearViscoelasticFrictionSpecies.h"
#include "Mercury3D.h"
#include "Boundaries/PeriodicBoundary.h"
#include "Particles/BaseParticle.h"
#include <Logger.h>

Go to the source code of this file.

Classes
class	PeriodicWallsWithSlidingFrictionUnitTest
	This test is a UnitTest for: Periodic Particles in combination with HGrid Tangential Springs over periodic contacts and during transitions from periodic to normal and reverse Removal of Particles while tangential spring information has to be conserved. More...

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

Function Documentation

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

Start off by solving the default problem

Definition at line 174 of file PeriodicWallsWithSlidingFrictionUnitTest.cpp.

References FATAL, BaseHandler< T >::getObject(), BaseInteractable::getPosition(), BaseInteractable::getVelocity(), Vec3D::isEqualTo(), logger, DPMBase::particleHandler, Files::setName(), and DPMBase::solve().

{    
    PeriodicWallsWithSlidingFrictionUnitTest periodicWallsWithSlidingFrictionUnitTest;
    periodicWallsWithSlidingFrictionUnitTest.solve();
    
    periodicWallsWithSlidingFrictionUnitTest.setName("PeriodicWallsWithSlidingFrictionUnitTest");

    BaseParticle* PNormal;
    BaseParticle* PCase1;
    BaseParticle* PCase2;
    Vec3D shift;
    PNormal = periodicWallsWithSlidingFrictionUnitTest.particleHandler.getObject(4);
    PCase1 = periodicWallsWithSlidingFrictionUnitTest.particleHandler.getObject(0);
    PCase2 = periodicWallsWithSlidingFrictionUnitTest.particleHandler.getObject(2);

    shift = Vec3D(-0.58, -0.4, 0.0);
    if (!PCase1->getPosition().isEqualTo(PNormal->getPosition() + shift, 1e-10))
    {
        logger(FATAL, "E0 The particle is in the wrong position. It is %, however is should be %", PCase1->getPosition(), PNormal->getPosition() + shift);
    }
    if (!PCase1->getVelocity().isEqualTo(PNormal->getVelocity(), 1e-10))
    {
        logger(FATAL, "E1 The particle has the wrong velocity. It is %, however is should be %", PCase1->getVelocity(), PNormal->getVelocity());
    }

    shift = Vec3D(0.0, -0.779, 0.0);
    if (!PCase2->getPosition().isEqualTo(PNormal->getPosition() + shift, 1e-10))
    {
        logger(FATAL, "E2 The particle is in the wrong position. It is %, however is should be %", PCase2->getPosition(), PNormal->getPosition() + shift);
    }
    if (!PCase2->getVelocity().isEqualTo(PNormal->getVelocity(), 1e-10))
    {
        logger(FATAL, "E3 The particle has the wrong velocity. It is %, however is should be %", PCase2->getVelocity(), PNormal->getVelocity());
    }

    PNormal = periodicWallsWithSlidingFrictionUnitTest.particleHandler.getObject(5);
    PCase1 = periodicWallsWithSlidingFrictionUnitTest.particleHandler.getObject(1);
    PCase2 = periodicWallsWithSlidingFrictionUnitTest.particleHandler.getObject(3);

    shift = Vec3D(1.0 - 0.58, -0.4, 0.0);
    if (!PCase1->getPosition().isEqualTo(PNormal->getPosition() + shift, 1e-10))
    {
        std::cout << PNormal->getPosition() - PCase1->getPosition() << " " << PNormal->getPosition() + PCase1->getPosition() << std::endl;
        logger(FATAL, "E4 The particle is in the wrong position. It is %, however is should be %", PCase1->getPosition(), PNormal->getPosition() + shift);
    }
    if (!PCase1->getVelocity().isEqualTo(PNormal->getVelocity(), 1e-10))
    {
        logger(FATAL, "E5 The particle has the wrong velocity. It is %, however is should be %", PCase1->getVelocity(), PNormal->getVelocity());
    }

    shift = Vec3D(0.0, 1.0 - 0.779, 0.0);
    if (!PCase2->getPosition().isEqualTo(PNormal->getPosition() + shift, 1e-10))
    {
        logger(FATAL, "E6 The particle is in the wrong position. It is %, however is should be %", PCase2->getPosition(), PNormal->getPosition() + shift);
    }
    if (!PCase2->getVelocity().isEqualTo(PNormal->getVelocity(), 1e-10))
    {
        logger(FATAL, "E7 The particle has the wrong velocity. It is %, however is should be %", PCase2->getVelocity(), PNormal->getVelocity());
    }
}