GitLab is used only for code review, issue tracking and project management. Canonical locations for source code are still https://gitweb.torproject.org/ https://git.torproject.org/ and git-rw.torproject.org.

v3bwfile.py 62.1 KB
Newer Older
juga  's avatar
juga committed
1 2 3 4
# -*- coding: utf-8 -*-
"""Classes and functions that create the bandwidth measurements document
(v3bw) used by bandwidth authorities."""

juga  's avatar
juga committed
5
import copy
juga  's avatar
juga committed
6
import logging
7
import math
juga  's avatar
juga committed
8
import os
9
from itertools import combinations
juga  's avatar
juga committed
10
from statistics import median, mean
11
from stem.descriptor import parse_file
juga  's avatar
juga committed
12

13
from sbws import __version__
14
from sbws.globals import (SPEC_VERSION, BW_LINE_SIZE, SBWS_SCALE_CONSTANT,
15 16
                          TORFLOW_SCALING, SBWS_SCALING, TORFLOW_BW_MARGIN,
                          TORFLOW_OBS_LAST, TORFLOW_OBS_MEAN,
17
                          PROP276_ROUND_DIG, MIN_REPORT, MAX_BW_DIFF_PERC)
juga  's avatar
juga committed
18
from sbws.lib.resultdump import ResultSuccess, _ResultType
19
from sbws.util.filelock import DirectoryLock
juga  's avatar
juga committed
20
from sbws.util.timestamp import (now_isodt_str, unixts_to_isodt_str,
21 22
                                 now_unixts, isostr_to_dt_obj,
                                 dt_obj_to_isodt_str)
23
from sbws.util.state import State
juga  's avatar
juga committed
24 25 26

log = logging.getLogger(__name__)

27
LINE_SEP = '\n'
28 29
KEYVALUE_SEP_V1 = '='
KEYVALUE_SEP_V2 = ' '
30 31 32 33 34 35

# NOTE: in a future refactor make make all the KeyValues be a dictionary
# with their type, so that it's more similar to stem parser.

# Header KeyValues
# =================
36 37
# KeyValues that need to be in a specific order in the Bandwidth File.
HEADER_KEYS_V1_1_ORDERED = ['version']
38 39 40 41 42 43 44 45 46
# KeyValues that are not initialized from the state file nor the measurements.
# They can also be pass as an argument to `Header` to overwrite default values,
# what is done in unit tests.
# `latest bandwidth` is special cause it gets its value from timestamp, which
# is not a KeyValue, but it's always pass as an agument.
# It could be separaed in other list, but so far there is no need, cause:
# 1. when it's pass to the Header to initialize it, it's just ignored.
# 2. when the file is created, it's took into account.
HEADER_KEYS_V1_1_SELF_INITIALIZED = [
47 48 49
    "software",
    "software_version",
    "file_created",
50 51 52 53
    "latest_bandwidth",
]
# KeyValues that are initialized from arguments.
HEADER_KEYS_V1_1_TO_INIT = [
54 55
    "earliest_bandwidth",
    "generator_started",
56 57
]

58 59 60 61
# number_eligible_relays is the number that ends in the bandwidth file
# ie, have not been excluded by one of the filters in 4. below
# They should be call recent_measurement_included_count to be congruent
# with the other KeyValues.
62 63 64 65 66 67 68
HEADER_KEYS_V1_2 = [
    "number_eligible_relays",
    "minimum_number_eligible_relays",
    "number_consensus_relays",
    "percent_eligible_relays",
    "minimum_percent_eligible_relays",
]
69

70 71 72 73 74 75
# KeyValues added in the Bandwidth File v1.3.0
HEADER_KEYS_V1_3 = [
    "scanner_country",
    "destinations_countries",
]

76 77 78
# KeyValues that count the number of relays that are in the bandwidth file,
# but ignored by Tor when voting, because they do not have a
# measured bandwidth.
79
HEADER_RECENT_MEASUREMENTS_EXCLUDED_KEYS = [
80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96
    # Number of relays that were measured but all the measurements failed
    # because of network failures or it was
    # not found a suitable helper relay
    'recent_measurements_excluded_error_count',
    # Number of relays that have successful measurements but the measurements
    # were not away from each other in X time (by default 1 day).
    'recent_measurements_excluded_near_count',
    # Number of relays that have successful measurements and they are away from
    # each other but they are not X time recent.
    # By default this is 5 days, which is the same time the older
    # the measurements can be by default.
    'recent_measurements_excluded_old_count',
    # Number of relays that have successful measurements and they are away from
    # each other and recent
    # but the number of measurements are less than X (by default 2).
    'recent_measurements_excluded_few_count',
]
97
# Added in #29591
98 99 100 101
# NOTE: recent_consensus_count, recent_priority_list_count,
# recent_measurement_attempt_count and recent_priority_relay_count
# are not reset when the scanner is stop.
# They will accumulate the values since the scanner was ever started.
102
HEADER_KEYS_V1_4 = [
103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120
    # 1.1 header: the number of different consensuses, that sbws has seen,
    # since the last 5 days
    'recent_consensus_count',
    # 2.4 Number of times a priority list has been created
    'recent_priority_list_count',
    # 2.5 Number of relays that there were in a priority list
    # [50, number of relays in the network * 0.05]
    'recent_priority_relay_count',
    # 3.6 header: the number of times that sbws has tried to measure any relay,
    # since the last 5 days
    # This would be the number of times a relays were in a priority list
    'recent_measurement_attempt_count',
    # 3.7 header: the number of times that sbws has tried to measure any relay,
    # since the last 5 days, but it didn't work
    # This should be the number of attempts - number of ResultSuccess -
    # something else we don't know yet
    # So far is the number of ResultError
    'recent_measurement_failure_count',
121 122
    # The time it took to report about half of the network.
    'time_to_report_half_network',
123
] + HEADER_RECENT_MEASUREMENTS_EXCLUDED_KEYS
124

