pi.py

#!/usr/bin/env python3
#[x-cmds]: UPDATE
"""Calculate the value of pi to a specified number of decimal places."""
from typing import Iterator
from xulbux.console import Throbber
from xulbux import FormatCodes, Console
import math
import time
import sys


ARGS = Console.get_args({"decimal_places": "before"})
REFERENCE_TIMES: dict[int, float] = {
    1000: 0.01,  # 1K DIGITS
    5000: 0.175,  # 5K DIGITS
    10000: 0.75,  # 10K DIGITS
    25000: 5.10,  # 25K DIGITS
    50000: 25,  # 50K DIGITS
    100000: 120,  # 100K DIGITS
    500000: 3000,  # 500K DIGITS
    1000000: 75000,  # 1M DIGITS
}


def get_hardware_score() -> float:
    try:
        import psutil

        cpu_freq = psutil.cpu_freq()
        max_freq = cpu_freq.max if cpu_freq else 3000
        cpu_count = psutil.cpu_count(logical=False) or 1
        memory = psutil.virtual_memory()
        memory_factor = 1 + (0.3 * (1 - memory.available / memory.total))
        return ((max_freq * math.sqrt(cpu_count)) / 4000) / memory_factor

    except (ImportError, AttributeError):
        # FALLBACK: USE DEFAULT HARDWARE SCORE (ASSUMES MODERN MID-RANGE SYSTEM)
        return 1.0


def estimate_runtime(precision: int) -> float:
    ref_points = sorted(REFERENCE_TIMES.keys())
    if precision <= 100:
        start_time = time.time()
        _ = pi(precision)
        return time.time() - start_time
    if precision >= max(ref_points):
        base_time = REFERENCE_TIMES[max(ref_points)]
        scaling = (precision / max(ref_points))**2.0
        if precision > 1000000:
            scaling *= 1.2
    else:
        upper_idx = next(i for i, x in enumerate(ref_points) if x >= precision)
        lower_idx = max(0, upper_idx - 1)
        lower_point = ref_points[lower_idx]
        upper_point = ref_points[upper_idx]
        lower_time = REFERENCE_TIMES[lower_point]
        upper_time = REFERENCE_TIMES[upper_point]
        if lower_point == upper_point:
            log_factor = 1
        else:
            raw_factor = precision / lower_point
            log_factor = (math.log(raw_factor)**2) / (math.log(upper_point / lower_point))
        base_time = lower_time * (upper_time / lower_time)**log_factor
        scaling = 1.0
    estimated_time = (base_time * scaling) / get_hardware_score()
    if estimated_time < 0.01:
        estimated_time = 0.01
    if precision <= 5000:
        correction_factor = 1.0
    elif precision <= 10000:
        correction_factor = 1.5
    elif precision <= 25000:
        correction_factor = 1.8
    elif precision <= 50000:
        correction_factor = 1.0
    elif precision <= 100000:
        correction_factor = 0.9
    else:
        correction_factor = 1.0
    estimated_time *= correction_factor
    return round(estimated_time, 2)


