hstop_demo.cpp File Reference

#include "SilbertPeriodic.h"

Classes
class	SilbertHstop

Functions
bool	PointIsAboveCurve (Mdouble h, Mdouble a, int study_num)
int	HstopCurve (string cmd UNUSED, int study_num, Mdouble h, Mdouble hMax, Mdouble a)
int	main (int argc UNUSED, char *argv[])

Variables
fstream	ReportFile

Function Documentation

HstopCurve()

int HstopCurve	(	string cmd	UNUSED,
		int	study_num,
		Mdouble	h,
		Mdouble	hMax,
		Mdouble	a
	)

{
    cout << "restart at study_num=" << study_num << ", h=" << h << ", a=" << a << endl;
    
    stringstream name;
    name << "Report" << study_num;
    ReportFile.open(name.str().c_str(),fstream::out | fstream::app);
    //~ ReportFile << "height\tangle\tlambda\tmu\tabovemuBottom\tabove\n";
    bool piac = PointIsAboveCurve(h,a,study_num);
    ReportFile.close();
 
    //now increase height gradually; decrease angle if flow starts
    if (h<hMax) {
        if (piac) a-=.5;
        else h+=2;
        stringstream command;
        // for msm1:
        command << "echo " << study_num << " " << h << " " << hMax << " " << a << " > arg";
        // for einder:
        //command << "name="<<cmd<<" && cd ${name/hstop_StudyHeightHmaxAngle.exe/} && ~/clusterscriptexecute $name " << study_num << " " << h << " " << hMax << " " << a << "&";
        cout << command.str() << endl;
        return system (command.str().c_str());
    } else {
        return 0;   
    }
}

References units::name, PointIsAboveCurve(), and ReportFile.

Referenced by main().

main()

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

{
    HstopCurve(argv[0],4,4.,60.,24.);   
}

References HstopCurve().

PointIsAboveCurve()

bool PointIsAboveCurve	(	Mdouble	h,
		Mdouble	a,
		int	study_num
	)

                                                            {
    cout << "PointIsAboveCurve(" << study_num << ", h=" << h << ", a=" << a << ")" << endl;
    SilbertHstop problem;
    problem.setInflowHeight(h) ;
    problem.setChuteAngle(a);
    problem.set_study(study_num);
    stringstream name;
    name << "H" << problem.getInflowHeight() 
         << "A" << problem.getChuteAngleDegrees() 
         << "L" << round(100.*problem.getFixedParticleRadius()*2.)/100.
         << "M" << problem.getSlidingFrictionCoefficient()
         << "B" << problem.getSlidingFrictionCoefficientBottom();
    problem.setName(name.str().c_str());
    cout << "starting " << name.str().c_str() << endl;
    // if you want more restart output, use these options
    //problem.restartFile.setFileType(FileType::MULTIPLE_FILES_PADDED);
    // for real runs:
    //problem.solve();
    // for test runs:
    problem.solve_analytic();
    stringstream com;
    ReportFile   << problem.getInflowHeight() 
         << "\t" << problem.getChuteAngleDegrees() 
         << "\t" << round(100.*problem.getFixedParticleRadius()*2.)/100.
         << "\t" << problem.getSlidingFrictionCoefficient()
         << "\t" << problem.getSlidingFrictionCoefficientBottom()
         << "\t" << problem.pointIsAboveCurve
         << endl;
    return problem.pointIsAboveCurve;
}

References Chute::getChuteAngleDegrees(), Chute::getFixedParticleRadius(), Chute::getInflowHeight(), SilbertPeriodic::getSlidingFrictionCoefficientBottom(), units::name, SilbertHstop::pointIsAboveCurve, ReportFile, helpers::round(), SilbertPeriodic::set_study(), Chute::setChuteAngle(), Chute::setInflowHeight(), DPMBase::setName(), and SilbertHstop::solve_analytic().

Referenced by HstopCurve().

Variable Documentation

ReportFile

fstream ReportFile

Referenced by HstopCurve(), and PointIsAboveCurve().