125 126 127 128 129 130 131 132
# KeyValues added in the Bandwidth File v1.5.0
# XXX: Change SPEC_VERSION when all the v1.5.0 keys are added, before a new
# sbws release.
# Tor version will be obtained from the state file, so it won't be pass as an
# argument, but will be self-initialized.
HEADER_KEYS_V1_5_TO_INIT = ['tor_version']
HEADER_KEYS_V1_5 = HEADER_KEYS_V1_5_TO_INIT

133 134 135 136 137
# KeyValues that are initialized from arguments, not self-initialized.
HEADER_INIT_KEYS = (
    HEADER_KEYS_V1_1_TO_INIT
    + HEADER_KEYS_V1_3
    + HEADER_KEYS_V1_2
138
    + HEADER_KEYS_V1_4
139
    + HEADER_KEYS_V1_5_TO_INIT
140
)
141

142
HEADER_INT_KEYS = HEADER_KEYS_V1_2 + HEADER_KEYS_V1_4
143
# List of all unordered KeyValues currently being used to generate the file
144 145 146 147 148 149
HEADER_UNORDERED_KEYS = (
    HEADER_KEYS_V1_1_SELF_INITIALIZED
    + HEADER_KEYS_V1_1_TO_INIT
    + HEADER_KEYS_V1_3
    + HEADER_KEYS_V1_2
    + HEADER_KEYS_V1_4
150
    + HEADER_KEYS_V1_5
151
)
152
# List of all the KeyValues currently being used to generate the file
153
HEADER_ALL_KEYS = HEADER_KEYS_V1_1_ORDERED + HEADER_UNORDERED_KEYS
154

155
TERMINATOR = '====='
156 157 158

# Bandwidth Lines KeyValues
# =========================
159
# Num header lines in v1.X.X using all the KeyValues
160
NUM_LINES_HEADER_V1 = len(HEADER_ALL_KEYS) + 2
161 162
LINE_TERMINATOR = TERMINATOR + LINE_SEP

juga  's avatar
juga committed
163
# KeyValue separator in Bandwidth Lines
164
BWLINE_KEYVALUES_SEP_V1 = ' '
165
# not inclding in the files the extra bws for now
166
BWLINE_KEYS_V0 = ['node_id', 'bw']
167 168 169 170 171 172 173 174 175 176 177 178 179
BWLINE_KEYS_V1_1 = [
    "master_key_ed25519",
    "nick",
    "rtt",
    "time",
    "success",
    "error_stream",
    "error_circ",
    "error_misc",
    # Added in #292951
    "error_second_relay",
    "error_destination",
]
180 181 182 183 184 185 186 187 188 189
BWLINE_KEYS_V1_2 = [
    "bw_median",
    "bw_mean",
    "desc_bw_avg",
    "desc_bw_bur",
    "desc_bw_obs_last",
    "desc_bw_obs_mean",
    "consensus_bandwidth",
    "consensus_bandwidth_is_unmeasured",
]
190

191 192
# There were no bandwidth lines key added in the specification version 1.3

193
# Added in #292951
194
BWLINE_KEYS_V1_4 = [
195 196 197
    # 1.2 relay: the number of different consensuses, that sbws has seen,
    # since the last 5 days, that have this relay
    'relay_in_recent_consensus_count',
198 199 200
    # 2.6 relay: the number of times a relay was "prioritized" to be measured
    # in the recent days (by default 5).
    'relay_recent_priority_list_count',
201 202 203 204 205 206 207 208 209 210 211 212
    # 3.8 relay:  the number of times that sbws has tried to measure
    # this relay, since the last 5 days
    # This would be the number of times a relay was in a priority list (2.6)
    # since once it gets measured, it either returns ResultError,
    # ResultSuccess or something else happened that we don't know yet
    'relay_recent_measurement_attempt_count',
    # 3.9 relay:  the number of times that sbws has tried to measure
    # this relay, since the last 5 days, but it didn't work
    # This should be the number of attempts - number of ResultSuccess -
    # something else we don't know yet
    # So far is the number of ResultError
    'relay_recent_measurement_failure_count',
213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229
    # Number of error results created in the last 5 days that are excluded.
    # This is the sum of all the errors.
    'relay_recent_measurements_excluded_error_count',
    # The number of successful results, created in the last 5 days,
    # that were excluded by a rule, for this relay.
    # 'relay_recent_measurements_excluded_error_count' would be the
    # sum of the following 3 + the number of error results.

    # The number of successful measurements that are not X time away
    # from each other (by default 1 day).
    'relay_recent_measurements_excluded_near_count',
    # The number of successful measurements that are away from each other
    # but not X time recent (by default 5 days).
    'relay_recent_measurements_excluded_old_count',
    # The number of measurements excluded because they are not at least X
    # (by default 2).
    'relay_recent_measurements_excluded_few_count',
230 231 232 233 234 235 236 237 238 239 240 241 242
    # `vote=0` is used for the relays that were excluded to
    # be reported in the bandwidth file and now they are
    # reported.
    # It tells Tor to do not vote on the relay.
    # `unmeasured=1` is used for the same relays and it is
    # added in case Tor would vote on them in future versions.
    # Maybe these keys should not be included for the relays
    # in which vote=1 and unmeasured=0.
    'vote', 'unmeasured',
    # When there not enough eligible relays (not excluded)
    # under_min_report is 1, `vote` is 0.
    # Added in #29853.
    'under_min_report',
243
]
244
BWLINE_KEYS_V1 = BWLINE_KEYS_V0 + BWLINE_KEYS_V1_1 + BWLINE_KEYS_V1_2 \
245
               + BWLINE_KEYS_V1_4
