demo updated

demo/init/vecellipt.dat.1d

@@ -9,9 +9,9 @@ vecellipt->components:                   2
 vecellipt->polynomial degree[0]:         1
 vecellipt->polynomial degree[1]:         1
 
-vecellipt->solver:                       bicgstab
+vecellipt->solver:                       gmres
 vecellipt->solver->max iteration:        1000
-vecellipt->solver->tolerance:            1.e-10
+vecellipt->solver->tolerance:            1.e-8
 vecellipt->solver->info:                 2
 vecellipt->solver->left precon:          diag
 vecellipt->solver->right precon:         no

demo/init/vecellipt.dat.2d

@@ -9,9 +9,9 @@ vecellipt->components:                   2
 vecellipt->polynomial degree[0]:         1
 vecellipt->polynomial degree[1]:         1
 
-vecellipt->solver:                       bicgstab 
+vecellipt->solver:                       gmres
 vecellipt->solver->max iteration:        1000
-vecellipt->solver->tolerance:            1.e-10
+vecellipt->solver->tolerance:            1.e-8
 vecellipt->solver->info:                 2
 vecellipt->solver->left precon:          diag
 vecellipt->solver->right precon:         no

demo/init/vecellipt.dat.3d

@@ -9,11 +9,11 @@ vecellipt->components:                   2
 vecellipt->polynomial degree[0]:         1
 vecellipt->polynomial degree[1]:         1
 
-vecellipt->solver:                       bicgstab
+vecellipt->solver:                       gmres
 vecellipt->solver->max iteration:        1000
-vecellipt->solver->tolerance:            1.e-10
+vecellipt->solver->tolerance:            1.e-8
 vecellipt->solver->info:                 2
-vecellipt->solver->left precon:          ilu
+vecellipt->solver->left precon:          diag
 vecellipt->solver->right precon:         no
 
 vecellipt->estimator[0]:                 residual

demo/src/vecellipt.cc

@@ -63,7 +63,7 @@ int main(int argc, char* argv[])
   vecellipt.addMatrixOperator(matrixOperator00, 0, 0);
 
   Operator matrixOperator10(vecellipt.getFeSpace(1), vecellipt.getFeSpace(0));
-  matrixOperator10.addTerm(new Simple_ZOT);
+  matrixOperator10.addTerm(new Simple_ZOT(-1.0));
   vecellipt.addMatrixOperator(matrixOperator10, 1, 0);
 
   Operator matrixOperator11(vecellipt.getFeSpace(0), vecellipt.getFeSpace(0));
@@ -78,7 +78,7 @@ int main(int argc, char* argv[])
   vecellipt.addVectorOperator(rhsOperator0, 0);
 
   // ===== add boundary conditions =====
-  vecellipt.addDirichletBC(1, 0, 0, new G);
+  vecellipt.addDirichletBC(1, 1, 1, new G);
 
 
   // ===== start adaption loop =====