Implementing tridiagonal matrix algorithm (TDMA) with NumPy
If its tridiagonal systems you want to solve there is solve_banded() in numpy.linalg. Not sure if that's what you're looking for.
Apparently, there is no way to do this in Python without using C or its pythonic variations.