def METHOD_NAME(
    rng: np.random.Generator,
) -> Tuple[pn.linops.LinearOperator, np.ndarray]:
    matrix = scipy.sparse.rand(
        10, 10, density=0.1, format="coo", dtype=np.double, random_state=rng
    )
    matrix.setdiag(2)
    matrix = matrix.tocsr()
    return pn.linops.Matrix(matrix), matrix.toarray()