From b3a6d16c3d76a5201e46e4a5595298c90bb8f8b5 Mon Sep 17 00:00:00 2001
From: atse <atse@m-labs.hk>
Date: Mon, 28 Aug 2023 14:52:09 +0800
Subject: [PATCH] steinhart_hart: Beta Parameter uom dimensions

The Beta Parameter is not dimensionless, and has unit kelvin.
Incorporate that into the type system.
---
 src/command_handler.rs |  4 +++-
 src/steinhart_hart.rs  | 14 ++++++++------
 2 files changed, 11 insertions(+), 7 deletions(-)

diff --git a/src/command_handler.rs b/src/command_handler.rs
index 3f30083..11e1595 100644
--- a/src/command_handler.rs
+++ b/src/command_handler.rs
@@ -34,11 +34,13 @@ use uom::{
             ElectricPotential,
             ElectricalResistance,
             ThermodynamicTemperature,
+            TemperatureInterval,
         },
         electric_current::ampere,
         electric_potential::volt,
         electrical_resistance::ohm,
         thermodynamic_temperature::degree_celsius,
+        temperature_interval::kelvin,
     },
 };
 
@@ -243,7 +245,7 @@ impl Handler {
         use super::command_parser::ShParameter::*;
         match parameter {
             T0 => sh.t0 = ThermodynamicTemperature::new::<degree_celsius>(value),
-            B => sh.b = value,
+            B => sh.b = TemperatureInterval::new::<kelvin>(value),
             R0 => sh.r0 = ElectricalResistance::new::<ohm>(value),
         }
         send_line(socket, b"{}");
diff --git a/src/steinhart_hart.rs b/src/steinhart_hart.rs
index 87149df..4e05220 100644
--- a/src/steinhart_hart.rs
+++ b/src/steinhart_hart.rs
@@ -3,10 +3,12 @@ use uom::si::{
     f64::{
         ElectricalResistance,
         ThermodynamicTemperature,
+        TemperatureInterval,
     },
     electrical_resistance::ohm,
     ratio::ratio,
     thermodynamic_temperature::{degree_celsius, kelvin},
+    temperature_interval::kelvin as kelvin_interval,
 };
 use serde::{Deserialize, Serialize};
 
@@ -15,17 +17,17 @@ use serde::{Deserialize, Serialize};
 pub struct Parameters {
     /// Base temperature
     pub t0: ThermodynamicTemperature,
-    /// Base resistance
+    /// Resistance at base temperature
     pub r0: ElectricalResistance,
     /// Beta
-    pub b: f64,
+    pub b: TemperatureInterval,
 }
 
 impl Parameters {
-    /// Perform the voltage to temperature conversion.
+    /// Perform the resistance to temperature conversion.
     pub fn get_temperature(&self, r: ElectricalResistance) -> ThermodynamicTemperature {
-        let inv_temp = 1.0 / self.t0.get::<kelvin>() + (r / self.r0).get::<ratio>().ln() / self.b;
-        ThermodynamicTemperature::new::<kelvin>(1.0 / inv_temp)
+        let temp = (self.t0.recip() + (r / self.r0).get::<ratio>().ln() / self.b).recip();
+        ThermodynamicTemperature::new::<kelvin>(temp.get::<kelvin_interval>())
     }
 }
 
@@ -34,7 +36,7 @@ impl Default for Parameters {
         Parameters {
             t0: ThermodynamicTemperature::new::<degree_celsius>(25.0),
             r0: ElectricalResistance::new::<ohm>(10_000.0),
-            b: 3800.0,
+            b: TemperatureInterval::new::<kelvin_interval>(3800.0),
         }
     }
 }