spectral_connectivity.transforms.tridisolve#

tridisolve(d: ndarray[tuple[int, ...], dtype[floating]], e: ndarray[tuple[int, ...], dtype[floating]], b: ndarray[tuple[int, ...], dtype[floating]], overwrite_b: bool = True) → ndarray[tuple[int, ...], dtype[floating]][source]#

Symmetric tridiagonal system solver, from Golub and Van Loan p157.

Note

Copied from NiTime.

Parameters:

  • d (ndarray) – main diagonal stored in d[:]

  • e (ndarray) – superdiagonal stored in e[:-1]

  • b (ndarray) – RHS vector

Returns:

x – Solution to Ax = b (if overwrite_b is False). Otherwise solution is stored in previous RHS vector b

Return type:

ndarray