diff --git a/STATUS.md b/STATUS.md
index 1ef99cd7..4e5acc7a 100644
--- a/STATUS.md
+++ b/STATUS.md
@@ -189,7 +189,7 @@ Floating point mathematical routines - not official, but well known and broadly
 | `$B794`   | `convert_A_to_FAC1`    | NOT DONE |                                                    |
 | `$B7B5`   | `STRVAL`               | NOT DONE | imports string to FAC1                             |
 | `$B7F7`   | `convert_FAC1_to_ADDR` | NOT DONE |                                                    |
-| `$B849`   | `add_HALF_FAC1`        | NOT DONE |                                                    |
+| `$B849`   | `add_HALF_FAC1`        | DONE     |                                                    |
 | `$B850`   | `sub_MEM_FAC1`         | DONE     |                                                    |
 | `$B853`   | `sub_FAC2_FAC1`        | DONE     |                                                    |
 | `$B862`   | `add_align_exponents`  | DONE     |                                                    |
@@ -203,7 +203,7 @@ Floating point mathematical routines - not official, but well known and broadly
 | `$BA8C`   | `mov_MEM_FAC2`         | DONE     |                                                    |
 | `$BA90`   | `get_FAC2_via_INDEX`   | NOT DONE |                                                    |
 | `$BAE2`   | `mul10_FAC1`           | NOT DONE |                                                    |
-| `$BAFE`   | `div10_FAC1_p`         | NOT DONE |                                                    |
+| `$BAFE`   | `div10_FAC1_p`         | PARTIAL  | needs `div_FAC2_FAC1`                              |
 | `$BB0F`   | `div_MEM_FAC1`         | PARTIAL  | needs `div_FAC2_FAC1`                              |
 | `$BB12`   | `div_FAC2_FAC1`        | NOT DONE | some work already started                          |
 | `$BBA2`   | `mov_MEM_FAC1`         | DONE     |                                                    |
@@ -216,7 +216,7 @@ Floating point mathematical routines - not official, but well known and broadly
 | `$BC0C`   | `mov_r_FAC1_FAC2`      | PARTIAL  | not fully tested yet                               |
 | `$BC0F`   | `mov_FAC1_FAC2`        | DONE     |                                                    |
 | `$BC1B`   | `round_FAC1`           | PARTIAL  | not fully tested yet                               |
-| `$BC2B`   | `sgn_FAC1_A`           | NOT DONE |                                                    |
+| `$BC2B`   | `sgn_FAC1_A`           | DONE     |                                                    |
 | `$BC39`   | `sgn_FAC1`             | NOT DONE |                                                    |
 | `$BC3C`   | `convert_A_to_FAC1`    | NOT DONE |                                                    |
 | `$BC44`   | `convert_i16_to_FAC1`  | NOT DONE |                                                    |
@@ -227,15 +227,15 @@ Floating point mathematical routines - not official, but well known and broadly
 | `$BCF3`   | `FIN`                  | NOT DONE | imports string to FAC1, ignores spaces             |
 | `$BD7E`   | `FINLOG`               | NOT DONE |                                                    |
 | `$BDDD`   | `FOUT`                 | NOT DONE | outputs FAC1 to string at $0100                    |
-| `$BF71`   | `sqr_FAC1`             | NOT DONE |                                                    |
+| `$BF71`   | `sqr_FAC1`             | PARTIAL  | needs `sqr_FAC2`                                   |
 | `$BF74`   | `sqr_FAC2`             | NOT DONE |                                                    |
-| `$BF78`   | `FPWR`                 | NOT DONE | set FAC1 as RAM raised to the power of FAC1        |
-| `$BF7B`   | `FPWRT`                | NOT DONE | set FAC1 as FAC2 raised to the power of FAC1       |
+| `$BF78`   | `pwr_FAC2_MEM`         | PARTIAL  | needs `pwr_FAC2_FAC1`                              |
+| `$BF7B`   | `pwr_FAC2_FAC1`        | NOT DONE |                                                    |
 | `$BFB4`   | `toggle_sign_FAC1`     | DONE     |                                                    |
 | `$BFED`   | `exp_FAC1`             | NOT DONE |                                                    |
 | `$E043`   | `poly1_FAC1`           | NOT DONE |                                                    |
 | `$E059`   | `poly2_FAC1`           | NOT DONE |                                                    |
