Commit c8cdf0c1 authored by Michael's avatar Michael
Browse files

Print descriptive names for dimtype and expansion

At startup, print a string describing the dimension
type and spectral expansion type, rather than the
internal enum value.
parent 6d7144d8
......@@ -48,6 +48,7 @@ namespace NSIntegrator {
using std::set;
using namespace Transformer;
// List supported physical dimension types (Fourier, free slip, no slip, mapped)
static const char* DIMTYPE_NAME[] = { "PERIODIC", "FREE_SLIP", "NO_SLIP" };
enum DIMTYPE {
PERIODIC,
FREE_SLIP, // Cosine expansion
......@@ -228,9 +229,9 @@ namespace NSIntegrator {
}
if (master()) {
fprintf(stderr,"Beginning Navier-Stokes timestepping, on a %d x %d x %d grid\n",szx,szy,szz);
fprintf(stderr, "X-dimension: type %d, expansion %d\n",tx,Sx);
fprintf(stderr, "Y-dimension: type %d, expansion %d\n",ty,Sy);
fprintf(stderr, "Z-dimension: type %d, expansion %d\n",tz,Sz);
fprintf(stderr, "X-dimension: type %s, expansion %s\n",DIMTYPE_NAME[tx],S_EXP_NAME[Sx]);
fprintf(stderr, "Y-dimension: type %s, expansion %s\n",DIMTYPE_NAME[ty],S_EXP_NAME[Sy]);
fprintf(stderr, "Z-dimension: type %s, expansion %s\n",DIMTYPE_NAME[tz],S_EXP_NAME[Sz]);
fprintf(stderr, "%d tracers\n",(int) tracers.size());
for (unsigned int k = 0; k < tracers.size(); k++) {
if (t_diffusivity[k] == 0) {
......
......@@ -17,6 +17,7 @@ namespace Transformer {
using blitz::TinyVector;
using blitz::GeneralArrayStorage;
static const char *S_EXP_NAME[] = { "COSINE", "SINE", "FOURIER", "CHEBY", "REAL", "COMPLEX", "NONE" };
enum S_EXP { /* Spectral expansions */
COSINE = 0, // DCT, for no-normal boundaries
SINE = 1, // DST, for derivatives of no-normal boundaries
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment