42 int main(
int argc,
char *argv[])
51 scalar fullCycle = 360.0;
56 while (runTime.userTimeValue() > ca0)
62 while (
mag(runTime.userTimeValue() - ca0) > eps)
64 scalar t0 = runTime.userTimeToTime(ca0 - runTime.userTimeValue());
65 runTime.setDeltaT(t0);
67 Info<<
"CA = " << runTime.userTimeValue() <<
endl;
71 scalar Vmax =
sum(mesh.V().field());
73 while (
mag(runTime.userTimeValue()-ca1) > eps)
75 scalar t1 = runTime.userTimeToTime(ca1 - runTime.userTimeValue());
76 runTime.setDeltaT(t1);
78 Info<<
"CA = " << runTime.userTimeValue() <<
endl;
82 scalar Vmin =
sum(mesh.V().field());
84 Info<<
"\nVmax = " << Vmax;
86 Info<<
"Compression ratio Vmax/Vmin = " << Vmax/Vmin <<
endl;
int main(int argc, char *argv[])
Ostream & endl(Ostream &os)
Add newline and flush stream.
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
dimensioned< scalar > mag(const dimensioned< Type > &)