TimePrintStep.c

Go to the documentation of this file.

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <simulation_object.h>
#include <timeintegration.h>

int TimePrintStep(void *ts )
{
  TimeIntegration* TS = (TimeIntegration*) ts;
  SimulationObject* sim = (SimulationObject*) TS->simulation;
  int ns, nsims = TS->nsims;

  if ((!TS->rank) && ((TS->iter+1)%sim[0].solver.screen_op_iter == 0)) {
    if (nsims > 1) {
      printf("--\n");
      printf("iter=%7d,  t=%1.3e\n", TS->iter+1, TS->waqt);
      if (TS->max_cfl >= 0) printf("  CFL=%1.3E\n", TS->max_cfl);
      if (TS->norm >= 0) printf("  norm=%1.4E\n", TS->norm);
      printf("  wctime=%1.1E (s)\n",TS->iter_wctime);
    } else {
      printf("iter=%7d  ",TS->iter+1 );
      printf("t=%1.3E  ",TS->waqt );
      if (TS->max_cfl >= 0) printf("CFL=%1.3E  ",TS->max_cfl );
      if (TS->norm >= 0) printf("norm=%1.4E  ",TS->norm );
      printf("wctime: %1.1E (s)  ",TS->iter_wctime);
    }

    /* calculate and print conservation error */
    if (!strcmp(sim[0].solver.ConservationCheck,"yes")) {

      double error = 0;
      for  (ns = 0; ns < nsims; ns++) {
        int v;
        for (v=0; v<sim[ns].solver.nvars; v++) {
          error += (sim[ns].solver.ConservationError[v] * sim[ns].solver.ConservationError[v]);
        }
      }
      error = sqrt(error);
      printf("  cons_err=%1.4E\n", error);

    } else {

      if (nsims == 1) printf("\n");

    }

    /* print physics-specific info, if available */
    for (ns = 0; ns < nsims; ns++) {
      if (sim[ns].solver.PrintStep) {
        if (nsims > 1) printf("Physics-specific output for domain %d:\n", ns);
        sim[ns].solver.PrintStep( &(sim[ns].solver),
                                  &(sim[ns].mpi),
                                  TS->waqt );
      }
    }
    if (nsims > 1) {
      printf("--\n");
      printf("\n");
    }
  }

  return(0);
}