Hi Guosheng,
Adding integration rules with the dx syntax was missing, I just added it.
You can now provide an integration rule for every element shape like that:
Code:
a += u*v*dx(intrules={TRIG:ir})
Actually, it is also possible with the current code by this workaround:
Code:
a.integrators[0].SetIntegrationRule(TRIG,ir)
By adding a dx-integral to a BilinearForm, it is internally converted to a SymbolicBFI.
best, Joachim