[aspect-devel] matrix assembly

Thomas Geenen geenen at gmail.com
Fri Aug 30 06:52:45 PDT 2013


looks good to me.

cheers
Thomas


On Fri, Aug 30, 2013 at 3:09 PM, Wolfgang Bangerth
<bangerth at math.tamu.edu>wrote:

>
>  The reasoning may have been that we have to assemble a term of the form
>>    u . grad phi_i * phi_j
>> with a total polynomial degree of
>>    stokes_degree + 2 * temp_degree - 1
>> which would suggest a Gauss formula of degree
>>    temp_degree + stokes_degree/2
>> rounded up if necessary. Does this sound reasonable? It would boil down to
>>    temp_degree + 1
>> in the most common case.
>>
>
> This would be the attached patch. Makes sense?
> W.
>
>
> --
> ------------------------------**------------------------------**
> ------------
> Wolfgang Bangerth               email:            bangerth at math.tamu.edu
>                                 www: http://www.math.tamu.edu/~**bangerth/<http://www.math.tamu.edu/~bangerth/>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://geodynamics.org/pipermail/aspect-devel/attachments/20130830/11142dc0/attachment-0001.html>


More information about the Aspect-devel mailing list