[aspect-devel] VectorTools::integrate_difference fails to return an error that decreases under grid refinement (Harsha Lokavarapu)
Harsha Lokavarapu
lokavarapuh at gmail.com
Thu Oct 27 11:13:18 PDT 2016
With "set use year instead of second = true," the convergence results
for the 2D spherical shell model using
VectorTools::integrate_difference() are:
TABLE 03
Refinement
Level | L1 Velocity | Rate | L2
Velocity | Rate
-----------------------------------------------------------------------------------------
2 | 2.570197e-03 | | 1.433416e-03 |
3 | 5.152276e-04 | 2.319 | 2.951068e-04 |
2.280
4 | 1.121537e-04 | 2.200 | 6.326037e-05 |
2.222
5 | 2.590115e-05 | 2.114 | 1.422224e-05 |
2.153
Learned that one needs to always be vary of dimensions.
Thanks for the comments.
- Harsha
On Wednesday 26 October 2016 12:00 PM,
aspect-devel-request at geodynamics.org wrote:
> Send Aspect-devel mailing list submissions to
> aspect-devel at geodynamics.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
> or, via email, send a message with subject or body 'help' to
> aspect-devel-request at geodynamics.org
>
> You can reach the person managing the list at
> aspect-devel-owner at geodynamics.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Aspect-devel digest..."
>
>
> Today's Topics:
>
> 1. Re: VectorTools::integrate_difference fails to return an
> error that decreases under grid refinement (Elbridge Gerry Puckett)
> 2. Re: VectorTools::integrate_difference fails to return an
> error that decreases under grid refinement (Wolfgang Bangerth)
>
>
> _______________________________________________
> 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/20161027/16863780/attachment.html>
More information about the Aspect-devel
mailing list