ufl team mailing list archive

[Harish Narayanan <harish.mlists@xxxxxxxxx>]

Anders Logg wrote:
> On Fri, Dec 04, 2009 at 11:45:54AM +0100, Anders Logg wrote:
>> On Fri, Dec 04, 2009 at 08:31:41PM +0100, Marie Rognes wrote:
>>> noreply@xxxxxxxxxxxxx wrote:
>>>> ------------------------------------------------------------
>>>> revno: 860
>>>> committer: Anders Logg <logg@xxxxxxxxx>
>>>> branch nick: ufl-dev
>>>> timestamp: Tue 2009-12-01 18:09:26 +0100
>>>> message:
>>>>  Redefine grad().T --> grad()
>>>
>>> Have you updated div() too?
>> No, I forgot that. I'll take a look.
>>
>> Are there any other operators that do something strange and should be fixed?
> 
> Strange... UFL seemed to already to compute the correct div:
> 
>   def div(self, o, a):
>         i = Index()
>         g = a[i] if a.rank() == 1 else a[...,i]
>         return g.dx(i)

This looks correct. Someone who uses curl should look at curl as well.

Harish