[aspect-devel] Benchmark

[I. van Zelst]

Dear Timo,

I ran the benchmark with only 1 core, and I got the same results as you!
Which makes me wonder, why it can only be run on 1 core. Do you know that?
However, besides running your test, I also ran a model with increased
resolution of ncellx = 32 or ref = 5. As soon as I did that, the errors and
maximum and minimum pressure were all wrong again, as you can see in the
enclosed figures. The reason for this is probably that the corners of the
domain show very extreme pressures, almost as if the model is not well
resolved there. Of course, these excesses in pressure lead automatically to
an increase in error. Could the reason for these badly resolved pressures
in the corners of the domain be that the solver can not handle the
resulting matrix very well, because of the many boundary conditions at the
corner elements?
Thanks in advance!

Kind regards,
Iris


2014-07-01 20:56 GMT+02:00 Timo Heister <heister at clemson.edu>:

> > 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/
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140702/939a11b2/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: errors.pdf
Type: application/pdf
Size: 15779 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140702/939a11b2/attachment-0002.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pressure.pdf
Type: application/pdf
Size: 8928 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20140702/939a11b2/attachment-0003.pdf>