Diffusion akantu (bugfixes/contact_mechanics_model) Browse

contact_mechan…

aka_common.hh
View Options

			/**
			* @file aka_common.hh
			*
			* @author Nicolas Richart <nicolas.richart@epfl.ch>
			*
			* @date creation: Mon Jun 14 2010
			* @date last modification: Thu Jan 21 2016
			*
			* @brief common type descriptions for akantu
			*
			* @section LICENSE
			*
			* Copyright (©) 2010-2012, 2014, 2015 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/>.
			*
			* @section DESCRIPTION
			*
			* All common things to be included in the projects files
			*
			*/

			/* -------------------------------------------------------------------------- */

			#ifndef __AKANTU_COMMON_HH__
			#define __AKANTU_COMMON_HH__

			/* -------------------------------------------------------------------------- */
			#include <list>
			#include <limits>

			/* -------------------------------------------------------------------------- */
			#define __BEGIN_AKANTU__ namespace akantu {
			#define __END_AKANTU__ }
			/* -------------------------------------------------------------------------- */
			#define __BEGIN_AKANTU_DUMPER__ namespace dumper {
			#define __END_AKANTU_DUMPER__ }
			/* -------------------------------------------------------------------------- */
			#if defined(WIN32)
			# define __attribute__(x)
			#endif

			/* -------------------------------------------------------------------------- */
			#include "aka_config.hh"
			#include "aka_error.hh"
			#include "aka_safe_enum.hh"
			/* -------------------------------------------------------------------------- */

			__BEGIN_AKANTU__

			/* -------------------------------------------------------------------------- */
			/* Common types */
			/* -------------------------------------------------------------------------- */
			typedef std::string ID;

			static const Real UINT_INIT_VALUE = Real(0.);
			#ifdef AKANTU_NDEBUG
			static const Real REAL_INIT_VALUE = Real(0.);
			#else
			static const Real REAL_INIT_VALUE = std::numeric_limits<Real>::quiet_NaN();
			#endif

			/* -------------------------------------------------------------------------- */
			/* Memory types */
			/* -------------------------------------------------------------------------- */

			typedef UInt MemoryID;


			typedef std::string Surface;
			typedef std::pair<Surface, Surface> SurfacePair;
			typedef std::list< SurfacePair > SurfacePairList;

			/* -------------------------------------------------------------------------- */
			extern const UInt _all_dimensions;


			/* -------------------------------------------------------------------------- */
			/* Mesh/FEM/Model types */
			/* -------------------------------------------------------------------------- */
			__END_AKANTU__

			#include "aka_element_classes_info.hh"

			__BEGIN_AKANTU__


			/// small help to use names for directions
			enum SpacialDirection {
			_x = 0,
			_y = 1,
			_z = 2
			};

			/// enum MeshIOType type of mesh reader/writer
			enum MeshIOType {
			_miot_auto, ///< Auto guess of the reader to use based on the extension
			_miot_gmsh, ///< Gmsh files
			_miot_gmsh_struct, ///< Gsmh reader with reintpretation of elements has structures elements
			_miot_diana, ///< TNO Diana mesh format
			_miot_abaqus ///< Abaqus mesh format
			};

			/// enum AnalysisMethod type of solving method used to solve the equation of motion
			enum AnalysisMethod {
			_static,
			_implicit_dynamic,
			_explicit_lumped_mass,
			_explicit_lumped_capacity,
			_explicit_consistent_mass
			};

			//! enum ContactResolutionMethod types of solving for the contact
			enum ContactResolutionMethod {
			_penalty,
			_lagrangian,
			_augmented_lagrangian,
			_nitsche,
			_mortar
			};

			//! enum ContactImplementationMethod types for different contact implementations
			enum ContactImplementationMethod {
			_none,
			_uzawa,
			_generalized_newton
			};

			/// enum SolveConvergenceMethod different resolution algorithms
			enum SolveConvergenceMethod {
			_scm_newton_raphson_tangent, ///< Newton-Raphson with tangent matrix
			_scm_newton_raphson_tangent_modified, ///< Newton-Raphson with constant tangent matrix
			_scm_newton_raphson_tangent_not_computed ///< Newton-Raphson with constant tangent matrix (user has to assemble it)
			};

			/// enum SolveConvergenceCriteria different convergence criteria
			enum SolveConvergenceCriteria {
			_scc_residual, ///< Use residual to test the convergence
			_scc_increment, ///< Use increment to test the convergence
			_scc_residual_mass_wgh ///< Use residual weighted by inv. nodal mass to testb
			};

			/// enum CohesiveMethod type of insertion of cohesive elements
			enum CohesiveMethod {
			_intrinsic,
			_extrinsic
			};

			/// myfunction(double * position, double * stress/force, double * normal, unsigned int material_id)
			typedef void (BoundaryFunction)(double ,double , double, unsigned int);

			/// @enum BoundaryFunctionType type of function passed for boundary conditions
			enum BoundaryFunctionType {
			_bft_stress,
			_bft_traction,
			_bft_traction_local
			};

			/// @enum SparseMatrixType type of sparse matrix used
			enum SparseMatrixType {
			_unsymmetric,
			_symmetric
			};

			/* -------------------------------------------------------------------------- */
			/* Contact */
			/* -------------------------------------------------------------------------- */
			typedef ID ContactID;
			typedef ID ContactSearchID;
			typedef ID ContactNeighborStructureID;

			enum ContactType {
			_ct_not_defined = 0,
			_ct_2d_expli = 1,
			_ct_3d_expli = 2,
			_ct_rigid = 3
			};

			enum ContactSearchType {
			_cst_not_defined = 0,
			_cst_2d_expli = 1,
			_cst_expli = 2
			};

			enum ContactNeighborStructureType {
			_cnst_not_defined = 0,
			_cnst_regular_grid = 1,
			_cnst_2d_grid = 2
			};

			/* -------------------------------------------------------------------------- */
			/* Friction */
			/* -------------------------------------------------------------------------- */
			typedef ID FrictionID;

			/* -------------------------------------------------------------------------- */
			/* Ghosts handling */
			/* -------------------------------------------------------------------------- */

			typedef ID SynchronizerID;

			/// @enum CommunicatorType type of communication method to use
			enum CommunicatorType {
			_communicator_mpi,
			_communicator_dummy
			};

			/// @enum SynchronizationTag type of synchronizations
			enum SynchronizationTag {
			//--- SolidMechanicsModel tags ---
			_gst_smm_mass, //< synchronization of the SolidMechanicsModel.mass
			_gst_smm_for_gradu, //< synchronization of the SolidMechanicsModel.displacement
			_gst_smm_boundary, //< synchronization of the boundary, forces, velocities and displacement
			_gst_smm_uv, //< synchronization of the nodal velocities and displacement
			_gst_smm_res, //< synchronization of the nodal residual
			_gst_smm_init_mat, //< synchronization of the data to initialize materials
			_gst_smm_stress, //< synchronization of the stresses to compute the internal forces
			_gst_smmc_facets, //< synchronization of facet data to setup facet synch
			_gst_smmc_facets_conn, //< synchronization of facet global connectivity
			_gst_smmc_facets_stress, //< synchronization of facets' stress to setup facet synch
			_gst_smmc_damage, //< synchronization of damage
			//--- GlobalIdsUpdater tags ---
			_gst_giu_global_conn, //< synchronization of global connectivities
			//--- CohesiveElementInserter tags ---
			_gst_ce_groups, //< synchronization of cohesive element insertion depending on facet groups
			//--- GroupManager tags ---
			_gst_gm_clusters, //< synchronization of clusters
			//--- HeatTransfer tags ---
			_gst_htm_capacity, //< synchronization of the nodal heat capacity
			_gst_htm_temperature, //< synchronization of the nodal temperature
			_gst_htm_gradient_temperature, //< synchronization of the element gradient temperature
			//--- LevelSet tags ---
			/// synchronization of the nodal level set value phi
			_gst_htm_phi,
			/// synchronization of the element gradient phi
			_gst_htm_gradient_phi,
			//--- Material non local ---
			_gst_mnl_for_average, //< synchronization of data to average in non local material
			_gst_mnl_weight, //< synchronization of data for the weight computations
			//--- NeighborhoodSynchronization tags ---
			_gst_nh_criterion,
			//--- General tags ---
			_gst_test, //< Test tag
			_gst_user_1, //< tag for user simulations
			_gst_user_2, //< tag for user simulations
			_gst_material_id, //< synchronization of the material ids
			_gst_for_dump, //< everything that needs to be synch before dump
			//--- Contact & Friction ---
			_gst_cf_nodal, //< synchronization of disp, velo, and current position
			_gst_cf_incr, //< synchronization of increment
			///--- Solver tags ---
			_gst_solver_solution //< synchronization of the solution obained with the PETSc solver
			};

			/// standard output stream operator for SynchronizationTag
			inline std::ostream & operator <<(std::ostream & stream, SynchronizationTag type);

			/// @enum GhostType type of ghost
			enum GhostType {
			_not_ghost,
			_ghost,
			_casper // not used but a real cute ghost
			};

			/* -------------------------------------------------------------------------- */
			struct GhostType_def {
			typedef GhostType type;
			static const type _begin_ = _not_ghost;
			static const type _end_ = _casper;
			};

			typedef safe_enum<GhostType_def> ghost_type_t;

			/// standard output stream operator for GhostType
			inline std::ostream & operator <<(std::ostream & stream, GhostType type);


			/// @enum SynchronizerOperation reduce operation that the synchronizer can perform
			enum SynchronizerOperation {
			_so_sum,
			_so_min,
			_so_max,
			_so_prod,
			_so_land,
			_so_band,
			_so_lor,
			_so_bor,
			_so_lxor,
			_so_bxor,
			_so_min_loc,
			_so_max_loc,
			_so_null
			};


			/* -------------------------------------------------------------------------- */
			/* Global defines */
			/* -------------------------------------------------------------------------- */
			#define AKANTU_MIN_ALLOCATION 2000

			#define AKANTU_INDENT " "
			#define AKANTU_INCLUDE_INLINE_IMPL

			/* -------------------------------------------------------------------------- */
			// int 2 type construct
			template <int d>
			struct Int2Type {
			enum { result = d };
			};

			// type 2 type construct
			template <class T>
			class Type2Type {
			typedef T OriginalType;
			};

			/* -------------------------------------------------------------------------- */
			template<class T>
			struct is_scalar {
			enum{ value = false };
			};

			#define AKANTU_SPECIFY_IS_SCALAR(type) \
			template<> \
			struct is_scalar<type> { \
			enum { value = true }; \
			}


			AKANTU_SPECIFY_IS_SCALAR(Real);
			AKANTU_SPECIFY_IS_SCALAR(UInt);
			AKANTU_SPECIFY_IS_SCALAR(Int);
			AKANTU_SPECIFY_IS_SCALAR(bool);

			template < typename T1, typename T2 >
			struct is_same {
			enum { value = false }; // is_same represents a bool.
			};

			template < typename T >
			struct is_same<T, T> {
			enum { value = true };
			};

			/* -------------------------------------------------------------------------- */
			#define AKANTU_SET_MACRO(name, variable, type) \
			inline void set##name (type variable) { \
			this->variable = variable; \
			}

			#define AKANTU_GET_MACRO(name, variable, type) \
			inline type get##name () const { \
			return variable; \
			}

			#define AKANTU_GET_MACRO_NOT_CONST(name, variable, type) \
			inline type get##name () { \
			return variable; \
			}

			#define AKANTU_GET_MACRO_BY_SUPPORT_TYPE(name, variable, type, \
			support, con) \
			inline con Array<type> & \
			get##name (const support & el_type, \
			const GhostType & ghost_type = _not_ghost) con { \
			return variable(el_type, ghost_type); \
			}

			#define AKANTU_GET_MACRO_BY_ELEMENT_TYPE(name, variable, type) \
			AKANTU_GET_MACRO_BY_SUPPORT_TYPE(name, variable, type, ElementType,)
			#define AKANTU_GET_MACRO_BY_ELEMENT_TYPE_CONST(name, variable, type) \
			AKANTU_GET_MACRO_BY_SUPPORT_TYPE(name, variable, type, ElementType, const)

			#define AKANTU_GET_MACRO_BY_GEOMETRIE_TYPE(name, variable, type) \
			AKANTU_GET_MACRO_BY_SUPPORT_TYPE(name, variable, type, GeometricalType,)
			#define AKANTU_GET_MACRO_BY_GEOMETRIE_TYPE_CONST(name, variable, type) \
			AKANTU_GET_MACRO_BY_SUPPORT_TYPE(name, variable, type, GeometricalType, const)

			/* -------------------------------------------------------------------------- */
			/// initialize the static part of akantu
			void initialize(int & argc, char ** & argv);
			/// initialize the static part of akantu and read the global input_file
			void initialize(const std::string & input_file, int & argc, char ** & argv);
			/* -------------------------------------------------------------------------- */
			/// finilize correctly akantu and clean the memory
			void finalize ();
			/* -------------------------------------------------------------------------- */
			/// Read an new input file
			void readInputFile(const std::string & input_file);
			/* -------------------------------------------------------------------------- */


			/*
			* For intel compiler annoying remark
			*/
			#if defined(__INTEL_COMPILER)
			/// remark #981: operands are evaluated in unspecified order
			#pragma warning ( disable : 981 )

			/// remark #383: value copied to temporary, reference to temporary used
			#pragma warning ( disable : 383 )

			#endif //defined(__INTEL_COMPILER)

			/* -------------------------------------------------------------------------- */
			/* string manipulation */
			/* -------------------------------------------------------------------------- */
			inline std::string to_lower(const std::string & str);
			/* -------------------------------------------------------------------------- */
			inline std::string trim(const std::string & to_trim);
			/* -------------------------------------------------------------------------- */

			/* -------------------------------------------------------------------------- */
			/// give a string representation of the a human readable size in bit
			template<typename T>
			std::string printMemorySize(UInt size);
			/* -------------------------------------------------------------------------- */


			__END_AKANTU__

			#include "aka_fwd.hh"

			__BEGIN_AKANTU__

			/// get access to the internal argument parser
			cppargparse::ArgumentParser & getStaticArgumentParser();

			/// get access to the internal input file parser
			Parser & getStaticParser();

			/// get access to the user part of the internal input file parser
			const ParserSection & getUserParser();

			__END_AKANTU__

			#include "aka_common_inline_impl.cc"

			/* -------------------------------------------------------------------------- */

			#if defined(AKANTU_UNORDERED_MAP_IS_CXX11)

			__BEGIN_AKANTU_UNORDERED_MAP__

			#if AKANTU_INTEGER_SIZE == 4
			#define AKANTU_HASH_COMBINE_MAGIC_NUMBER 0x9e3779b9
			#elif AKANTU_INTEGER_SIZE == 8
			#define AKANTU_HASH_COMBINE_MAGIC_NUMBER 0x9e3779b97f4a7c13LL
			#endif

			/**
			* Hashing function for pairs based on hash_combine from boost The magic number
			* is coming from the golden number @f[\phi = \frac{1 + \sqrt5}{2}@f]
			* @f[\frac{2^32}{\phi} = 0x9e3779b9@f]
			* http://stackoverflow.com/questions/4948780/magic-number-in-boosthash-combine
			* http://burtleburtle.net/bob/hash/doobs.html
			*/
			template <typename a, typename b> struct hash< std::pair<a, b> > {
			public:
			hash() : ah(), bh() {}
			size_t operator()(const std::pair<a, b> & p) const {
			size_t seed = ah(p.first);
			return bh(p.second) + AKANTU_HASH_COMBINE_MAGIC_NUMBER + (seed << 6) +
			(seed >> 2);
			}

			private:
			const hash<a> ah;
			const hash<b> bh;
			};

			__END_AKANTU_UNORDERED_MAP__

			#endif


			#endif /* __AKANTU_COMMON_HH__ */

rAKA/src/common/aka_common.hh8b77901c9c9fbugfixes/contact_mechan…

aka_common.hhView Options

rAKA/src/common/aka_common.hh
8b77901c9c9fbugfixes/contact_mechan…

aka_common.hh
View Options