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.
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)