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. -- Martin _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
