ufl team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Wed, Mar 04, 2009 at 02:27:09PM +0100, Martin Sandve Alnæs wrote:
> On Wed, Mar 4, 2009 at 2:12 PM, Anders Logg <logg@xxxxxxxxx> wrote:
> > I get an error when trying purge_list_tensors.
> >
> > Here's a simple example:
> >
> >  element = VectorElement("Lagrange", "triangle", 1)
> >
> >  v = TestFunction(element)
> >  u = TrialFunction(element)
> >
> >  a = inner(as_tensor([v[0], v[1]]), as_tensor([u[0], u[1]]))
> >
> > This seems to work if I first call expand_compounds, but not without.
> > Is expand_compounds a requirement before calling purge_list_tensors?
> 
> Yes. Actually, I assumed expand_derivatives had been applied.

ok. What about the second example, should that work?

  b = inner(as_tensor([v[0], v[1]]), as_tensor([u[0], u[1].dx(0)]))

-- 
Anders

Attachment: signature.asc
Description: Digital signature