PeriodicClumpSelfTest.cpp File Reference

#include "Mercury3D.h"
#include "Walls/InfiniteWall.h"
#include "Species/LinearViscoelasticFrictionSpecies.h"
#include "Particles/ClumpParticle.h"
#include "../../Clump/ClumpHeaders/ClumpInput.h"
#include "../../Clump/ClumpHeaders/Mercury3DClump.h"
#include <stdlib.h>
#include "Boundaries/PeriodicBoundary.h"

Classes
class	clumpTest

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

Variables
Mdouble	f_min = -4
Mdouble	f_max = 4

Function Documentation

main()

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

{
    clumpTest problem;
    auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
    species->setDensity(1.0); // sets the species type-0 density
    //species->setConstantRestitution(0);
    std::cout<<species->getConstantRestitution()<<std::endl;
    species->setDissipation(0.0);
    species->setStiffness(1e6);
    const Mdouble collisionTime = species->getCollisionTime(problem.getClumpMass());
    problem.setClumpDamping(0);
    problem.setTimeStep(collisionTime / 50.0);
 
    // Quick demonstration
    problem.setSaveCount(100);
    problem.setTimeMax(1);
 
 
 
    problem.removeOldFiles();
    problem.solve();
 
    ClumpParticle* p = dynamic_cast<ClumpParticle*>(problem.particleHandler.getLastObject());
    Vec3D angVel = p->getAngularVelocity();
    Vec3D known_angVel = Vec3D(-1.30004, 19.9932, 0.944699);
 
    Vec3D pos = p->getPosition();
    Vec3D known_pos = Vec3D(-0.201397, -2.37542, -0.284277);
 
 
    helpers::check(angVel,known_angVel,1e-4, "Angular velocity check");
    helpers::check(pos,known_pos,1e-4, "Pebble position");
    
    return 0;
}

References helpers::check(), e(), p, and problem.

Variable Documentation

f_max

Mdouble f_max = 4

Referenced by clumpTest::clumpTest().

f_min

Mdouble f_min = -4

Referenced by clumpTest::clumpTest().