used solve_banded for SCIPYsolve