multicalc 0.7.0

Rust scientific computing for single and multi-variable calculus
Documentation 
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

"""Generate Gaussian quadrature node/weight tables for orders 1..=MAX_ORDER.

Requires numpy. Emits one module file per family into src/gaussian_tables/,
each holding `(weight, abscissa)` arrays plus a flat index table used by
`gaussian_tables::nodes`.

Run from anywhere:  python3 scripts/build_gaussian_integration_tables.py
"""

import os

from numpy.polynomial.hermite import hermgauss
from numpy.polynomial.laguerre import laggauss
from numpy.polynomial.legendre import leggauss

MAX_ORDER = 30
OUT_DIR = os.path.join(os.path.dirname(__file__), "..", "src", "gaussian_tables")

FAMILIES = {
    "legendre": ("LEGENDRE", leggauss),
    "hermite": ("HERMITE", hermgauss),
    "laguerre": ("LAGUERRE", laggauss),
}


def render(name, gauss):
    out = [
        "// Auto-generated by scripts/build_gaussian_integration_tables.py",
        "// Do not edit manually.",
        "",
    ]
    for order in range(1, MAX_ORDER + 1):
        abscissae, weights = gauss(order)
        out.append("pub static {}_{}: [(f64, f64); {}] = [".format(name, order, order))
        for w, x in zip(weights, abscissae):
            out.append("    ({!r}, {!r}),".format(float(w), float(x)))
        out.append("];")
    out.append("")
    out.append("pub static {}: [&[(f64, f64)]; {}] = [".format(name, MAX_ORDER + 1))
    out.append("    &[],")
    for i in range(1, MAX_ORDER + 1):
        out.append("    &{}_{},".format(name, i))
    out.append("];")
    out.append("")
    return "\n".join(out)


def main():
    os.makedirs(OUT_DIR, exist_ok=True)
    for stem, (name, gauss) in FAMILIES.items():
        path = os.path.join(OUT_DIR, stem + ".rs")
        with open(path, "w") as f:
            f.write(render(name, gauss))
        print("wrote", path)


if __name__ == "__main__":
    main()