Optimization and extension of pentadiagonal band(J) solver to multiregion systems containing interior boundaries