specify intrule for edge-wise calculations

Main Forum Old-Forum - do not post here
1

Hello,

I am having trouble using the intrule flag. I want to calculate the edge-wise integral using a higher order quadrature rule. The last line of the following code gives me a segmentation fault:

V = L2(mesh, order=1,dgjumps=True)
u,v = V.TnT()
a = BilinearForm(V)
ir = IntegrationRule(SEGM, 5)
a += SymbolicBFI(u*v, skeleton=True, intrule = ir)

Best,
Guosheng

2

you can use the other (now recommended) version:

a += SymbolicBFI(u*v, skeleton=True).SetIntegrationRule(SEGM, ir)

the faulting version is also fixed

3

Thanks!