LeesEdwardsSelfTest.cpp File Reference

#include "Mercury2D.h"
#include "Boundaries/LeesEdwardsBoundary.h"
#include "Species/LinearViscoelasticSlidingFrictionSpecies.h"

Classes
class	LeesEdwardsSelfTest
	[Lees:headers] More...

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

Function Documentation

main()

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

[Lees:class] [Lees:main]

{
    //instantiate the class
    LeesEdwardsSelfTest problem;
    //set output and time stepping properties
    problem.setXBallsAdditionalArguments("-w0 -v0 -solidf -cmode 5");
    //solve
    problem.solve(argc, argv);
 
    logger(INFO,"Execute 'source LeesEdwardsSelfTest.sh' to get coarse-grained statistics");
    helpers::writeToFile("LeesEdwardsSelfTest.sh","../MercuryCG/fstatistics LeesEdwardsSelfTest -stattype XY -w 0.25 -h 0.25 -tmin 2\n"
     "../MercuryCG/fstatistics LeesEdwardsSelfTest -stattype Y -w 0.25 -h 0.1 -tmin 2 -o LeesEdwardsSelfTest.Y.stat");
 
    logger(INFO,"Run 'LeesEdwardsSelfTest.m' in MATLAB/octave to visualise the statistical output");
    helpers::writeToFile("LeesEdwardsSelfTest.m","%% 2D velocity field v_x(x,y)\n"
     "addpath('../MercuryCG/')\n"
     "data = loadstatistics('LeesEdwardsSelfTest.stat');\n"
     "colormap(jet)\n"
     "contourf(data.x,data.y,data.VelocityX,20,'EdgeColor','none')\n"
     "c = colorbar\n"
     "c.Label.String = '\\rho';\n"
     "title('Velocity')\n"
     "xlabel('x')\n"
     "ylabel('z');\n"
     "axis equal\n"
     "%% 1D velocity field v_x(y)\n"
     "dataY = loadstatistics('LeesEdwardsSelfTest.Y.stat');\n"
     "plot(dataY.y,dataY.VelocityX)\n"
     "xlabel('y')\n"
     "ylabel('v_x');\n"
     "axis equal");
 
    return 0;
}

References INFO, logger, problem, and helpers::writeToFile().