246 247
# NOTE: tech-debt: assign boolean type to vote and unmeasured,
# when the attributes are defined with a type, as stem does.
248 249 250 251 252 253 254 255 256 257 258 259
BWLINE_INT_KEYS = (
    [
        "bw",
        "rtt",
        "success",
        "error_stream",
        "error_circ",
        "error_misc",
    ]
    + BWLINE_KEYS_V1_2
    + BWLINE_KEYS_V1_4
)
juga  's avatar
juga committed
260 261


262
def round_sig_dig(n, digits=PROP276_ROUND_DIG):
263 264 265 266 267 268 269
    """Round n to 'digits' significant digits in front of the decimal point.
       Results less than or equal to 1 are rounded to 1.
       Returns an integer.

       digits must be greater than 0.
       n must be less than or equal to 2**73, to avoid floating point errors.
       """
270
    digits = int(digits)
271 272 273 274 275 276 277 278 279
    assert digits >= 1
    if n <= 1:
        return 1
    digits_in_n = int(math.log10(n)) + 1
    round_digits = max(digits_in_n - digits, 0)
    rounded_n = round(n, -round_digits)
    return int(rounded_n)


280
def kb_round_x_sig_dig(bw_bs, digits=PROP276_ROUND_DIG):
281 282 283 284 285 286 287 288 289 290 291
    """Convert bw_bs from bytes to kilobytes, and round the result to
       'digits' significant digits.
       Results less than or equal to 1 are rounded up to 1.
       Returns an integer.

       digits must be greater than 0.
       n must be less than or equal to 2**82, to avoid floating point errors.
       """
    # avoid double-rounding by using floating-point
    bw_kb = bw_bs / 1000.0
    return round_sig_dig(bw_kb, digits=digits)
292 293


juga  's avatar
juga committed
294 295 296 297 298 299 300 301 302
def num_results_of_type(results, type_str):
    return len([r for r in results if r.type == type_str])


# Better way to use enums?
def result_type_to_key(type_str):
    return type_str.replace('-', '_')


303
class V3BWHeader(object):
juga  's avatar
juga committed
304 305
    """
    Create a bandwidth measurements (V3bw) header
306
    following bandwidth measurements document spec version 1.X.X.
juga  's avatar
juga committed
307

308
    :param str timestamp: timestamp in Unix Epoch seconds of the most recent
309
        generator result.
juga  's avatar
juga committed
310 311 312
    :param str version: the spec version
    :param str software: the name of the software that generates this
    :param str software_version: the version of the software
313
    :param dict kwargs: extra headers. Currently supported:
juga  's avatar
juga committed
314

315 316 317 318
        - earliest_bandwidth: str, ISO 8601 timestamp in UTC time zone
          when the first bandwidth was obtained
        - generator_started: str, ISO 8601 timestamp in UTC time zone
          when the generator started
juga  's avatar
juga committed
319
    """
320
    def __init__(self, timestamp, **kwargs):
juga  's avatar
juga committed
321 322 323 324
        assert isinstance(timestamp, str)
        for v in kwargs.values():
            assert isinstance(v, str)
        self.timestamp = timestamp
325 326 327 328
        # KeyValues with default value when not given by kwargs
        self.version = kwargs.get('version', SPEC_VERSION)
        self.software = kwargs.get('software', 'sbws')
        self.software_version = kwargs.get('software_version', __version__)
329
        self.file_created = kwargs.get('file_created', now_isodt_str())
juga  's avatar
juga committed
330
        # latest_bandwidth should not be in kwargs, since it MUST be the
331
        # same as timestamp
juga  's avatar
juga committed
332
        self.latest_bandwidth = unixts_to_isodt_str(timestamp)
333
        [setattr(self, k, v) for k, v in kwargs.items()
334
         if k in HEADER_INIT_KEYS]
335

336
    def __str__(self):
337
        if self.version.startswith('1.'):
338 339
            return self.strv1
        return self.strv2
340

juga  's avatar
juga committed
341
    @classmethod
342 343
    def from_results(cls, results, scanner_country=None,
                     destinations_countries=None, state_fpath=''):
juga  's avatar
juga committed
344 345 346
        kwargs = dict()
        latest_bandwidth = cls.latest_bandwidth_from_results(results)
        earliest_bandwidth = cls.earliest_bandwidth_from_results(results)
347
        # NOTE: Blocking, reads file
348
        generator_started = cls.generator_started_from_file(state_fpath)
349
        recent_consensus_count = cls.consensus_count_from_file(state_fpath)
juga  's avatar
juga committed
350
        timestamp = str(latest_bandwidth)
351 352 353 354 355 356 357 358 359

        # XXX: tech-debt: obtain the other values from the state file using
        # this state variable.
        # Store the state as an attribute of the object?
        state = State(state_fpath)
        tor_version = state.get('tor_version', None)
        if tor_version:
            kwargs['tor_version'] = tor_version

juga  's avatar
juga committed
360 361 362 363
        kwargs['latest_bandwidth'] = unixts_to_isodt_str(latest_bandwidth)
        kwargs['earliest_bandwidth'] = unixts_to_isodt_str(earliest_bandwidth)
        if generator_started is not None:
            kwargs['generator_started'] = generator_started
364 365 366
        # To be compatible with older bandwidth files, do not require it.
        if scanner_country is not None:
            kwargs['scanner_country'] = scanner_country
367 368
        if destinations_countries is not None:
            kwargs['destinations_countries'] = destinations_countries
369 370
        if recent_consensus_count is not None:
            kwargs['recent_consensus_count'] = str(recent_consensus_count)
