def METHOD_NAME():
    mode = get_default_mode()
    mode = mode.including("specialize", "local_sampling_dot_csr")
    for sp_format in ["csr"]:  # Not implemented for other format
        inputs = [
            matrix(),
            matrix(),
            getattr(aesara.sparse, sp_format + "_matrix")(),
        ]
        f = aesara.function(inputs, sparse.sampling_dot(*inputs), mode=mode)
        if aesara.config.blas__ldflags:
            assert not any(
                isinstance(node.op, sparse.SamplingDot)
                for node in f.maker.fgraph.toposort()
            )
        else:
            # SamplingDotCSR's C implementation needs blas, so it should not
            # be inserted
            assert not any(
                isinstance(node.op, SamplingDotCSR)
                for node in f.maker.fgraph.toposort()
            )