Dears, I need your help with the inner boundary condition setting.

For simplicity, we take a 2-D Poisson equation for example, which correspond to the “s_Poisson_2D_f.f90” in the MKL library.

Case A

In case A, it is clear that we can assign the following array to set the boundary condition.

For example in s_Poisson_2D_f.f90

bd_ax(iy) = 1.0E0,

bd_bx(iy) = 1.0E0 ,

bd_ay(ix) = -2.0*pi*sin(2*pi*(ix-1)/nx),

bd_by(ix) = 2.0*pi*sin(2*pi*(ix-1)/nx)

Case B

In case B , in the analysis volume, there are some other charge sources in the analysis volume. The Potential of the black box is set to 0. And the line charge denity of the long wire is given.

So how to solve the problem of case B with “s_Poisson_2D_f.f90” in the MKL library?