← Back to team overview

ffc team mailing list archive

Re: define a "diagonal" form

 

I think this is already handled. The code generated by FFC is optimized so that things that are known to be zero a priori are never computed.

Try the simplest possible example you can think of and look at the generated code (the function tabulate_tensor) and see if this is correct. (And let us know what you find.)

/Anders


Alessio Quaglino wrote:
I'm wondering if it's possible to tell FFC that he doesn't have to test
all the basis functions against each other, but only against themself
once, so that I get only the diagonal terms of the matrix. I guess this
would speedup *a lot* the assembly in the case of piecewise elements
having support on only one tetrahedra, or the case when you need only
diagonal elements. Am I right or this special case is already handled? Can
I use a special notation to achieve this aim? Thanks.

Alessio

_______________________________________________
FFC-dev mailing list
FFC-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/ffc-dev


References