← Back to team overview

ffc team mailing list archive

Re: outer products

 

> 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
>
> /Anders
> _______________________________________________
> FFC-dev mailing list
> FFC-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/ffc-dev
>




Follow ups

References