wip, f16 not working

Author: tgross35

src/math/generic/scalbn.rs

@@ -54,8 +54,8 @@ where
     // 2 ^ sig_total_bits, representation of what can be accounted for with subnormals
     let f_exp_subnorm = F::from_parts(false, sig_total_bits + F::EXP_BIAS, zero);
 
-    dbg!((Hexf(f_exp_max), Hexf(f_exp_min), Hexf(f_exp_subnorm)));
-
+    dbg!(exp_max, exp_min, sig_total_bits, sig_total_bits + F::EXP_BIAS);
+    dbg!(Hexf(f_exp_max), Hexf(f_exp_min), Hexf(f_exp_subnorm));
     dbg!(Hexf(x), n);
 
     // The goal is to multiply `x` by a scale factor that applies `n`. However, there are cases
@@ -86,13 +86,16 @@ where
         let add = -exp_min - sig_total_bits as i32;
         dbg!(Hexf(mul), add);
 
-        if n < exp_min + sig_total_bits as i32 {
+        if n < exp_min - sig_total_bits as i32 {
+            dbg!("block 1");
             x *= f_exp_min;
             n += -exp_min;
+            dbg!(Hexf(x), n);
 
-            if n < exp_min + sig_total_bits as i32 {
+            if n < exp_min - sig_total_bits as i32 {
                 x *= f_exp_min;
                 n += -exp_min;
+                dbg!(Hexf(x), n);
 
                 // if n < exp_min + sig_total_bits as i32 {
                 //     x *= f_exp_min;
@@ -105,9 +108,33 @@ where
 
                 if n < exp_min {
                     n = exp_min;
+                    dbg!(Hexf(x), n);
                 }
                 // }
+            } else if n < exp_min {
+                // only for f16
+
+                x *= mul;
+                n += add;
+
+                dbg!(Hexf(x), n);
+                if n < exp_min {
+                    n = exp_min;
+                    dbg!(Hexf(x), n);
+                }
             }
+        } else if n < exp_min {
+            // dbg!("block 2");
+            // // only for f16
+
+            // x *= mul;
+            // n += add;
+
+            // dbg!(Hexf(x), n);
+            // if n < exp_min {
+            //     n = exp_min;
+            //     dbg!(Hexf(x), n);
+            // }
         }
 
         // if n < exp_min {
@@ -262,28 +289,4 @@ mod tests {
     fn spec_test_f128() {
         spec_test::<f128>();
     }
-
-    // #[test]
-    // fn foobar32() {
-    //     let x = hf32!("0x1.fffffep+127");
-    //     let n = -2147483639;
-    //     scalbn(x, n);
-    //     std::eprintln!();
-    //     let x = hf32!("0x1.fffffep-126");
-    //     let n = 2147483639;
-    //     scalbn(x, n);
-    //     panic!();
-    // }
-
-    // #[test]
-    // fn foobar16() {
-    //     let x = hf16!("0x1.ffp+15");
-    //     let n = -2147483639;
-    //     scalbn(x, n);
-    //     std::eprintln!();
-    //     let x = hf16!("0x1.ffp-15");
-    //     let n = 2147483639;
-    //     scalbn(x, n);
-    //     panic!();
-    // }
 }