ffc team mailing list archive
-
ffc team
-
Mailing list archive
-
Message #00680
Re: outer products
On Tue, Sep 26, 2006 at 07:33:27PM +0200, Dag Lindbo wrote:
> > On Tue, Sep 26, 2006 at 06:00:46PM +0200, Johan Jansson wrote:
> >> On Tue, Sep 26, 2006 at 05:53:51PM +0200, Anders Logg wrote:
> >>
> >> ...
> >>
> >> > Good point. The following example seems to work:
> >> >
> >> > element = FiniteElement("Vector Lagrange", "triangle", 1)
> >> > v = BasisFunction(element)
> >> > print mult(transp([vec(v)]), [vec(v)])
> >> > print mult([vec(v)], transp([vec(v)]))
> >> >
> >> > Should we add a new operator mat() that returns [vec()] or should we
> >> > make vec() return this directly so it works like a column vector?
> >> >
> >> >
> >> > /Anders
> >>
> >> The second alternative is probably best. I think mat() should be
> >> reserved for matrix-valued functions, to perform the equivalent of
> >> vec().
> >>
> >> Johan
> >
> > Sounds good, but this will require some work. A number of other
> > operators in operators.py will need to be updated correspondingly.
> >
> > If anyone is willing to try, you're more than welcome. If not, I'll
> > wait until the reimplementation (and extension) of the form language.
>
> Tomrrow morning I'll write a simple operator 'outer(vec(n),vec(n))' that
> might (if I'm sucessful) be enough until the extension of the language is
> complete.
>
> /Dag
ok.
/Anders
Follow ups
References