← Back to team overview

ufl team mailing list archive

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