savvy 0.10.0

A simple R extension interface
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
76
77
78
79
80
81
82
83
84

// NOTE: No implementation is provided for bool because R's bool is tricky.
// https://cpp11.r-lib.org/articles/cpp11.html#na

pub trait NotAvailableValue {
    fn is_na(&self) -> bool;
    fn na() -> Self;
}

impl NotAvailableValue for f64 {
    fn is_na(&self) -> bool {
        unsafe { savvy_ffi::R_IsNA(*self) != 0 }
    }

    fn na() -> Self {
        unsafe { savvy_ffi::R_NaReal }
    }
}

impl NotAvailableValue for i32 {
    fn is_na(&self) -> bool {
        unsafe { *self == savvy_ffi::R_NaInt }
    }

    fn na() -> Self {
        unsafe { savvy_ffi::R_NaInt }
    }
}

#[cfg(feature = "complex")]
impl NotAvailableValue for num_complex::Complex64 {
    fn is_na(&self) -> bool {
        unsafe { self.re == savvy_ffi::R_NaReal }
    }

    fn na() -> Self {
        unsafe {
            num_complex::Complex64 {
                re: savvy_ffi::R_NaReal,
                im: savvy_ffi::R_NaReal,
            }
        }
    }
}

use std::sync::OnceLock;

pub(crate) static NA_CHAR_PTR: OnceLock<&str> = OnceLock::new();

impl NotAvailableValue for &str {
    fn is_na(&self) -> bool {
        self.as_ptr() == Self::na().as_ptr()
    }

    // I use the underlying "NA" string of R_NaString directry here, but this
    // wasn't possible on extendr due to some unobvious reason related to
    // concurrency.
    //
    // cf., https://github.com/extendr/extendr/issues/483#issuecomment-1435499525
    fn na() -> Self {
        NA_CHAR_PTR.get_or_init(|| unsafe {
            let c_ptr = savvy_ffi::R_CHAR(savvy_ffi::R_NaString) as _;
            std::str::from_utf8_unchecked(std::slice::from_raw_parts(c_ptr, 2))
        })
    }
}

/// Return true if the SEXP is a length-1 of vector containing NA.
pub(crate) unsafe fn is_scalar_na(x: savvy_ffi::SEXP) -> bool {
    if unsafe { savvy_ffi::Rf_xlength(x) } != 1 {
        return false;
    }

    let ty = unsafe { savvy_ffi::TYPEOF(x) };
    match ty {
        savvy_ffi::INTSXP => unsafe { savvy_ffi::INTEGER_ELT(x, 0) }.is_na(),
        savvy_ffi::REALSXP => unsafe { savvy_ffi::REAL_ELT(x, 0) }.is_na(),
        #[cfg(feature = "complex")]
        savvy_ffi::CPLXSXP => unsafe { savvy_ffi::COMPLEX_ELT(x, 0) }.is_na(),
        savvy_ffi::LGLSXP => unsafe { savvy_ffi::LOGICAL_ELT(x, 0) }.is_na(),
        savvy_ffi::RAWSXP => false, // raw doesn't have NA
        savvy_ffi::STRSXP => unsafe { savvy_ffi::STRING_ELT(x, 0) == savvy_ffi::R_NaString },
        _ => false,
    }
}