Boundary Integral Equation Method For Electrostatic Field Prediction In Piecewise-Homogeneous Electrolytes