← Back to team overview

ffc team mailing list archive

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