diff --git a/demo/src/ball.cc b/demo/src/ball.cc
index 516337f90cd812f937ddcf6d3f3edfeb0f6ac46d..486d6310973725aae85054b61f40511de0084071 100644
--- a/demo/src/ball.cc
+++ b/demo/src/ball.cc
@@ -70,13 +70,13 @@ int main(int argc, char* argv[])
 					       adaptInfo);
   
   // ===== create matrix operator =====
-  Operator matrixOperator(Operator::MATRIX_OPERATOR, ball.getFESpace());
+  Operator matrixOperator(Operator::MATRIX_OPERATOR, ball.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   ball.addMatrixOperator(&matrixOperator);
 
   // ===== create rhs operator =====
-  int degree = ball.getFESpace()->getBasisFcts()->getDegree();
-  Operator rhsOperator(Operator::VECTOR_OPERATOR, ball.getFESpace());
+  int degree = ball.getFeSpace()->getBasisFcts()->getDegree();
+  Operator rhsOperator(Operator::VECTOR_OPERATOR, ball.getFeSpace());
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   ball.addVectorOperator(&rhsOperator);
 
diff --git a/demo/src/bunny.cc b/demo/src/bunny.cc
index 50733f3560bb9ea97b69dcda48ab8bf0a470ca54..45d108e03fc24aeb3a55c8ee0449ac3ac444bdad 100644
--- a/demo/src/bunny.cc
+++ b/demo/src/bunny.cc
@@ -68,15 +68,15 @@ int main(int argc, char* argv[])
   
   // ===== create matrix operator =====
   Operator matrixOperator(Operator::MATRIX_OPERATOR,
-			  bunny.getFESpace());
+			  bunny.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   bunny.addMatrixOperator(&matrixOperator);
 
   // ===== create rhs operator =====
   Operator rhsOperator(Operator::VECTOR_OPERATOR,
-		       bunny.getFESpace());
+		       bunny.getFeSpace());
 
-  int degree = bunny.getFESpace()->getBasisFcts()->getDegree();
+  int degree = bunny.getFeSpace()->getBasisFcts()->getDegree();
 
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   bunny.addVectorOperator(&rhsOperator);
diff --git a/demo/src/couple.cc b/demo/src/couple.cc
index f1af8da1d3828553e9818d65623e78e68d5a1d90..0bf3783197a684d242b537d8b9adb8902cce4ea0 100644
--- a/demo/src/couple.cc
+++ b/demo/src/couple.cc
@@ -142,21 +142,21 @@ int main(int argc, char* argv[])
 		      adoptFlag);
 
   // ===== create operators for problem1 =====
-  Operator matrixOperator1(Operator::MATRIX_OPERATOR, problem1.getFESpace());
+  Operator matrixOperator1(Operator::MATRIX_OPERATOR, problem1.getFeSpace());
   matrixOperator1.addSecondOrderTerm(new Laplace_SOT);
   problem1.addMatrixOperator(&matrixOperator1);
 
-  int degree = problem1.getFESpace()->getBasisFcts()->getDegree();
-  Operator rhsOperator1(Operator::VECTOR_OPERATOR, problem1.getFESpace());
+  int degree = problem1.getFeSpace()->getBasisFcts()->getDegree();
+  Operator rhsOperator1(Operator::VECTOR_OPERATOR, problem1.getFeSpace());
   rhsOperator1.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   problem1.addVectorOperator(&rhsOperator1);
 
   // ===== create operators for problem2 =====
-  Operator matrixOperator2(Operator::MATRIX_OPERATOR, problem2.getFESpace());
+  Operator matrixOperator2(Operator::MATRIX_OPERATOR, problem2.getFeSpace());
   matrixOperator2.addZeroOrderTerm(new Simple_ZOT);
   problem2.addMatrixOperator(&matrixOperator2);
   
-  Operator rhsOperator2(Operator::VECTOR_OPERATOR, problem2.getFESpace());
+  Operator rhsOperator2(Operator::VECTOR_OPERATOR, problem2.getFeSpace());
   rhsOperator2.addZeroOrderTerm(new VecAtQP_ZOT(problem1.getSolution(), 
 						new Identity(degree)));
   problem2.addVectorOperator(&rhsOperator2);
