def to_netflux(flux):
    r"""Compute the netflux.

    f_ij^{+}=max{0, f_ij-f_ji}
    for all pairs i,j

    Parameters
    ----------
    flux : (M, M) scipy.sparse matrix
        Matrix of flux values between pairs of states.

    Returns
    -------
    netflux : (M, M) scipy.sparse matrix
        Matrix of netflux values between pairs of states.

    """
    netflux = flux - flux.T

    """Set negative entries to zero"""
    netflux = remove_negative_entries(netflux)
    return netflux
r"""Compute the netflux.

    f_ij^{+}=max{0, f_ij-f_ji}
    for all pairs i,j

    Parameters
    ----------
    flux : (M, M) scipy.sparse matrix
        Matrix of flux values between pairs of states.

    Returns
    -------
    netflux : (M, M) scipy.sparse matrix
        Matrix of netflux values between pairs of states.