dolfin team mailing list archive

["Martin Sandve Alnæs" <martinal@xxxxxxxxx>]

2007/9/14, Martin Sandve Alnæs <martinal@xxxxxxxxx>:
> I only needed to do
>   coefficients = ArrayFunctionPtr()
> instead of
>   coefficients =[]
>
> We should have a typemap for this, and more assemble(...) signatures
> taking ufc::form so there's no need to do the whole
>
> asm.assemble(A, form, coefficients, cell_domains,
> exterior_facet_domains, interior_facet_domains, reset_tensor)
>
> But the resulting matrix is completely wrong:
> | (0, 0, 5.56e-02) (0, 1, 0.00e+00) |
> | (1, 0, 0.00e+00) (1, 1, 0.00e+00) (1, 2, 0.00e+00) |
> | (2, 1, 0.00e+00) (2, 2, 3.47e-18) (2, 3, 0.00e+00) |
> | (3, 2, 0.00e+00) (3, 3, -5.56e-02) |
>
> Guess I need to do some debugging...


The bugs where all mine, now it works!

We plan to release a new version of SyFi including SFC (the SyFi Form
Compiler ) soon, and will do some benchmarking of FFC vs SFC as well.

-- 
Martin