diff --git a/demo/src/ellipt.cc b/demo/src/ellipt.cc
index 1eaf426ce411234ea8f5a74f8bdae11fb7beb075..4498e05e530a180edd1a21928483f0035c696597 100644
--- a/demo/src/ellipt.cc
+++ b/demo/src/ellipt.cc
@@ -66,14 +66,14 @@ int main(int argc, char* argv[])
 
   
   // ===== create matrix operator =====
-  Operator matrixOperator(Operator::MATRIX_OPERATOR, ellipt.getFESpace());
+  Operator matrixOperator(Operator::MATRIX_OPERATOR, ellipt.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   ellipt.addMatrixOperator(matrixOperator);
 
 
   // ===== create rhs operator =====
-  int degree = ellipt.getFESpace()->getBasisFcts()->getDegree();
-  Operator rhsOperator(Operator::VECTOR_OPERATOR, ellipt.getFESpace());
+  int degree = ellipt.getFeSpace()->getBasisFcts()->getDegree();
+  Operator rhsOperator(Operator::VECTOR_OPERATOR, ellipt.getFeSpace());
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   ellipt.addVectorOperator(rhsOperator);
 
diff --git a/demo/src/heat.cc b/demo/src/heat.cc
index bb414fee1bc18db79bbf2074caa1d502e06749e4..a605ae8eaab985a27046b2c1dd057dfa8d2bc2d1 100644
--- a/demo/src/heat.cc
+++ b/demo/src/heat.cc
@@ -194,7 +194,7 @@ int main(int argc, char** argv)
 
 
   // ===== create rhs functions =====
-  int degree = heatSpace.getFESpace()->getBasisFcts()->getDegree();
+  int degree = heatSpace.getFeSpace()->getBasisFcts()->getDegree();
   F *rhsFct = new F(degree);
   rhsFct->setTimePtr(heat.getRHSTimePtr());
 
@@ -205,7 +205,7 @@ int main(int argc, char** argv)
 
   // create laplace
   Operator A(Operator::MATRIX_OPERATOR | Operator::VECTOR_OPERATOR, 
-	     heatSpace.getFESpace());
+	     heatSpace.getFeSpace());
   A.addSecondOrderTerm(new Laplace_SOT);
   A.setUhOld(heat.getOldSolution());
   if (*(heat.getThetaPtr()) != 0.0)
@@ -215,14 +215,14 @@ int main(int argc, char** argv)
 
   // create zero order operator
   Operator C(Operator::MATRIX_OPERATOR | Operator::VECTOR_OPERATOR,
-	     heatSpace.getFESpace());
+	     heatSpace.getFeSpace());
   C.addZeroOrderTerm(new Simple_ZOT);
   C.setUhOld(heat.getOldSolution());
   heatSpace.addMatrixOperator(C, heat.getTau1Ptr(), heat.getTau1Ptr());
   heatSpace.addVectorOperator(C, heat.getTau1Ptr(), heat.getTau1Ptr());
 
   // create RHS operator
-  Operator F(Operator::VECTOR_OPERATOR, heatSpace.getFESpace());
+  Operator F(Operator::VECTOR_OPERATOR, heatSpace.getFeSpace());
   F.addZeroOrderTerm(new CoordsAtQP_ZOT(rhsFct));
   heatSpace.addVectorOperator(F);
 
diff --git a/demo/src/neumann.cc b/demo/src/neumann.cc
index cd825be11a386c99e42185abb82950a3fb9744fe..dc134bddade7961901e57719926037b3fefc37d1 100644
--- a/demo/src/neumann.cc
+++ b/demo/src/neumann.cc
@@ -63,13 +63,13 @@ int main(int argc, char* argv[])
 					       adaptInfo);
   
   // ===== create matrix operator =====
