orchid/src/parse/numeric.rs

use std::num::IntErrorKind;
use std::ops::Range;
use std::rc::Rc;

use ordered_float::NotNan;

use super::context::Context;
use super::errors::NaNLiteral;
use super::lexer::split_filter;
use crate::error::{ProjectError, ProjectResult};
use crate::foreign::Atom;
use crate::systems::stl::Numeric;

#[derive(Debug, Clone, PartialEq, Eq)]
pub enum NumErrorKind {
  NaN,
  Overflow,
  InvalidDigit,
}
impl NumErrorKind {
  fn from_int(kind: &IntErrorKind) -> Self {
    match kind {
      IntErrorKind::InvalidDigit => Self::InvalidDigit,
      IntErrorKind::NegOverflow | IntErrorKind::PosOverflow => Self::Overflow,
      _ => panic!("Impossible error condition"),
    }
  }
}

#[derive(Debug, Clone, PartialEq, Eq)]
pub struct NumError {
  pub range: Range<usize>,
  pub kind: NumErrorKind,
}

impl NumError {
  pub fn into_proj(
    self,
    len: usize,
    tail: &str,
    ctx: &(impl Context + ?Sized),
  ) -> Rc<dyn ProjectError> {
    let start = ctx.source().len() - tail.len() - len + self.range.start;
    let location = ctx.range_loc(start..start + self.range.len());
    match self.kind {
      NumErrorKind::NaN => NaNLiteral(location).rc(),
      _ => panic!(),
      // NumErrorKind::Int(iek) => IntError(location, iek).rc(),
    }
  }
}

pub fn parse_num(string: &str) -> Result<Numeric, NumError> {
  let overflow_err =
    NumError { range: 0..string.len(), kind: NumErrorKind::Overflow };
  let (radix, noprefix, pos) =
    (string.strip_prefix("0x").map(|s| (16u8, s, 2)))
      .or_else(|| string.strip_prefix("0b").map(|s| (2u8, s, 2)))
      .or_else(|| string.strip_prefix("0o").map(|s| (8u8, s, 2)))
      .unwrap_or((10u8, string, 0));
  // identity
  let (base, exponent) = match noprefix.split_once('p') {
    Some((b, e)) => {
      let (s, d, len) = e.strip_prefix('-').map_or((1, e, 0), |ue| (-1, ue, 1));
      (b, s * int_parse(d, radix, pos + b.len() + 1 + len)? as i32)
    },
    None => (noprefix, 0),
  };
  match base.split_once('.') {
    None => {
      let base_usize = int_parse(base, radix, pos)?;
      if let Ok(pos_exp) = u32::try_from(exponent) {
        if let Some(radical) = usize::from(radix).checked_pow(pos_exp) {
          let number = base_usize.checked_mul(radical).ok_or(overflow_err)?;
          return Ok(Numeric::Uint(number));
        }
      }
      let f = (base_usize as f64) * (radix as f64).powi(exponent);
      let err = NumError { range: 0..string.len(), kind: NumErrorKind::NaN };
      Ok(Numeric::Float(NotNan::new(f).map_err(|_| err)?))
    },
    Some((whole, part)) => {
      let whole_n = int_parse(whole, radix, pos)? as f64;
      let part_n = int_parse(part, radix, pos + whole.len() + 1)? as f64;
      let real_val = whole_n + (part_n / radix.pow(part.len() as u32) as f64);
      let f = real_val * (radix as f64).powi(exponent);
      Ok(Numeric::Float(NotNan::new(f).expect("None of the inputs are NaN")))
    },
  }
}

fn int_parse(s: &str, radix: u8, start: usize) -> Result<usize, NumError> {
  let s = s.chars().filter(|c| *c != '_').collect::<String>();
  let range = start..(start + s.len());
  usize::from_str_radix(&s, radix as u32)
    .map_err(|e| NumError { range, kind: NumErrorKind::from_int(e.kind()) })
}

pub fn numchar(c: char) -> bool { c.is_alphanumeric() | "._-".contains(c) }
pub fn numstart(c: char) -> bool { c.is_ascii_digit() }

pub fn lex_numeric<'a>(
  data: &'a str,
  ctx: &dyn Context,
) -> Option<ProjectResult<(Atom, &'a str)>> {
  data.chars().next().filter(|c| numstart(*c)).map(|_| {
    let (num_str, tail) = split_filter(data, numchar);
    match parse_num(num_str) {
      Ok(Numeric::Float(f)) => Ok((Atom::new(f), tail)),
      Ok(Numeric::Uint(i)) => Ok((Atom::new(i), tail)),
      Err(e) => Err(e.into_proj(num_str.len(), tail, ctx)),
    }
  })
}

#[cfg(test)]
mod test {
  use crate::parse::numeric::parse_num;
  use crate::systems::stl::Numeric;

  #[test]
  fn just_ints() {
    let test = |s, n| assert_eq!(parse_num(s), Ok(Numeric::Uint(n)));
    test("12345", 12345);
    test("0xcafebabe", 0xcafebabe);
    test("0o751", 0o751);
    test("0b111000111", 0b111000111);
  }

  #[test]
  fn decimals() {
    let test = |s, n| assert_eq!(parse_num(s).map(|n| n.as_f64()), Ok(n));
    test("3.1417", 3.1417);
    test("3.1417", 3_f64 + 1417_f64 / 10000_f64);
    test("0xf.cafe", 0xf as f64 + 0xcafe as f64 / 0x10000 as f64);
    test("34p3", 34000f64);
    test("0x2p3", (0x2 * 0x1000) as f64);
    test("1.5p3", 1500f64);
    test("0x2.5p3", (0x25 * 0x100) as f64);
  }
}

#[must_use]
pub fn print_nat16(num: NotNan<f64>) -> String {
  let exp = num.log(16.0).floor();
  let man = num / 16_f64.powf(exp);
  format!("{man}p{exp:.0}")
}