ufl team mailing list archive
-
ufl team
-
Mailing list archive
-
Message #00907
Re: Bug in purge_list_tensors?
Yes. Now it does.
Martin
On Wed, Mar 4, 2009 at 2:36 PM, Anders Logg <logg@xxxxxxxxx> wrote:
> 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
>
> -----BEGIN PGP SIGNATURE-----
> Version: GnuPG v1.4.9 (GNU/Linux)
>
> iEYEARECAAYFAkmug9MACgkQTuwUCDsYZdGqhwCePAZefk8o9kR01fFlGNlXscnW
> wTkAnj16YbCStRHfDqhXxk1VcJ146eEW
> =8iKN
> -----END PGP SIGNATURE-----
>
>
Follow ups
References