-  Operator matrixOperator(Operator::MATRIX_OPERATOR, neumann.getFESpace());
+  Operator matrixOperator(Operator::MATRIX_OPERATOR, neumann.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   neumann.addMatrixOperator(&matrixOperator);
 
   // ===== create rhs operator =====
-  int degree = neumann.getFESpace()->getBasisFcts()->getDegree();
-  Operator rhsOperator(Operator::VECTOR_OPERATOR, neumann.getFESpace());
+  int degree = neumann.getFeSpace()->getBasisFcts()->getDegree();
+  Operator rhsOperator(Operator::VECTOR_OPERATOR, neumann.getFeSpace());
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   neumann.addVectorOperator(&rhsOperator);
 
diff --git a/demo/src/parallelellipt.cc b/demo/src/parallelellipt.cc
index fe06109906d09c5b18159af3b6e85a0cccc2101a..8b19ff431f85e3f6127edc30e85c5f33bb5b0847 100644
--- a/demo/src/parallelellipt.cc
+++ b/demo/src/parallelellipt.cc
@@ -72,13 +72,13 @@ int main(int argc, char* argv[])
 					       adaptInfo);
   
   // ===== create matrix operator =====
-  Operator matrixOperator(Operator::MATRIX_OPERATOR, ellipt.getFESpace());
+  Operator matrixOperator(Operator::MATRIX_OPERATOR, ellipt.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   ellipt.addMatrixOperator(&matrixOperator);
 
   // ===== create rhs operator =====
-  int degree = ellipt.getFESpace()->getBasisFcts()->getDegree();
-  Operator rhsOperator(Operator::VECTOR_OPERATOR, ellipt.getFESpace());
+  int degree = ellipt.getFeSpace()->getBasisFcts()->getDegree();
+  Operator rhsOperator(Operator::VECTOR_OPERATOR, ellipt.getFeSpace());
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   ellipt.addVectorOperator(&rhsOperator);
 
diff --git a/demo/src/parallelheat.cc b/demo/src/parallelheat.cc
index 6c774dcee46ad3b0c5d4746d3145519fac33bf3a..9927c532da24cfeb654e27f0808bdf45b26addac 100644
--- a/demo/src/parallelheat.cc
+++ b/demo/src/parallelheat.cc
@@ -175,7 +175,7 @@ int main(int argc, char** argv)
 			  adaptInfoInitial);
 
   // ===== create rhs functions =====
-  int degree = heatSpace->getFESpace()->getBasisFcts()->getDegree();
+  int degree = heatSpace->getFeSpace()->getBasisFcts()->getDegree();
   F *rhsFct = new F(degree);
   rhsFct->setTimePtr(heat->getTimePtr());
 
@@ -186,7 +186,7 @@ int main(int argc, char** argv)
   // create laplace
   Operator *A = new Operator(Operator::MATRIX_OPERATOR | 
 			     Operator::VECTOR_OPERATOR,
-			     heatSpace->getFESpace());
+			     heatSpace->getFeSpace());
 
   A->addSecondOrderTerm(new Laplace_SOT);
 
@@ -197,7 +197,7 @@ int main(int argc, char** argv)
   // create zero order operator
   Operator *C = new Operator(Operator::MATRIX_OPERATOR | 
 			     Operator::VECTOR_OPERATOR,
-			     heatSpace->getFESpace());
+			     heatSpace->getFeSpace());
 
   C->addZeroOrderTerm(new Simple_ZOT);
 
@@ -208,7 +208,7 @@ int main(int argc, char** argv)
 
   // create RHS operator
   Operator *F = new Operator(Operator::VECTOR_OPERATOR,
-			     heatSpace->getFESpace());
+			     heatSpace->getFeSpace());
 
   F->addZeroOrderTerm(new CoordsAtQP_ZOT(rhsFct));
 
diff --git a/demo/src/parametric.cc b/demo/src/parametric.cc
index e5c56040f828675747c2b8cb42f32052996e866b..c5097634fb0c866cc97a8e8dec88853827fb64f5 100644
--- a/demo/src/parametric.cc
+++ b/demo/src/parametric.cc
@@ -60,7 +60,7 @@ public:
 
     // ===== create coords vector =====
     for (int i = 0; i < Global::getGeo(WORLD); i++)