371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399

        recent_measurement_attempt_count = \
            cls.recent_measurement_attempt_count_from_file(state_fpath)
        if recent_measurement_attempt_count is not None:
            kwargs['recent_measurement_attempt_count'] = \
                str(recent_measurement_attempt_count)

        # If it is a failure that is not a ResultError, then
        # failures = attempts - all mesaurements
        # Works only in the case that old measurements files already had
        # measurements count
        if recent_measurement_attempt_count is not None:
            all_measurements = 0
            for result_list in results.values():
                all_measurements += len(result_list)
            measurement_failures = (recent_measurement_attempt_count
                                    - all_measurements)
            kwargs['recent_measurement_failure_count'] = \
                str(measurement_failures)

        priority_lists = cls.recent_priority_list_count_from_file(state_fpath)
        if priority_lists is not None:
            kwargs['recent_priority_list_count'] = str(priority_lists)

        priority_relays = \
            cls.recent_priority_relay_count_from_file(state_fpath)
        if priority_relays is not None:
            kwargs['recent_priority_relay_count'] = str(priority_relays)

juga  's avatar
juga committed
400 401 402
        h = cls(timestamp, **kwargs)
        return h

403
    @classmethod
404
    def from_lines_v1(cls, lines):
405 406
        """
        :param list lines: list of lines to parse
407
        :returns: tuple of V3BWHeader object and non-header lines
408 409 410 411
        """
        assert isinstance(lines, list)
        try:
            index_terminator = lines.index(TERMINATOR)
Matt Traudt's avatar
Matt Traudt committed
412
        except ValueError:
413 414 415
            # is not a bw file or is v100
            log.warn('Terminator is not in lines')
            return None
416
        ts = lines[0]
417
        kwargs = dict([l.split(KEYVALUE_SEP_V1)
418
                       for l in lines[:index_terminator]
419
                       if l.split(KEYVALUE_SEP_V1)[0] in HEADER_ALL_KEYS])
420
        h = cls(ts, **kwargs)
juga  's avatar
juga committed
421
        # last line is new line
422
        return h, lines[index_terminator + 1:-1]
423 424

    @classmethod
425
    def from_text_v1(self, text):
426
        """
427
        :param str text: text to parse
428
        :returns: tuple of V3BWHeader object and non-header lines
429 430
        """
        assert isinstance(text, str)
431
        return self.from_lines_v1(text.split(LINE_SEP))
juga  's avatar
juga committed
432

juga  's avatar
juga committed
433 434 435 436 437 438 439 440 441 442 443
    @classmethod
    def from_lines_v100(cls, lines):
        """
        :param list lines: list of lines to parse
        :returns: tuple of V3BWHeader object and non-header lines
        """
        assert isinstance(lines, list)
        h = cls(lines[0])
        # last line is new line
        return h, lines[1:-1]

juga  's avatar
juga committed
444
    @staticmethod
445
    def generator_started_from_file(state_fpath):
446 447 448 449
        '''
        ISO formatted timestamp for the time when the scanner process most
        recently started.
        '''
450
        state = State(state_fpath)
451
        if 'scanner_started' in state:
452 453
            # From v1.1.0-dev `state` is capable of converting strs to datetime
            return dt_obj_to_isodt_str(state['scanner_started'])
454 455
        else:
            return None
juga  's avatar
juga committed
456

457 458 459
    @staticmethod
    def consensus_count_from_file(state_fpath):
        state = State(state_fpath)
460 461
        if 'recent_consensus_count' in state:
            return state['recent_consensus_count']
462 463 464
        else:
            return None

465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494
    # NOTE: in future refactor store state in the class
    @staticmethod
    def recent_measurement_attempt_count_from_file(state_fpath):
        """
        Returns the number of times any relay was queued to be measured
        in the recent (by default 5) days from the state file.
        """
        state = State(state_fpath)
        return state.get('recent_measurement_attempt_count', None)

    @staticmethod
    def recent_priority_list_count_from_file(state_fpath):
        """
        Returns the number of times
        :meth:`~sbws.lib.relayprioritizer.RelayPrioritizer.best_priority`
        was run
        in the recent (by default 5) days from the state file.
        """
        state = State(state_fpath)
        return state.get('recent_priority_list_count', None)

    @staticmethod
    def recent_priority_relay_count_from_file(state_fpath):
        """
        Returns the number of times any relay was "prioritized" to be measured
        in the recent (by default 5) days from the state file.
        """
        state = State(state_fpath)
        return state.get('recent_priority_relay_count', None)

juga  's avatar
juga committed
495
    @staticmethod
juga  's avatar
juga committed
496
    def latest_bandwidth_from_results(results):
497
        return round(max([r.time for fp in results for r in results[fp]]))
juga  's avatar
juga committed
498 499 500

    @staticmethod
    def earliest_bandwidth_from_results(results):
501
        return round(min([r.time for fp in results for r in results[fp]]))
juga  's avatar
juga committed
502

juga  's avatar
juga committed
503 504 505 506 507
    @property
    def keyvalue_unordered_tuple_ls(self):
        """Return list of KeyValue tuples that do not have specific order."""
        # sort the list to generate determinist headers
        keyvalue_tuple_ls = sorted([(k, v) for k, v in self.__dict__.items()
508
                                    if k in HEADER_UNORDERED_KEYS])
juga  's avatar
juga committed
509 510 511 512 513 514 515 516
        return keyvalue_tuple_ls

    @property
    def keyvalue_tuple_ls(self):
        """Return list of all KeyValue tuples"""
        return [('version', self.version)] + self.keyvalue_unordered_tuple_ls

    @property
517 518 519
    def keyvalue_v1str_ls(self):
        """Return KeyValue list of strings following spec v1.X.X."""
        keyvalues = [self.timestamp] + [KEYVALUE_SEP_V1.join([k, v])
juga  's avatar
juga committed
520 521 522 523
                                        for k, v in self.keyvalue_tuple_ls]
        return keyvalues

    @property
524 525 526
    def strv1(self):
        """Return header string following spec v1.X.X."""
        header_str = LINE_SEP.join(self.keyvalue_v1str_ls) + LINE_SEP + \
juga  's avatar
juga committed
527 528 529 530
            LINE_TERMINATOR
        return header_str

    @property
