sparse_matrix_aij.hh
No OneTemporary
Actions

Subscribers

None

File Metadata

Created: Tue, May 14, 15:21

sparse_matrix_aij.hh
View Options

	/**
	* @file sparse_matrix_aij.hh
	*
	* @author Nicolas Richart <nicolas.richart@epfl.ch>
	*
	* @date Mon Aug 17 21:22:57 2015
	*
	* @brief AIJ implementation of the SparseMatrix (this the format used by Mumps)
	*
	* @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/>.
	*
	*/

	/* -------------------------------------------------------------------------- */
	#include "aka_common.hh"
	#include "aka_array.hh"
	#include "sparse_matrix.hh"
	/* -------------------------------------------------------------------------- */

	#ifndef __AKANTU_SPARSE_MATRIX_AIJ_HH__
	#define __AKANTU_SPARSE_MATRIX_AIJ_HH__

	namespace akantu {
	class DOFManagerDefault;
	}

	__BEGIN_AKANTU__

	class SparseMatrixAIJ : public SparseMatrix {
	/* ------------------------------------------------------------------------ */
	/* Constructors/Destructors */
	/* ------------------------------------------------------------------------ */
	public:
	SparseMatrixAIJ(DOFManagerDefault & dof_manager,
	const MatrixType & matrix_type,
	const ID & id = "sparse_matrix_aij");

	SparseMatrixAIJ(const SparseMatrixAIJ & matrix,
	const ID & id = "sparse_matrix_aij");

	virtual ~SparseMatrixAIJ();

	/* ------------------------------------------------------------------------ */
	/* Methods */
	/* ------------------------------------------------------------------------ */
	public:
	/// remove the existing profile
	inline void clearProfile();

	/// add a non-zero element
	virtual UInt addToProfile(UInt i, UInt j);

	/// set the matrix to 0
	virtual void clear();

	/// assemble a local matrix in the sparse one
	inline void addToMatrix(UInt i, UInt j, Real value);

	/// set the size of the matrix
	void resize(UInt size) { this->size = size; }

	/// modify the matrix to "remove" the blocked dof
	virtual void applyBoundary(Real block_val = 1.);

	/// save the profil in a file using the MatrixMarket file format
	virtual void saveProfile(const std::string & filename) const;

	/// save the matrix in a file using the MatrixMarket file format
	virtual void saveMatrix(const std::string & filename) const;

	/// copy assuming the profile are the same
	virtual void copyContent(const SparseMatrix & matrix);

	/// multiply the matrix by a scalar
	void mul(Real alpha);

	/// add matrix *this += B
	//virtual void add(const SparseMatrix & matrix, Real alpha);

	/// Equivalent of *gemv in blas
	virtual void matVecMul(const Array<Real> & x, Array<Real> & y,
	Real alpha = 1., Real beta = 0.) const;

	/* ------------------------------------------------------------------------ */
	/// accessor to A_{ij} - if (i, j) not present it returns 0
	inline Real operator()(UInt i, UInt j) const;

	/// accessor to A_{ij} - if (i, j) not present it fails, (i, j) should be
	/// first added to the profile
	inline Real & operator()(UInt i, UInt j);

	protected:
	/// This is the revert of add B += \alpha * *this;
	virtual void addMeTo(SparseMatrix & B, Real alpha) const;

	private:
	/// This is just to inline the addToMatrix function
	template<class MatrixType>
	void addMeToTemplated(MatrixType & B, Real alpha) const;

	/* ------------------------------------------------------------------------ */
	/* Accessors */
	/* ------------------------------------------------------------------------ */
	public:
	AKANTU_GET_MACRO(IRN, irn, const Array<Int> &);

	AKANTU_GET_MACRO(JCN, jcn, const Array<Int> &);

	AKANTU_GET_MACRO(A, a, const Array<Real> &);

	/// The release changes at each call of a function that changes the profile,
	/// it in increasing but could overflow so it should be checked as
	/// (my_release != release) and not as (my_release < release)
	AKANTU_GET_MACRO(ProfileRelease, profile_release, UInt);
	AKANTU_GET_MACRO(ValueRelease, value_release, UInt);

	protected:
	typedef std::pair<UInt, UInt> KeyCOO;
	typedef unordered_map<KeyCOO, UInt>::type coordinate_list_map;

	/// get the pair corresponding to (i, j)
	inline KeyCOO key(UInt i, UInt j) const {
	if (this->matrix_type == _symmetric && (i > j))
	return std::make_pair(j, i);
	return std::make_pair(i, j);
	}

	/* ------------------------------------------------------------------------ */
	/* Class Members */
	/* ------------------------------------------------------------------------ */
	private:
	DOFManagerDefault & dof_manager;

	/// row indexes
	Array<Int> irn;

	/// column indexes
	Array<Int> jcn;

	/// values : A[k] = Matrix[irn[k]][jcn[k]]
	Array<Real> a;

	/// Profile release
	UInt profile_release;

	/// Value release
	UInt value_release;

	/*
	* map for (i, j) -> k correspondence
	*/
	coordinate_list_map irn_jcn_k;
	};

	__END_AKANTU__

	/* -------------------------------------------------------------------------- */
	/* inline functions */
	/* -------------------------------------------------------------------------- */
	#include "sparse_matrix_aij_inline_impl.cc"

	#endif /* __AKANTU_SPARSE_MATRIX_AIJ_HH__ */

sparse_matrix_aij.hhNo OneTemporaryActions

File Metadata

sparse_matrix_aij.hhView Options

Event Timeline

sparse_matrix_aij.hh
No OneTemporary
Actions

sparse_matrix_aij.hh
View Options