-      parametricCoords[i] = new DOFVector<double>(this->getFESpace(), 
+      parametricCoords[i] = new DOFVector<double>(this->getFeSpace(), 
 						  "parametric coords");
 
     // ===== create parametric object =====
@@ -76,7 +76,7 @@ public:
   void buildBeforeCoarsen(AdaptInfo *adaptInfo, Flag flag) {
     FUNCNAME("ParametricSphere::buildAfterCoarsen()");
     MSG("calculation of parametric coordinates\n");
-    int preDOFs = this->getFESpace()->getAdmin()->getNumberOfPreDOFs(VERTEX);
+    int preDOFs = this->getFeSpace()->getAdmin()->getNumberOfPreDOFs(VERTEX);
     int dim = this->getMesh()->getDim();
     int dow = Global::getGeo(WORLD);
     WorldVector<double> vertexCoords;
@@ -142,15 +142,15 @@ int main(int argc, char* argv[])
   
   // ===== create matrix operator =====
   Operator matrixOperator(Operator::MATRIX_OPERATOR,
-			  parametric.getFESpace());
+			  parametric.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   parametric.addMatrixOperator(&matrixOperator);
 
   // ===== create rhs operator =====
   Operator rhsOperator(Operator::VECTOR_OPERATOR,
-		       parametric.getFESpace());
+		       parametric.getFeSpace());
 
-  int degree = parametric.getFESpace()->getBasisFcts()->getDegree();
+  int degree = parametric.getFeSpace()->getBasisFcts()->getDegree();
 
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   parametric.addVectorOperator(&rhsOperator);
diff --git a/demo/src/periodic.cc b/demo/src/periodic.cc
index 628a4568696889b4309fa8fb95c5ed6d313b73b4..6330155f50a022567c75cf7dfb5f315ec8a8190d 100644
--- a/demo/src/periodic.cc
+++ b/demo/src/periodic.cc
@@ -53,13 +53,13 @@ int main(int argc, char* argv[])
 					       adaptInfo);
   
   // ===== create matrix operator =====
-  Operator matrixOperator(Operator::MATRIX_OPERATOR, periodic.getFESpace());
+  Operator matrixOperator(Operator::MATRIX_OPERATOR, periodic.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   periodic.addMatrixOperator(&matrixOperator);
 
   // ===== create rhs operator =====
-  int degree = periodic.getFESpace()->getBasisFcts()->getDegree();
-  Operator rhsOperator(Operator::VECTOR_OPERATOR, periodic.getFESpace());
+  int degree = periodic.getFeSpace()->getBasisFcts()->getDegree();
+  Operator rhsOperator(Operator::VECTOR_OPERATOR, periodic.getFeSpace());
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   periodic.addVectorOperator(&rhsOperator);
 
diff --git a/demo/src/sphere.cc b/demo/src/sphere.cc
index cc70c3d5f2486ee3809a311d51c7fabaac2b6e60..9b83b208d56af98e2901fd4263e5f2f1367f6a9f 100644
--- a/demo/src/sphere.cc
+++ b/demo/src/sphere.cc
@@ -56,15 +56,15 @@ int main(int argc, char* argv[])
   
   // ===== create matrix operator =====
   Operator matrixOperator(Operator::MATRIX_OPERATOR,
-			  sphere.getFESpace());
+			  sphere.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   sphere.addMatrixOperator(&matrixOperator);
 
   // ===== create rhs operator =====
   Operator rhsOperator(Operator::VECTOR_OPERATOR,
-		       sphere.getFESpace());
+		       sphere.getFeSpace());
 
