map2fig 0.7.8

Fast, publication-quality HEALPix sky map visualization in Rust
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

use map2fig::healpix::{ang2pix_ring, pix2ang_ring};
use std::f64::consts::PI;

fn main() {
    divan::main();
}

#[divan::bench]
fn bench_pix2ang_ring_bulk() {
    let nside = divan::black_box(512i64);
    let mut sum = 0.0f64;

    for pix in 0..1000 {
        let pix = pix % (12 * nside * nside);
        let (theta, phi) = pix2ang_ring(nside, pix);
        sum += theta + phi;
    }

    divan::black_box(sum);
}

#[divan::bench]
fn bench_ang2pix_ring_bulk() {
    let nside = divan::black_box(512i64);
    let mut sum = 0i64;

    for i in 0..1000 {
        let theta = PI * i as f64 / 1000.0;
        let phi = 2.0 * PI * i as f64 / 1000.0;
        sum += ang2pix_ring(nside, theta, phi);
    }

    divan::black_box(sum);
}

#[divan::bench]
fn bench_pix2ang_ring_scattered() {
    let nside = divan::black_box(8192i64);
    let mut sum = 0.0f64;

    // Simulate scattered access pattern (cache-hostile)
    for i in 0..100 {
        let pix = (i * 31337) % (12 * nside * nside);
        let (theta, phi) = pix2ang_ring(nside, pix);
        sum += theta + phi;
    }

    divan::black_box(sum);
}

#[divan::bench]
fn bench_ang2pix_ring_scattered() {
    let nside = divan::black_box(8192i64);
    let mut sum = 0i64;

    // Wide range of angles to avoid branch prediction bias
    for i in 0..100 {
        let angle_varied = (i as f64 * 0.0432) % 1.0; // Irrational multiplier
        let theta = PI * angle_varied;
        let phi = 2.0 * PI * angle_varied.sin().abs();
        sum += ang2pix_ring(nside, theta, phi);
    }

    divan::black_box(sum);
}