Hello SfePy,

May be I am wrong, but it seems that there is an issue with the minus sign before volume integral in Poisson problem definition.

I am trying to solve Poisson problem for the simple analytical function $u(x) = x^2$, so the $\Delta u = 2$ and the week form is $\int_\Omega u v \, d\Omega = \int_\Omega 2 v \, d\Omega$.

In SfePy I used the following definition:
dw_laplace.i1.Omega(m.val, v, u) = dw_volume_integrate.i1.Omega(f.val, v)
but I got a quite big relative error: 1e-1

But in case I add the minus sign:
dw_laplace.i1.Omega(m.val, v, u) = -dw_volume_integrate.i1.Omega(f.val, v)
the accuracy of the solution becomes very good, the relative error: 1e-3

So my question, why it is necessary to add minus sign before volume integral? Yes, I know that the classical Poisson problem is $u(x) = -f(x)$. May be we assume minus sign implicitly?

The test script demonstrated a problem is in the attachment.

Sincerely,
Alec