obtain integrand from integral

Question

I try to generalize some methods. In one step I need to modify the integrand in a integral.

Can I obtain integrand from a integral Form in FEniCS?

v = (rho**2) *dx
integrand = v.??
new_v = integrand*sin(rho)*dx

Thx!

Answer 1

Hi, consider

from dolfin import *

mesh = UnitSquareMesh(10, 10)

f = Constant(10)
v = f*ds(domain=mesh) + (f**2)*ds(domain=mesh) + (f**3)*ds(domain=mesh)

for integral in v.integrals():
    integrand = integral.integrand()
    form = integrand*dx(domain=mesh)

    print form, assemble(form)