From 22315a82dca63167d9e7bab8036efbce1dfbcea3 Mon Sep 17 00:00:00 2001
From: Thomas Witkowski <thomas.witkowski@gmx.de>
Date: Mon, 20 Jul 2009 11:23:01 +0000
Subject: [PATCH] Small bug in ellipt example.
---
demo/src/ellipt.cc | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/demo/src/ellipt.cc b/demo/src/ellipt.cc
index d924ab4c..6989f421 100644
--- a/demo/src/ellipt.cc
+++ b/demo/src/ellipt.cc
@@ -24,6 +24,8 @@ class F : public AbstractFunction<double, WorldVector<double> >
{
public:
+ F(int degree) : AbstractFunction<double, WorldVector<double> >(degree) {}
+
/// Implementation of AbstractFunction::operator().
double operator()(const WorldVector<double>& x) const
{
@@ -69,8 +71,9 @@ int main(int argc, char* argv[])
ellipt.addMatrixOperator(&matrixOperator);
// ===== create rhs operator =====
+ int degree = ellipt.getFESpace()->getBasisFcts()->getDegree();
Operator rhsOperator(Operator::VECTOR_OPERATOR, ellipt.getFESpace());
- rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F()));
+ rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
ellipt.addVectorOperator(&rhsOperator);
// ===== start adaption loop =====
--
GitLab