531 532 533
    def keyvalue_v2_ls(self):
        """Return KeyValue list of strings following spec v2.X.X."""
        keyvalue = [self.timestamp] + [KEYVALUE_SEP_V2.join([k, v])
juga  's avatar
juga committed
534 535 536 537
                                       for k, v in self.keyvalue_tuple_ls]
        return keyvalue

    @property
538 539 540
    def strv2(self):
        """Return header string following spec v2.X.X."""
        header_str = LINE_SEP.join(self.keyvalue_v2_ls) + LINE_SEP + \
juga  's avatar
juga committed
541 542 543 544 545 546
            LINE_TERMINATOR
        return header_str

    @property
    def num_lines(self):
        return len(self.__str__().split(LINE_SEP))
juga  's avatar
juga committed
547

juga  's avatar
juga committed
548 549 550
    def add_stats(self, **kwargs):
        # Using kwargs because attributes might chage.
        [setattr(self, k, str(v)) for k, v in kwargs.items()
551
         if k in HEADER_KEYS_V1_2]
juga  's avatar
juga committed
552

553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605
    def add_time_report_half_network(self):
        """Add to the header the time it took to measure half of the network.

        It is not the time the scanner actually takes on measuring all the
        network, but the ``number_eligible_relays`` that are reported in the
        bandwidth file and directory authorities will vote on.

        This is calculated for half of the network, so that failed or not
        reported relays do not affect too much.

        For instance, if there are 6500 relays in the network, half of the
        network would be 3250. And if there were 4000 eligible relays
        measured in an interval of 3 days, the time to measure half of the
        network would be 3 days * 3250 / 4000.

        Since the elapsed time is calculated from the earliest and the
        latest measurement and a relay might have more than 2 measurements,
        this would give an estimate on how long it would take to measure
        the network including all the valid measurements.

        Log also an estimated on how long it would take with the current
        number of relays included in the bandwidth file.
        """
        # NOTE: in future refactor do not convert attributes to str until
        # writing to the file, so that they do not need to be converted back
        # to do some calculations.
        elapsed_time = (
            (isostr_to_dt_obj(self.latest_bandwidth)
             - isostr_to_dt_obj(self.earliest_bandwidth))
            .total_seconds())

        # This attributes were added later and some tests that
        # do not initialize them would fail.
        eligible_relays = int(getattr(self, 'number_eligible_relays', 0))
        consensus_relays = int(getattr(self, 'number_consensus_relays', 0))
        if not(eligible_relays and consensus_relays):
            return

        half_network = consensus_relays / 2
        # Calculate the time it would take to measure half of the network
        if eligible_relays >= half_network:
            time_half_network = round(
                elapsed_time * half_network / eligible_relays
            )
            self.time_to_report_half_network = str(time_half_network)

        # In any case log an estimated on the time to measure all the network.
        estimated_time = round(
            elapsed_time * consensus_relays / eligible_relays
        )
        log.info("Estimated time to measure the network: %s hours.",
                 round(estimated_time / 60 / 60))

606 607 608 609 610 611 612 613 614
    def add_relays_excluded_counters(self, exclusion_dict):
        """
        Add the monitoring KeyValues to the header about the number of
        relays not included because they were not ``eligible``.
        """
        log.debug("Adding relays excluded counters.")
        for k, v in exclusion_dict.items():
            setattr(self, k, str(v))

juga  's avatar
juga committed
615

juga  's avatar
juga committed
616
class V3BWLine(object):
juga  's avatar
juga committed
617
    """
618
    Create a Bandwidth List line following the spec version 1.X.X.
juga  's avatar
juga committed
619

620 621 622 623
    :param str node_id: the relay fingerprint
    :param int bw: the bandwidth value that directory authorities will include
        in their votes.
    :param dict kwargs: extra headers.
juga  's avatar
juga committed
624

625 626
    .. note:: tech-debt: move node_id and bw to kwargs and just ensure that
       the required values are in **kwargs
juga  's avatar
juga committed
627
    """
juga  's avatar
juga committed
628 629
    def __init__(self, node_id, bw, **kwargs):
        assert isinstance(node_id, str)
630
        assert node_id.startswith('$')
juga  's avatar
juga committed
631 632 633
        self.node_id = node_id
        self.bw = bw
        [setattr(self, k, v) for k, v in kwargs.items()
634
         if k in BWLINE_KEYS_V1]
juga  's avatar
juga committed
635

juga  's avatar
juga committed
636
    def __str__(self):
637
        return self.bw_strv1
juga  's avatar
juga committed
638

juga  's avatar
juga committed
639
    @classmethod
juga  's avatar
juga committed
640 641
    def from_results(cls, results, secs_recent=None, secs_away=None,
                     min_num=0):
642 643 644
        """Convert sbws results to relays' Bandwidth Lines

        ``bs`` stands for Bytes/seconds
645
        ``bw_mean`` means the bw is obtained from the mean of the all the
646 647 648 649 650
        downloads' bandwidth.
        Downloads' bandwidth are calculated as the amount of data received
        divided by the the time it took to received.
        bw = data (Bytes) / time (seconds)
        """
651
        # log.debug("Len success_results %s", len(success_results))
juga  's avatar
juga committed
652 653 654 655 656 657 658
        node_id = '$' + results[0].fingerprint
        kwargs = dict()
        kwargs['nick'] = results[0].nickname
        if getattr(results[0], 'master_key_ed25519'):
            kwargs['master_key_ed25519'] = results[0].master_key_ed25519
        kwargs['time'] = cls.last_time_from_results(results)
        kwargs.update(cls.result_types_from_results(results))
659 660 661
        consensuses_count = \
            [r.relay_in_recent_consensus_count for r in results
             if getattr(r, 'relay_in_recent_consensus_count', None)]
662 663
        if consensuses_count:
            consensus_count = max(consensuses_count)
664
            kwargs['relay_in_recent_consensus_count'] = str(consensus_count)
