ffc team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

On Tuesday 26 January 2010 07:28:33 Kristian Oelgaard wrote:
> 2010/1/26 Garth N. Wells <gnw20@xxxxxxxxx>:
> > Kristian Oelgaard wrote:
> >> 2010/1/26 Anders Logg <logg@xxxxxxxxx>:
> >>> I'm just about to take off for the Depeche Mode concert in Göteborg
> >>> but let's discuss it tomorrow. Or feel free to figure out what
> >>> tabulate should do... :-)
> >>
> >> OK enjoy, I might have  a look although I have an important Handball
> >> match to watch tonight,
> >
> > Are you joking?
> 
> Hmm, I guess handball is still not very popular (or even known?) in Great
>  Britain. Denmark is playing Norway and only the winner has a chance of
>  qualifying for the semi-finals in the mens European championships.

*laugh*

God kamp!

Johan

> Kristian
> 
> > Garth
> >
> >> but I guess you don't follow the championships
> >> anymore in Sweden? :)
> >>
> >> Kristian
> >>
> >>> --
> >>> Anders
> >>>
> >>> On Tue, Jan 26, 2010 at 04:08:07PM +0100, k.b.oelgaard@xxxxxxxxx wrote:
> >>>> MixedElement.tabulate() produces some mysterious results. Try the
> >>>> following in a ufl file:
> >>>>
> >>>> from ffc.fiatinterface import create_element
> >>>> elem = FiniteElement("Lagrange", triangle, 1)
> >>>> print create_element(elem).tabulate(1, [(0.333, 0.333)])
> >>>>
> >>>> This gives:
> >>>>
> >>>> {(0, 1): array([[-1.0],
> >>>>       [0.0],
> >>>>       [1.0]], dtype=object), (1, 0): array([[-1.0],
> >>>>       [1.0],
> >>>>       [0.0]], dtype=object), (0, 0): array([[0.334],
> >>>>       [0.333],
> >>>>       [0.333]], dtype=object)}
> >>>>
> >>>> which is correct.
> >>>>
> >>>> For a vector (mixed) element
> >>>>
> >>>> elem = VectorElement("Lagrange", triangle, 1)
> >>>>
> >>>> the result is
> >>>>
> >>>> {(0, 1): array([[[ 0.334],                                  [ 0.
> >>>> ]],
> >>>>
> >>>>       [[ 0.333],
> >>>>        [ 0.   ]],
> >>>>
> >>>>       [[ 0.333],
> >>>>        [ 0.   ]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.334]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.333]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.333]]]), (1, 0): array([[[ 0.334],
> >>>>        [ 0.   ]],
> >>>>
> >>>>       [[ 0.333],
> >>>>        [ 0.   ]],
> >>>>
> >>>>       [[ 0.333],
> >>>>        [ 0.   ]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.334]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.333]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.333]]]), (0, 0): array([[[ 0.334],
> >>>>        [ 0.   ]],
> >>>>
> >>>>       [[ 0.333],
> >>>>        [ 0.   ]],
> >>>>
> >>>>       [[ 0.333],
> >>>>        [ 0.   ]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.334]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.333]],
> >>>>
> >>>>       [[ 0.   ],
> >>>>        [ 0.333]]])}
> >>>>
> >>>> which is very wrong. I had a look at the code in
> >>>> MixedElement.tabulate(), but couldn't figure out what was going on.
> >>>>
> >>>> Kristian
> >>>>
> >>>>
> >>>>
> >>>>
> >>>> _______________________________________________
> >>>> Mailing list: https://launchpad.net/~ffc
> >>>> Post to     : ffc@xxxxxxxxxxxxxxxxxxx
> >>>> Unsubscribe : https://launchpad.net/~ffc
> >>>> More help   : https://help.launchpad.net/ListHelp
> >>>
> >>> -----BEGIN PGP SIGNATURE-----
> >>> Version: GnuPG v1.4.9 (GNU/Linux)
> >>>
> >>> iEYEARECAAYFAktfBtEACgkQTuwUCDsYZdE/DACfYRm5wWiLHtVoRDMmw9syr9vA
> >>> JnYAnjJHHAtl0oST5d9aMjJ0acbtb7rE
> >>> =qBbe
> >>> -----END PGP SIGNATURE-----
> >>
> >> ------------------------------------------------------------------------
> >>
> >> _______________________________________________
> >> Mailing list: https://launchpad.net/~ffc
> >> Post to     : ffc@xxxxxxxxxxxxxxxxxxx
> >> Unsubscribe : https://launchpad.net/~ffc
> >> More help   : https://help.launchpad.net/ListHelp
>