From b2f37389aa2dd9f91ccbd645becb0e3e5749d4b6 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Robert=20Kl=C3=B6fkorn?= <robertk@dune-project.org>
Date: Mon, 2 Feb 2004 15:43:33 +0000
Subject: [PATCH] Minor changes.
[[Imported from SVN: r330]]
---
fem/common/combinedoperator.hh | 5 +++--
fem/common/discretefunction.cc | 1 +
fem/common/quadrature.hh | 8 +++++++-
3 files changed, 11 insertions(+), 3 deletions(-)
diff --git a/fem/common/combinedoperator.hh b/fem/common/combinedoperator.hh
index d628bff2c..8c5f8cca6 100644
--- a/fem/common/combinedoperator.hh
+++ b/fem/common/combinedoperator.hh
@@ -131,7 +131,7 @@ namespace Dune
RangeFieldType _fB;
//! for temporary use
- Range *_tmp;
+ mutable Range *_tmp;
private:
//! Barton Nackman
@@ -167,7 +167,7 @@ namespace Dune
}
template <class GridIteratorType>
- void applyLocal ( GridIteratorType &it , const Domain & arg , Range & dest )
+ void applyLocal ( GridIteratorType &it , const Domain & arg , Range & dest ) const
{
std::cerr << "Combined Operator::applyLocal: No Default Implemenation is provieded!\n";
abort();
@@ -200,6 +200,7 @@ namespace Dune
template <class GridIteratorType>
void applyLocal ( GridIteratorType &it , const Domain & arg , Range & dest ) const
{
+ std::cout << "Warning: CombinedOperator<ADD>::applyLocal not correct! \n";
_b.applyLocal ( it , arg , dest );
_a.applyLocal ( it , arg , dest );
diff --git a/fem/common/discretefunction.cc b/fem/common/discretefunction.cc
index 1facb1267..f19a852b0 100644
--- a/fem/common/discretefunction.cc
+++ b/fem/common/discretefunction.cc
@@ -267,6 +267,7 @@ namespace Dune
file << d << " " << r << " ";
file << n << " " << m << "\n";
file << myId_ << " " << ftype << "\n";
+ //file << functionSpace_.continuous() << " " << functionSpace_.polynomOrder() << "\n";
file.close();
}
diff --git a/fem/common/quadrature.hh b/fem/common/quadrature.hh
index d3290fbb5..648bf644c 100644
--- a/fem/common/quadrature.hh
+++ b/fem/common/quadrature.hh
@@ -72,7 +72,13 @@ namespace Dune {
};
-
+ //***********************************************************************
+ //!
+ //! Note: The sum over all weigths for the quadrature is the volume of the
+ //! reference element, for example the sum over three quadrature point for
+ //! a triangle is 0.5 which is the volume on reference triangle
+ //!
+ //***********************************************************************
template <class Domain, class RangeField, ElementType ElType, int polOrd>
struct QuadraturePoints
{
--
GitLab