dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Fri, Dec 14, 2007 at 12:40:13PM +0100, Anders Logg wrote:
> On Thu, Dec 13, 2007 at 05:42:49PM +0100, Anne Voigt wrote:
> > I have the following problem. I defined a trilinearform with the help of
> > FFC. Now I wanted to write the results into a 3D tensor which is
> > assembled with the help of assemble.h . What I already did I implemented
> > a subclass of GenericTensor called Tensor. But unfortunately I am still
> > not able to assemble such a tenor. The reason is
> > SparsityPatternBuilder.cpp which is not implemented for problems >2D:
> > 
> > void SparsityPatternBuilder::build(SparsityPattern& sparsity_pattern,
> > Mesh& mesh,
> >                                    UFC& ufc)
> > {
> >   if (ufc.form.rank() == 0)
> >     scalarBuild(sparsity_pattern);
> >   else if (ufc.form.rank() == 1)
> >     vectorBuild(sparsity_pattern, ufc);
> >   else if (ufc.form.rank() == 2)
> >     matrixBuild(sparsity_pattern, mesh, ufc);
> >   else
> >     error("Cannot compute sparsity patterm for size > 2.");
> > }
> > 
> > Is there anybody who can help me with this problem???
> > I would be in big trouble if not so please try to help!!!
> > 
> > Anne

Excellent question... It seems we need to fix this. Ola and I are
working on it.

-- 
Anders