dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Fri, Jan 18, 2008 at 03:56:46PM +0100, Anne Voigt wrote:
> I just wanted to ask if there is any progress with the assembling
> process of a trilinearform.
> As a little reminder, here is what I wrote in December:
> 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.");
>  }
>  
> I know that you all are very busy that's why I would be already pleased
> if you would be able to tell how much longer I may have to wait?!
> 
> Thanks for your help
> Anne
> 
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev

We're working on this. The hangup right now is to solve a problem
related to circular dependencies between matrix classes and matrix
factories (in particular for uBLAS where things are templated).

We may have found a solution today (in particular Martin who has a
black belt in templates). Ola, did it work?

-- 
Anders