ufl team mailing list archive
-
ufl team
-
Mailing list archive
-
Message #01559
Re: [Branch ~ufl-core/ufl/main] Rev 860: Redefine grad().T --> grad()
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)
But I spotted what I think is a bug in the Div class. I pushed a fix
for that but I don't think it has any consequences.
--
Anders
Attachment:
signature.asc
Description: Digital signature
Follow ups
References