dolfin team mailing list archive

["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...

-- 
Martin