-  int degree = sphere.getFESpace()->getBasisFcts()->getDegree();
+  int degree = sphere.getFeSpace()->getBasisFcts()->getDegree();
 
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   sphere.addVectorOperator(&rhsOperator);
diff --git a/demo/src/torus.cc b/demo/src/torus.cc
index 3dd603f58a4b07bf98896cffc68cb19ea6bfce28..98bb884b64718a8f9fe332831a76ca3ef1cbef66 100644
--- a/demo/src/torus.cc
+++ b/demo/src/torus.cc
@@ -117,15 +117,15 @@ int main(int argc, char* argv[])
   
   // ===== create matrix operator =====
   Operator matrixOperator(Operator::MATRIX_OPERATOR,
-			  torus.getFESpace());
+			  torus.getFeSpace());
   matrixOperator.addSecondOrderTerm(new Laplace_SOT);
   torus.addMatrixOperator(&matrixOperator);
 
   // ===== create rhs operator =====
   Operator rhsOperator(Operator::VECTOR_OPERATOR,
-		       torus.getFESpace());
+		       torus.getFeSpace());
 
-  int degree = torus.getFESpace()->getBasisFcts()->getDegree();
+  int degree = torus.getFeSpace()->getBasisFcts()->getDegree();
 
   F f(degree);
   rhsOperator.addZeroOrderTerm(new CoordsAtQP_ZOT(&f));
@@ -144,7 +144,7 @@ int main(int argc, char* argv[])
   DegreeOfFreedom dof;
   WorldVector<double> x;
  
-  const FiniteElemSpace *feSpace = torus.getFESpace();
+  const FiniteElemSpace *feSpace = torus.getFeSpace();
   const BasisFunction *basFcts = feSpace->getBasisFcts();
   int numBasFcts = basFcts->getNumber();
   DegreeOfFreedom *localIndices = new DegreeOfFreedom[numBasFcts];
diff --git a/demo/src/vecellipt.cc b/demo/src/vecellipt.cc
index 36c2f3fe46f9c908bd04de15b21003133a1ace52..605738eda6d0fb19251f712a5524d45572072339 100644
--- a/demo/src/vecellipt.cc
+++ b/demo/src/vecellipt.cc
@@ -65,24 +65,24 @@ int main(int argc, char* argv[])
   
   // ===== create matrix operators =====
   Operator matrixOperator00(Operator::MATRIX_OPERATOR, 
-			    vecellipt.getFESpace(0), vecellipt.getFESpace(0));
+			    vecellipt.getFeSpace(0), vecellipt.getFeSpace(0));
   matrixOperator00.addSecondOrderTerm(new Laplace_SOT);
   vecellipt.addMatrixOperator(matrixOperator00, 0, 0);
 
   Operator matrixOperator10(Operator::MATRIX_OPERATOR, 
-			    vecellipt.getFESpace(1), vecellipt.getFESpace(0));
+			    vecellipt.getFeSpace(1), vecellipt.getFeSpace(0));
   matrixOperator10.addZeroOrderTerm(new Simple_ZOT);
   vecellipt.addMatrixOperator(matrixOperator10, 1, 0);
 
   Operator matrixOperator11(Operator::MATRIX_OPERATOR,
-			    vecellipt.getFESpace(1), vecellipt.getFESpace(1));
+			    vecellipt.getFeSpace(1), vecellipt.getFeSpace(1));
   matrixOperator11.addZeroOrderTerm(new Simple_ZOT(-1.0));
   vecellipt.addMatrixOperator(matrixOperator11, 1, 1);
 
 
   // ===== create rhs operator =====
-  int degree = vecellipt.getFESpace(0)->getBasisFcts()->getDegree();
-  Operator rhsOperator0(Operator::VECTOR_OPERATOR, vecellipt.getFESpace(0));
+  int degree = vecellipt.getFeSpace(0)->getBasisFcts()->getDegree();
+  Operator rhsOperator0(Operator::VECTOR_OPERATOR, vecellipt.getFeSpace(0));
   rhsOperator0.addZeroOrderTerm(new CoordsAtQP_ZOT(new F(degree)));
   vecellipt.addVectorOperator(rhsOperator0, 0);
 
