dolfin team mailing list archive

[logg@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.

Excellent!

/Anders