dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #06013
Re: assemebling process for a trilinearform
Anders Logg skrev den 18/01-2008 følgende:
> 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?
Sure, it did. I have to compile DOLFIN using PETSc for all the part to
fit together. I'll probably commit early next week.
Ola
> --
> Anders
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
Follow ups
References