astroapers 0.1.1

Rust-backed aperture overlap and summation utilities
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

//! Core aperture overlap and summation kernels.
//!
//! Pixel coordinates follow the SEP/Photutils convention: pixel `(i, j)` is
//! centered at `(i, j)` and covers `[i - 0.5, i + 0.5]` in x and y.

use numpy::{
    Element, PyArray1, PyArrayMethods, PyReadonlyArray1, PyReadonlyArray2, PyReadonlyArrayDyn,
    PyUntypedArrayMethods,
};
use pyo3::exceptions::{PyMemoryError, PyValueError};
use pyo3::prelude::*;
use rayon::prelude::*;
use std::env;
use std::f64::consts::PI;
use std::sync::atomic::{AtomicUsize, Ordering};

use crate::geometry::{
    bbox_in_frame, circ_bbox, circ_pixel_area, ellip_bbox, ellip_inner_outer_radius2,
    ellip_pixel_area, extent_bbox, rect_pixel_area, Ellip, SQRT_HALF,
};

const DEFAULT_PARALLEL_THRESHOLD: usize = 64;
static PARALLEL_THRESHOLD: AtomicUsize = AtomicUsize::new(0);

type BBoxColumns = (Vec<isize>, Vec<isize>, Vec<isize>, Vec<isize>);
type WeightManyResult = (
    Vec<Py<PyAny>>,
    Vec<isize>,
    Vec<isize>,
    Vec<isize>,
    Vec<isize>,
);

mod apsum;
mod bbox;
mod mask;
mod npix;
mod parallel_api;
mod path;
mod raster;
mod summation;
mod support;
mod weights;

pub(crate) fn register(m: &Bound<'_, PyModule>) -> PyResult<()> {
    apsum::register(m)?;
    npix::register(m)?;
    bbox::register(m)?;
    mask::register(m)?;
    parallel_api::register(m)?;
    path::register(m)?;
    weights::register(m)?;
    Ok(())
}