“…Neumann BCs can be enforced both at the stencil level, i.e., when the local interpolants exactly satisfy the BCs at the boundary nodes [16,17,18,10,19], or at the assembly level, i.e., when the equations for the BCs at the boundary nodes appear explicitly in the final sparse matrix [20,21]. In both cases the RBF-FD collocation method, however, may lead to large errors and/or ill-conditioning problems.…”