def LU(data):
        """LU-decomposition using either Cython or Numba

        Parameters
        ----------
        data : 2D-array
            Storage for dia-matrix on entry and L and U matrices
            on exit.
        """
        raise NotImplementedError