changed heat demo to have some working example

demo/init/heat.dat.2d

@@ -15,16 +15,17 @@ heat->space->solver->left precon:     no
 heat->space->solver->right precon:    no
 
 heat->space->estimator:              residual
-heat->space->estimator->C0:          1.0
-heat->space->estimator->C1:          1.0
+heat->space->estimator->C0:          0.1
+heat->space->estimator->C1:          0.1
 
 heat->theta:                         1.0
 
-heat->adapt->timestep:               0.1
+heat->adapt->timestep:               0.01
+heat->adapt->min timestep:           0.01
 heat->adapt->start time:             0.0
-heat->adapt->end time:               1.0
+heat->adapt->end time:               10.0
 
-heat->adapt->tolerance:              0.1
+heat->adapt->tolerance:              0.0001
 heat->adapt->rel space error:        0.5
 heat->adapt->rel time error:         0.5
 heat->adapt->time theta 1:           1.0
@@ -36,13 +37,14 @@ heat->adapt->time delta 2:           1.4142
 heat->adapt->max iteration:          100
 
 heat->adapt->coarsen allowed: 1   % 0|1
-heat->adapt->refine bisections: 2
-heat->adapt->coarsen bisections: 2
+heat->adapt->refine bisections: 1
+heat->adapt->coarsen bisections: 1
 
-heat->space->marker->strategy:       2    % 0=none, 1=GR, 2=MS, 3=ES, 4=GERS
+heat->space->marker->strategy:      4   % 0=none, 1=GR, 2=MS, 3=ES, 4=GERS
 heat->space->marker->MSGamma:        0.5
 heat->space->marker->MSGammaC:       0.1
 
+heat->space->output->ParaView animation: 1
 heat->space->output->filename:       output/heat
 heat->space->output->ParaView format: 1
 heat->space->output->write every i-th timestep: 1

demo/src/heat.cc

@@ -230,7 +230,7 @@ int main(int argc, char** argv)
   boundaryFct->setTimePtr(heat.getBoundaryTimePtr());
   heat.setExactSolution(boundaryFct);
 
-  heatSpace.addDirichletBC(DIRICHLET, boundaryFct);
+  // heatSpace.addDirichletBC(DIRICHLET, boundaryFct);
 
 
   // ===== start adaption loop =====