def differencing_by_DD(aflx, bflx):
    """Differencing scheme based on unweighted diamond difference; aflx and
    bflx are the cell-averaged and the cell-edge fluxes, respectively."""
    return 2 * aflx - bflx