ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

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.

/Anders