semiflow 0.15.0-beta

no_std Rust PDE solver: evolution equations via Chernoff approximation of operator semigroups (heat, diffusion, Schrodinger)
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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75

//! Gauss-Laguerre 32-pt quadrature data (Abramowitz & Stegun Table 25.9).
//!
//! Shared by `resolvent.rs`. All values are `f64` constants.

/// 32-point Gauss-Laguerre quadrature nodes `s_k` for `∫₀^∞ e^{-s} h(s) ds`.
pub(crate) const GL32_NODES: [f64; 32] = [
    0.044_489_365_833_267,
    0.234_526_109_519_619,
    0.576_884_629_301_886,
    1.072_448_753_817_82,
    1.722_408_776_444_65,
    2.528_336_706_425_79,
    3.492_213_273_021_99,
    4.616_456_769_749_77,
    5.903_958_504_174_24,
    7.358_126_733_186_24,
    8.982_940_924_212_59,
    10.783_018_632_54,
    12.763_697_986_742_8,
    14.931_139_755_522_6,
    17.292_454_336_715_3,
    19.855_860_940_336_1,
    22.630_889_013_196_8,
    25.628_636_022_459_2,
    28.862_101_816_323_5,
    32.346_629_153_964_8,
    36.100_494_805_752,
    40.145_719_771_539_4,
    44.509_207_995_754_9,
    49.224_394_987_308_6,
    54.333_721_333_397,
    59.892_509_162_134,
    65.975_377_287_935_2,
    72.687_628_090_662_7,
    80.187_446_977_913_5,
    88.735_340_417_892_4,
    98.829_542_868_283_9,
    111.751_398_097_938,
];

/// 32-point Gauss-Laguerre quadrature weights `w_k`.
pub(crate) const GL32_WEIGHTS: [f64; 32] = [
    0.109_218_341_952_385,
    0.210_443_107_938_813,
    0.235_213_229_669_848,
    0.195_903_335_972_881,
    0.129_983_786_286_071,
    0.070_578_623_865_717_4,
    0.031_760_912_509_175_1,
    0.011_918_214_834_838_5,
    0.003_738_816_294_611_53,
    0.000_980_803_306_614_955,
    0.000_214_864_918_801_364,
    3.920_341_967_987_95e-5,
    5.934_541_612_868_63e-6,
    7.416_404_578_667_55e-7,
    7.604_567_879_120_78e-8,
    6.350_602_226_625_81e-9,
    4.281_382_971_040_93e-10,
    2.305_899_491_891_34e-11,
    9.799_379_288_727_09e-13,
    3.237_801_657_729_27e-14,
    8.171_823_443_420_72e-16,
    1.542_133_833_393_82e-17,
    2.119_792_290_163_62e-19,
    2.054_429_673_788_05e-21,
    1.346_982_586_637_39e-23,
    5.661_294_130_397_35e-26,
    1.418_560_545_463_04e-28,
    1.913_375_494_454_22e-31,
    1.192_248_760_098_54e-34,
    2.671_511_219_240_13e-38,
    1.338_616_942_106_25e-42,
    4.510_536_193_898_97e-48,
];