differential_evolution/population_initialization.py at master · vrbaj/differential_evolution · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
"""Backward-compatible population initialization helpers."""

from __future__ import annotations

import random
from collections.abc import Callable, Sequence

from differential_evolution.initializers import (
    OppositionInitializer,
    QuasiOppositionInitializer,
    RandomInitializer,
    SobolInitializer,
    TentInitializer,
)


def random_initialization(
    population_size: int,
    bounds: Sequence[Sequence[float]],
    seed: int | float | str | bytes | bytearray | None = None,
) -> list[list[float]]:
    initializer = RandomInitializer()
    return initializer(
        population_size,
        [(lower, upper) for lower, upper in bounds],
        lambda _: 0.0,
        random.Random(seed),
    )


def tent_initialization(
    population_size: int,
    bounds: Sequence[Sequence[float]],
    seed: int | float | str | bytes | bytearray | None = None,
) -> list[list[float]]:
    initializer = TentInitializer()
    return initializer(
        population_size,
        [(lower, upper) for lower, upper in bounds],
        lambda _: 0.0,
        random.Random(seed),
    )


def obl_initialization(
    population_size: int,
    bounds: Sequence[Sequence[float]],
    cost_function: Callable[[list[float]], float],
    seed: int | float | str | bytes | bytearray | None = None,
) -> list[list[float]]:
    initializer = OppositionInitializer()
    return initializer(
        population_size,
        [(lower, upper) for lower, upper in bounds],
        cost_function,
        random.Random(seed),
    )


def qobl_initialization(
    population_size: int,
    bounds: Sequence[Sequence[float]],
    cost_function: Callable[[list[float]], float],
    seed: int | float | str | bytes | bytearray | None = None,
) -> list[list[float]]:
    initializer = QuasiOppositionInitializer()
    return initializer(
        population_size,
        [(lower, upper) for lower, upper in bounds],
        cost_function,
        random.Random(seed),
    )


def sobol_initialization(
    population_size: int,
    bounds: Sequence[Sequence[float]],
) -> list[list[float]]:
    initializer = SobolInitializer()
    return initializer(
        population_size,
        [(lower, upper) for lower, upper in bounds],
        lambda _: 0.0,
        random.Random(0),
    )