Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F75501551
model.hh
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Fri, Aug 2, 15:04
Size
9 KB
Mime Type
text/x-c++
Expires
Sun, Aug 4, 15:04 (2 d)
Engine
blob
Format
Raw Data
Handle
19559425
Attached To
rAKA akantu
model.hh
View Options
/**
* @file model.hh
*
* @author Guillaume Anciaux <guillaume.anciaux@epfl.ch>
* @author David Simon Kammer <david.kammer@epfl.ch>
* @author Nicolas Richart <nicolas.richart@epfl.ch>
*
* @date Tue Jul 27 18:15:37 2010
*
* @brief Interface of a model
*
* @section LICENSE
*
* Copyright (©) 2010-2011 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* Akantu is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Akantu is distributed in the hope that it will be useful, but WITHOUT ANY
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
* A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Akantu. If not, see <http://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#ifndef __AKANTU_MODEL_HH__
#define __AKANTU_MODEL_HH__
/* -------------------------------------------------------------------------- */
#include "aka_common.hh"
#include "aka_memory.hh"
#include "mesh.hh"
#include "fem.hh"
#include "mesh_utils.hh"
#include "synchronizer_registry.hh"
#include "distributed_synchronizer.hh"
#include "static_communicator.hh"
#include "mesh_partition.hh"
#include "dof_synchronizer.hh"
#include "pbc_synchronizer.hh"
/* -------------------------------------------------------------------------- */
__BEGIN_AKANTU__
class DumperIOHelper;
class Model : public Memory {
/* ------------------------------------------------------------------------ */
/* Constructors/Destructors */
/* ------------------------------------------------------------------------ */
public:
typedef Mesh mesh_type;
Model(Mesh& mesh, UInt spatial_dimension = _all_dimensions,
const ID & id = "model",
const MemoryID & memory_id = 0);
virtual ~Model();
typedef std::map<std::string, FEM *> FEMMap;
/* ------------------------------------------------------------------------ */
/* Methods */
/* ------------------------------------------------------------------------ */
public:
virtual void initModel() = 0;
/// create the synchronizer registry object
void createSynchronizerRegistry(DataAccessor * data_accessor);
/// create a parallel synchronizer and distribute the mesh
Synchronizer & createParallelSynch(MeshPartition * partition,
DataAccessor * data_accessor);
/// change local equation number so that PBC is assembled properly
void changeLocalEquationNumberforPBC(std::map<UInt,UInt> & pbc_pair,UInt dimension);
/// function to print the containt of the class
virtual void printself(std::ostream & stream, int indent = 0) const = 0;
/// initialize the model for PBC
void setPBC(UInt x, UInt y, UInt z);
void setPBC(SurfacePairList & surface_pairs,
ElementType surface_e_type);
virtual void initPBC();
/* ------------------------------------------------------------------------ */
/* Access to the dumpable interface of the boundaries */
/* ------------------------------------------------------------------------ */
/// Dump the data for a given boundary
void dumpBoundary(const std::string & boundary_name);
void dumpBoundary(const std::string & boundary_name,
const std::string & dumper_name);
/// Dump the data for all boundaries
void dumpBoundary();
/// Set the directory for a given boundary
void setBoundaryDirectory(const std::string & directory,
const std::string & boundary_name);
/// Set the directory for all boundaries
void setBoundaryDirectory(const std::string & directory);
/// Set the base name for a given boundary
void setBoundaryBaseName(const std::string & basename,
const std::string & boundary_name);
/// Get the internal dumper of a given boundary
DumperIOHelper & getBoundaryDumper(const std::string & boundary_name);
/* ------------------------------------------------------------------------ */
/* Accessors */
/* ------------------------------------------------------------------------ */
public:
/// get id of model
AKANTU_GET_MACRO(ID, id, const ID)
/// get the number of surfaces
AKANTU_GET_MACRO(Mesh, mesh, Mesh&);
/// return the object hadling equation numbers
AKANTU_GET_MACRO(DOFSynchronizer, *dof_synchronizer, const DOFSynchronizer &)
/// synchronize the boundary in case of parallel run
virtual void synchronizeBoundaries() {};
/// return the fem object associated with a provided name
inline FEM & getFEM(const ID & name = "") const;
/// return the fem boundary object associated with a provided name
virtual FEM & getFEMBoundary(const ID & name = "");
/// register a fem object associated with name
template <typename FEMClass> inline void registerFEMObject(const std::string & name,
Mesh & mesh,
UInt spatial_dimension);
/// unregister a fem object associated with name
inline void unRegisterFEMObject(const std::string & name);
/// return the synchronizer registry
SynchronizerRegistry & getSynchronizerRegistry();
public:
/// return the fem object associated with a provided name
template <typename FEMClass>
inline FEMClass & getFEMClass(std::string name = "") const;
/// return the fem boundary object associated with a provided name
template <typename FEMClass>
inline FEMClass & getFEMClassBoundary(std::string name = "");
/// get the pbc pairs
std::map<UInt,UInt> & getPBCPairs(){return pbc_pair;};
/* ------------------------------------------------------------------------ */
/* Pack and unpack helper functions */
/* ------------------------------------------------------------------------ */
public:
inline UInt getNbQuadraturePoints(const Array<Element> & elements,
const ID & fem_id = ID()) const;
protected:
/// returns if node is slave in pbc
inline bool getIsPBCSlaveNode(const UInt node);
public:
template<typename T>
inline void packElementalDataHelper(const ByElementTypeArray<T> & data_to_pack,
CommunicationBuffer & buffer,
const Array<Element> & elements,
bool per_quadrature_point = true,
const ID & fem_id = ID()) const;
template<typename T>
inline void unpackElementalDataHelper(ByElementTypeArray<T> & data_to_unpack,
CommunicationBuffer & buffer,
const Array<Element> & elements,
bool per_quadrature_point = true,
const ID & fem_id = ID()) const;
protected:
template<typename T, bool pack_helper>
inline void packUnpackElementalDataHelper(ByElementTypeArray<T> & data_to_pack,
CommunicationBuffer & buffer,
const Array<Element> & element,
bool per_quadrature_point,
const ID & fem_id) const;
/* -------------------------------------------------------------------------- */
template<typename T>
inline void packNodalDataHelper(Array<T> & data_to_pack,
CommunicationBuffer & buffer,
const Array<Element> & element) const;
template<typename T>
inline void unpackNodalDataHelper(Array<T> & data_to_unpack,
CommunicationBuffer & buffer,
const Array<Element> & element) const;
template<typename T, bool pack_helper>
inline void packUnpackNodalDataHelper(Array<T> & data,
CommunicationBuffer & buffer,
const Array<Element> & elements) const;
/* ------------------------------------------------------------------------ */
/* Class Members */
/* ------------------------------------------------------------------------ */
protected:
/// Mesh
Mesh & mesh;
/// id
ID id;
/// Spatial dimension of the problem
UInt spatial_dimension;
/// the main fem object present in all models
FEMMap fems;
/// the fem object present in all models for boundaries
FEMMap fems_boundary;
/// default fem object
std::string default_fem;
/// synchronizer registry
SynchronizerRegistry * synch_registry;
/// handle the equation number things
DOFSynchronizer * dof_synchronizer;
/// pbc pairs
std::map<UInt,UInt> pbc_pair;
/// flag per node to know is pbc slave
Array<bool> is_pbc_slave_node;
};
/* -------------------------------------------------------------------------- */
/* inline functions */
/* -------------------------------------------------------------------------- */
#if defined (AKANTU_INCLUDE_INLINE_IMPL)
# include "model_inline_impl.cc"
#endif
/// standard output stream operator
inline std::ostream & operator <<(std::ostream & stream, const Model & _this)
{
_this.printself(stream);
return stream;
}
__END_AKANTU__
#endif /* __AKANTU_MODEL_HH__ */
Event Timeline
Log In to Comment