1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
11 OpenFOAM is free software; you can redistribute it and/or modify it
12 under the terms of the GNU General Public License as published by the
13 Free Software Foundation; either version 2 of the License, or (at your
14 option) any later version.
16 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM; if not, write to the Free Software Foundation,
23 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
25 \*---------------------------------------------------------------------------*/
27 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
29 template<class Form, class Type>
30 inline Foam::Matrix<Form, Type>::Matrix()
38 template<class Form, class Type>
39 inline Foam::autoPtr<Foam::Matrix<Form, Type> > Foam::Matrix<Form, Type>::clone() const
41 return autoPtr<Matrix<Form, Type> >(new Matrix<Form, Type>(*this));
45 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
47 template<class Form, class Type>
48 inline const Foam::Matrix<Form, Type>& Foam::Matrix<Form, Type>::null()
50 return *reinterpret_cast< Matrix<Form, Type>* >(0);
54 //- Return the number of rows
55 template<class Form, class Type>
56 inline Foam::label Foam::Matrix<Form, Type>::n() const
62 template<class Form, class Type>
63 inline Foam::label Foam::Matrix<Form, Type>::m() const
69 template<class Form, class Type>
70 inline Foam::label Foam::Matrix<Form, Type>::size() const
76 template<class Form, class Type>
77 inline void Foam::Matrix<Form, Type>::checki(const label i) const
81 FatalErrorIn("Matrix<Form, Type>::checki(const label)")
82 << "attempt to access element from zero sized row"
85 else if (i<0 || i>=n_)
87 FatalErrorIn("Matrix<Form, Type>::checki(const label)")
88 << "index " << i << " out of range 0 ... " << n_-1
94 template<class Form, class Type>
95 inline void Foam::Matrix<Form, Type>::checkj(const label j) const
99 FatalErrorIn("Matrix<Form, Type>::checkj(const label)")
100 << "attempt to access element from zero sized column"
101 << abort(FatalError);
103 else if (j<0 || j>=m_)
105 FatalErrorIn("Matrix<Form, Type>::checkj(const label)")
106 << "index " << j << " out of range 0 ... " << m_-1
107 << abort(FatalError);
112 // * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
114 template<class Form, class Type>
115 inline Type* Foam::Matrix<Form, Type>::operator[](const label i)
124 template<class Form, class Type>
125 inline const Type* Foam::Matrix<Form, Type>::operator[](const label i) const
134 // ************************************************************************* //