def Solve(u, data, axis=0):
        """Fast solve using either Cython or Numba

        Parameters
        ----------
        u : array
            rhs on entry, solution on exit
        data : 2D-array
            Storage for dia-matrix containing L and U matrices
        axis : int, optional
            The axis we are solving over
        """
        raise NotImplementedError