vibratedBed.cpp File Reference

#include <iomanip>
#include "SilbertPeriodic.h"

Classes
class	vibratedBed

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

Function Documentation

main()

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

{
    vibratedBed problem;
    //set case, height, angle to given or default values
    std::vector<Mdouble> studyNumber;
    studyNumber.resize(3);
    const unsigned int numArgs = 5;
 
    if (argc>numArgs)
    {
        studyNumber[0]=atof(argv[1]);
        studyNumber[1]=atof(argv[2]);
        studyNumber[2]=atof(argv[3]);
        problem.setAmplitude(atof(argv[4]));
        problem.setFrequency(atof(argv[5]));
        problem.setTimeStep(problem.getTimeStep()*2.0);
        problem.setSaveCount(50.0/problem.getTimeStep());
        problem.eneFile.setSaveCount(0.5/problem.getTimeStep());
        problem.setTimeMax(2000.0);
        problem.dataFile.setFileType(FileType::ONE_FILE);
        problem.restartFile.setFileType(FileType::MULTIPLE_FILES_PADDED);
        problem.run(studyNumber,argc-numArgs,argv+numArgs);
    }
    else
    {
        logger(ERROR, "Not enough input arguments given (./vibratedBed $study $height $angle $amplitude $frequency);\n"
                      "using demo values (equivalent to ./vibratedBed 4 5 10 0.02 1 -tmax 1000 -dt 0.001)");
        studyNumber[0] = 4;
        studyNumber[1] = 0.1;
        studyNumber[2] = 10;
        problem.setAmplitude(0.01);
        problem.setFrequency(1.0);
        problem.setRoughBottomType(MONOLAYER_DISORDERED);
        
        problem.setTimeStep(problem.getTimeStep() * 10.0);
        problem.setTimeMax(2.0);
        problem.setSaveCount(1e1);
        problem.eneFile.setSaveCount(1e2);
        problem.dataFile.setFileType(FileType::ONE_FILE);
        problem.run(studyNumber, 1, argv);
    }
}

References ERROR, logger, MONOLAYER_DISORDERED, MULTIPLE_FILES_PADDED, ONE_FILE, and problem.