diff --git a/demo/src/vecheat.cc b/demo/src/vecheat.cc
index a2dfddaa4f7fea3f3d35dfdd177291a04fcdf3de..3936632e67543a7816dfdbebf8f9992e0f156bed 100644
--- a/demo/src/vecheat.cc
+++ b/demo/src/vecheat.cc
@@ -249,9 +249,9 @@ int main(int argc, char** argv)
   
 
   // ===== create rhs functions =====
-  F *rhsFct0 = new F(vecheatSpace.getFESpace(0)->getBasisFcts()->getDegree());
+  F *rhsFct0 = new F(vecheatSpace.getFeSpace(0)->getBasisFcts()->getDegree());
   rhsFct0->setTimePtr(vecheat.getRHSTimePtr());
-  F *rhsFct1 = new F(vecheatSpace.getFESpace(1)->getBasisFcts()->getDegree());
+  F *rhsFct1 = new F(vecheatSpace.getFeSpace(1)->getBasisFcts()->getDegree());
   rhsFct1->setTimePtr(vecheat.getRHSTimePtr());
 
   // ===== create operators =====
@@ -260,12 +260,12 @@ int main(int argc, char** argv)
 
   // create laplace
   Operator A00 (Operator::MATRIX_OPERATOR | Operator::VECTOR_OPERATOR,
-		vecheatSpace.getFESpace(0), vecheatSpace.getFESpace(0));
+		vecheatSpace.getFeSpace(0), vecheatSpace.getFeSpace(0));
   A00.addSecondOrderTerm(new Laplace_SOT);
   A00.setUhOld(vecheat.getOldSolution()->getDOFVector(0));
 
   Operator A11(Operator::MATRIX_OPERATOR | Operator::VECTOR_OPERATOR,
-	       vecheatSpace.getFESpace(1), vecheatSpace.getFESpace(1));
+	       vecheatSpace.getFeSpace(1), vecheatSpace.getFeSpace(1));
   A11.addSecondOrderTerm(new Laplace_SOT);
   A11.setUhOld(vecheat.getOldSolution()->getDOFVector(1));
 
@@ -281,7 +281,7 @@ int main(int argc, char** argv)
 
   // create zero order operator
   Operator C00(Operator::MATRIX_OPERATOR | Operator::VECTOR_OPERATOR,
-	       vecheatSpace.getFESpace(0), vecheatSpace.getFESpace(0));
+	       vecheatSpace.getFeSpace(0), vecheatSpace.getFeSpace(0));
   C00.addZeroOrderTerm(new Simple_ZOT);
   C00.setUhOld(vecheat.getOldSolution()->getDOFVector(0));
   vecheatSpace.addMatrixOperator(C00, 0, 0, 
@@ -290,7 +290,7 @@ int main(int argc, char** argv)
   
   
   Operator C11(Operator::MATRIX_OPERATOR | Operator::VECTOR_OPERATOR,
-	       vecheatSpace.getFESpace(1), vecheatSpace.getFESpace(1));
+	       vecheatSpace.getFeSpace(1), vecheatSpace.getFeSpace(1));
   C11.addZeroOrderTerm(new Simple_ZOT);
   C11.setUhOld(vecheat.getOldSolution()->getDOFVector(1));
   vecheatSpace.addMatrixOperator(C11, 1, 1, 			    
@@ -299,11 +299,11 @@ int main(int argc, char** argv)
 				 vecheat.getTau1Ptr(), vecheat.getTau1Ptr());
 
   // create RHS operator
-  Operator F0(Operator::VECTOR_OPERATOR, vecheatSpace.getFESpace(0));
+  Operator F0(Operator::VECTOR_OPERATOR, vecheatSpace.getFeSpace(0));
   F0.addZeroOrderTerm(new CoordsAtQP_ZOT(rhsFct0));
   vecheatSpace.addVectorOperator(F0, 0);
 
-  Operator F1 = Operator(Operator::VECTOR_OPERATOR, vecheatSpace.getFESpace(1));
+  Operator F1 = Operator(Operator::VECTOR_OPERATOR, vecheatSpace.getFeSpace(1));
   F1.addZeroOrderTerm(new CoordsAtQP_ZOT(rhsFct1));
   vecheatSpace.addVectorOperator(F1, 1);