OutputROMSolution.cpp File Reference

Write out the solution to file. More...

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <basic.h>
#include <common_cpp.h>
#include <arrayfunctions.h>
#include <io_cpp.h>
#include <plotfuncs_cpp.h>
#include <mpivars_cpp.h>
#include <simulation_object.h>

Go to the source code of this file.

Functions
void	IncrementFilenameIndex (char *, int)
int	OutputROMSolution (void *s, int nsims, double a_time)

Detailed Description

Write out the solution to file.

Author

Debojyoti Ghosh

Definition in file OutputROMSolution.cpp.

Function Documentation

IncrementFilenameIndex()

void IncrementFilenameIndex	(	char *	f,
		int	len
	)

Increment the output filename of the form op_nnnnn.dat by 1, i.e., the number represented by the string nnnnn is incremented by 1.

Increment a character string representing an integer by 1. For example: "00002" -> "00003"; "3421934" -> "3421935"; "999" -> "000". The string can be of arbitrary length.

Parameters

f	Character string representing the integer
len	Length of the string

Definition at line 46 of file IncrementFilename.c.

{
  int i;
  for (i=len-1; i>=0; i--) {
    if (f[i] == '9') {
      f[i] = '0';
      if (!i) fprintf(stderr,"Warning: file increment hit max limit. Resetting to zero.\n");
    } else {
      f[i]++;
      break;
    }
  }
}

OutputROMSolution()

int OutputROMSolution	(	void *	s,
		int	nsims,
		double	a_time
	)

Write out the ROM solution to file

Parameters

s	Array of simulation objects of type SimulationObject
nsims	Number of simulation objects
a_time	Current simulation time

Definition at line 24 of file OutputROMSolution.cpp.

{
  SimulationObject* simobj = (SimulationObject*) s;
  int ns;

  for (ns = 0; ns < nsims; ns++) {

    HyPar*        solver = &(simobj[ns].solver);
    MPIVariables* mpi    = &(simobj[ns].mpi);
    
    if ((!solver->WriteOutput) && (strcmp(solver->plot_solution,"yes"))) continue;
  
    char fname_root[_MAX_STRING_SIZE_];
    strcpy(fname_root, solver->op_rom_fname_root);

    if (nsims > 1) {
      char index[_MAX_STRING_SIZE_];
      GetStringFromInteger(ns, index, (int)log10(nsims)+1);
      strcat(fname_root, "_");
      strcat(fname_root, index);
    }
  
    WriteArray( solver->ndims,
                solver->nvars,
                solver->dim_global,
                solver->dim_local,
                solver->ghosts,
                solver->x,
                solver->u_rom_predicted,
                solver,
                mpi,
                fname_root );

    if (!strcmp(solver->plot_solution, "yes")) {
      PlotArray(   solver->ndims,
                   solver->nvars,
                   solver->dim_global,
                   solver->dim_local,
                   solver->ghosts,
                   solver->x,
                   solver->u,
                   a_time,
                   solver,
                   mpi,
                   fname_root );
    }

    /* increment the index string, if required */
    if ((!strcmp(solver->output_mode,"serial")) && (!strcmp(solver->op_overwrite,"no"))) {
      IncrementFilenameIndex(solver->filename_index,solver->index_length);
    }

  }
  
  return 0;
}