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