метод фиктивных областей для решения уравнения Пуассона

[Tuzembobel]

Здравствуйте! Хотел задаться вопросом: допустим у нас есть трехмерная задача Пуассона:
$\Delta u = -f$ в области $V$, $u\bigl|_{\partial V}=0$.
Хотим применить быстрое дискретное преобразование Фурье для его решения. Но дискретное преобразование Фурье можно выполнять в кубе (везде где встречал, он выполнялся в прямоугольных областях). Пусть $V_1$ - дополнение нашей области $V$ до куба $K$.
Интересует численное (приближенное решение), тогда исходная задача

$$ \begin{cases} 1\cdot \frac{\partial^2u}{\partial x^2} + 1\cdot\frac{\partial^2u}{\partial y^2} + 1\cdot\frac{\partial^2u}{\partial z^2} = -f, & x\in V \\ u(x,y,z)=0, & x\in \partial V \end{cases} $$

сведется к

$$ \begin{cases} a\cdot \frac{\partial^2u}{\partial x^2} + a\cdot\frac{\partial^2u}{\partial y^2} + a\cdot\frac{\partial^2u}{\partial z^2} = -f, & x\in K \\ u(x,y,z)=0, & x\in \partial K \\ \end{cases} $$

плюс добавляется условие, что функция

$$ \sum_{i=1}^{3}a\cdot \cos (\vec{n}, x_i)\frac{\partial u}{\partial x_i} $$

имеет нулевой скачок на $S$ - общей части границ $V$ и $V_1$.

Вопрос в следующем: как численно задать условие того, что у вышеуказанной функции нулевой скачок на $S$?