[aspect-devel] Benchmark

[Timo Heister]

> the updated figures that are attached. I also enclose the new .cc file. Do
> you have any other ideas for this behaviour?

Maybe:
./ This warning highlights a problem in your code:
burstedde.cc: In member function ‘dealii::Tensor<1, dim>
aspect::BursteddeBenchmark::BursteddeGravity<dim>::gravity_vector(const
dealii::Point<dim, double>&) const [with int dim = 2]’:
burstedde.cc:512:69: warning: ‘beta’ is used uninitialized in this
function [-Wuninitialized]
       const double mu=exp(1. - beta * (x*(1.-x)+y*(1.-y) + z*(1.-z)));
./ Your error computations won't work correctly as is when using more
than one MPI process. Are you?
./ This doesn't look too bad:
     Errors u_L1, p_L1, u_L2, p_L2: 2.781543e-04, 1.430527e-03,
3.116033e-04, 2.626537e-03 ref=2
     Errors u_L1, p_L1, u_L2, p_L2: 3.434312e-05, 3.123564e-04,
3.892204e-05, 6.297030e-04 ref=3
     Errors u_L1, p_L1, u_L2, p_L2: 4.505512e-06, 8.855139e-05,
4.880252e-06, 1.904190e-04 ref=4


-- 
Timo Heister
http://www.math.clemson.edu/~heister/