dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

This is not the ordering used by FFC. FFC uses another ordering where
we first take all dofs for one component, then all dofs for the next
and so on. (This gives a block diagonal matrix for a system of
Poisson's equation.)

FFC uses the ordering

      j*no_scalar_dofs + i

Ideally, DOLFIN should not need to worry about this but instead call
the proper function from the element generated by FFC. If there is no
such function, just assume that there is (add it to the file generated
by FFC) and I can put the functionality into FFC later.

/Anders

On Tue, Feb 15, 2005 at 11:30:49AM +0100, dolfin@xxxxxxxxxx wrote:
> Commit from hoffman (2005-02-15 11:30 CET)
> 
> Added vector function interface to NewFunction. The internal ordering of
> in particular project() needs to be compatbile with that of FFC. The
> present ordering follows the rule: value for vector component j at dof i,
> is stored at i*no_comp+j (with no_comp the number of vector components).
> 
>   dolfin  src/demo/la/Makefile                      1.193
>   dolfin  src/demo/solvers/ode/stiff/Makefile       1.148
>   dolfin  src/kernel/fem/Makefile                   1.200
>   dolfin  src/kernel/function/NewFunction.cpp       1.9
>   dolfin  src/kernel/function/dolfin/NewFunction.h  1.7
> 
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/cgi-bin/mailman/listinfo/dolfin-dev
>