665

666 667 668 669 670 671 672 673 674 675 676 677 678 679
        measurements_attempts = \
            [r.relay_recent_measurement_attempt_count for r in results
             if getattr(r, 'relay_recent_measurement_attempt_count', None)]
        if measurements_attempts:
            kwargs['relay_recent_measurement_attempt_count'] = \
                str(max(measurements_attempts))

        relay_recent_priority_list_counts = \
            [r.relay_recent_priority_list_count for r in results
             if getattr(r, 'relay_recent_priority_list_count', None)]
        if relay_recent_priority_list_counts:
            kwargs['relay_recent_priority_list_count'] = \
                str(max(relay_recent_priority_list_counts))

680
        success_results = [r for r in results if isinstance(r, ResultSuccess)]
681 682 683 684

        # NOTE: The following 4 conditions exclude relays from the bandwidth
        # file when the measurements does not satisfy some rules, what makes
        # the relay non-`eligible`.
685 686
        # In BWLINE_KEYS_V1_4 it is explained what they mean.
        # In HEADER_RECENT_MEASUREMENTS_EXCLUDED_KEYS it is also
687 688 689 690 691 692
        # explained the what it means the strings returned.
        # They rules were introduced in #28061 and #27338
        # In #28565 we introduce the KeyValues to know why they're excluded.
        # In #28563 we report these relays, but make Tor ignore them.
        # This might confirm #28042.

693 694 695 696 697 698 699 700 701 702 703 704 705 706
        # If the relay is non-`eligible`:
        # Create a bandwidth line with the relay, but set ``vote=0`` so that
        # Tor versions with patch #29806 does not vote on the relay.
        # Set ``bw=1`` so that Tor versions without the patch,
        # will give the relay low bandwidth.
        # Include ``unmeasured=1`` in case Tor would vote on unmeasured relays
        # in future versions.
        # And return because there are not bandwidth values.
        # NOTE: the bandwidth values could still be obtained if:
        # 1. ``ResultError`` will store them
        # 2. assign ``results_recent = results`` when there is a ``exclusion
        # reason.
        # This could be done in a better way as part of a refactor #28684.

juga  's avatar
juga committed
707 708
        kwargs['vote'] = 0
        kwargs['unmeasured'] = 1
709 710

        exclusion_reason = None
711

712 713 714 715 716
        number_excluded_error = len(results) - len(success_results)
        if number_excluded_error > 0:
            # then the number of error results is the number of results
            kwargs['relay_recent_measurements_excluded_error_count'] = \
                number_excluded_error
717
        if not success_results:
718
            exclusion_reason = 'recent_measurements_excluded_error_count'
719
            return (cls(node_id, 1, **kwargs), exclusion_reason)
720

721 722
        results_away = \
            cls.results_away_each_other(success_results, secs_away)
723 724 725
        number_excluded_near = len(success_results) - len(results_away)
        if number_excluded_near > 0:
            kwargs['relay_recent_measurements_excluded_near_count'] = \
726
                number_excluded_near
727
        if not results_away:
728
            exclusion_reason = \
729
                'recent_measurements_excluded_near_count'
730
            return (cls(node_id, 1, **kwargs), exclusion_reason)
731 732
        # log.debug("Results away from each other: %s",
        #           [unixts_to_isodt_str(r.time) for r in results_away])
733

734
        results_recent = cls.results_recent_than(results_away, secs_recent)
735 736 737 738
        number_excluded_old = len(results_away) - len(results_recent)
        if number_excluded_old > 0:
            kwargs['relay_recent_measurements_excluded_old_count'] = \
                number_excluded_old
739
        if not results_recent:
740
            exclusion_reason = \
741
                'recent_measurements_excluded_old_count'
742
            return (cls(node_id, 1, **kwargs), exclusion_reason)
743

744
        if not len(results_recent) >= min_num:
745 746
            kwargs['relay_recent_measurements_excluded_few_count'] = \
                len(results_recent)
747
            # log.debug('The number of results is less than %s', min_num)
748
            exclusion_reason = \
749
                'recent_measurements_excluded_few_count'
750 751 752 753 754 755
            return (cls(node_id, 1, **kwargs), exclusion_reason)

        # For any line not excluded, do not include vote and unmeasured
        # KeyValues
        del kwargs['vote']
        del kwargs['unmeasured']
756

757 758 759 760 761 762 763 764 765 766 767 768 769 770 771
        rtt = cls.rtt_from_results(results_recent)
        if rtt:
            kwargs['rtt'] = rtt
        bw = cls.bw_median_from_results(results_recent)
        kwargs['bw_mean'] = cls.bw_mean_from_results(results_recent)
        kwargs['bw_median'] = cls.bw_median_from_results(
            results_recent)
        kwargs['desc_bw_avg'] = \
            cls.desc_bw_avg_from_results(results_recent)
        kwargs['desc_bw_bur'] = \
            cls.desc_bw_bur_from_results(results_recent)
        kwargs['consensus_bandwidth'] = \
            cls.consensus_bandwidth_from_results(results_recent)
        kwargs['consensus_bandwidth_is_unmeasured'] = \
            cls.consensus_bandwidth_is_unmeasured_from_results(
772
                results_recent)
773 774 775 776 777
        kwargs['desc_bw_obs_last'] = \
            cls.desc_bw_obs_last_from_results(results_recent)
        kwargs['desc_bw_obs_mean'] = \
            cls.desc_bw_obs_mean_from_results(results_recent)
        bwl = cls(node_id, bw, **kwargs)
778
        return bwl, None
juga  's avatar
juga committed
779 780 781 782 783 784

    @classmethod
    def from_data(cls, data, fingerprint):
        assert fingerprint in data
        return cls.from_results(data[fingerprint])

juga  's avatar
juga committed
785
    @classmethod
786
    def from_bw_line_v1(cls, line):
