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

Classes
class	ChangingTOIParticle

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

Variables
Mdouble	f_min = -4

Mdouble	f_max = 4

Mdouble	z_shift = 2

Mdouble	margin = 2

Mdouble	z_margin = 2.8

Mdouble	av_min = 0

Mdouble	av_max = 0

Mdouble	tv_min = 0

Mdouble	tv_max = 0

int	N_att = 1000

int	SAVECOUNT = 400

Function Documentation

◆ main()

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

 {
     ChangingTOIParticle problem;
     auto species = problem.speciesHandler.copyAndAddObject(LinearViscoelasticFrictionSpecies());
     species->setDensity(1.0); // sets the species type-0 density
     species->setDissipation(50.0);
     species->setStiffness(1e6);
  
     species->setSlidingFrictionCoefficient(0.6);
     species->setSlidingStiffness(5e5);
     species->setRollingFrictionCoefficient(0.0);
     species->setRollingStiffness(5e5);
  
     const Mdouble collisionTime = species->getCollisionTime(problem.getClumpMass());
     problem.setClumpDamping(0);
     problem.setTimeStep(collisionTime / 50.0);
     problem.setSaveCount(SAVECOUNT);
     problem.setHGridMaxLevels(1);
     //problem.setTimeMax(48.0);
     problem.setTimeMax(0.08);
  
     problem.removeOldFiles();
     problem.solve();
     return 0;
 }