Add MS BASIC equality test in similarity.c

Plus fixing mistakes

src/basic/math/b947.inv_FAC1.s

@@ -39,6 +39,7 @@ inv_FAC1_mantissa:
     lda FACOV
     eor #$FF
     sta FACOV
+    inc FACOV
     bne inc_FAC1_ret
     
     ; Fallthrough to inc_FAC1

src/basic/math/b983.shiftr_FAC.s

@@ -45,7 +45,7 @@ shiftr_FAC:
     sbc #8              ; Go back to previous value
     tay
     lda FACOV
-    bcs shiftr_FAC
+    bcs shftrt
 .sext:
     asl 1,X
     bcc @2

src/basic/math/bc1b.round_FAC1.jmp.s

@@ -5,5 +5,6 @@
 
 ; Just a compatibility jump - our routine is longer than original one
 
-	jmp round_FAC1
+round_FAC1:
+	jmp round_FAC1_impl
 	

src/basic/math/round_FAC1.s

@@ -12,7 +12,7 @@
 
 ; XXX test this
 
-round_FAC1:
+round_FAC1_impl:
 
 	; If number is 0, do not do anything
 

src/basic/math_mov/bbd0.mov_r_FAC1_VAR.s

@@ -20,7 +20,7 @@
 
 mov_r_FAC1_VAR:
 
-	ldx #FORPNT+0
-	ldy #FORPNT+1
+	ldx FORPNT+0
+	ldy FORPNT+1
 
 	; FALLTROUGH to mov_r_FAC1_MEM

src/basic/math_mov/bbfc.mov_FAC2_FAC1.s

@@ -32,7 +32,7 @@ mov_FAC2_FAC1:
     lda FAC2_sign
 
 mov_FAC2_FAC1_sign:
-    sta FAC2_sign
+    sta FAC1_sign
 
     ; Copy the mantissa and the exponent
 

tools/similarity.c

@@ -92,7 +92,11 @@ int main(int argc,char **argv)
   if (basic) {
     for (int i = 0; i < NUM_CHUNKS; i++) {
       for (int addr = microsoft_basic_chunks[i * 2]; addr <= microsoft_basic_chunks[i * 2 + 1]; addr++) {
-        f1[addr - 0xA000] = 0;
+        short index = addr - 0xA000;
+        if (f1[index] != f2[index]) {
+            printf("BASIC does not match MS BASIC at %X (%X != %X)\n", addr, f1[index], f2[index]);
+        }
+        f1[index] = 0;
       }
     }
   }