implemented another batch of small floating point routiines

2020-02-27 17:03:55 +01:00
1 parent 4a510d78fd
commit 6f7e099004
@@ -0,0 +1,17 @@
+// #LAYOUT# STD *       #TAKE
+// #LAYOUT# X16 BASIC_0 #TAKE-OFFSET 2000
+// #LAYOUT# *   BASIC_0 #TAKE
+// #LAYOUT# *   *       #IGNORE
+
+//
+// Math package - add 0.5 to FAC1
+//
+// See also:
+// - [CM64] Computes Mapping the Commodore 64 - page 112
+
+add_HALF_FAC1:
+
+	ldy #>const_HALF
+	lda #<const_HALF
+
+	jmp add_MEM_FAC1
@@ -0,0 +1,4 @@
+This is an important floating point math routine, documented in many places.
+
+Example code utilizing this routine:
+- https://github.com/irmen/prog8/blob/master/compiler/res/prog8lib/c64flt.p8
@@ -0,0 +1,29 @@
+// #LAYOUT# STD *       #TAKE
+// #LAYOUT# X16 BASIC_0 #TAKE-OFFSET 2000
+// #LAYOUT# *   BASIC_0 #TAKE
+// #LAYOUT# *   *       #IGNORE
+
+//
+// Math package - divide FAC1 by 10, result is always positive
+//
+// Notes:
+// - shifts FAC1 to FAC2, loads FAC1 with the constant 10, falls through to the division routine
+
+//
+// See also:
+// - [CM64] Computes Mapping the Commodore 64 - page 114
+// - https://www.c64-wiki.com/wiki/Floating_point_arithmetic
+// - https://codebase64.org/doku.php?id=base:kernal_floating_point_mathematics
+//
+
+// XXX test this
+
+div10_FAC1_p:
+
+	lda #$00
+	sta FAC1_sign
+
+	ldy #>const_TEN
+	lda #<const_TEN
+
+	jmp div_MEM_FAC1
@@ -0,0 +1,5 @@
+This is an important floating point math routine, documented in many places.
+
+Example code utilizing this routine:
+- https://github.com/irmen/prog8/blob/master/compiler/res/prog8lib/c64flt.p8
+
@@ -0,0 +1,4 @@
+This is an important floating point math routine, documented in many places.
+
+Example code utilizing this routine:
+- https://github.com/irmen/prog8/blob/master/compiler/res/prog8lib/c64flt.p8
@@ -0,0 +1,26 @@
+// #LAYOUT# STD *       #TAKE
+// #LAYOUT# X16 BASIC_0 #TAKE-OFFSET 2000
+// #LAYOUT# *   BASIC_0 #TAKE
+// #LAYOUT# *   *       #IGNORE
+
+//
+// Math package - evaluate FAC1 sign, to .A
+//
+// Output:
+// - .A - 0 for FAC1 equal 0, 1 for FAC1 > 0, -1 for FAC1 < 0
+//
+// See also:
+// - [CM64] Computes Mapping the Commodore 64 - page 115
+// - https://www.c64-wiki.com/wiki/Floating_point_arithmetic
+// - https://codebase64.org/doku.php?id=base:kernal_floating_point_mathematics
+//
+
+sgn_FAC1_A:
+
+	lda FAC1_exponent
+	beq !+
+	lda FAC1_sign
+	bne !+
+	lda #$01
+!:
+	rts
@@ -0,0 +1,4 @@
+This is an important floating point math routine, documented in many places.
+
+Example code utilizing this routine:
+- https://github.com/irmen/prog8/blob/master/compiler/res/prog8lib/c64flt.p8
@@ -0,0 +1,22 @@
+// #LAYOUT# STD *       #TAKE
+// #LAYOUT# X16 BASIC_0 #TAKE-OFFSET 2000
+// #LAYOUT# *   BASIC_0 #TAKE
+// #LAYOUT# *   *       #IGNORE
+
+//
+// Math package - square root of FAC1
+//
+// Note:
+// - transfers FAC1 to FAC2 and falls through to SQR2
+//
+// See also:
+// - [CM64] Computes Mapping the Commodore 64 - page 117
+// - https://www.c64-wiki.com/wiki/Floating_point_arithmetic
+// - https://codebase64.org/doku.php?id=base:kernal_floating_point_mathematics
+//
+
+sqr_FAC1:
+
+	jsr mov_FAC1_FAC2
+	
+	// FALLLTROUGH to sqr_FAC2
@@ -0,0 +1,4 @@
+This is an important floating point math routine, documented in many places.
+
+Example code utilizing this routine:
+- https://github.com/irmen/prog8/blob/master/compiler/res/prog8lib/c64flt.p8
@@ -0,0 +1,24 @@
+// #LAYOUT# STD *       #TAKE
+// #LAYOUT# X16 BASIC_0 #TAKE-OFFSET 2000
+// #LAYOUT# *   BASIC_0 #TAKE
+// #LAYOUT# *   *       #IGNORE
+
+//
+// Math package - set FAC1 as RAM raised to the power of FAC2
+//
+// Input:
+// - .A - address low byte
+// - .Y - address high byte
+//
+// 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
+//
+
+pwr_FAC2_MEM:
+
+	jsr mov_MEM_FAC1
+
+	// FALLTROUGH to pwr_FAC2_FAC1