-| `$E097`   | `rnd_FAC1`             | PARTIAL  | requires `sgn_FAC1_A`                              |
+| `$E097`   | `rnd_FAC1`             | DONE     |                                                    |
 | `$E09A`   | `rnd_A`                | DONE     |                                                    |
 | `$E0BE`   | `rnd_generate`         | DONE     |                                                    |
 | `$E264`   | `cos_FAC1`             | NOT DONE |                                                    |
diff --git a/c64/basic/,stubs_math/bf7b.FPWRT.s b/c64/basic/,stubs_math/bf7b.pwr_FAC2_FAC1.s
similarity index 97%
rename from c64/basic/,stubs_math/bf7b.FPWRT.s
rename to c64/basic/,stubs_math/bf7b.pwr_FAC2_FAC1.s
index 1bc6264a..6bb1a7f6 100644
--- a/c64/basic/,stubs_math/bf7b.FPWRT.s
+++ b/c64/basic/,stubs_math/bf7b.pwr_FAC2_FAC1.s
@@ -21,5 +21,5 @@
 
 // XXX provide implementation
 
-FPWRT:
+pwr_FAC2_FAC1:
 	STUB_IMPLEMENTATION()
diff --git a/c64/basic/,stubs_math/b849.add_HALF_FAC1.s b/c64/basic/math/b849.add_HALF_FAC1.s
similarity index 82%
rename from c64/basic/,stubs_math/b849.add_HALF_FAC1.s
rename to c64/basic/math/b849.add_HALF_FAC1.s
index deac2795..76f21ff0 100644
--- a/c64/basic/,stubs_math/b849.add_HALF_FAC1.s
+++ b/c64/basic/math/b849.add_HALF_FAC1.s
@@ -9,8 +9,9 @@
 // See also:
 // - [CM64] Computes Mapping the Commodore 64 - page 112
 
-// XXX implement
-
 add_HALF_FAC1:
 
-	STUB_IMPLEMENTATION()
+	ldy #>const_HALF
+	lda #<const_HALF
+
+	jmp add_MEM_FAC1
diff --git a/c64/basic/,stubs_math/b849.interop b/c64/basic/math/b849.interop
similarity index 100%
rename from c64/basic/,stubs_math/b849.interop
rename to c64/basic/math/b849.interop
diff --git a/c64/basic/,stubs_math/bafe.div10_FAC1_p.s b/c64/basic/math/bafe.div10_FAC1_p.s
similarity index 84%
rename from c64/basic/,stubs_math/bafe.div10_FAC1_p.s
rename to c64/basic/math/bafe.div10_FAC1_p.s
index 0be04ac9..0c1af809 100644
--- a/c64/basic/,stubs_math/bafe.div10_FAC1_p.s
+++ b/c64/basic/math/bafe.div10_FAC1_p.s
@@ -16,7 +16,14 @@
 // - https://codebase64.org/doku.php?id=base:kernal_floating_point_mathematics
 //
 
-// XXX provide implementation
+// XXX test this
 
 div10_FAC1_p:
-	STUB_IMPLEMENTATION()
+
+	lda #$00
+	sta FAC1_sign
+
+	ldy #>const_TEN
+	lda #<const_TEN
+
+	jmp div_MEM_FAC1
diff --git a/c64/basic/,stubs_math/bafe.interop b/c64/basic/math/bafe.interop
similarity index 100%
rename from c64/basic/,stubs_math/bafe.interop
rename to c64/basic/math/bafe.interop
diff --git a/c64/basic/,stubs_math/bc2b.interop b/c64/basic/math/bc2b.interop
similarity index 100%
rename from c64/basic/,stubs_math/bc2b.interop
rename to c64/basic/math/bc2b.interop
diff --git a/c64/basic/,stubs_math/bc2b.sgn_FAC1_A.s b/c64/basic/math/bc2b.sgn_FAC1_A.s
similarity index 87%
rename from c64/basic/,stubs_math/bc2b.sgn_FAC1_A.s
rename to c64/basic/math/bc2b.sgn_FAC1_A.s
index 866ac209..66f449e6 100644
--- a/c64/basic/,stubs_math/bc2b.sgn_FAC1_A.s
+++ b/c64/basic/math/bc2b.sgn_FAC1_A.s
@@ -15,7 +15,12 @@
 // - https://codebase64.org/doku.php?id=base:kernal_floating_point_mathematics
 //
 
-// XXX provide implementation
-
 sgn_FAC1_A:
-	STUB_IMPLEMENTATION()
+
+	lda FAC1_exponent
+	beq !+
+	lda FAC1_sign
+	bne !+
+	lda #$01
+!:
+	rts
diff --git a/c64/basic/,stubs_math/bf71.interop b/c64/basic/math/bf71.interop
similarity index 100%
rename from c64/basic/,stubs_math/bf71.interop
rename to c64/basic/math/bf71.interop
diff --git a/c64/basic/,stubs_math/bf71.sqr_FAC1.s b/c64/basic/math/bf71.sqr_FAC1.s
similarity index 89%
rename from c64/basic/,stubs_math/bf71.sqr_FAC1.s
rename to c64/basic/math/bf71.sqr_FAC1.s
index ac9a07c3..a670c052 100644
--- a/c64/basic/,stubs_math/bf71.sqr_FAC1.s
+++ b/c64/basic/math/bf71.sqr_FAC1.s
@@ -15,7 +15,8 @@
 // - https://codebase64.org/doku.php?id=base:kernal_floating_point_mathematics
 //
 
-// XXX provide implementation
-
 sqr_FAC1:
-	STUB_IMPLEMENTATION()
+
+	jsr mov_FAC1_FAC2
+	
+	// FALLLTROUGH to sqr_FAC2
diff --git a/c64/basic/,stubs_math/bf78.interop b/c64/basic/math/bf78.interop
similarity index 100%
rename from c64/basic/,stubs_math/bf78.interop
rename to c64/basic/math/bf78.interop
diff --git a/c64/basic/,stubs_math/bf78.FPWR.s b/c64/basic/math/bf78.pwr_FAC2_MEM.s
similarity index 71%
rename from c64/basic/,stubs_math/bf78.FPWR.s
rename to c64/basic/math/bf78.pwr_FAC2_MEM.s
index 3225ae9c..8b785204 100644
--- a/c64/basic/,stubs_math/bf78.FPWR.s
+++ b/c64/basic/math/bf78.pwr_FAC2_MEM.s
@@ -4,7 +4,7 @@
 // #LAYOUT# *   *       #IGNORE
 
 //
-// Math package - set FAC1 as RAM raised to the power of FAC1
+// Math package - set FAC1 as RAM raised to the power of FAC2
 //
 // Input:
 // - .A - address low byte
@@ -12,11 +12,13 @@
 //
 // See also:
 // - [CM64] Computes Mapping the Commodore 64 - page 117 (FPWRT)
+// - https://www.c64-wiki.com/wiki/BASIC-ROM
 // - https://www.c64-wiki.com/wiki/Floating_point_arithmetic
 // - https://codebase64.org/doku.php?id=base:kernal_floating_point_mathematics
 //
 
-// XXX provide implementation
+pwr_FAC2_MEM:
 
-FPWR:
-	STUB_IMPLEMENTATION()
+	jsr mov_MEM_FAC1
+
+	// FALLTROUGH to pwr_FAC2_FAC1