juga  's avatar
juga committed
787
        assert isinstance(line, str)
788
        kwargs = dict([kv.split(KEYVALUE_SEP_V1)
789
                       for kv in line.split(BWLINE_KEYVALUES_SEP_V1)
790
                       if kv.split(KEYVALUE_SEP_V1)[0] in BWLINE_KEYS_V1])
juga  's avatar
juga committed
791
        for k, v in kwargs.items():
792
            if k in BWLINE_INT_KEYS:
juga  's avatar
juga committed
793
                kwargs[k] = int(v)
794 795 796 797 798
        node_id = kwargs['node_id']
        bw = kwargs['bw']
        del kwargs['node_id']
        del kwargs['bw']
        bw_line = cls(node_id, bw, **kwargs)
juga  's avatar
juga committed
799 800
        return bw_line

juga  's avatar
juga committed
801 802
    @staticmethod
    def results_away_each_other(results, secs_away=None):
juga  's avatar
juga committed
803 804
        # log.debug("Checking whether results are away from each other in %s "
        #           "secs.", secs_away)
juga  's avatar
juga committed
805 806
        if secs_away is None or len(results) < 2:
            return results
807 808 809 810 811
        for a, b in combinations(results, 2):
            if abs(a.time - b.time) > secs_away:
                return results
        # log.debug("Results are NOT away from each other in at least %ss: %s",
        #           secs_away, [unixts_to_isodt_str(r.time) for r in results])
812
        return []
juga  's avatar
juga committed
813 814 815 816 817

    @staticmethod
    def results_recent_than(results, secs_recent=None):
        if secs_recent is None:
            return results
818
        results_recent = list(filter(
juga  's avatar
juga committed
819
                            lambda x: (now_unixts() - x.time) < secs_recent,
820 821 822 823 824 825
                            results))
        # if not results_recent:
        #     log.debug("Results are NOT more recent than %ss: %s",
        #               secs_recent,
        #               [unixts_to_isodt_str(r.time) for r in results])
        return results_recent
juga  's avatar
juga committed
826

827
    @staticmethod
828
    def bw_median_from_results(results):
829 830 831 832
        return max(round(median([dl['amount'] / dl['duration']
                                 for r in results for dl in r.downloads])), 1)

    @staticmethod
833
    def bw_mean_from_results(results):
834 835 836
        return max(round(mean([dl['amount'] / dl['duration']
                               for r in results for dl in r.downloads])), 1)

juga  's avatar
juga committed
837 838 839 840 841 842 843 844
    @staticmethod
    def last_time_from_results(results):
        return unixts_to_isodt_str(round(max([r.time for r in results])))

    @staticmethod
    def rtt_from_results(results):
        # convert from miliseconds to seconds
        rtts = [(round(rtt * 1000)) for r in results for rtt in r.rtts]
juga  's avatar
juga committed
845
        rtt = round(median(rtts)) if rtts else None
juga  's avatar
juga committed
846 847 848 849 850 851 852 853 854
        return rtt

    @staticmethod
    def result_types_from_results(results):
        rt_dict = dict([(result_type_to_key(rt.value),
                         num_results_of_type(results, rt.value))
                        for rt in _ResultType])
        return rt_dict

855 856 857 858 859 860 861 862
    @staticmethod
    def desc_bw_avg_from_results(results):
        """Obtain the last descriptor bandwidth average from the results."""
        for r in reversed(results):
            if r.relay_average_bandwidth is not None:
                return r.relay_average_bandwidth
        return None

863 864 865 866 867 868 869 870
    @staticmethod
    def desc_bw_bur_from_results(results):
        """Obtain the last descriptor bandwidth burst from the results."""
        for r in reversed(results):
            if r.relay_burst_bandwidth is not None:
                return r.relay_burst_bandwidth
        return None

871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886
    @staticmethod
    def consensus_bandwidth_from_results(results):
        """Obtain the last consensus bandwidth from the results."""
        for r in reversed(results):
            if r.consensus_bandwidth is not None:
                return r.consensus_bandwidth
        return None

    @staticmethod
    def consensus_bandwidth_is_unmeasured_from_results(results):
        """Obtain the last consensus unmeasured flag from the results."""
        for r in reversed(results):
            if r.consensus_bandwidth_is_unmeasured is not None:
                return r.consensus_bandwidth_is_unmeasured
        return None

887
    @staticmethod
888 889
    def desc_bw_obs_mean_from_results(results):
        desc_bw_obs_ls = []
890 891
        for r in results:
            if r.relay_observed_bandwidth is not None:
892 893 894
                desc_bw_obs_ls.append(r.relay_observed_bandwidth)
        if desc_bw_obs_ls:
            return max(round(mean(desc_bw_obs_ls)), 1)
895 896 897
        return None

    @staticmethod
898
    def desc_bw_obs_last_from_results(results):
899 900 901 902 903 904
        # the last is at the end of the list
        for r in reversed(results):
            if r.relay_observed_bandwidth is not None:
                return r.relay_observed_bandwidth
        return None

juga  's avatar
juga committed
905 906 907 908 909
    @property
    def bw_keyvalue_tuple_ls(self):
        """Return list of KeyValue Bandwidth Line tuples."""
        # sort the list to generate determinist headers
        keyvalue_tuple_ls = sorted([(k, v) for k, v in self.__dict__.items()
910
                                    if k in BWLINE_KEYS_V1])
juga  's avatar
juga committed
911
        return keyvalue_tuple_ls
juga  's avatar
juga committed
912

juga  's avatar
juga committed
913
    @property
914
    def bw_keyvalue_v1str_ls(self):
juga  's avatar
juga committed
915
        """Return list of KeyValue Bandwidth Line strings following
916
        spec v1.X.X.
juga  's avatar
juga committed
917
        """
918
        bw_keyvalue_str = [KEYVALUE_SEP_V1 .join([k, str(v)])
juga  's avatar
juga committed
919 920
                           for k, v in self.bw_keyvalue_tuple_ls]
        return bw_keyvalue_str
