* be prepended by the user time used
* by the running program so far.
*/
- void log_execution_time(bool flag);
+ void log_execution_time(bool flag);
/**
* Output a constant something through
template <class T>
inline
LogStream &
-LogStream::
-operator<< (const T& t)
+LogStream::operator<< (const T& t)
{
// if the previous command was an
// #endl#, print the topmost prefix
template <int rank_, int dim>
inline
-Tensor<rank_,dim>::Tensor () {
+Tensor<rank_,dim>::Tensor ()
+{
// default constructor. not specifying an initializer list calls
// the default constructor of the subobjects, which initialize them
// selves. therefore, the tensor is set to zero this way
template <int rank_, int dim>
inline
-Tensor<rank_,dim>::Tensor (const array_type &initializer) {
+Tensor<rank_,dim>::Tensor (const array_type &initializer)
+{
for (unsigned int i=0; i<dim; ++i)
subtensor[i] = Tensor<rank_-1,dim>(initializer[i]);
};
template <int rank_, int dim>
inline
Tensor<rank_-1,dim> &
-Tensor<rank_,dim>::
-operator[] (const unsigned int i) {
+Tensor<rank_,dim>::operator[] (const unsigned int i)
+{
Assert (i<dim, ExcIndexRange(i, 0, dim));
return subtensor[i];
template <int rank_, int dim>
inline
const Tensor<rank_-1,dim> &
-Tensor<rank_,dim>::
-operator[] (const unsigned int i) const {
+Tensor<rank_,dim>::operator[] (const unsigned int i) const
+{
Assert (i<dim, ExcIndexRange(i, 0, dim));
return subtensor[i];
template <int rank_, int dim>
inline
Tensor<rank_,dim> &
-Tensor<rank_,dim>::
-operator = (const Tensor<rank_,dim> &t) {
+Tensor<rank_,dim>::operator = (const Tensor<rank_,dim> &t)
+{
for (unsigned int i=0; i<dim; ++i)
subtensor[i] = t.subtensor[i];
return *this;
template <int rank_, int dim>
inline
bool
-Tensor<rank_,dim>::
-operator == (const Tensor<rank_,dim> &p) const {
+Tensor<rank_,dim>::operator == (const Tensor<rank_,dim> &p) const
+{
for (unsigned int i=0; i<dim; ++i)
if (subtensor[i] != p.subtensor[i]) return false;
return true;
template <int rank_, int dim>
inline
bool
-Tensor<rank_,dim>::
-operator != (const Tensor<rank_,dim> &p) const {
+Tensor<rank_,dim>::operator != (const Tensor<rank_,dim> &p) const
+{
return !((*this) == p);
};
template <int rank_, int dim>
inline
Tensor<rank_,dim> &
-Tensor<rank_,dim>::
-operator += (const Tensor<rank_,dim> &p) {
+Tensor<rank_,dim>::operator += (const Tensor<rank_,dim> &p)
+{
for (unsigned int i=0; i<dim; ++i)
subtensor[i] += p.subtensor[i];
return *this;
template <int rank_, int dim>
inline
Tensor<rank_,dim> &
-Tensor<rank_,dim>::
-operator -= (const Tensor<rank_,dim> &p) {
+Tensor<rank_,dim>::operator -= (const Tensor<rank_,dim> &p)
+{
for (unsigned int i=0; i<dim; ++i)
subtensor[i] -= p.subtensor[i];
return *this;
template <int rank_, int dim>
inline
Tensor<rank_,dim> &
-Tensor<rank_,dim>::
-operator *= (const double &s) {
+Tensor<rank_,dim>::operator *= (const double &s)
+{
for (unsigned int i=0; i<dim; ++i)
subtensor[i] *= s;
return *this;
template <int rank_, int dim>
inline
Tensor<rank_,dim> &
-Tensor<rank_,dim>::
-operator /= (const double &s) {
+Tensor<rank_,dim>::operator /= (const double &s)
+{
for (unsigned int i=0; i<dim; ++i)
subtensor[i] /= s;
return *this;
template <int rank_, int dim>
inline
Tensor<rank_,dim>
-Tensor<rank_,dim>::
-operator + (const Tensor<rank_,dim> &t) const {
+Tensor<rank_,dim>::operator + (const Tensor<rank_,dim> &t) const
+{
Tensor<rank_,dim> tmp(*this);
for (unsigned int i=0; i<dim; ++i)
template <int rank_, int dim>
inline
Tensor<rank_,dim>
-Tensor<rank_,dim>::
-operator - (const Tensor<rank_,dim> &t) const {
+Tensor<rank_,dim>::operator - (const Tensor<rank_,dim> &t) const
+{
Tensor<rank_,dim> tmp(*this);
for (unsigned int i=0; i<dim; ++i)
template <int rank_, int dim>
inline
-void Tensor<rank_,dim>::clear () {
+void Tensor<rank_,dim>::clear ()
+{
for (unsigned int i=0; i<dim; ++i)
subtensor[i].clear();
};
template <int dim>
inline
-Tensor<1,dim>::Tensor (const bool initialize) {
+Tensor<1,dim>::Tensor (const bool initialize)
+{
Assert (dim>0, ExcDimTooSmall(dim));
if (initialize)
template <int dim>
inline
-Tensor<1,dim>::Tensor (const array_type &initializer) {
+Tensor<1,dim>::Tensor (const array_type &initializer)
+{
for (unsigned int i=0; i<dim; ++i)
values[i] = initializer[i];
};
template <int dim>
inline
-Tensor<1,dim>::Tensor (const Tensor<1,dim> &p) {
+Tensor<1,dim>::Tensor (const Tensor<1,dim> &p)
+{
for (unsigned int i=0; i<dim; ++i)
values[i] = p.values[i];
};
template <int dim>
inline
-double Tensor<1,dim>::
-operator [] (const unsigned int index) const {
+double Tensor<1,dim>::operator [] (const unsigned int index) const
+{
Assert (index<dim, ExcIndexRange (index, 0, dim));
return values[index];
};
template <int dim>
inline
-double & Tensor<1,dim>::
-operator [] (const unsigned int index) {
+double & Tensor<1,dim>::operator [] (const unsigned int index)
+{
Assert (index<dim, ExcIndexRange (index, 0, dim));
return values[index];
};
template <int dim>
inline
-Tensor<1,dim> & Tensor<1,dim>::
-operator = (const Tensor<1,dim> &p) {
+Tensor<1,dim> & Tensor<1,dim>::operator = (const Tensor<1,dim> &p)
+{
for (unsigned int i=0; i<dim; ++i)
values[i] = p.values[i];
return *this;
template <int dim>
inline
-bool Tensor<1,dim>::
-operator == (const Tensor<1,dim> &p) const {
+bool Tensor<1,dim>::operator == (const Tensor<1,dim> &p) const
+{
for (unsigned int i=0; i<dim; ++i)
if (values[i] != p.values[i]) return false;
return true;
template <int dim>
inline
-bool Tensor<1,dim>::
-operator != (const Tensor<1,dim> &p) const {
+bool Tensor<1,dim>::operator != (const Tensor<1,dim> &p) const
+{
return !((*this) == p);
};
template <int dim>
inline
-Tensor<1,dim> & Tensor<1,dim>::
-operator += (const Tensor<1,dim> &p) {
+Tensor<1,dim> & Tensor<1,dim>::operator += (const Tensor<1,dim> &p)
+{
for (unsigned int i=0; i<dim; ++i)
values[i] += p.values[i];
return *this;
template <int dim>
inline
-Tensor<1,dim> & Tensor<1,dim>::
-operator -= (const Tensor<1,dim> &p) {
+Tensor<1,dim> & Tensor<1,dim>::operator -= (const Tensor<1,dim> &p)
+{
for (unsigned int i=0; i<dim; ++i)
values[i] -= p.values[i];
return *this;
template <int dim>
inline
-Tensor<1,dim> & Tensor<1,dim>::
-operator *= (const double &s) {
+Tensor<1,dim> & Tensor<1,dim>::operator *= (const double &s)
+{
for (unsigned int i=0; i<dim; ++i)
values[i] *= s;
return *this;
template <int dim>
inline
-Tensor<1,dim> & Tensor<1,dim>::
-operator /= (const double &s) {
+Tensor<1,dim> & Tensor<1,dim>::operator /= (const double &s)
+{
for (unsigned int i=0; i<dim; ++i)
values[i] /= s;
return *this;
template <int dim>
inline
-double Tensor<1,dim>::
-operator * (const Tensor<1,dim> &p) const {
+double Tensor<1,dim>::operator * (const Tensor<1,dim> &p) const
+{
double q=0;
for (unsigned int i=0; i<dim; ++i)
q += values[i] * p.values[i];
template <int dim>
inline
-Tensor<1,dim> Tensor<1,dim>::
-operator + (const Tensor<1,dim> &p) const {
+Tensor<1,dim> Tensor<1,dim>::operator + (const Tensor<1,dim> &p) const
+{
return (Tensor<1,dim>(*this) += p);
};
template <int dim>
inline
-Tensor<1,dim> Tensor<1,dim>::
-operator - (const Tensor<1,dim> &p) const {
+Tensor<1,dim> Tensor<1,dim>::operator - (const Tensor<1,dim> &p) const
+{
return (Tensor<1,dim>(*this) -= p);
};
template <int dim>
inline
-void Tensor<1,dim>::clear () {
+void Tensor<1,dim>::clear ()
+{
for (unsigned int i=0; i<dim; ++i)
values[i] = 0;
};
template <int dim>
inline
-ostream & operator << (ostream &out, const Tensor<1,dim> &p) {
+ostream & operator << (ostream &out, const Tensor<1,dim> &p)
+{
for (unsigned int i=0; i<dim-1; ++i)
out << p[i] << ' ';
out << p[dim-1];
template <>
inline
-ostream & operator << (ostream &out, const Tensor<1,1> &p) {
+ostream & operator << (ostream &out, const Tensor<1,1> &p)
+{
out << p[0];
return out;
template <int dim, typename Accessor>
inline
bool
-TriaRawIterator<dim,Accessor>::
-operator < (const TriaRawIterator &i) const {
+TriaRawIterator<dim,Accessor>::operator < (const TriaRawIterator &i) const {
Assert (state() != invalid, ExcDereferenceInvalidObject());
Assert (i.state() != invalid, ExcDereferenceInvalidObject());
Assert (&accessor.get_triangulation() == &i.accessor.get_triangulation(),
template <int dim, typename Accessor>
inline
-TriaIterator<dim,Accessor>::
-TriaIterator (const TriaIterator<dim,Accessor> &i) :
+TriaIterator<dim,Accessor>::TriaIterator (const TriaIterator<dim,Accessor> &i) :
TriaRawIterator<dim,Accessor> (static_cast<TriaRawIterator<dim,Accessor> >(i)) {};
template <int dim, typename Accessor>
inline
-TriaIterator<dim,Accessor>::
-TriaIterator (const TriaRawIterator<dim,Accessor> &i) :
+TriaIterator<dim,Accessor>::TriaIterator (const TriaRawIterator<dim,Accessor> &i) :
TriaRawIterator<dim,Accessor> (i)
{
#ifdef DEBUG
template <int dim, typename Accessor>
inline
-TriaActiveIterator<dim,Accessor>::
-TriaActiveIterator (const TriaActiveIterator<dim,Accessor> &i) :
+TriaActiveIterator<dim,Accessor>::TriaActiveIterator (const TriaActiveIterator<dim,Accessor> &i) :
TriaIterator<dim,Accessor> (static_cast<TriaIterator<dim,Accessor> >(i)) {};
template <int dim, typename Accessor>
inline
-TriaActiveIterator<dim,Accessor>::
-TriaActiveIterator (const TriaRawIterator<dim,Accessor> &i) :
+TriaActiveIterator<dim,Accessor>::TriaActiveIterator (const TriaRawIterator<dim,Accessor> &i) :
TriaIterator<dim,Accessor> (i)
{
#ifdef DEBUG
template <int dim, typename Accessor>
inline
-TriaActiveIterator<dim,Accessor>::
-TriaActiveIterator (const TriaIterator<dim,Accessor> &i) :
+TriaActiveIterator<dim,Accessor>::TriaActiveIterator (const TriaIterator<dim,Accessor> &i) :
TriaIterator<dim,Accessor> (i)
{
#ifdef DEBUG
template<class SOLVER, class MATRIX, class PRECOND>
-MGCoarseGridLACIteration<SOLVER,MATRIX,PRECOND>::
-MGCoarseGridLACIteration(SOLVER &s,
+MGCoarseGridLACIteration<SOLVER,MATRIX,PRECOND>::MGCoarseGridLACIteration(SOLVER &s,
const MATRIX &m,
const PRECOND &p)
:
template<class SOLVER, class MATRIX, class PRECOND>
void
-MGCoarseGridLACIteration<SOLVER,MATRIX,PRECOND>::
-operator() (const unsigned int /* level */,
+MGCoarseGridLACIteration<SOLVER,MATRIX,PRECOND>::operator() (const unsigned int /* level */,
Vector<double> &dst,
const Vector<double> &src) const
{
template<typename number>
template<int dim>
-MGSmootherRelaxation<number>::
-MGSmootherRelaxation (const MGDoFHandler<dim> &mg_dof,
+MGSmootherRelaxation<number>::MGSmootherRelaxation (const MGDoFHandler<dim> &mg_dof,
const MGLevelObject<SparseMatrix<number> > &matrix,
const function_ptr relaxation,
const unsigned int steps,