Homec4science
Diffusion LibMultiScale (features/new_parralelism) Browse

libmultiscale/src/continuum/domain_continuum.hh
b05d5f2aba78features/new_parralelism

Last Change
Raw File

domain_continuum.hh
View Options

/**
* @file domain_continuum.hh
*
* @author Guillaume Anciaux <guillaume.anciaux@epfl.ch>
* @author Till Junge <till.junge@epfl.ch>
*
* @date Tue Oct 29 22:09:05 2013
*
* @brief Mother of all continuum domains
*
* @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_DOMAIN_CONTINUUM_HH__
#define __LIBMULTISCALE_DOMAIN_CONTINUUM_HH__
/* -------------------------------------------------------------------------- */
#include "domain_continuum_interface.hh"
#include "container_mesh.hh"
/* -------------------------------------------------------------------------- */
__BEGIN_LIBMULTISCALE__
/* -------------------------------------------------------------------------- */
template <typename Ref> class Filter;
/* -------------------------------------------------------------------------- */
/**
Class DomainContinu
Template class from DomainContinuumInterface
that helps developper to implement with optimisation
processes by template a Continuous domain object
*/
template <typename ContElems, typename ContNodes, UInt Dim>
class DomainContinuum :
public DomainContinuumInterface,
public ComponentLMOutTyped<ContainerMesh<ContNodes,ContElems>,OUTPUT_REF>
{
/* ------------------------------------------------------------------------ */
/* Constructors/Destructors */
/* ------------------------------------------------------------------------ */
public:
DomainContinuum(DomainID ID,CommGroup GID);
virtual ~DomainContinuum();
/* ------------------------------------------------------------------------ */
/* Typedefs */
/* ------------------------------------------------------------------------ */
//! generic container type to course over nodes
typedef ContNodes ContainerNodes;
//! generic iterator type to course over nodes
typedef typename ContainerNodes::iterator IteratorNodes;
//! generic reference to nodes
typedef typename ContainerNodes::Ref RefNode;
//! generic reference to nodes
typedef typename ContainerNodes::Ref RefPoint;
//! generic reference to container of points (contains also the mesh)
typedef ContainerMesh<ContNodes,ContElems> ContainerPoints;
//! generic container type to course over elems
typedef ContElems ContainerElems;
//! generic iterator type to course over elems
typedef typename ContainerElems::iterator IteratorElems;
//! generic reference to nodes
typedef typename ContainerElems::Ref RefElem;
//! generic reference to container of points
typedef typename ContainerPoints::ContainerSubset ContainerSubset;
/* ------------------------------------------------------------------------ */
/* Accessors */
/* ------------------------------------------------------------------------ */
//! return the geometric definition of the domain
Geometry & getGeom();
//! get internal timestep to a given value
virtual Real getTimeStep(){return timeStep;};
//! return the container of elements
ContainerElems & getContainerElems();
//! return the container of nodes
ContainerNodes & getContainerNodes();
//! return the container of points
ContainerPoints & getContainer();
ContainerPoints & getOutput();
//! return PBC pairs in a vector
std::vector< std::pair<UInt ,UInt > > & getPBCpairs();
//! return value of arlequin flag
bool arlequin(){return arlequin_flag;}
//! returns whether an md domain is perdiodic in direction dir
virtual bool isPeriodic(UInt dir){LM_TOIMPLEMENT;}
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
void build(){};
//! implementation of the parsing method
virtual void declareParams();
//!reload xml file of a saved situation
void readXMLFile(const std::string & filename);
//! connect and build the adapted object
void connect(FilterManager & f);
//! connect and build the adapted object
void connect(ActionManager & f);
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
protected:
//! flag to active arlequin energy weigthing
bool arlequin_flag;
//! PBC pairs
std::vector< std::pair<UInt ,UInt > > pbc_pairs;
//! slave node initial displacements;
std::vector<Real> slave_displacements;
//! container for elements
ContElems elems;
//! container for nodes
ContNodes nodes;
//! container of mesh
ContainerPoints mesh_container;
//! whether the initial displacements have been computed
bool slave_displacements_computed;
//! mapping from zone id to geometries
std::vector<GeomID> geometries;
//! mapping between zone id and BC types
// std::vector<BoundaryType *> boundaries;
//! contrained geom
GeomID geom_constrained;
//! timestep value
Quantity<Time> timeStep;
//! flag to activate correct handling of prestressed meshes in pbc
bool prestressed_periodicity_flag;
};
/* -------------------------------------------------------------------------- */
__END_LIBMULTISCALE__
#endif /* __LIBMULTISCALE_DOMAIN_CONTINUUM_HH__ */
c4science · Help