numrs2 0.3.3

A Rust implementation inspired by NumPy for numerical computing (NumRS2)
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

// Test complex seed scenarios that might reveal issues
use numrs2::random::advanced_distributions::vonmises;
use numrs2::random::distributions::{normal, set_seed};

#[test]
#[ignore = "Seeding behavior changed during SciRS2 migration - requires seeding implementation fix"]
fn test_complex_seed_scenarios() {
    println!("Testing complex seed scenarios...");

    // Scenario 1: Mixed function calls with same seed
    println!("\nScenario 1: Mixed functions, same seed");
    set_seed(42);
    let normal_a = normal(0.0, 1.0, &[3]).unwrap();
    let vonmises_a = vonmises(0.0, 1.0, &[3]).unwrap();

    set_seed(42);
    let normal_b = normal(0.0, 1.0, &[3]).unwrap();
    let vonmises_b = vonmises(0.0, 1.0, &[3]).unwrap();

    println!("Normal A: {:?}", normal_a.to_vec());
    println!("Normal B: {:?}", normal_b.to_vec());
    println!("Vonmises A: {:?}", vonmises_a.to_vec());
    println!("Vonmises B: {:?}", vonmises_b.to_vec());

    assert_eq!(
        normal_a.to_vec(),
        normal_b.to_vec(),
        "Normal sequences should be identical"
    );
    assert_eq!(
        vonmises_a.to_vec(),
        vonmises_b.to_vec(),
        "Vonmises sequences should be identical"
    );

    // Scenario 2: Reset seed multiple times
    println!("\nScenario 2: Multiple seed resets");
    set_seed(123);
    let v1 = vonmises(0.0, 1.0, &[2]).unwrap();
    set_seed(456);
    let v2 = vonmises(0.0, 1.0, &[2]).unwrap();
    set_seed(123);
    let v3 = vonmises(0.0, 1.0, &[2]).unwrap();

    println!("Vonmises 123: {:?}", v1.to_vec());
    println!("Vonmises 456: {:?}", v2.to_vec());
    println!("Vonmises 123 again: {:?}", v3.to_vec());

    assert_eq!(
        v1.to_vec(),
        v3.to_vec(),
        "Same seed should produce identical sequences"
    );
    assert_ne!(
        v1.to_vec(),
        v2.to_vec(),
        "Different seeds should produce different sequences"
    );

    // Scenario 3: vonmises with special parameters
    println!("\nScenario 3: Vonmises with special parameters");
    set_seed(42);
    let vm_uniform = vonmises(0.0, 0.0, &[3]).unwrap(); // Should be uniform-like
    set_seed(42);
    let vm_uniform2 = vonmises(0.0, 0.0, &[3]).unwrap();

    println!("Vonmises uniform 1: {:?}", vm_uniform.to_vec());
    println!("Vonmises uniform 2: {:?}", vm_uniform2.to_vec());

    assert_eq!(
        vm_uniform.to_vec(),
        vm_uniform2.to_vec(),
        "Even uniform vonmises should be reproducible"
    );
}