From 55163c2d1f5487e4d5213b256652527843ce002b Mon Sep 17 00:00:00 2001
From: Luke Kenneth Casson Leighton <lkcl@lkcl.net>
Date: Sat, 31 Jul 2021 19:38:07 +0100
Subject: [PATCH] add SINGLE function to double2single, to replace
 manually-created version

---
 openpower/isafunctions/double2single.mdwn | 32 +++++++++++++++++++++++
 1 file changed, 32 insertions(+)

diff --git a/openpower/isafunctions/double2single.mdwn b/openpower/isafunctions/double2single.mdwn
index 31686e4d..28e48e29 100644
--- a/openpower/isafunctions/double2single.mdwn
+++ b/openpower/isafunctions/double2single.mdwn
@@ -139,3 +139,35 @@ Round to Single-Precision instruction.
 
         return result
 
+<!-- Power ISA v3.0B p140 section 4.6.2 -->
+
+    def SINGLE(WORD):
+        exp <- [0] * 11
+        frac <- [0] * 53
+        sign <- 0b0
+        # Normalized Operand
+        if (WORD[1:8] >u 0) & (WORD[1:8] <u 255) then
+            FRT[0:1] <- WORD[0:1]
+            FRT[2] <- ¬WORD[1]
+            FRT[3] <- ¬WORD[1]
+            FRT[4] <- ¬WORD[1]
+            FRT[5:63] <- WORD[2:31] || [0]*29
+        # Denormalized Operand
+        if (WORD[1:8] = 0) & (WORD[9:31] != 0) then
+            sign <- WORD[0]
+            exp <- -126
+            frac[0:52] <- 0b0 || WORD[9:31] || [0]*29
+            #normalize the operand
+            do while frac[0] = 0
+                frac[0:52] <- frac[1:52] || 0b0
+                exp <- exp - 1
+            FRT[0] <- sign
+            FRT[1:11] <- exp + 1023
+            FRT[12:63] <- frac[1:52]
+        # Zero / Infinity / NaN
+        if (WORD[1:8] = 255) | (WORD[1:31] = 0) then
+            FRT[0:1] <- WORD[0:1]
+            FRT[2] <- WORD[1]
+            FRT[3] <- WORD[1]
+            FRT[4] <- WORD[1]
+            FRT[5:63] <- WORD[2:31] || [0]*29
-- 
2.30.2