OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Added: Simulation driver for flexible structure in FSI-simulation (forces not coupled yet)

Browse Source 
git-svn-id: http://svn.sintef.no/trondheim/IFEM/trunk@2369 e10b68d5-8a6e-419e-a041-bce267b0401d
This commit is contained in:
kmo
2013-05-07 17:02:06 +00:00
committed by Knut Morten Okstad
parent 5a07d21e08
commit 57f4329df0
1 changed files with 69 additions and 60 deletions
Download Patch File Download Diff File Expand all files Collapse all files

129
Apps/Common/AppCommon.h
View File

@@ -1,3 +1,4 @@
// $Id$
//==============================================================================
//!
//! \file AppCommon.h
@@ -6,76 +7,84 @@
//!
//! \author Arne Morten Kvarving / SINTEF
//!
//! \brief Common helper templates for applications
//! \brief Common helper templates for applications.
//!
//==============================================================================
#ifndef APP_COMMON_H_
#define APP_COMMON_H_
#ifndef _APP_COMMON_H_
#define _APP_COMMON_H_
#include "HDF5Writer.h"
#include "XMLWriter.h"
namespace SIM {
//! \brief Handle application restarts
//! \param[template] Simulator The top SIM of your application
//! \param[template] Solver The SIMSolver of your application
//! \param[in] Simulator The top SIM instance of your application
//! \param[in] Simulator The SIMSolver instance of your application
//! \param[in] interval The stride in the input file
//! \param[in] steps The number of timesteps to load
//! \param[in] restartfile The file to read from
template<class Simulator, class Solver>
void handleRestart(Simulator& simulator, Solver& solver,
int interval, int steps, const char* restartfile)
{
DataExporter reader(true, interval, steps);
simulator.registerFields(reader);
XMLWriter* xml = new XMLWriter(restartfile);
HDF5Writer* hdf = new HDF5Writer(restartfile,true);
reader.registerWriter(xml);
reader.registerWriter(hdf);
int max = reader.getTimeLevel();
if (steps == 2 && reader.getTimeLevel() > -1) {
reader.loadTimeLevel(reader.getTimeLevel()-1,xml,hdf);
solver.advanceStep();
}
if (max > -1 && reader.loadTimeLevel(-1,xml,hdf)) {
xml->writeTimeInfo(0, interval, steps, solver.getTimePrm());
solver.fastForward(reader.realTimeLevel(max));
}
}
//! \brief Handle application data output
//! \param[template] Simulator The top SIM of your application
//! \param[template] Solver The SIMSolver of your application
//! \param[in] Simulator The top SIM instance of your application
//! \param[in] Simulator The SIMSolver instance of your application
//! \param[in] append Whether or not to append to file
//! \param[in] hdf5file The file to save to
//! \param[in] interval The stride in the input file
//! \param[in] steps The number of timesteps to load
template<class Simulator, class Solver>
DataExporter* handleDataOutput(Simulator& simulator, Solver& solver,
bool append,
const std::string& hdf5file,
int interval,
int steps)
namespace SIM
{
DataExporter* writer = new DataExporter(true, interval, steps);
simulator.registerFields(*writer);
XMLWriter* xml = new XMLWriter(hdf5file);
HDF5Writer* hdf = new HDF5Writer(hdf5file, append);
writer->registerWriter(xml);
writer->registerWriter(hdf);
if (!append)
writer->dumpTimeLevel(&solver.getTimePrm()); // initial state
//! \brief Handles application restarts.
//! \param[in] simulator The top SIMbase instance of your application
//! \param[in] solver The SIMSolver instance of your application
//! \param[in] interval The stride in the input file
//! \param[in] steps The number of timesteps to load
//! \param[in] restartfile The file to read from
template<class Simulator, class Solver>
void handleRestart(Simulator& simulator, Solver& solver,
const std::string& restartfile, int interval, int steps)
{
DataExporter reader(true,interval,steps);
simulator.registerFields(reader);
XMLWriter* xml = new XMLWriter(restartfile);
HDF5Writer* hdf = new HDF5Writer(restartfile,true);
reader.registerWriter(xml);
reader.registerWriter(hdf);
int max = reader.getTimeLevel();
if (steps == 2 && reader.getTimeLevel() > -1) {
reader.loadTimeLevel(reader.getTimeLevel()-1,xml,hdf);
solver.advanceStep();
}
if (max > -1 && reader.loadTimeLevel(-1,xml,hdf)) {
xml->writeTimeInfo(0, interval, steps, solver.getTimePrm());
solver.fastForward(reader.realTimeLevel(max));
}
}
return writer;
//! \brief Handles application restarts.
//TODO: Remove this later. Obsolete wrapper.
template<class Simulator, class Solver>
void handleRestart(Simulator& simulator, Solver& solver,
int interval, int steps, const char* restartfile)
{ return handleRestart(simulator,solver,restartfile,interval,steps); }
//! \brief Handles application data output.
//! \param[in] simulator The top SIMbase instance of your application
//! \param[in] solver The SIMSolver instance of your application
//! \param[in] hdf5file The file to save to
//! \param[in] append Whether or not to append to file
//! \param[in] interval The stride in the input file
//! \param[in] steps The number of timesteps to load
template<class Simulator, class Solver>
DataExporter* handleDataOutput(Simulator& simulator, Solver& solver,
const std::string& hdf5file,
bool append, int interval, int steps)
{
DataExporter* writer = new DataExporter(true,interval,steps);
simulator.registerFields(*writer);
XMLWriter* xml = new XMLWriter(hdf5file);
HDF5Writer* hdf = new HDF5Writer(hdf5file,append);
writer->registerWriter(xml);
writer->registerWriter(hdf);
if (!append)
writer->dumpTimeLevel(&solver.getTimePrm()); // initial state
return writer;
}
//! \brief Handles application data output.
//TODO: Remove this later. Obsolete wrapper.
template<class Simulator, class Solver>
DataExporter* handleDataOutput(Simulator& simulator, Solver& solver,
bool append, const std::string& hdf5file,
int interval, int steps)
{ return handleDataOutput(simulator,solver,hdf5file,append,interval,steps); }
}
}
#endif