146 FILE *counter_file=fopen(
"record",
"w");
148 for (
Mdouble a=0; a<71; a+=8) {
150 sprintf (buffer,
"rail%fR", a); md.setName(buffer);
151 md.species->setTorsionFrictionCoefficient(0);
152 md.solve(argc, argv);
153 fprintf(counter_file,
"%f %.10e %.10e\n",a,md.wallHandler.getObject(0)->getForce().X,md.wallHandler.getObject(0)->getForce().Y);
158 sprintf(buffer,
"rail%f", a);
161 md.solve(argc, argv);
162 fprintf(counter_file,
"%f %.10e %.10e\n", a, md.wallHandler.getObject(0)->getForce().X,
163 md.wallHandler.getObject(0)->getForce().Y);
165 fclose(counter_file);
double Mdouble
Definition: GeneralDefine.h:34
LL< Log::INFO > INFO
Info log level.
Definition: Logger.cc:55
Logger< MERCURYDPM_LOGLEVEL > logger("MercuryKernel")
Definition of different loggers with certain modules. A user can define its own custom logger here.
Definition: SiegenRail.cpp:31