statistics_while_running_XZ.cpp File Reference

#include "scr/Chute.h"
#include "scr/StatisticsVector.h"

Classes
class	statistics_while_running< T >

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

Function Documentation

main()

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

{
    //select file name to restart from
    string name;
    if (argc>1) {name=argv[1]; argv++; argc--;} else name="../ini_statistics/H10A22L1M0.5B0.5";
    //select time interval
    Mdouble tint=.2; 
    if (argc>1) {tint=atof(argv[1]); argv++; argc--; cout << "tint=" << tint << endl;}
    argv--; argc++;
    //load restart data
    cout << "restart data: " << name << ".restart" << endl;
    statistics_while_running<XZ> problem;
    problem.setName(name.c_str());
    problem.load_restart_data();
    problem.restartFile.setFileType(FileType::ONE_FILE);
    problem.writeRestartFile();
    problem.auto_set_z();
    //keep file name but create files in the local directory, i.e. remove folder
    cout << "old name: " << problem.getName() << endl;
    size_t found=name.find_last_of("/\\");
    problem.setName(name.substr(found+1).c_str());
    cout << "new name: " << problem.getName() << endl;
    //set output to minimum
    problem.dataFile.setFileType(FileType::NO_FILE);
    problem.restartFile.setFileType(FileType::NO_FILE);
    problem.fStatFile.setFileType(FileType::NO_FILE);
    problem.eneFile.setFileType(FileType::ONE_FILE);
    //set statistical parameters
    problem.setDoPeriodicWalls(false);
    //~ problem.setZMinStat(-1);
    problem.setNX(50);
    problem.setNZ(10);
    problem.setCGWidth(.1);
    problem.setSaveCount(25);
    problem.set_infiniteStressForFixedParticles(true);
    problem.setCGTimeMin(problem.getTime());
    problem.setTimeMax(problem.getTime()+tint);
    //solve and create live statistics
    problem.readStatArguments(argc, argv);
    problem.write(std::cout,false);
    problem.solve();
}