#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 <CMakeDefinitions.h>

Classes
class	ChangingTOIParticle

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

Variables
Mdouble	f_min = -20

Mdouble	f_max = 20

Mdouble	margin = 0

Mdouble	av_min = -0

Mdouble	av_max = 0

Mdouble	tv_min = -100

Mdouble	tv_max = 100

int	N_att = 2000

std::vector< int >	D_h

Function Documentation

◆ main()

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

 {
     ChangingTOIParticle problem;
     problem.setNumberOfOMPThreads(helpers::readFromCommandLine(argc, argv, "-omp",4));
     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);
     std::cout <<"coll time"<<collisionTime<<std::endl;
     problem.setTimeStep(collisionTime / 50.0);
  
     // Quick demonstration
     problem.setSaveCount(50);
     problem.setTimeMax(0.1);
  
  
     problem.removeOldFiles();
     problem.solve();
  
  
     // Return the log of the Dzh particles
     std::ofstream D_hyst;  D_hyst.open ("Dzh_hystory.txt");
     for (int i = 0; i<D_h.size(); i+=50){D_hyst <<D_h[i]<<std::endl;}
     D_hyst.close();
  
     return 0;
 }