Small bug in ellipt example.

22315a82 · Thomas Witkowski · ba416bd2 · 22315a82
Commit 22315a82 authored Jul 20, 2009 by Thomas Witkowski
Hide whitespace changes
Inline Side-by-side

Showing with 4 additions and 1 deletion

demo/src/ellipt.cc demo/src/ellipt.cc +4 -1

No files found.
--- 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 =====