[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