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