def METHOD_NAME(dtype=dtype):
    def I(Lmax,m,s):
        n = spin2Jacobi(Lmax,m,s)[0]
        N = np.ones(n,dtype=dtype)
        return infinite_csr(banded((N,[0]),(max(n,0),max(n,0))))
    return Operator(I,SphereCodomain(0,0,0,0))