Diffusion akantu (bugfixes/bug_fix_in_synchronization_of_mat_vec_mul_product_when_assembled_to_residual) Browse

bug_fix_in_syn…

dof_manager_default.hh
View Options

			/**
			* @file dof_manager_default.hh
			*
			* @author Nicolas Richart <nicolas.richart@epfl.ch>
			*
			* @date creation: Tue Aug 18 2015
			* @date last modification: Wed Jan 31 2018
			*
			* @brief Default implementation of the dof manager
			*
			*
			* Copyright (©) 2015-2018 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/>.
			*
			*/

			/* -------------------------------------------------------------------------- */
			#include "dof_manager.hh"
			/* -------------------------------------------------------------------------- */
			#include <functional>
			#include <unordered_map>
			/* -------------------------------------------------------------------------- */

			#ifndef AKANTU_DOF_MANAGER_DEFAULT_HH_
			#define AKANTU_DOF_MANAGER_DEFAULT_HH_

			namespace akantu {
			class SparseMatrixAIJ;
			class NonLinearSolverDefault;
			class TimeStepSolverDefault;
			class DOFSynchronizer;
			} // namespace akantu

			namespace akantu {

			class DOFManagerDefault : public DOFManager {
			/* ------------------------------------------------------------------------ */
			/* Constructors/Destructors */
			/* ------------------------------------------------------------------------ */
			public:
			DOFManagerDefault(const ID & id = "dof_manager_default");
			DOFManagerDefault(Mesh & mesh, const ID & id = "dof_manager_default");
			~DOFManagerDefault() override;

			protected:
			struct DOFDataDefault : public DOFData {
			explicit DOFDataDefault(const ID & dof_id);
			};

			/* ------------------------------------------------------------------------ */
			/* Methods */
			/* ------------------------------------------------------------------------ */
			public:
			// /// register an array of degree of freedom
			// void registerDOFs(const ID & dof_id, Array<Real> & dofs_array,
			// const DOFSupportType & support_type) override;

			// /// the dof as an implied type of _dst_nodal and is defined only on a
			// subset
			// /// of nodes
			// void registerDOFs(const ID & dof_id, Array<Real> & dofs_array,
			// const ID & group_support) override;

			/**
			* Assemble elementary values to the global matrix. The dof number is
			* implicitly considered as conn(el, n) * nb_nodes_per_element + d.
			* With 0 < n < nb_nodes_per_element and 0 < d < nb_dof_per_node
			**/
			void assembleElementalMatricesToMatrix(
			const ID & matrix_id, const ID & dof_id,
			const Array<Real> & elementary_mat, ElementType type,
			GhostType ghost_type, const MatrixType & elemental_matrix_type,
			const Array<UInt> & filter_elements) override;

			void assembleMatMulVectToArray(const ID & dof_id, const ID & A_id,
			const Array<Real> & x, Array<Real> & array,
			Real scale_factor = 1.) override;

			/// multiply a vector by a lumped matrix and assemble the result to the
			/// residual
			void assembleLumpedMatMulVectToResidual(const ID & dof_id, const ID & A_id,
			const Array<Real> & x,
			Real scale_factor = 1) override;

			/// assemble coupling terms between to dofs
			void assemblePreassembledMatrix(const ID & dof_id_m, const ID & dof_id_n,
			const ID & matrix_id,
			const TermsToAssemble & terms) override;

			protected:
			void assembleToGlobalArray(const ID & dof_id,
			const Array<Real> & array_to_assemble,
			SolverVector & global_array,
			Real scale_factor) override;

			template <typename T>
			void assembleToGlobalArray(const ID & dof_id,
			const Array<T> & array_to_assemble,
			Array<T> & global_array, T scale_factor);

			void getArrayPerDOFs(const ID & dof_id, const SolverVector & global,
			Array<Real> & local) override;

			template <typename T>
			void getArrayPerDOFs(const ID & dof_id, const Array<T> & global_array,
			Array<T> & local_array) const;
			void makeConsistentForPeriodicity(const ID & dof_id,
			SolverVector & array) override;

			public:
			/// update the global dofs vector
			void updateGlobalBlockedDofs() override;

			// /// apply boundary conditions to jacobian matrix
			// void applyBoundary(const ID & matrix_id = "J") override;

			private:
			/// Add a symmetric matrices to a symmetric sparse matrix
			void addSymmetricElementalMatrixToSymmetric(
			SparseMatrixAIJ & matrix, const Matrix<Real> & element_mat,
			const Vector<Int> & equation_numbers, UInt max_size);

			/// Add a unsymmetric matrices to a symmetric sparse matrix (i.e. cohesive
			/// elements)
			void addUnsymmetricElementalMatrixToSymmetric(
			SparseMatrixAIJ & matrix, const Matrix<Real> & element_mat,
			const Vector<Int> & equation_numbers, UInt max_size);

			/// Add a matrices to a unsymmetric sparse matrix
			void addElementalMatrixToUnsymmetric(SparseMatrixAIJ & matrix,
			const Matrix<Real> & element_mat,
			const Vector<Int> & equation_numbers,
			UInt max_size);

			void addToProfile(const ID & matrix_id, const ID & dof_id,
			ElementType type, GhostType ghost_type);

			/* ------------------------------------------------------------------------ */
			/* MeshEventHandler interface */
			/* ------------------------------------------------------------------------ */
			protected:
			std::tuple<UInt, UInt, UInt>
			registerDOFsInternal(const ID & dof_id, Array<Real> & dofs_array) override;

			// std::pair<UInt, UInt>
			// updateNodalDOFs(const ID & dof_id, const Array<UInt> & nodes_list)
			// override;

			void resizeGlobalArrays() override;

			public:
			/// function to implement to react on akantu::NewNodesEvent
			void onNodesAdded(const Array<UInt> & nodes_list,
			const NewNodesEvent & event) override;

			/* ------------------------------------------------------------------------ */
			/* Accessors */
			/* ------------------------------------------------------------------------ */
			public:
			/// Get an instance of a new SparseMatrix
			SparseMatrix & getNewMatrix(const ID & matrix_id,
			const MatrixType & matrix_type) override;

			/// Get an instance of a new SparseMatrix as a copy of the SparseMatrix
			/// matrix_to_copy_id
			SparseMatrix & getNewMatrix(const ID & matrix_id,
			const ID & matrix_to_copy_id) override;

			/// Get the reference of an existing matrix
			SparseMatrixAIJ & getMatrix(const ID & matrix_id);

			/// Get an instance of a new lumped matrix
			SolverVector & getNewLumpedMatrix(const ID & matrix_id) override;

			/* ------------------------------------------------------------------------ */
			/* Non Linear Solver */
			/* ------------------------------------------------------------------------ */
			/// Get instance of a non linear solver
			NonLinearSolver & getNewNonLinearSolver(
			const ID & nls_solver_id,
			const NonLinearSolverType & _non_linear_solver_type) override;

			/* ------------------------------------------------------------------------ */
			/* Time-Step Solver */
			/* ------------------------------------------------------------------------ */
			/// Get instance of a time step solver
			TimeStepSolver &
			getNewTimeStepSolver(const ID & id, const TimeStepSolverType & type,
			NonLinearSolver & non_linear_solver,
			SolverCallback & solver_callback) override;

			/* ------------------------------------------------------------------------ */
			private:
			/// Get the solution array
			Array<Real> & getSolutionArray();

			/// Get the residual array
			const Array<Real> & getResidualArray() const;

			/// Get the residual array
			Array<Real> & getResidualArray();

			public:
			/// access the internal dof_synchronizer
			AKANTU_GET_MACRO_NOT_CONST(Synchronizer, *synchronizer, DOFSynchronizer &);

			/// access the internal dof_synchronizer
			bool hasSynchronizer() const { return synchronizer != nullptr; }

			Array<bool> & getBlockedDOFs();
			const Array<bool> & getBlockedDOFs() const;

			protected:
			std::unique_ptr<DOFData> getNewDOFData(const ID & dof_id) override;

			/* ------------------------------------------------------------------------ */
			/* Class Members */
			/* ------------------------------------------------------------------------ */
			protected:
			using DOFToMatrixProfile =
			std::map<std::pair<ID, ID>,
			std::vector<std::pair<ElementType, GhostType>>>;

			/// contains the the dofs that where added to the profile of a given matrix.
			DOFToMatrixProfile matrix_profiled_dofs;

			/// synchronizer to maintain coherency in dof fields
			std::unique_ptr<DOFSynchronizer> synchronizer;

			friend class DOFSynchronizer;

			/// Array containing the true or false if the node is in global_blocked_dofs
			Array<bool> global_blocked_dofs_uint;
			};

			} // namespace akantu

			#include "dof_manager_default_inline_impl.hh"

			#endif /* AKANTU_DOF_MANAGER_DEFAULT_HH_ */

rAKA/src/model/common/dof_manager/dof_manager_default.hh33327d9f46e4bugfixes/bug_fix_in_syn…

dof_manager_default.hhView Options

rAKA/src/model/common/dof_manager/dof_manager_default.hh
33327d9f46e4bugfixes/bug_fix_in_syn…

dof_manager_default.hh
View Options