def format_time(seconds: float, short: bool = False, pretty_print: bool = False) -> str:
    units = (
        (
            ("SMBH", 1e106 * 365.25 * 24 * 60 * 60),
            ("HD", 1e100 * 365.25 * 24 * 60 * 60),
            ("BH", 1e40 * 365.25 * 24 * 60 * 60),
            ("DE", 1e14 * 365.25 * 24 * 60 * 60),
            ("SE", 1e12 * 365.25 * 24 * 60 * 60),
            ("GY", 225e6 * 365.25 * 24 * 60 * 60),
            ("HT", 13.8e9 * 365.25 * 24 * 60 * 60),
            ("Q", 1e30 * 365.25 * 24 * 60 * 60),
            ("R", 1e27 * 365.25 * 24 * 60 * 60),
            ("Y", 1e24 * 365.25 * 24 * 60 * 60),
            ("Z", 1e21 * 365.25 * 24 * 60 * 60),
            ("E", 1e18 * 365.25 * 24 * 60 * 60),
            ("P", 1e15 * 365.25 * 24 * 60 * 60),
            ("T", 1e12 * 365.25 * 24 * 60 * 60),
            ("G", 1e9 * 365.25 * 24 * 60 * 60),
            ("M", 1e6 * 365.25 * 24 * 60 * 60),
            ("k", 1e3 * 365.25 * 24 * 60 * 60),
            ("y", 365.25 * 24 * 60 * 60),
            ("mo", 30 * 24 * 60 * 60),
            ("w", 7 * 24 * 60 * 60),
            ("d", 24 * 60 * 60),
            ("h", 60 * 60),
            ("m", 60),
            ("s", 1),
        ),
        (
            ("supermassive black hole lifespan", 1e106 * 365.25 * 24 * 60 * 60),
            ("universe heat death", 1e100 * 365.25 * 24 * 60 * 60),
            ("black hole era", 1e40 * 365.25 * 24 * 60 * 60),
            ("degenerate era", 1e14 * 365.25 * 24 * 60 * 60),
            ("stelliferous era", 1e12 * 365.25 * 24 * 60 * 60),
            ("galactic year", 225e6 * 365.25 * 24 * 60 * 60),
            ("Hubble time", 13.8e9 * 365.25 * 24 * 60 * 60),
            ("quetta-year", 1e30 * 365.25 * 24 * 60 * 60),
            ("ronna-year", 1e27 * 365.25 * 24 * 60 * 60),
            ("yotta-year", 1e24 * 365.25 * 24 * 60 * 60),
            ("zetta-year", 1e21 * 365.25 * 24 * 60 * 60),
            ("exa-year", 1e18 * 365.25 * 24 * 60 * 60),
            ("peta-year", 1e15 * 365.25 * 24 * 60 * 60),
            ("tera-year", 1e12 * 365.25 * 24 * 60 * 60),
            ("giga-year", 1e9 * 365.25 * 24 * 60 * 60),
            ("mega-year", 1e6 * 365.25 * 24 * 60 * 60),
            ("kilo-year", 1e3 * 365.25 * 24 * 60 * 60),
            ("year", 365.25 * 24 * 60 * 60),
            ("month", 30 * 24 * 60 * 60),
            ("week", 7 * 24 * 60 * 60),
            ("day", 24 * 60 * 60),
            ("hour", 60 * 60),
            ("minute", 60),
            ("second", 1),
        ),
    )
    parts: list[str] = []
    b_val, val_name, a_name = (
        "[b|br:cyan]" if pretty_print else "",
        f"{'' if short else ' '}{'[_b|i|cyan]' if pretty_print else ''}",
        "[_i|br:cyan]" if pretty_print else "",
    )
    for name, formula in units[0 if short else 1]:
        if (val := int(seconds // formula)) > 0:
            if not short:
                val = f"{val:,}"
            parts.append(f"{b_val}{val}{val_name}{name if val == '1' or short else f'{name}s'}{a_name}")
            seconds %= formula
    if not parts:
        formatted_seconds = f"{f'{seconds:.3f}'.rstrip('0').rstrip('.')}"
        parts.append(
            f"{b_val}{formatted_seconds}{val_name}{units[0 if short else 1][-1][0] if seconds == 1 or short else f'{units[0 if short else 1][-1][0]}s'}{a_name}"
        )
    if short:
        return ("[dim](:)" if pretty_print else ":").join(parts)
    if len(parts) > 1:
        return (("[dim] + [_dim]" if pretty_print else ", ").join(parts[:-1]) +
                ("[dim] + [_dim]" if pretty_print else " and ") + parts[-1])
    return parts[0]


def p() -> Iterator[int]:
    q, r, t, j = 1, 180, 60, 2
    while True:
        u, y = 3 * (3 * j + 1) * (3 * j + 2), (q * (27 * j - 12) + 5 * r) // (5 * t)
        yield y
        q, r, t, j = (10 * q * j * (2 * j - 1), 10 * u * (q * (5 * j - 2) + r - y * t), t * u, j + 1)


def pi(decimals: int = 10) -> str:
    _p = p()
    return "3." + "".join(str(next(_p)) for _ in range(decimals + 1))[1:]


def main() -> None:
    global CALC_DONE
    input_k = (
        int(ARGS.decimal_places.values[0]) \
        if ARGS.decimal_places.values and ARGS.decimal_places.values[0].replace("_", "").isdigit()
        else 10
    )

    if (estimated_secs := estimate_runtime(input_k)) >= 604800:
        FormatCodes.print(
            f"\n[b|bg:black]( π [in]( CALCULATION WOULD TAKE TOO LONG ))\n\n{format_time(estimated_secs, pretty_print=True)}[_]\n"
        )
    else:
        FormatCodes.print(
            f"\n[dim](Will take about [b]{format_time(estimated_secs)}[_|dim] to calculate:)" if estimated_secs > 1 else ""
        )
        result = None

        try:
            with Throbber().context():
                result = pi(input_k)
        except MemoryError:
            FormatCodes.print(
                "\r[b|br:yellow](Your computer doesn't have enough memory for this calculation!)",
                "[b|bg:black]( π [in]( CALCULATION WOULD TAKE THIS LONG IF YOU HAD ENOUGH MEMORY ))\n",
                f"{format_time(estimated_secs, pretty_print=True)}[_]",
                end="\n\n",
                sep="\n",
            )
        except KeyboardInterrupt:
            FormatCodes.print("\r[b|br:red](✗)  \n")
            sys.exit(0)

        if result:
            FormatCodes.print(f"\r[br:cyan]({result})\n")
        else:
            FormatCodes.print("\r[b|br:red](✗)  \n")


if __name__ == "__main__":
    main()