castep-cell-io 0.6.0

use castep_cell_fmt::{Cell, CellValue, ToCell, ToCellValue, FromKeyValue, CResult};
use castep_cell_fmt::query::*;
use castep_cell_fmt::parse::FromCellValue;
use crate::units::EnergyUnit;

/// Controls the tolerance for accepting convergence of the Fermi-energy.
///
/// Keyword type: Real
///
/// Default: 0.1 × ELEC_EIGENVALUE_TOL
///
/// Example:
/// EFERMI_TOL : 0.0000007 eV
/// EFERMI_TOL : 0.0000007 (uses default unit, likely eV)
#[derive(Debug, Clone, Copy, PartialEq, PartialOrd)]
pub struct EFermiTol {
    /// The Fermi energy tolerance value.
    pub value: f64,
    /// The optional unit of the energy value.
    /// If None, a default unit (likely eV) is implied.
    pub unit: Option<EnergyUnit>,
}

impl FromKeyValue for EFermiTol {
    const KEY_NAME: &'static str = "EFERMI_TOL";

    fn from_cell_value_kv(value: &CellValue<'_>) -> CResult<Self> {
        match value {
            CellValue::Array(arr) if !arr.is_empty() => {
                let val = value_as_f64(&arr[0])?;
                let unit = if arr.len() > 1 {
                    Some(EnergyUnit::from_cell_value(&arr[1])?)
                } else {
                    None
                };
                Ok(Self { value: val, unit })
            }
            _ => {
                let val = value_as_f64(value)?;
                Ok(Self { value: val, unit: None })
            }
        }
    }
}

impl ToCell for EFermiTol {
    fn to_cell(&self) -> Cell<'_> {
        Cell::KeyValue("EFERMI_TOL", self.to_cell_value())
    }
}

impl ToCellValue for EFermiTol {
    fn to_cell_value(&self) -> CellValue<'_> {
        CellValue::Array(
            [
                CellValue::Float(self.value),
                self.unit
                    .as_ref()
                    .map(|u| u.to_cell_value())
                    .unwrap_or(CellValue::Null),
            ]
            .to_vec(),
        )
    }
}