> 2007/9/14, Martin Sandve Alnæs <[EMAIL PROTECTED]>: >> 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 _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