juga  's avatar
juga committed
921

juga  's avatar
juga committed
922
    @property
923 924
    def bw_strv1(self):
        """Return Bandwidth Line string following spec v1.X.X."""
925
        bw_line_str = BWLINE_KEYVALUES_SEP_V1.join(
926
                        self.bw_keyvalue_v1str_ls) + LINE_SEP
juga  's avatar
juga committed
927 928 929 930 931 932
        if len(bw_line_str) > BW_LINE_SIZE:
            # if this is the case, probably there are too many KeyValues,
            # or the limit needs to be changed in Tor
            log.warn("The bandwidth line %s is longer than %s",
                     len(bw_line_str), BW_LINE_SIZE)
        return bw_line_str
juga  's avatar
juga committed
933

juga  's avatar
juga committed
934

935
class V3BWFile(object):
juga  's avatar
juga committed
936
    """
937
    Create a Bandwidth List file following spec version 1.X.X
juga  's avatar
juga committed
938 939 940 941

    :param V3BWHeader v3bwheader: header
    :param list v3bwlines: V3BWLines
    """
juga  's avatar
juga committed
942 943 944 945 946
    def __init__(self, v3bwheader, v3bwlines):
        self.header = v3bwheader
        self.bw_lines = v3bwlines

    def __str__(self):
947
        return str(self.header) + ''.join([str(bw_line) or ''
juga  's avatar
juga committed
948 949
                                           for bw_line in self.bw_lines])

juga  's avatar
juga committed
950
    @classmethod
951 952
    def from_results(cls, results, scanner_country=None,
                     destinations_countries=None, state_fpath='',
953
                     scale_constant=SBWS_SCALE_CONSTANT,
954 955
                     scaling_method=TORFLOW_SCALING,
                     torflow_obs=TORFLOW_OBS_LAST,
juga  's avatar
juga committed
956
                     torflow_cap=TORFLOW_BW_MARGIN,
957
                     round_digs=PROP276_ROUND_DIG,
958
                     secs_recent=None, secs_away=None, min_num=0,
959 960
                     consensus_path=None, max_bw_diff_perc=MAX_BW_DIFF_PERC,
                     reverse=False):
961 962 963 964 965 966
        """Create V3BWFile class from sbws Results.

        :param dict results: see below
        :param str state_fpath: path to the state file
        :param int scaling_method:
            Scaling method to obtain the bandwidth
967
            Possible values: {None, SBWS_SCALING, TORFLOW_SCALING} = {0, 1, 2}
968 969 970 971 972 973 974 975 976 977 978
        :param int scale_constant: sbws scaling constant
        :param int torflow_obs: method to choose descriptor observed bandwidth
        :param bool reverse: whether to sort the bw lines descending or not

        Results are in the form::

            {'relay_fp1': [Result1, Result2, ...],
             'relay_fp2': [Result1, Result2, ...]}

        """
        log.info('Processing results to generate a bandwidth list file.')
979 980
        header = V3BWHeader.from_results(results, scanner_country,
                                         destinations_countries, state_fpath)
981
        bw_lines_raw = []
982
        bw_lines_excluded = []
juga  's avatar
juga committed
983 984
        number_consensus_relays = cls.read_number_consensus_relays(
            consensus_path)
985
        state = State(state_fpath)
986 987 988

        # Create a dictionary with the number of relays excluded by any of the
        # of the filtering rules that makes relays non-`eligible`.
989
        # NOTE: In HEADER_RECENT_MEASUREMENTS_EXCLUDED_KEYS it is
990 991 992
        # explained what are the KeyValues.
        # See also the comments in `from_results`.
        exclusion_dict = dict(
993
            [(k, 0) for k in HEADER_RECENT_MEASUREMENTS_EXCLUDED_KEYS]
994
            )
juga  's avatar
juga committed
995
        for fp, values in results.items():
996
            # log.debug("Relay fp %s", fp)
997 998
            line, reason = V3BWLine.from_results(values, secs_recent,
                                                 secs_away, min_num)
999 1000
            # If there is no reason it means the line will not be excluded.
            if not reason:
juga  's avatar
juga committed
1001
                bw_lines_raw.append(line)
1002
            else:
1003 1004 1005
                # Store the excluded lines to include them in the bandwidth
                # file.
                bw_lines_excluded.append(line)
1006 1007 1008 1009
                exclusion_dict[reason] = exclusion_dict.get(reason, 0) + 1
        # Add the headers with the number of excluded relays by reason
        header.add_relays_excluded_counters(exclusion_dict)

1010
        if not bw_lines_raw:
1011 1012
            # It could be possible to scale the lines that were successful
            # even if excluded, but is not done here.
1013 1014
            log.info("After applying restrictions to the raw results, "
                     "there is not any. Scaling can not be applied.")
1015
            # Update the header and log the progress.
juga  's avatar
juga committed
1016
            cls.update_progress(
1017 1018 1019 1020 1021
                cls, 0, header, number_consensus_relays, state)
            # Set the lines that would be excluded anyway (`vote=0`) with
            # `under_min_report=1`
            cls.set_under_min_report(bw_lines_excluded)
            # Create the bandwidth file with the lines that would be excluded.
1022
            return cls(header, bw_lines_excluded)
1023 1024
        if scaling_method == SBWS_SCALING:
            bw_lines = cls.bw_sbws_scale(bw_lines_raw, scale_constant)
1025
            cls.warn_if_not_accurate_enough(bw_lines, scale_constant)
1026 1027
            # log.debug(bw_lines[-1])
        elif scaling_method == TORFLOW_SCALING:
juga  's avatar
juga committed
1028
            bw_lines = cls.bw_torflow_scale(bw_lines_raw, torflow_obs,
1029
                                            torflow_cap, round_digs)