I am now trying to build a generic assembling template to integrate a specific source term in a boundary. It doesn't seem to be complicated but Iam missing something..
Here are the variables:
U: 2d vector field to be pushed as data, (mf_u : FEM_PK(2,2), mf_u.qdim()=2). N: unit vector normal to the boundary (extern normal) (a 2d vector?) dp: scalar field (mf_p : FEM_PK(2,2), mf_p.qdim()=1).
The source term to be integrated over the boundaries:
f_i=(U.N).dp_i
And here is the piece of code I am trying to make:
============================================ | An error has been detected !!! | ============================================ Error in getfem_mat_elem.cc, line 121 : Exact integration not allowed in this context
And then, changing the integration method (mim) from "IM_EXACT_SIMPLEX(2)" to "IM_NC(2,6)" (newton-cotes of degree 6 on triangles) I started havin this new message: