Homec4science
Diffusion LibMultiScale (features/coupler_component) Browse

libmultiscale/src/dumper/dumper_paraview.hh
b05d5f2aba78features/coupler_compon…

Last Change
Raw File

dumper_paraview.hh
View Options

/**
* @file dumper_paraview.hh
*
* @author Guillaume Anciaux <guillaume.anciaux@epfl.ch>
*
* @date Fri Jul 04 21:55:04 2014
*
* @brief This dumper allows to generate paraview files for vizualization
*
* @section LICENSE
*
* Copyright INRIA and CEA
*
* The LibMultiScale is a C++ parallel framework for the multiscale
* coupling methods dedicated to material simulations. This framework
* provides an API which makes it possible to program coupled simulations
* and integration of already existing codes.
*
* This Project was initiated in a collaboration between INRIA Futurs Bordeaux
* within ScAlApplix team and CEA/DPTA Ile de France.
* The project is now continued at the Ecole Polytechnique Fédérale de Lausanne
* within the LSMS/ENAC laboratory.
*
* This software is governed by the CeCILL-C license under French law and
* abiding by the rules of distribution of free software. You can use,
* modify and/ or redistribute the software under the terms of the CeCILL-C
* license as circulated by CEA, CNRS and INRIA at the following URL
* "http://www.cecill.info".
*
* As a counterpart to the access to the source code and rights to copy,
* modify and redistribute granted by the license, users are provided only
* with a limited warranty and the software's author, the holder of the
* economic rights, and the successive licensors have only limited
* liability.
*
* In this respect, the user's attention is drawn to the risks associated
* with loading, using, modifying and/or developing or reproducing the
* software by the user in light of its specific status of free software,
* that may mean that it is complicated to manipulate, and that also
* therefore means that it is reserved for developers and experienced
* professionals having in-depth computer knowledge. Users are therefore
* encouraged to load and test the software's suitability as regards their
* requirements in conditions enabling the security of their systems and/or
* data to be ensured and, more generally, to use and operate it in the
* same conditions as regards security.
*
* The fact that you are presently reading this means that you have had
* knowledge of the CeCILL-C license and that you accept its terms.
*
*/
#ifndef __LIBMULTISCALE_DUMPER_PARAVIEW_HH__
#define __LIBMULTISCALE_DUMPER_PARAVIEW_HH__
/* -------------------------------------------------------------------------- */
#include "dumper.hh"
#include "paraview_helper.hh"
/* -------------------------------------------------------------------------- */
__BEGIN_LIBMULTISCALE__
/** Class DumperParaview
* Implementation of a dumper to paraview vtu files
*/
template <typename _Input>
class DumperParaview : public Dumper<_Input>{
/* ------------------------------------------------------------------------ */
/* Constructors/Destructors */
/* ------------------------------------------------------------------------ */
public:
DumperParaview(const std::string & name, ComponentLMInterface & d);
~DumperParaview();
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
//! most generic dump function
void dump(_Input & cont);
template <typename ContField>
void dumpField(const std::string & name_field, _Input & cont,
ContField & contField, const UInt Dim);
template <typename ContField, UInt Dim>
void dumpField(const std::string & name_field,_Input & cont,ContField & contField);
//! initialisation of the dumper (do nothing)
void init();
//! set the parameters to dumper
/** such as binary, text mode or
serie that need to be dumped */
void declareParams();
void openPVTU();
void openVTU();
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
private:
//! flag for displacement serie
bool dep_serie;
//! flag for initial position serie
bool p0_serie;
//! flag for velocity serie
bool vel_serie;
//! flag for force/acceleration serie
bool force_serie;
//! flag for applied force serie
bool applied_force_serie;
//! flag for distribution serie
bool proc_serie;
//! flag for constraUInt serie
bool stress_serie;
//! flag for deformation serie
bool strain_serie;
//! flag for potential energy serie
bool epot_serie;
//! flag for mass serie
bool mass_serie;
//! flag for charge serie
bool charge_serie;
//! flag for id serie
bool id_serie;
//! flag to produce zipped files
bool flag_compressed;
//! flag to produce text files
bool flag_text;
//! flag to produce base64 files
bool flag_base64;
//! flag for weight serie
bool alpha_serie;
//! flag for electronic density
bool e_density_serie;
//! flag for barycenter computation
bool barycenters_serie;
//! flag for temperature serie
bool temperature_serie;
//! flag for temperature variation serie
bool temperatureVar_serie;
//! flag for heat capacity serie
bool heatCapacity_serie;
//! flag for heat capacity serie
bool heatRate_serie;
//! flag for external heat rate serie
bool external_heatRate_serie;
//! flag for tag serie
bool tag_serie;
//! additional field
std::vector<FilterID> additional_fields;
//! pointer to paraview helper facilities
ParaviewHelper paraHelper;
//! the file where to write
LMFile file;
//! flag for output of nodal blockage
bool flag_boundary;
};
/* -------------------------------------------------------------------------- */
/* -------------------------------------------------------------------------- */
template <typename _Input>
template <typename ContField, UInt Dim>
void DumperParaview<_Input>::dumpField(const std::string & name_field, _Input & cont,
ContField & contField){
if (contField.nbElem()/Dim != cont.nbElem())
LM_FATAL("compute/filter " << name_field << " contains " << contField.nbElem()/Dim << " entries"
<< " while container has " << cont.nbElem() << " : incompatible match");
paraHelper.startData(name_field,Dim);
UInt cpt = 0;
typename ContField::iterator it = contField.getIterator();
for (Real value = it.getFirst();!it.end();value = it.getNext(), ++cpt) {
paraHelper.pushReal(value);
}
paraHelper.endData();
}
/* -------------------------------------------------------------------------- */
template <typename _Input>
template <typename ContField>
void DumperParaview<_Input>::dumpField(const std::string & name_field, _Input & cont,
ContField & contField, const UInt Dim){
switch (Dim){
case 1 : dumpField<ContField,1>(name_field,cont,contField);break;
case 2 : dumpField<ContField,2>(name_field,cont,contField);break;
case 3 : dumpField<ContField,3>(name_field,cont,contField);break;
case 4 : dumpField<ContField,3>(name_field,cont,contField);break;
case 5 : dumpField<ContField,3>(name_field,cont,contField);break;
case 6 : dumpField<ContField,3>(name_field,cont,contField);break;
case 7 : dumpField<ContField,3>(name_field,cont,contField);break;
case 8 : dumpField<ContField,3>(name_field,cont,contField);break;
case 9 : dumpField<ContField,3>(name_field,cont,contField);break;
case 10: dumpField<ContField,3>(name_field,cont,contField);break;
default: LM_FATAL("internal error " << Dim << " " << name_field);
}
}
/* -------------------------------------------------------------------------- */
__END_LIBMULTISCALE__
#endif /* __LIBMULTISCALE_DUMPER_PARAVIEW_HH__ */
c4science · Help