Re: [firedrake] Integral is missing an integration domain.

On 06/04/16 17:14, Fryderyk Wilczynski wrote:

Thank you, that conditional worked.

I have realized that I should in fact have trial functions, (equations are linear in (n+1) terms), I have been mistakenly following example scripts for nonlinear problems.

I have changed

v, p, phi, mu = split(u)

rho = Function(P1)

to

v, p, phi, mu = TrialFunctions(W)

rho = TrialFunction(P1)

You have a term in a3: inner(chi, phi*grad(mu))*dx but phi and mu are both trial functions. Did you mean phi0? Lawrence