Blame | Last modification | View Log | Download | RSS feed
Tue Jul 17 11:00:22 2018 Page 12500 A.D. 65816 Macro Assembler #26960 - Version 5.02g-----------------------------------------------------Input Filename : src\FE\float.asmOutput Filename : obj\FE\float.objListing Has Been Relocated2585 .LIST on25862587 ;.INCLUDE inc\dirp01.inc2588 F8FFB1 .INCLUDE inc\dpcbm.inc2589 ;; cbm emulation25902591 ;;CBMP0 .EQU $7F0025922593 ; sezione COMMON -- questo permette di includere il file in piu' file25942595 DPCBM: .SECTION page0, ref_only, common25962597 000000 .ABSOLUTE ;; inizia sempre da $002598 000000 .ORG 0x0025992600 ;;MAXLF .EQU $10 ; max. logical files2601 00000B MAXLF .EQU 11 ; max. logical files2602 000008 MAXVIEC .EQU 8 ; max. virtual iec device #2603 000000 .DS 162604 0000102605 000010 00 subflg .DB ; ($10)2606 000011 00 .DB ; ($11)2607 000012 00 .DB ; ($12)2608 000013 00 .DB ; ($13)2609 000014 0000 lineno .DW ; ($14-$15) line number used by editor2610 000016 00 tmppt .DB ; ($16)2611 000017 00 tcnt .DB ; ($17) temp. counter/index2612 000018 0000 .DW ; ($18-$19)2613 00001A 0000 .DW ; ($1A-$1B)2614 00001C 0000 .DW ; ($1C-$1D)2615 00001E 0000 .DW ; ($1E-$1F)2616 000020 00 .DB ; ($20)26172618 ; basic/p-code virtual machine work area2619 000021 0000 .DW ; ($21-$22)2620 000023 0000 .DW ; ($23-$24)2621 000025 0000 .DW ; ($25-$26)2622 000027 0000 .DW ; ($27-$28)2623 000029 0000 .DW ; ($29-$2A)2624 00002B 0000 .DW ; ($2B-$2C)2625 00002D 0000 .DW ; ($2D-$2E)2626 00002F 0000 arytab .DW ; ($2F-$30)2627 000031 0000 strend .DW ; ($31-$32)2628 000033 0000 freetop .DW ; ($33-$34)2629 000035 0000 freespc .DW ; ($35-$36)2630 000037 0000 memsiz .DW ; ($37-$38)Tue Jul 17 11:00:22 2018 Page 22631 000039 0000 curline .DW ; ($39-$3A) current line number in execution2632 00003B 0000 oldline .DW ; ($3B-$3C) previous line number in execution2633 00003D 0000 oldtxt .DW ; ($3D-$3E) saved pointer to text in basic execution2634 00003F 0000 datlin .DW ; ($3F-$40)2635 000041 0000 datptr .DW ; ($41-$42)2636 000043 0000 inpptr .DW ; ($43-$44)2637 000045 0000 varnam .DW ; ($45-$46)2638 000047 0000 varpnt .DW ; ($47-$48)2639 000049 0000 forpnt .DW ; ($49-$4A) basic2640 00004B 0000 opptr .DW ; ($4B-$4C) basic2641 00004D 00 opmask .DB ; ($4D)2642 00004E 0000 defpnt .DW ; ($4E-$4F)2643 000050 0000 dscpnt .DW ; ($50-$51)2644 000052 00 dsclen .DB ; ($52)2645 000053 00 prompt .DB ; ($53) prompt flag/CMD channel - editor/basic ($13)26462647 ; equates for p-code virtual machine2648 00003B prtrap .EQU oldline ; ($3B) p-code trap errors flag2649 00003C pcstack .EQU oldline+1 ; ($3C-$3D) p-code save stack2650 00003E prterr .EQU oldtxt+1 ; ($3E) p-code run time error number2651 00003F psrline .EQU datlin ; ($3F-$40) pascal source line number2652 000041 prseed .EQU datptr ; ($41-$42) p-code random generator2653 000043 pbrkcnt .EQU inpptr ; ($43) p-code counter for break check2654 000044 tmpy .EQU inpptr+1 ; ($44) p-code temp. for save Y2655 000045 inqidx .EQU varnam ; ($45) p-code index to input buffer2656 000046 inqcnt .EQU varnam+1 ; ($46) p-code input buffer counter2657 000047 tmpv1 .EQU varpnt ; ($47) p-code temp.2658 000048 tmpv2 .EQU varpnt+1 ; ($48) p-code temp.2659 000049 tmpa .EQU forpnt ; ($49) p-code temp. save A2660 00004A tmpwa .EQU forpnt+1 ; ($4A) p-code tmp. save byte to write2661 00004B tmpx .EQU opptr ; ($4B) p-code tmp. save X when print2662 00004C pcterr .EQU opptr+1 ; ($4C) p-code compiler error flag2663 00004D tmpdsk .EQU opmask ; ($4D) p-code default temp. disk unit2664 00004E defprn .EQU defpnt ; ($4E) p-code default printer #device2665 00004F pcdbg .EQU defpnt+1 ; ($4F) p-code debug run after error2666 000050 pcptr .EQU dscpnt ; ($50-$51) pointer to p-code (debug)2667 000052 autolf .EQU dsclen ; ($52) p-code auto line feed flag26682669 000054 tmpwrk .DS 4 ; ($54-$57) temp. area editor/virtual iec2670 000058 00 tmpb0 .DB ; ($58)2671 000059 00 argbits .DB ; ($59) FAC #1 guard bits2672 00005A 00 tmpb2 .DB ; ($5A) tmp. byte - used by fpu2673 00005B 00 txtlen .DB ; ($5B) tmp. word length of text to insert in buffer2674 00005C 00 tmpb1 .DB ; ($5C) tmp. byte used by editor and float point unit2675 00005D 0000 tmpval .DW ; ($5D-$5E) temp. word used by editor/fpu2676 00005F 0000 tlnptr .DW ; ($5F-$60) temp. line text pointer used by editor/fpu2677 000061 00 facexp .DB ; ($61) FAC #1 exponent2678 000062 facm .DS 4 ; ($62-$65) FAC #1 mantissa2679 000066 00 facsgn .DB ; ($66) FAC #1 sign2680 000067 00 sercnt .DB ; ($67) FAC #1 series evaluation2681 000068 00 facov .DB ; ($68) FAC #1 overflow2682 000069 00 argexp .DB ; ($69) FAC #2 exponent2683 00006A argm .DS 4 ; ($6A-$6D) FAC #2 mantissa2684 00006E 00 argsgn .DB ; ($6E) FAC #2 sign2685 00006F 00 arisgn .DB ; ($6F) FAC #1 & #2 sign comparison2686 000070 00 facbits .DB ; ($70) FAC #1 guard bits2687 000071 0000 fbufpt .DW ; ($71-$72) FAC series evaluation pointerTue Jul 17 11:00:22 2018 Page 32688 000073 00 fbufpt2 .DB ; ($73) FAC series evaluation bank (progr. bank)26892690 ; tmp. variables (equates) used by float point unit2691 000054 factmp .EQU tmpwrk ; ($54-$57) used by fpu mult/div2692 00005A cprmsk .EQU tmpb2 ; ($5A) used by trigs. funcs.2693 00005B faclsb .EQU txtlen ; ($5B) used by float point unit2694 00005D dccnt .EQU tmpval ; ($5D) digit count after decimal dot2695 00005E expval .EQU tmpval+1 ; ($5E) exponent value when convert2696 00005F dotfg .EQU tlnptr ; ($5F) flag decimal dot when convert2697 000060 expsgn .EQU tlnptr+1 ; ($60) flag exponent sign when conver2698 000067 sgnfg .EQU sercnt ; ($67) flag: minus sign when convert2699 000071 decidx .EQU fbufpt ; ($71) string index when convert2700 000072 s2fer .EQU fbufpt+1 ; ($72) string to float flag error2701 000073 decidx2 .EQU fbufpt+2 ; ($73) used when convert float to str27022703 000074 0000 s2fstk .DW ; ($74-$75) save stack when string to float conversion2704 000076 longp LP ; ($76-$77-$78) temp. long pointer2705 000079 00 s2fptr .DB ; ($79) flag trap overflow err. string to float conv.27062707 00007A 0000 txtp .DW ; ($7A-$7B) current pointer to text2708 00007C 00 pcdsk .DB ; ($7C) flag for pascal compiler2709 ; <7> : compile in memory27102711 00007D iecrp LP ; ($7D) iec: pointer to read buffer for include file2712 000080 00 iecalt .DB ; ($80) iec: flag read from include file2713 000081 iectp LP ; ($81) iec: pointer to read main text buffer2714 000084 00 iecndx .DB ; ($84) iec: index main text buffer2715 000085 00 ieceof .DB ; ($85) iec: end of file flag main text buffer2716 000086 iecwp LP ; ($86) iec: pointer to write buffer2717 000089 00 iecfg .DB ; ($89) iec: flag main text buffer2718 00008A 0000 iecndx2 .DW ; ($8A) iec: index buffer include file2719 00008C 00 iecbuf .DB ; ($8C) iec: byte buffer2720 00008D 00 svbnk .DB ; ($8D) save vbnk when run pascal compiler2721 00008E 00 ieceof2 .DB ; ($8E) iec: flag eof include file2722 00008F 00 vxflag .DB ; ($8F) flag valid text in x-bank2723 000090 00 kstat .DB ; ($90) CBM emulated IEC status2724 000091 00 errno .DB ; ($91) saved error number2725 000092 00 keyfg .DB ; ($92) keyboard flag: KbdFlag2726 000093 00 keyctl .DB ; ($93) flag control key from get key2727 000094 00 keysh .DB ; ($94) flag keyboard: KbdShift2728 000095 00 keytg .DB ; ($95) flag keyboard: KbdToggle2729 000096 00 basic .DB ; ($96) flag current 'basic' mode2730 ; <7>: flag basic mode2731 ; <6>: flag pascal compiler running2732 000097 00 runf .DB ; ($97) flag current run mode2733 000098 00 ldtnd .DB ; ($98) number of open files2734 000099 00 dfltn .DB ; ($99) default input device2735 00009A 00 dflto .DB ; ($9A) default output device2736 00009B 00 skpbk .DB ; ($9B) skip blank flag2737 00009C 00 w2sf .DB ; ($9C) integer to string flag conversion2738 00009D 00 msgflg .DB ; ($9D) flag message2739 00009E 00 autof .DB ; ($9E) flag auto-increment line numbering2740 00009F 0000 wauto .DW ; ($9F-$A0) auto-increment for line numbering2741 0000A1 00 tbnk .DB ; ($A1) bank that hold basic/pascal text2742 0000A2 00 vbnk .DB ; ($A2) bank that hold basic var/pascal pcode-stack2743 0000A3 00 cbnk .DB ; ($A3) compiler bank/temp. bank for exchange2744 0000A4 00 jbnk .DB ; ($A4) bank for obj, when compiling (3 bank max.)Tue Jul 17 11:00:22 2018 Page 42745 0000A5 00 xbnk .DB ; ($A5) exchange bank for text buffer2746 0000A6 00 dbnk .DB ; ($A6) bank for ram disk (jbnk + 3)2747 0000A7 0000 botmem .DW ; ($A7-$A8) start of memory on the text bank2748 0000A9 0000 topmem .DW ; ($A9-$AA) top of memory on the text bank2749 0000AB 0000 txtbeg .DW ; ($AB-$AC) begin of text buffer2750 0000AD 0000 txtend .DW ; ($AD-$AE) end of text buffer + 12751 0000AF 0000 txttop .DW ; ($AF-$B0) top of memory + 1 reserved to text buffer2752 0000B1 0000 varbeg .DW ; ($B1-$B2) start of basic var/pcode2753 0000B3 0000 vartop .DW ; ($B3-$B4) top of basic var/pcode + 12754 0000B5 0000 xtop .DW ; ($B5-$B6) top of text + 1 in x-bank2755 0000B7 00 fnlen .DB ; ($B7) length of current filename2756 0000B8 00 lacur .DB ; ($B8) current logical file number2757 0000B9 00 sacur .DB ; ($B9) current secondary address2758 0000BA 00 facur .DB ; ($BA) current device number2759 0000BB 0000 fnadr .DW ; ($BB-$BC) pointer to current filename2760 0000BD 00 fnadr2 .DB ; ($BD) - bank that hold file name2761 0000BE 00 objdsk .DB ; ($BE) compiler: option '1' specified2762 0000BF 00 truedsk .DB ; ($BF) true disk for load pascal soource file2763 0000C0 0000 curstk .DW ; ($C0-$C1) current stack when in edit mode2764 0000C2 00 pcodef .DB ; ($C2) flag p-code present for execution27652766 0000C3 latb .DS MAXLF ; logical files table2767 0000CE fatb .DS MAXLF ; devices numbers table2768 0000D9 satb .DS MAXLF ; secondary addresses table27692770 ;=============================================================27712772 0000E4 0000 cmdlin .DW ; command line ptr for pscal program2773 0000E6 00 fpcflag .DB ; fpc flags (new style)2774 ; <7>: new style compiler2775 ; <6>: no save exec. file2776 ; <5>: save obj. file27772778 0000E6 .RELATIVE2779 .ENDS27802781 ; basic work area starting at $FF00 of bank that hold text buffer2782 CBMHIM: .SECTION ref_only, common2783 000000 .ABSOLUTE2784 00FF00 .ORG $FF002785 ; .ORG 027862787 0000A2 TBUFMAX .EQU $A2 ; size of input buffer2788 00AA55 VALIDSGN .EQU $AA55 ; valid signature27892790 00FF00 cbmwork .DS 4 ; first 4 bytes free for insert text line2791 00FF04 TXTBUF .DS TBUFMAX ; input buffer for text insertion2792 00FFA6 NUMSTR .DS 32 ; room for numeric string27932794 00FFC6 0000 TMPPTR .DW ; temp. pointer2795 00FFC8 0000 .DW2796 00FFCA 0000 .DW2797 00FFCC 0000 .DW2798 00FFCE 0000 .DW2799 00FFD0 GFNAME .DS 32 ; global file name2800 00FFF0 .DS 82801 00FFF8Tue Jul 17 11:00:22 2018 Page 52802 ;;.ORG $FFF82803 00FFF8 0000 TXTEND2 .DW ; saved text end + 12804 00FFFA 0000 TXTVSGN .DW ; valid text buffer signature2805 00FFFC 0000 XTOP2 .DW ; saved text end + 1 in exchange bank2806 00FFFE 0000 XBVSGN .DW ; valid x-bank buffer signature28072808 00FFFE .RELATIVE2809 .ENDS28102811 .COMMENT @2812 cbmwork .EQU $FF00 ; first 4 bytes free for insert text line2813 TXTBUF .EQU $FF04 ; input buffer for text insertion2814 NUMSTR .EQU $FFA6 ; room for numeric string28152816 TMPPTR .EQU $FFC6 ; temp. pointer28172818 TXTEND2 .EQU $FFF8 ; saved text end + 12819 TXTVSGN .EQU $FFFA ; valid text buffer signature2820 XTOP2 .EQU $FFFC ; saved text end + 1 in exchange bank2821 XBVSGN .EQU $FFFE ; valid x-bank buffer signature2822 @28232824 ; basic/pascal work area starting at $FF00 of bank that hold variables/p-code2825 CBMHIV: .SECTION common, ref_only2826 000000 .ABSOLUTE2827 00FF00 .ORG $FF002828 ; .ORG 028292830 00005A QINSIZ .EQU $5A ; input buffer max. 90 chars.2831 000020 NUMSIZ .EQU $20 ; size of numeric buffer for conversion2832 000018 MAX2LF .EQU ((MAXLF + 1) * 2) ; size of p-code logical files table2833 000013 MAXFNLEN .EQU 19 ; max. length of file name '@x:yyyyyyyyyyyyyyyy'28342835 00FF00 MINBUF .DS QINSIZ2836 00FF5A MNUMSTR .DS NUMSIZ2837 00FF7A 0000 MTMPW0 .DW2838 00FF7C 0000 MTMPW1 .DW2839 00FF7E 00 MTMPB0 .DB2840 00FF7F MFNAME .DS 16 ; file name buffer2841 00FF8F MTMPNAM .DS 15 ; tmp. file name2842 00FF9E MPOW10 .DS 10 ; power of 10 table2843 00FFA8 MTRUE .DS 62844 00FFAE MFALSE .DS 62845 00FFB4 MSTDF0 .DS 4 ; descriptor for file #0 (consolle)2846 00FFB8 MSTDF1 .DS 4 ; descriptor for file #1 (disk command chann.)2847 00FFBC MSTDF2 .DS 4 ; descriptor for file #2 (used by compiler)2848 00FFC0 MMSET .DS 16 ; set variables2849 00FFD0 MFILET .DS MAX2LF ; table of logical files: #0 consolle, #1 disk2850 00FFE8 FAC1TMP .DS 5 ; these 3 used by float routines2851 00FFED FAC2TMP .DS 52852 00FFF2 FAC3TMP .DS 52853 00FFF7 FAC4TMP .DS 5 ; free to use2854 00FFFC 0000 MTPTR .DW ; used by new/dispose2855 00FFFE 0000 MTMPW2 .DW28562857 00FFFE .RELATIVE2858 .ENDSTue Jul 17 11:00:22 2018 Page 628592860 .COMMENT @2861 MINBUF .EQU $FF002862 MNUMSTR .EQU $FF5A2863 MTMPW0 .EQU $FF7A2864 MTMPW1 .EQU $FF7C2865 MTMPB0 .EQU $FF7E2866 MFNAME .EQU $FF7F ; file name buffer2867 MTMPNAM .EQU $FF8F ; tmp. file name2868 MPOW10 .EQU $FF9E ; power of 10 table2869 MTRUE .EQU $FFA82870 MFALSE .EQU $FFAE2871 MSTDF0 .EQU $FFB4 ; descriptor for file #0 (consolle)2872 MSTDF1 .EQU $FFB8 ; descriptor for file #1 (disk command chann.)2873 MSTDF2 .EQU $FFBC ; descriptor for file #2 (used by compiler)2874 MMSET .EQU $FFC0 ; set variables2875 MFILET .EQU $FFD0 ; table of logical files: #0 consolle, #1 disk2876 FAC1TMP .EQU $FFE8 ; these 3 used by float routines2877 FAC2TMP .EQU $FFED2878 FAC3TMP .EQU $FFF22879 FAC4TMP .EQU $FFF7 ; free to use2880 MTPTR .EQU $FFFC ; used by new/dispose2881 MTMPW2 .EQU $FFFE2882 @28832884 ; temp. file names in x-bank2885 00FF00 PASSRC .EQU $FF002886 00FF20 OBJDST .EQU $FF202887 00FF40 EXEDST .EQU $FF402888 00FF60 OBJDEL .EQU $FF60288928902891 00FF5A TTNUMSTR .EQU MNUMSTR28922893 .CODEFE2894 .LONGA off2895 .LONGI off28962897 .EXTERN cbm_chrget, cbm_chrgot, cbm_err, vpmerr2898 .GLOBAL XBBD4, XBBA2, XBC5B, XBA8C, XBB12, XBA2B, XB853, XB86A2899 .GLOBAL XBAFE, XBAE2, XBCF3P, XBCF3, XBBA2T,XBDDD, XE26B, XE2642900 .GLOBAL XE2B4, XE30E, XBFED, XBF7B, XBF71, XB9EA, XBC49, XBB0F2901 .GLOBAL XB9BC, XBBC7, XBBA2T, XB9E0, FCON05, XB867, FCON102902 .GLOBAL XBBA2L, XBA8CK, XBBA2K, XLN10, XBC0C, XB9BC, XBA282903 FE938D2904 ;---------------------------------------------------------------------------2905 ; move FAC #2 routines2906 ;---------------------------------------------------------------------------29072908 ; move float number at address AY to FAC #22909 ; address is in current data bank2910 FE938D XBA8C:2911 FE938D 85 76 sta <longp2912 FE938F 84 77 sty <longp+12913 FE9391 A0 04 ldy #42914 FE9393 B1 76 lda (longp),y2915 FE9395 85 6D sta <argm+3Tue Jul 17 11:00:22 2018 Page 72916 FE9397 88 dey2917 FE9398 B1 76 lda (longp),y2918 FE939A 85 6C sta <argm+22919 FE939C 88 dey2920 FE939D B1 76 lda (longp),y2921 FE939F 85 6B sta <argm+12922 FE93A1 88 dey2923 FE93A2 B1 76 lda (longp),y2924 FE93A4 85 6E sta <argsgn2925 FE93A6 45 66 eor <facsgn2926 FE93A8 85 6F sta <arisgn2927 FE93AA A5 6E lda <argsgn2928 FE93AC 09 80 ora #$802929 FE93AE 85 6A sta <argm2930 FE93B0 88 dey2931 FE93B1 B1 76 lda (longp),y2932 FE93B3 85 69 sta <argexp2933 FE93B5 A5 61 lda <facexp2934 FE93B7 60 rts29352936 ; move float number at address AY in current program bank to FAC #22937 FE93B8 XBA8CK:2938 FE93B8 85 76 sta longp2939 FE93BA 84 77 sty longp+12940 FE93BC A2 FE ldx #^XBA8CK2941 FE93BE 80 04 bra XBA8CK229422943 ; move float number at address AYX to FAC #22944 FE93C0 XBA8CL:2945 FE93C0 85 76 sta <longp2946 FE93C2 84 77 sty <longp+129472948 FE93C4 XBA8CK2:2949 FE93C4 86 78 stx <longp+22950 FE93C6 A0 04 ldy #42951 FE93C8 B7 76 lda [longp],y2952 FE93CA 85 6D sta <argm+32953 FE93CC 88 dey2954 FE93CD B7 76 lda [longp],y2955 FE93CF 85 6C sta <argm+22956 FE93D1 88 dey2957 FE93D2 B7 76 lda [longp],y2958 FE93D4 85 6B sta <argm+12959 FE93D6 88 dey2960 FE93D7 B7 76 lda [longp],y2961 FE93D9 85 6E sta <argsgn2962 FE93DB 45 66 eor <facsgn2963 FE93DD 85 6F sta <arisgn2964 FE93DF A5 6E lda <argsgn2965 FE93E1 09 80 ora #$802966 FE93E3 85 6A sta <argm2967 FE93E5 88 dey2968 FE93E6 B7 76 lda [longp],y2969 FE93E8 85 69 sta <argexp2970 FE93EA A5 61 lda <facexp2971 FE93EC XB848:2972 FE93EC 60 rtsTue Jul 17 11:00:22 2018 Page 829732974 ;---------------------------------------------------------------------------2975 ; subctraction implementation2976 ;---------------------------------------------------------------------------29772978 XB849: ; compute FAC #1 + 0.52979 FE93ED A9 AC lda #<FCON052980 FE93EF A0 9C ldy #>FCON052981 FE93F1 80 16 bra XB867 ; add costant stored in current program bank2982 FE93F32983 XB850: ; compute [float at address AYK] - FAC #12984 FE93F3 20 B8 93 jsr XBA8CK ; move float at AYK address to FAC #229852986 XB853: ; compute FAC #2 - FAC #12987 FE93F6 A5 66 lda <facsgn2988 FE93F8 49 FF eor #$FF2989 FE93FA 85 66 sta <facsgn2990 FE93FC 45 6E eor <argsgn2991 FE93FE 85 6F sta <arisgn2992 FE9400 A5 61 lda <facexp2993 FE9402 80 08 bra XB86A29942995 ;---------------------------------------------------------------------------2996 ; addition implementation2997 ;---------------------------------------------------------------------------29982999 FE9404 XB862:3000 FE9404 20 6B 95 jsr XB999 ; shift right mantissa3001 FE9407 90 3C bcc XB8A3 ; always return CF = 030023003 ; add float at address AYK to FAC #13004 FE9409 XB867:3005 FE9409 20 B8 93 jsr XBA8CK ; move float in mem. to FAC #230063007 ; add FAC #1 to FAC #2 - in input A = FAC #1 exponent3008 FE940C XB86A:3009 FE940C D0 03 bne ?023010 FE940E 4C BE 97 jmp XBBFC ; move FAC #2 to FAC #13011 FE9411 ?02:3012 FE9411 A6 70 ldx <facbits3013 FE9413 86 59 stx <argbits3014 FE9415 A2 69 ldx #argexp3015 FE9417 A5 69 lda <argexp3016 FE9419 XB877:3017 FE9419 A8 tay3018 FE941A F0 D0 beq XB848 ; exit (FAC #2 = 0)3019 FE941C 38 sec3020 FE941D E5 61 sbc <facexp3021 FE941F F0 24 beq XB8A3 ; same exponent3022 FE9421 90 12 bcc ?023023 FE9423 84 61 sty <facexp3024 FE9425 A4 6E ldy <argsgn3025 FE9427 84 66 sty <facsgn3026 FE9429 49 FF eor #$FF3027 FE942B 69 00 adc #$003028 FE942D A0 00 ldy #$003029 FE942F 84 59 sty <argbitsTue Jul 17 11:00:22 2018 Page 93030 FE9431 A2 61 ldx #facexp3031 FE9433 D0 04 bne ?043032 FE9435 A0 00 ?02: ldy #$003033 FE9437 84 70 sty <facbits3034 FE9439 C9 F9 ?04: cmp #$F93035 FE943B 30 C7 bmi XB8623036 FE943D A8 tay3037 FE943E A5 70 lda <facbits3038 FE9440 56 01 lsr <$01,x3039 FE9442 20 82 95 jsr XB9B03040 FE9445 XB8A3:3041 FE9445 24 6F bit <arisgn3042 FE9447 10 5D bpl XB8FE ; add mantissa3043 FE9449 A0 61 ldy #facexp3044 FE944B E0 69 cpx #argexp3045 FE944D F0 02 beq ?023046 FE944F A0 69 ldy #argexp3047 FE9451 38 ?02: sec3048 FE9452 49 FF eor #$FF3049 FE9454 65 59 adc <argbits3050 FE9456 85 70 sta <facbits3051 FE9458 8B phb ; save data bank3052 FE9459 A9 00 lda #0 ; set data bank = 03053 FE945B 48 pha ; need this because don't exist DP,y3054 FE945C AB plb3055 FE945D B9 04 7F lda !CBMP0+$04,y3056 FE9460 F5 04 sbc <$04,x3057 FE9462 85 65 sta <facm+33058 FE9464 B9 03 7F lda !CBMP0+$03,y3059 FE9467 F5 03 sbc <$03,x3060 FE9469 85 64 sta <facm+23061 FE946B B9 02 7F lda !CBMP0+$02,y3062 FE946E F5 02 sbc <$02,x3063 FE9470 85 63 sta <facm+13064 FE9472 B9 01 7F lda !CBMP0+$01,y3065 FE9475 F5 01 sbc <$01,x3066 FE9477 85 62 sta <facm3067 FE9479 AB plb ; restore data bank30683069 FE947A XB8D2:3070 FE947A B0 03 bcs XB8D73071 FE947C 20 EE 94 jsr XB9473072 FE947F3073 FE947F XB8D7:3074 FE947F A0 00 ldy #$003075 FE9481 98 tya3076 FE9482 18 clc3077 FE9483 A6 62 ?02: ldx <facm3078 FE9485 D0 49 bne XB9293079 FE9487 A6 63 ldx <facm+13080 FE9489 86 62 stx <facm3081 FE948B A6 64 ldx <facm+23082 FE948D 86 63 stx <facm+13083 FE948F A6 65 ldx <facm+33084 FE9491 86 64 stx <facm+23085 FE9493 A6 70 ldx <facbits3086 FE9495 86 65 stx <facm+3Tue Jul 17 11:00:22 2018 Page 103087 FE9497 84 70 sty <facbits3088 FE9499 69 08 adc #$083089 FE949B C9 20 cmp #$203090 FE949D D0 E4 bne ?0230913092 XB8F7: ; clear FAC #13093 FE949F A9 00 lda #$003094 XB8F9: ; set FAC #1 exponent3095 FE94A1 85 61 sta <facexp3096 XB8FB: ; set FAC #1 sign3097 FE94A3 85 66 sta <facsgn3098 FE94A5 60 rts30993100 XB8FE: ; add mantissa of FAC #2 to mantissa of FAC #13101 FE94A6 65 59 adc <argbits3102 FE94A8 85 70 sta <facbits3103 FE94AA A5 65 lda <facm+33104 FE94AC 65 6D adc <argm+33105 FE94AE 85 65 sta <facm+33106 FE94B0 A5 64 lda <facm+23107 FE94B2 65 6C adc <argm+23108 FE94B4 85 64 sta <facm+23109 FE94B6 A5 63 lda <facm+13110 FE94B8 65 6B adc <argm+13111 FE94BA 85 63 sta <facm+13112 FE94BC A5 62 lda <facm3113 FE94BE 65 6A adc <argm3114 FE94C0 85 62 sta <facm3115 FE94C2 80 19 bra XB936 ; normalize31163117 XB91D: ; shift left mantissa bit by bit3118 FE94C4 69 01 adc #$013119 FE94C6 06 70 asl <facbits3120 FE94C8 26 65 rol <facm+33121 FE94CA 26 64 rol <facm+23122 FE94CC 26 63 rol <facm+13123 FE94CE 26 62 rol <facm31243125 XB929: ; shift left mantissa if need -- A = exponent3126 FE94D0 10 F2 bpl XB91D3127 FE94D2 38 sec3128 FE94D3 E5 61 sbc <facexp3129 FE94D5 B0 C8 bcs XB8F7 ; underflow -> clear FAC #13130 FE94D7 49 FF eor #$FF ; 2's complement3131 FE94D9 69 01 adc #$013132 FE94DB 85 61 sta facexp31333134 XB936: ; normalize FAC #1 if carry3135 FE94DD 90 0E bcc XB946 ; no carry: exit31363137 XB938: ; increment FAC #1 exponent3138 FE94DF E6 61 inc <facexp3139 FE94E1 F0 42 beq XB97E ;OVERFLOW3140 FE94E3 66 62 ror <facm3141 FE94E5 66 63 ror <facm+13142 FE94E7 66 64 ror <facm+23143 FE94E9 66 65 ror <facm+3Tue Jul 17 11:00:22 2018 Page 113144 FE94EB 66 70 ror <facbits3145 FE94ED 60 XB946: rts31463147 XB947: ; FAC #1 2's complement3148 FE94EE A5 66 lda <facsgn3149 FE94F0 49 FF eor #$FF3150 FE94F2 85 66 sta <facsgn31513152 XB94D: ; mantissa 2's complement3153 FE94F4 A5 62 lda <facm3154 FE94F6 49 FF eor #$FF3155 FE94F8 85 62 sta <facm3156 FE94FA A5 63 lda <facm+13157 FE94FC 49 FF eor #$FF3158 FE94FE 85 63 sta <facm+13159 FE9500 A5 64 lda <facm+23160 FE9502 49 FF eor #$FF3161 FE9504 85 64 sta <facm+23162 FE9506 A5 65 lda <facm+33163 FE9508 49 FF eor #$FF3164 FE950A 85 65 sta <facm+33165 FE950C A5 70 lda <facbits3166 FE950E 49 FF eor #$FF3167 FE9510 85 70 sta <facbits3168 FE9512 E6 70 inc <facbits3169 FE9514 D0 0E bne XB97D ; exit31703171 XB96F: ; increment mantissa3172 FE9516 E6 65 inc <facm+33173 FE9518 D0 0A bne XB97D ; exit3174 FE951A E6 64 inc <facm+23175 FE951C D0 06 bne XB97D ; exit3176 FE951E E6 63 inc <facm+13177 FE9520 D0 02 bne XB97D ; exit3178 FE9522 E6 62 inc <facm3179 FE9524 60 XB97D: rts31803181 XB97E: ; overflow error3182 FE9525 24 96 bit <basic3183 FE9527 30 27 bmi ?103184 FE9529 A6 97 ldx <runf3185 FE952B F0 23 beq ?103186 FE952D A5 79 lda <s2fptr ; error from conversion routine ?3187 FE952F 10 1A bpl ?08 ; no - handle runtime error3188 FE9531 85 72 sta <s2fer ; setup conversion error flag3189 FE9533 A0 00 ldy #$003190 FE9535 INDEX163191 FE9535 C2 10 rep #PXFLAG3192 .LONGI on3193 .MNLIST3194 FE9537 A6 74 ldx <s2fstk ; recover stack pointer3195 FE9539 9A txs3196 FE953A A5 A2 lda <vbnk ; setup right data bank3197 FE953C 48 pha3198 FE953D AB plb3199 FE953E 88 dey3200 FE953F 8C E8 FF sty FAC1TMPTue Jul 17 11:00:22 2018 Page 123201 FE9542 8C EA FF sty FAC1TMP+23202 FE9545 INDEX083203 FE9545 E2 10 sep #PXFLAG3204 .LONGI off3205 .MNLIST3206 FE9547 8C EC FF sty FAC1TMP+43207 FE954A 60 rts3208 FE954B A2 09 ?08: ldx #$09 ; pascal: floating point overflow3209 FE954D 4C 50 1B jmp vpmerr3210 FE9550 A2 0F ?10: ldx #$0F ; overflow3211 FE9552 4C 47 01 jmp cbm_err32123213 XB983: ; shift right temp. mantissa by 8 bits3214 FE9555 A2 53 ldx #factmp-132153216 XB985: ; shift right 8 bits mantissa pointed by X3217 ; A = negative of shift count * 83218 FE9557 B4 04 ldy <$04,x3219 FE9559 84 70 sty <facbits3220 FE955B B4 03 ldy <$03,x3221 FE955D 94 04 sty <$04,x3222 FE955F B4 02 ldy <$02,x3223 FE9561 94 03 sty <$03,x3224 FE9563 B4 01 ldy <$01,x3225 FE9565 94 02 sty <$02,x3226 FE9567 A4 68 ldy facov3227 FE9569 94 01 sty <$01,x32283229 XB999: ; shift right mantissa3230 FE956B 69 08 adc #$083231 FE956D 30 E8 bmi XB985 ; shift 8 bits3232 FE956F F0 E6 beq XB985 ; shift 8 bits3233 FE9571 E9 08 sbc #$083234 FE9573 A8 tay ; Y = number of negative right shift3235 FE9574 A5 70 lda <facbits3236 FE9576 B0 14 bcs XB9BA ; no shift32373238 FE9578 XB9A6:3239 FE9578 16 01 asl <$01,x3240 FE957A 90 02 bcc ?023241 FE957C F6 01 inc <$01,x3242 FE957E 76 01 ?02: ror <$01,x3243 FE9580 76 01 ror <$01,x32443245 FE9582 XB9B0:3246 FE9582 76 02 ror <$02,x3247 FE9584 76 03 ror <$03,x3248 FE9586 76 04 ror <$04,x3249 FE9588 6A ror a3250 FE9589 C8 iny3251 FE958A D0 EC bne XB9A63252 FE958C XB9BA:3253 FE958C 18 clc3254 FE958D 60 rts32553256 ;---------------------------------------------------------------------------3257 ; LN() implementationTue Jul 17 11:00:22 2018 Page 133258 ;---------------------------------------------------------------------------32593260 ; compute LN(FAC #1)3261 FE958E XB9EA:3262 FE958E 20 ED 97 jsr XBC2B ; check if FAC #1 is zero or negative3263 FE9591 F0 02 beq ?02 ; error3264 FE9593 10 12 bpl XB9F43265 FE9595 24 96 ?02: bit <basic3266 FE9597 30 09 bmi ?043267 FE9599 A6 97 ldx <runf3268 FE959B F0 05 beq ?043269 FE959D A2 0C ldx #$0C ; pascal error: complex log or square root3270 FE959F 4C 50 1B jmp vpmerr3271 FE95A2 A2 0E ?04: LDX #$0E ; error: illegal quantity3272 FE95A4 4C 47 01 JMP cbm_err3273 FE95A7 XB9F4:3274 FE95A7 A5 61 lda <facexp ; here always CF = 03275 FE95A9 E9 7F sbc #$7F ; depolarize3276 FE95AB 48 pha3277 FE95AC A9 80 lda #$80 ; reduce FAC #13278 FE95AE 85 61 sta <facexp3279 FE95B0 A9 C6 lda #<XB9D6 ; 0.5 * SQRT(2)3280 FE95B2 A0 9C ldy #>XB9D63281 FE95B4 20 09 94 jsr XB867 ; add to FAC #13282 FE95B7 A9 CB lda #<XB9DB ; SQRT(2)3283 FE95B9 A0 9C ldy #>XB9DB3284 FE95BB A2 FE ldx #^XB9F4 ; current program bank3285 FE95BD 20 91 96 jsr XBB0F ; divide by SQRT(2)3286 FE95C0 A9 A2 lda #<XB9BC ; 1.03287 FE95C2 A0 9C ldy #>XB9BC3288 FE95C4 20 F3 93 jsr XB850 ; compute 1.0 - [FAC #1]3289 FE95C7 A9 B1 lda #<XB9C1 ; polynomial table3290 FE95C9 A0 9C ldy #>XB9C13291 FE95CB 20 A7 9B jsr XE043 ; compute polynomial3292 FE95CE A9 D0 lda #<XB9E03293 FE95D0 A0 9C ldy #>XB9E03294 FE95D2 20 09 94 jsr XB867 ; add -0.5 to FAC #13295 FE95D5 68 pla3296 FE95D6 20 92 99 jsr XBD7E ; add exponent to FAC #13297 FE95D9 A9 D5 lda #<XB9E5 ; LN(2)3298 FE95DB A0 9C ldy #>XB9E53299 FE95DD A2 FE ldx #^XB9F4 ; current program bank3300 ; mult LN(2)3301 FE95DF3302 ;---------------------------------------------------------------------------3303 ; multiplication implementation3304 ;---------------------------------------------------------------------------33053306 ; compute [float at address AYX] * [FAC #1]3307 FE95DF XBA28:3308 FE95DF 20 C0 93 jsr XBA8CL ; move float from long address to FAC #233093310 ; compute [FAC #2] * [FAC #1]3311 XBA2B: ; A = FAC #1 exponent3312 FE95E2 F0 59 beq XBA8B ; exit if FAC #1 = 03313 FE95E4 20 3E 96 jsr XBAB7 ; add the 2 exponents3314 FE95E7 64 54 stz <factmp ; clear temp. mantissaTue Jul 17 11:00:22 2018 Page 143315 FE95E9 64 55 stz <factmp+13316 FE95EB 64 56 stz <factmp+23317 FE95ED 64 57 stz <factmp+33318 FE95EF A5 70 lda <facbits ; mult. byte by byte3319 FE95F1 20 0B 96 jsr XBA593320 FE95F4 A5 65 lda <facm+33321 FE95F6 20 0B 96 jsr XBA593322 FE95F9 A5 64 lda <facm+23323 FE95FB 20 0B 96 jsr XBA593324 FE95FE A5 63 lda <facm+13325 FE9600 20 0B 96 jsr XBA593326 FE9603 A5 62 lda <facm3327 FE9605 20 10 96 jsr XBA5E3328 FE9608 4C 20 97 jmp XBB8F ; copy temp. mantissa to FAC #133293330 XBA59: ; mult byte A with FAC #2 mantissa and put result in temp. mantissa LA07B3331 FE960B D0 03 bne XBA5E3332 FE960D 4C 55 95 jmp XB983 ; shift right temp. mantissa3333 FE9610 XBA5E:3334 FE9610 4A lsr a3335 FE9611 09 80 ora #$803336 FE9613 A8 ?02: tay3337 FE9614 90 19 bcc ?043338 FE9616 18 clc3339 FE9617 A5 57 lda <factmp+33340 FE9619 65 6D adc <argm+33341 FE961B 85 57 sta <factmp+33342 FE961D A5 56 lda <factmp+23343 FE961F 65 6C adc <argm+23344 FE9621 85 56 sta <factmp+23345 FE9623 A5 55 lda <factmp+13346 FE9625 65 6B adc <argm+13347 FE9627 85 55 sta <factmp+13348 FE9629 A5 54 lda <factmp3349 FE962B 65 6A adc <argm3350 FE962D 85 54 sta <factmp3351 FE962F 66 54 ?04: ror <factmp3352 FE9631 66 55 ror <factmp+13353 FE9633 66 56 ror <factmp+23354 FE9635 66 57 ror <factmp+33355 FE9637 66 70 ror <facbits3356 FE9639 98 tya3357 FE963A 4A lsr a3358 FE963B D0 D6 bne ?023359 FE963D 60 XBA8B: rts33603361 XBAB7: ; add exponents3362 FE963E A5 69 lda <argexp33633364 FE9640 XBAB9:3365 FE9640 F0 20 beq XBADA ; clear FAC #1 and exit3366 FE9642 18 clc3367 FE9643 65 61 adc <facexp3368 FE9645 90 05 bcc ?023369 FE9647 30 1E bmi XBADF3370 FE9649 18 clc3371 FE964A 80 02 bra ?04Tue Jul 17 11:00:22 2018 Page 153372 FE964C 10 14 ?02: bpl XBADA3373 FE964E 69 80 ?04: adc #$803374 FE9650 85 61 sta <facexp3375 FE9652 D0 03 bne ?063376 FE9654 4C A3 94 jmp XB8FB3377 FE9657 A5 6F ?06: lda <arisgn3378 FE9659 85 66 sta <facsgn3379 FE965B 60 rts33803381 FE965C XBAD4:3382 FE965C A5 66 lda <facsgn3383 FE965E 49 FF eor #$FF3384 FE9660 30 05 bmi XBADF3385 XBADA: ; skip return address and clear FAC #13386 FE9662 68 pla3387 FE9663 68 pla3388 FE9664 4C 9F 94 jmp XB8F73389 FE96673390 XBADF: ; overflow3391 FE9667 4C 25 95 jmp XB97E33923393 XBAE2: ; compute (FAC #1) * 103394 FE966A 20 CE 97 jsr XBC0C ; round FAC #1 and move it to FAC #23395 FE966D AA tax3396 FE966E F0 10 beq XBAF83397 FE9670 18 clc3398 FE9671 69 02 adc #$02 ; mult * 43399 FE9673 B0 F2 bcs XBADF ; overflow3400 FE9675 A2 00 XBAED: ldx #$003401 FE9677 86 6F stx <arisgn3402 FE9679 20 19 94 jsr XB877 ; add itself3403 FE967C E6 61 inc <facexp ; mult * 23404 FE967E F0 E7 beq XBADF ; overflow3405 FE9680 60 XBAF8: rts34063407 ;---------------------------------------------------------------------------3408 ; division implementation3409 ;---------------------------------------------------------------------------34103411 ; compute [FAC #1] / 103412 FE9681 XBAFE:3413 FE9681 20 CE 97 jsr XBC0C ; round FAC #1 and move it to FAC #23414 FE9684 A9 A7 lda #<FCON103415 FE9686 A0 9C ldy #>FCON103416 FE9688 A2 00 ldx #0 ; comparison sign34173418 ; compute [FAC #2] / [float at AY address (in current program bank)]3419 FE968A XBB07:3420 FE968A 86 6F stx <arisgn ; X = comparison sign3421 FE968C 20 49 97 jsr XBBA2K ; move from addr. AY (program bank) to FAC #13422 FE968F 80 03 bra XBB123423 FE96913424 ; compute [float at AYX address] / [FAC #1]3425 FE9691 XBB0F:3426 FE9691 20 C0 93 jsr XBA8CL ; move from address AYX to FAC #23427 FE96943428 ; compute [FAC #2] / [FAC #1]Tue Jul 17 11:00:22 2018 Page 163429 XBB12: ; A = FAC #1 exponent3430 FE9694 F0 78 beq ?16 ; division by zero error3431 FE9696 20 DD 97 jsr XBC1B ; round FAC #13432 FE9699 A9 00 lda #$003433 FE969B 38 sec3434 FE969C E5 61 sbc <facexp ; 2's complement3435 FE969E 85 61 sta <facexp3436 FE96A0 20 3E 96 jsr XBAB7 ; add exponents3437 FE96A3 E6 61 inc <facexp3438 FE96A5 F0 C0 beq XBADF ; overflow3439 ;ldx #$04 ; here different !3440 FE96A7 A2 FF ldx #$FF ; here different !3441 FE96A9 A9 01 lda #$013442 FE96AB A4 6A ?02: ldy <argm3443 FE96AD C4 62 cpy <facm3444 FE96AF D0 10 bne ?043445 FE96B1 A4 6B ldy <argm+13446 FE96B3 C4 63 cpy <facm+13447 FE96B5 D0 0A bne ?043448 FE96B7 A4 6C ldy <argm+23449 FE96B9 C4 64 cpy <facm+23450 FE96BB D0 04 bne ?043451 FE96BD A4 6D ldy <argm+33452 FE96BF C4 65 cpy <facm+33453 FE96C1 08 ?04: php3454 FE96C2 2A rol a3455 FE96C3 90 0D bcc ?063456 FE96C5 E8 inx3457 FE96C6 E0 04 cpx #$043458 FE96C8 B0 39 bcs ?143459 ;dex ; here different !3460 ;bmi ?14 ; no store !3461 FE96CA 95 54 sta <factmp,x3462 FE96CC E0 03 cpx #$033463 FE96CE F0 2F beq ?123464 FE96D0 A9 01 lda #$013465 FE96D2 28 ?06: plp3466 FE96D3 B0 0E bcs ?103467 FE96D5 06 6D ?08: asl <argm+33468 FE96D7 26 6C rol <argm+23469 FE96D9 26 6B rol <argm+13470 FE96DB 26 6A rol <argm3471 FE96DD B0 E2 bcs ?043472 FE96DF 30 CA bmi ?023473 FE96E1 10 DE bpl ?043474 FE96E3 A8 ?10: tay3475 FE96E4 A5 6D lda <argm+33476 FE96E6 E5 65 sbc <facm+33477 FE96E8 85 6D sta <argm+33478 FE96EA A5 6C lda <argm+23479 FE96EC E5 64 sbc <facm+23480 FE96EE 85 6C sta <argm+23481 FE96F0 A5 6B lda <argm+13482 FE96F2 E5 63 sbc <facm+13483 FE96F4 85 6B sta <argm+13484 FE96F6 A5 6A lda <argm3485 FE96F8 E5 62 sbc <facmTue Jul 17 11:00:22 2018 Page 173486 FE96FA 85 6A sta <argm3487 FE96FC 98 tya3488 FE96FD 80 D6 bra ?083489 FE96FF A9 40 ?12: lda #$403490 FE9701 D0 CF bne ?063491 FE9703 0A ?14: asl a3492 FE9704 0A asl a3493 FE9705 0A asl a3494 FE9706 0A asl a3495 FE9707 0A asl a3496 FE9708 0A asl a3497 FE9709 85 70 sta <facbits3498 FE970B 28 plp3499 FE970C 80 12 bra XBB8F3500 FE970E 24 96 ?16: bit <basic3501 FE9710 30 09 bmi ?183502 FE9712 A6 97 ldx <runf3503 FE9714 F0 05 beq ?183504 FE9716 A2 12 ldx #$12 ; pascal: division by zero3505 FE9718 4C 50 1B jmp vpmerr3506 FE971B A2 14 ?18: ldx #$14 ; division by zero3507 FE971D 4C 47 01 jmp cbm_err35083509 XBB8F: ; copy mantissa3510 FE9720 A5 54 lda <factmp3511 FE9722 85 62 sta <facm3512 FE9724 A5 55 lda <factmp+13513 FE9726 85 63 sta <facm+13514 FE9728 A5 56 lda <factmp+23515 FE972A 85 64 sta <facm+23516 FE972C A5 57 lda <factmp+33517 FE972E 85 65 sta <facm+33518 FE9730 4C 7F 94 jmp XB8D7 ; normalize FAC #135193520 ;---------------------------------------------------------------------------3521 ; move FAC #1 routines3522 ;---------------------------------------------------------------------------35233524 ; move float number at address AY in current data bank to FAC #13525 FE9733 XBBA2:3526 FE9733 85 76 sta <longp3527 FE9735 84 77 sty <longp+13528 FE9737 8B phb3529 FE9738 FA plx3530 FE9739 80 14 bra XBBA2X3531 FE973B3532 ; move float number at address AYX to FAC #13533 FE973B XBBA2L:3534 FE973B 85 76 sta <longp3535 FE973D 84 77 sty <longp+13536 FE973F 80 0E bra XBBA2X35373538 ; move float number at address AY in data bank that hold FAC1TMP to FAC #13539 FE9741 XBBA2T:3540 FE9741 85 76 sta <longp3541 FE9743 84 77 sty <longp+13542 FE9745 A6 A2 ldx <vbnkTue Jul 17 11:00:22 2018 Page 183543 FE9747 80 06 bra XBBA2X35443545 ; move float number at address AY in current program bank to FAC #13546 FE9749 XBBA2K:3547 FE9749 85 76 sta <longp3548 FE974B 84 77 sty <longp+13549 FE974D A2 FE ldx #^XBBA2K3550 FE974F XBBA2X:3551 FE974F 86 78 stx <longp+23552 FE9751 A0 04 ldy #43553 FE9753 B7 76 lda [longp],y3554 FE9755 85 65 sta <facm+33555 FE9757 88 dey3556 FE9758 B7 76 lda [longp],y3557 FE975A 85 64 sta <facm+23558 FE975C 88 dey3559 FE975D B7 76 lda [longp],y3560 FE975F 85 63 sta <facm+13561 FE9761 88 dey3562 FE9762 B7 76 lda [longp],y3563 FE9764 85 66 sta <facsgn3564 FE9766 09 80 ora #$803565 FE9768 85 62 sta <facm3566 FE976A 88 dey3567 FE976B B7 76 lda [longp],y3568 FE976D 85 61 sta <facexp3569 FE976F 84 70 sty <facbits3570 FE9771 60 rts35713572 ; round FAC #1 and move to FAC1TMP3573 FE9772 XBBC7:3574 FE9772 A2 E8 ldx #<FAC1TMP3575 FE9774 A0 FF ldy #>FAC1TMP3576 FE9776 A5 A2 lda <vbnk ; bank that hold FAC1TMP/FAC2TMP3577 FE9778 80 1A bra XBBD4X35783579 ; round FAC #1 and move to FAC2TMP3580 FE977A XBBCA:3581 FE977A A2 ED ldx #<FAC2TMP3582 FE977C A0 FF ldy #>FAC2TMP3583 FE977E A5 A2 lda <vbnk ; bank that hold FAC1TMP/FAC2TMP3584 FE9780 80 12 bra XBBD4X35853586 ; round FAC #1 and move to FAC3TMP3587 FE9782 XBBX3:3588 FE9782 A2 F2 ldx #<FAC3TMP3589 FE9784 A0 FF ldy #>FAC3TMP3590 FE9786 A5 A2 lda <vbnk ; bank that hold FAC1TMP/FAC2TMP3591 FE9788 80 0A bra XBBD4X35923593 ; round FAC #1 and move to forpnt and bank that hold basic variables3594 FE978A XBBD0:3595 FE978A A6 49 ldx <forpnt3596 FE978C A4 4A ldy <forpnt+135973598 ; round FAC #1 and move to XY addr, in data bank that hold basic variables3599 FE978E XBBD4V:Tue Jul 17 11:00:22 2018 Page 193600 FE978E A5 A2 lda <vbnk3601 FE9790 80 02 bra XBBD4X3602 FE97923603 ; round FAC #1 and move to XY address in current data bank3604 FE9792 XBBD4:3605 FE9792 8B phb ; A = current data bank3606 FE9793 68 pla36073608 FE9794 XBBD4L:3609 FE9794 XBBD4X:3610 FE9794 86 76 stx <longp3611 FE9796 84 77 sty <longp+13612 FE9798 85 78 sta <longp+23613 FE979A 20 DD 97 jsr XBC1B ; round FAC #13614 FE979D A0 04 ldy #$043615 FE979F A5 65 lda <facm+33616 FE97A1 97 76 sta [longp],y3617 FE97A3 88 dey3618 FE97A4 A5 64 lda <facm+23619 FE97A6 97 76 sta [longp],y3620 FE97A8 88 dey3621 FE97A9 A5 63 lda <facm+13622 FE97AB 97 76 sta [longp],y3623 FE97AD 88 dey3624 FE97AE A5 66 lda <facsgn3625 FE97B0 09 7F ora #$7F3626 FE97B2 25 62 and <facm3627 FE97B4 97 76 sta [longp],y3628 FE97B6 88 dey3629 FE97B7 A5 61 lda <facexp3630 FE97B9 97 76 sta [longp],y3631 FE97BB 84 70 sty <facbits3632 FE97BD 60 rts36333634 XBBFC: ; move FAC #2 to FAC #13635 FE97BE A5 6E lda <argsgn3636 XBBFE: ; set FAC #1 sign and move abs FAC #2 to FAC #13637 FE97C0 85 66 sta <facsgn3638 FE97C2 A2 05 ldx #$053639 FE97C4 B5 68 ?02: lda <argexp-1,x3640 FE97C6 95 60 sta <facexp-1,x3641 FE97C8 CA dex3642 FE97C9 D0 F9 bne ?023643 FE97CB 86 70 stx <facbits3644 FE97CD 60 rts36453646 XBC0C: ; round FAC #1 and move it to FAC #23647 FE97CE 20 DD 97 jsr XBC1B36483649 XBC0F: ; move FAC #1 to FAC #2 without before rounding FAC #13650 FE97D1 A2 06 ldx #$063651 FE97D3 B5 60 ?02: lda <facexp-1,x3652 FE97D5 95 68 sta <argexp-1,x3653 FE97D7 CA dex3654 FE97D8 D0 F9 bne ?023655 FE97DA 86 70 stx <facbits3656 FE97DC 60 XBC1A: rtsTue Jul 17 11:00:22 2018 Page 2036573658 XBC1B: ; round FAC #13659 FE97DD A5 61 lda <facexp3660 FE97DF F0 FB beq XBC1A3661 FE97E1 06 70 asl <facbits3662 FE97E3 90 F7 bcc XBC1A36633664 XBC23: ; increment mantissa and normalize3665 FE97E5 20 16 95 jsr XB96F3666 FE97E8 D0 F2 bne XBC1A3667 FE97EA 4C DF 94 jmp XB93836683669 ;---------------------------------------------------------------------------3670 ; sgn/abs function's3671 ;---------------------------------------------------------------------------36723673 XBC2B: ; get sign of FAC #13674 FE97ED A5 61 lda <facexp3675 FE97EF F0 09 beq XBC383676 FE97F1 A5 66 XBC2F: lda <facsgn3677 FE97F3 2A XBC31: rol a3678 FE97F4 A9 FF lda #$FF3679 FE97F6 B0 02 bcs XBC383680 FE97F8 A9 01 lda #$013681 FE97FA 60 XBC38: rts36823683 XBC39: ; sign function (return a float)3684 FE97FB 20 ED 97 jsr XBC2B36853686 XBC3C: ; convert a signed byte in A to float format3687 FE97FE 85 62 sta <facm3688 FE9800 64 63 stz <facm+13689 FE9802 A2 88 ldx #$88 ; exponent for a byte3690 FE9804 A5 62 XBC44: lda <facm3691 FE9806 49 FF eor #$FF3692 FE9808 2A rol a3693 FE9809 A9 00 XBC49: lda #$00 ; CF = 1 positive, CF = 0 negative3694 FE980B 85 65 sta <facm+33695 FE980D 85 64 sta <facm+23696 FE980F 86 61 XBC4F: stx <facexp3697 FE9811 85 70 sta <facbits3698 FE9813 85 66 sta <facsgn3699 FE9815 4C 7A 94 jmp XB8D2 ; normalize FAC #137003701 XBC58: ; absolute value3702 FE9818 64 66 stz <facsgn3703 FE981A 60 rts37043705 ;---------------------------------------------------------------------------3706 ; comparison function's3707 ;---------------------------------------------------------------------------37083709 ; compare float at address AY in current data bank with FAC #13710 FE981B XBC5B:3711 FE981B 85 76 sta <longp3712 FE981D XBC5D:3713 FE981D 84 77 sty <longp+1Tue Jul 17 11:00:22 2018 Page 213714 FE981F 8B phb3715 FE9820 FA plx3716 FE9821 86 78 XBC5DX: stx <longp+23717 FE9823 A0 00 ldy #$003718 FE9825 B7 76 lda [longp],y3719 FE9827 C8 iny3720 FE9828 AA tax3721 FE9829 F0 C2 beq XBC2B3722 FE982B B7 76 lda [longp],y3723 FE982D 45 66 eor <facsgn3724 FE982F 30 C0 bmi XBC2F3725 FE9831 E4 61 cpx <facexp3726 FE9833 D0 21 bne XBC923727 FE9835 B7 76 lda [longp],y3728 FE9837 09 80 ora #$803729 FE9839 C5 62 cmp <facm3730 FE983B D0 19 bne XBC923731 FE983D C8 iny3732 FE983E B7 76 lda [longp],y3733 FE9840 C5 63 cmp <facm+13734 FE9842 D0 12 bne XBC923735 FE9844 C8 iny3736 FE9845 B7 76 lda [longp],y3737 FE9847 C5 64 cmp <facm+23738 FE9849 D0 0B bne XBC923739 FE984B C8 iny3740 FE984C A9 7F lda #$7F3741 FE984E C5 70 cmp <facbits3742 FE9850 B7 76 lda [longp],y3743 FE9852 E5 65 sbc <facm+33744 FE9854 F0 37 beq XBCBA3745 FE9856 A5 66 XBC92: lda <facsgn3746 FE9858 90 99 bcc XBC313747 FE985A 49 FF eor #$FF3748 FE985C 80 95 bra XBC3137493750 ; compare float at address AYK with FAC #13751 FE985E XBC5BK:3752 FE985E 85 76 sta <longp3753 FE9860 84 77 sty <longp+13754 FE9862 A2 FE ldx #^XBC5BK3755 FE9864 80 BB bra XBC5DX37563757 ; compare float at FAC1/2/3TMP with FAC #13758 FE9866 XBC5BT:3759 FE9866 85 76 sta <longp3760 FE9868 84 77 sty <longp+13761 FE986A A6 A2 ldx <vbnk3762 FE986C 80 B3 bra XBC5DX3763 FE986E3764 ;---------------------------------------------------------------------------3765 ; extraction of the integer part of FAC #13766 ;---------------------------------------------------------------------------37673768 ; convert float FAC #1 in 4-bytes signed integer3769 FE986E XBC9B:3770 FE986E A5 61 lda <facexpTue Jul 17 11:00:22 2018 Page 223771 FE9870 F0 4A beq XBCE93772 FE9872 38 sec3773 FE9873 E9 A0 sbc #$A03774 FE9875 24 66 bit <facsgn3775 FE9877 10 09 bpl ?023776 FE9879 AA tax3777 FE987A A9 FF lda #$FF3778 FE987C 85 68 sta <facov3779 FE987E 20 F4 94 jsr XB94D3780 FE9881 8A txa3781 FE9882 A2 61 ?02: ldx #facexp3782 FE9884 C9 F9 cmp #$F93783 FE9886 10 06 bpl XBCBB3784 FE9888 20 6B 95 jsr XB9993785 FE988B 84 68 sty <facov3786 ;stz <facov3787 FE988D 60 XBCBA: rts3788 FE988E A8 XBCBB: tay3789 FE988F A5 66 lda <facsgn3790 FE9891 29 80 and #$803791 FE9893 46 62 lsr <facm3792 FE9895 05 62 ora <facm3793 FE9897 85 62 sta <facm3794 FE9899 20 82 95 jsr XB9B03795 FE989C 84 68 sty <facov3796 ;stz <facov3797 FE989E 60 rts37983799 XBCCC: ; basic int() function3800 FE989F A5 61 lda <facexp3801 FE98A1 C9 A0 cmp #$A03802 FE98A3 B0 20 bcs XBCF23803 FE98A5 20 6E 98 jsr XBC9B3804 FE98A8 84 70 sty <facbits3805 FE98AA A5 66 lda <facsgn3806 FE98AC 84 66 sty <facsgn3807 FE98AE 49 80 eor #$803808 FE98B0 2A rol a3809 FE98B1 A9 A0 lda #$A03810 FE98B3 85 61 sta <facexp3811 FE98B5 A5 65 lda <facm+33812 FE98B7 85 5B sta <faclsb3813 FE98B9 4C 7A 94 jmp XB8D2 ; normalize38143815 XBCE9: ; clear FAC #1 mantissa3816 FE98BC 85 62 sta <facm3817 FE98BE 85 63 sta <facm+13818 FE98C0 85 64 sta <facm+23819 FE98C2 85 65 sta <facm+33820 FE98C4 A8 tay3821 FE98C5 60 XBCF2: rts38223823 ;---------------------------------------------------------------------------3824 ; conversion decimal string to float3825 ;---------------------------------------------------------------------------38263827 ; conversion routine in pascal mode - return float in FAC1TMP and FAC #1Tue Jul 17 11:00:22 2018 Page 233828 FE98C6 XBCF3P:3829 FE98C6 A2 80 ldx #$80 ; trap overflow run time error3830 FE98C8 86 79 stx <s2fptr ; this avoid run time error in conversion3831 FE98CA 20 D8 98 jsr XBCF3 ; call conversion routine3832 ; here return after an overflow error too3833 FE98CD 64 9B stz <skpbk ; from now skip blank while scan string3834 FE98CF 64 79 stz <s2fptr ; re-enable overflow error3835 FE98D1 A5 72 lda <s2fer ; if not zero error pending exit with A <> 03836 FE98D3 D0 F0 bne XBCF23837 FE98D5 4C 04 03 jmp cbm_chrgot ; return last char scanned (must be zero)3838 FE98D83839 ; conversion routine different from the one of cbm basic3840 FE98D8 XBCF3:3841 FE98D8 INDEX163842 FE98D8 C2 10 rep #PXFLAG3843 .LONGI on3844 .MNLIST3845 FE98DA BA tsx ; save current stack pointer3846 FE98DB 86 74 stx <s2fstk ; for error recover if pascal mode3847 FE98DD INDEX083848 FE98DD E2 10 sep #PXFLAG3849 .LONGI off3850 .MNLIST3851 FE98DF A2 1A ldx #(fbufpt2-argbits) ; clear working area3852 FE98E1 74 59 ?02: stz <argbits,x3853 FE98E3 CA dex3854 FE98E4 10 FB bpl ?023855 FE98E6 64 9B stz <skpbk ; skip blank at string beginning3856 FE98E8 20 04 03 jsr cbm_chrgot ; get first char. not blank3857 FE98EB 90 11 bcc ?14 ; first char is digit: after don't skip blank3858 FE98ED C9 2B cmp #'+' ; at beginning only '+' or '-' allowed3859 FE98EF F0 08 beq ?10 ; next char (still skip blanks)3860 FE98F1 C9 2D cmp #'-'3861 FE98F3 D0 64 bne ?36 ; end conversion: exit with float = 03862 FE98F5 A2 80 ldx #$803863 FE98F7 86 67 stx <sgnfg ; set flag minus sign3864 FE98F9 20 FE 02 ?10: jsr cbm_chrget ; next char after sign: expected a digit3865 FE98FC B0 57 bcs ?35a ; error: no digit or end of string3866 FE98FE A2 80 ?14: ldx #$803867 FE9900 86 9B stx <skpbk ; from now don't skip blank3868 FE9902 20 84 99 ?16: jsr XBD6A ; add digit to partial FAC #13869 FE9905 20 FE 02 ?18: jsr cbm_chrget ; next char3870 FE9908 90 F8 bcc ?16 ; is a digit3871 FE990A C9 2E cmp #'.' ; is decimal dot ?3872 FE990C F0 3C beq ?35 ; set dot flag3873 FE990E C9 45 cmp #'E' ; exponential form ?3874 FE9910 F0 04 beq ?223875 FE9912 C9 65 cmp #'e'3876 FE9914 D0 43 bne ?36 ; end conversion: setup float (no exponent)3877 FE9916 20 FE 02 ?22: jsr cbm_chrget ; next char: expected a digit or a sign3878 FE9919 90 1F bcc ?32 ; if digit add to partial exponent value3879 FE991B C9 2B cmp #'+'3880 FE991D F0 16 beq ?30 ; get next exponent digit3881 FE991F 24 96 bit <basic3882 FE9921 10 04 bpl ?24 ; token only in basic mode3883 FE9923 C9 AA cmp #$AA ; token '+' in basic mode3884 FE9925 F0 0E beq ?30 ; get next exponent digitTue Jul 17 11:00:22 2018 Page 243885 FE9927 C9 2D ?24: cmp #'-'3886 FE9929 F0 08 beq ?28 ; set exponent sign flag3887 FE992B 24 96 bit <basic3888 FE992D 10 10 bpl ?34 ; end: setup float with exponent3889 FE992F C9 AB cmp #$AB ; token '-' in basic mode3890 FE9931 D0 0C bne ?34 ; end: setup float with exponent3891 FE9933 66 60 ?28: ror <expsgn ; here CF = 1 so bit 7 = 13892 FE9935 20 FE 02 ?30: jsr cbm_chrget ; next exponent digit3893 FE9938 B0 05 bcs ?34 ; no digit -- end of exponent scan3894 FE993A 20 A5 99 ?32: jsr XBD91 ; evaluate partial exponent3895 FE993D 80 F6 bra ?30 ; next digit3896 FE993F 24 60 ?34: bit <expsgn ; evaluate full exponent value3897 FE9941 10 16 bpl ?36 ; positive exponent3898 FE9943 A9 00 lda #03899 FE9945 38 sec3900 FE9946 E5 5E sbc <expval ; negative exponent3901 FE9948 80 11 bra ?383902 FE994A 66 5F ?35: ror <dotfg ; here CF = 1 - set dot flag3903 FE994C 24 5F bit <dotfg ; check if already dot flag3904 FE994E 70 05 bvs ?35a ; more than one only decimal dot3905 FE9950 20 FE 02 jsr cbm_chrget ; after decimal dot expected one digit3906 FE9953 90 AD bcc ?16 ; ok, add digit to float and continue scan3907 FE9955 A2 80 ?35a: ldx #$80 ; 2 dot or no digit after dot or after sign3908 FE9957 86 72 stx <s2fer ; set conversion error3909 FE9959 A5 5E ?36: lda <expval ; exponent value3910 FE995B 38 ?38: sec ; compute true exponent3911 FE995C E5 5D sbc <dccnt ; subtract decimal place number3912 FE995E 85 5E sta <expval3913 FE9960 F0 12 beq ?44 ; no adjustment3914 FE9962 10 09 bpl ?42 ; mult * 10 loop3915 FE9964 20 81 96 ?40: jsr XBAFE ; div / 10 loop3916 FE9967 E6 5E inc <expval3917 FE9969 D0 F9 bne ?403918 FE996B F0 07 beq ?443919 FE996D 20 6A 96 ?42: jsr XBAE2 ; mult * 10 loop3920 FE9970 C6 5E dec <expval3921 FE9972 D0 F9 bne ?423922 FE9974 A5 67 ?44: lda <sgnfg3923 FE9976 10 03 bpl ?46 ; positive float3924 FE9978 20 4A 9B jsr XBFB4 ; negate float3925 FE997B 64 9B ?46: stz <skpbk ; from now skip blank while scan string3926 FE997D 24 96 bit <basic3927 FE997F 30 42 bmi XBDB0 ; finish here in basic mode3928 FE9981 4C 72 97 jmp XBBC7 ; round FAC #1 and move to FAC1TMP39293930 XBD6A: ; accumulate current digit in A into FAC #13931 FE9984 48 pha3932 FE9985 24 5F bit <dotfg3933 FE9987 10 02 bpl ?023934 FE9989 E6 5D inc <dccnt ; increment decimal digit count3935 FE998B 20 6A 96 ?02: jsr XBAE2 ; [FAC #1] * 103936 FE998E 68 pla3937 FE998F 38 sec3938 FE9990 E9 30 sbc #$3039393940 XBD7E: ; add a byte value to FAC #13941 FE9992 48 phaTue Jul 17 11:00:22 2018 Page 253942 FE9993 20 CE 97 jsr XBC0C ; round FAC #1 and move to FAC #23943 FE9996 68 pla3944 FE9997 20 FE 97 jsr XBC3C ; convert byte A in float in FAC #13945 FE999A A5 6E lda <argsgn3946 FE999C 45 66 eor <facsgn3947 FE999E 85 6F sta <arisgn3948 FE99A0 A6 61 ldx <facexp3949 FE99A2 4C 0C 94 jmp XB86A ; add FAC #1 and FAC #239503951 XBD91: ; compute partial exponent3952 FE99A5 E9 2F sbc #$2F ; here come with CF = 03953 FE99A7 85 5C sta <tmpb1 ; save digit3954 FE99A9 A5 5E lda <expval ; current exponent value3955 FE99AB C9 0A cmp #$0A3956 FE99AD 90 09 bcc ?02 ; exponent < 10 -> can mult *103957 FE99AF A9 64 lda #$64 ; max. exp. = +/-1003958 FE99B1 24 60 bit <expsgn3959 FE99B3 30 0C bmi ?04 ; min. negative exponent = -1003960 FE99B5 4C 25 95 jmp XB97E ; overflow if exponent positive3961 FE99B8 0A ?02: asl a ; mult. partial exponent *103962 FE99B9 0A asl a3963 FE99BA 18 clc3964 FE99BB 65 5E adc <expval3965 FE99BD 0A asl a3966 FE99BE 18 clc3967 FE99BF 65 5C adc <tmpb1 ; add current digit3968 FE99C1 85 5E ?04: sta <expval ; store exponent3969 FE99C3 60 XBDB0: rts39703971 ;---------------------------------------------------------------------------3972 ; conversion float to decimal string3973 ;---------------------------------------------------------------------------39743975 ; final string is TTNUMSTR in variables/p-code bank3976 FE99C4 XBDDD:3977 FE99C4 A0 01 ldy #$01 ; index to final string3978 FE99C6 XBDDF:3979 FE99C6 8B phb ; save current data bank3980 FE99C7 A5 A2 lda <vbnk3981 FE99C9 48 pha3982 FE99CA AB plb3983 FE99CB A9 20 lda #' '3984 FE99CD 24 66 bit <facsgn3985 FE99CF 10 02 bpl ?023986 FE99D1 A9 2D lda #'-'3987 FE99D3 99 59 FF ?02: sta !TTNUMSTR-1,y3988 FE99D6 64 66 stz <facsgn3989 FE99D8 84 71 sty <decidx3990 FE99DA C8 iny3991 FE99DB A9 30 lda #$303992 FE99DD A6 61 ldx <facexp3993 FE99DF D0 03 bne ?043994 FE99E1 4C FC 9A jmp ?583995 FE99E4 A9 00 ?04: lda #$003996 FE99E6 E0 80 cpx #$803997 FE99E8 F0 02 beq ?063998 FE99EA B0 0B bcs ?08Tue Jul 17 11:00:22 2018 Page 263999 FE99EC A9 EE ?06: lda #<XBDBD4000 FE99EE A0 9C ldy #>XBDBD4001 FE99F0 A2 FE ldx #^XBDDD4002 FE99F2 20 DF 95 jsr XBA28 ; mult4003 FE99F5 A9 F7 lda #$F74004 FE99F7 85 5D ?08: sta <dccnt4005 FE99F9 A9 E9 ?10: lda #<XBDB8 ; manage float so 1e8 <= x <= 1e94006 FE99FB A0 9C ldy #>XBDB84007 FE99FD 20 5E 98 jsr XBC5BK4008 FE9A00 F0 1E beq ?204009 FE9A02 10 12 bpl ?164010 FE9A04 A9 E4 ?12: lda #<XBDB34011 FE9A06 A0 9C ldy #>XBDB34012 FE9A08 20 5E 98 jsr XBC5BK4013 FE9A0B F0 02 beq ?144014 FE9A0D 10 0E bpl ?184015 FE9A0F 20 6A 96 ?14: jsr XBAE2 ; mult * 104016 FE9A12 C6 5D dec <dccnt4017 FE9A14 D0 EE bne ?124018 FE9A16 20 81 96 ?16: jsr XBAFE ; div by 104019 FE9A19 E6 5D inc <dccnt4020 FE9A1B D0 DC bne ?104021 FE9A1D 20 ED 93 ?18: jsr XB849 ; round to +0.54022 FE9A20 20 6E 98 ?20: jsr XBC9B ; convert to integer 32 bits4023 FE9A23 A2 01 ldx #$014024 FE9A25 A5 5D lda <dccnt4025 FE9A27 18 clc4026 FE9A28 69 0A adc #$0A4027 FE9A2A 30 09 bmi ?224028 FE9A2C C9 0B cmp #$0B4029 FE9A2E B0 06 bcs ?244030 FE9A30 69 FF adc #$FF4031 FE9A32 AA tax4032 FE9A33 A9 02 lda #$024033 FE9A35 38 ?22: sec4034 FE9A36 E9 02 ?24: sbc #$024035 FE9A38 85 5E sta <expval4036 FE9A3A 86 5D stx <dccnt4037 FE9A3C 8A txa4038 FE9A3D F0 02 beq ?26 ; decimal form4039 FE9A3F 10 19 bpl ?30 ; exponential form4040 FE9A41 A4 71 ?26: ldy <decidx4041 FE9A43 A9 30 lda #'0' ; store '0.' or '0.0'4042 FE9A45 C8 iny4043 FE9A46 99 59 FF sta !TTNUMSTR-1,y4044 FE9A49 A9 2E lda #$2E4045 FE9A4B C8 iny4046 FE9A4C 99 59 FF sta !TTNUMSTR-1,y4047 FE9A4F 8A txa4048 FE9A50 F0 06 beq ?284049 FE9A52 A9 30 lda #$304050 FE9A54 C8 iny4051 FE9A55 99 59 FF sta !TTNUMSTR-1,y4052 FE9A58 84 71 ?28: sty <decidx4053 FE9A5A A0 00 ?30: ldy #$004054 FE9A5C A2 80 ldx #$804055 FE9A5E 8B ?32: phbTue Jul 17 11:00:22 2018 Page 274056 FE9A5F 4B phk4057 FE9A60 AB plb4058 FE9A61 A5 65 ?34: lda <facm+34059 FE9A63 18 clc4060 FE9A64 79 F6 9C adc !XBF16+3,y4061 FE9A67 85 65 sta <facm+34062 FE9A69 A5 64 lda <facm+24063 FE9A6B 79 F5 9C adc !XBF16+2,y4064 FE9A6E 85 64 sta <facm+24065 FE9A70 A5 63 lda <facm+14066 FE9A72 79 F4 9C adc !XBF16+1,y4067 FE9A75 85 63 sta <facm+14068 FE9A77 A5 62 lda <facm4069 FE9A79 79 F3 9C adc !XBF16,y4070 FE9A7C 85 62 sta <facm4071 FE9A7E E8 inx4072 FE9A7F B0 04 bcs ?364073 FE9A81 10 DE bpl ?344074 FE9A83 30 02 bmi ?384075 FE9A85 30 DA ?36: bmi ?344076 FE9A87 8A ?38: txa4077 FE9A88 90 04 bcc ?404078 FE9A8A 49 FF eor #$FF4079 FE9A8C 69 0A adc #$0A4080 FE9A8E 69 2F ?40: adc #$2F4081 FE9A90 C8 iny4082 FE9A91 C8 iny4083 FE9A92 C8 iny4084 FE9A93 C8 iny4085 FE9A94 84 73 sty <decidx24086 FE9A96 A4 71 ldy <decidx4087 FE9A98 C8 iny4088 FE9A99 AA tax4089 FE9A9A 29 7F and #$7F4090 FE9A9C AB plb4091 FE9A9D 99 59 FF sta !TTNUMSTR-1,y4092 FE9AA0 C6 5D dec <dccnt4093 FE9AA2 D0 06 bne ?424094 FE9AA4 A9 2E lda #$2E4095 FE9AA6 C8 iny4096 FE9AA7 99 59 FF sta !TTNUMSTR-1,y4097 FE9AAA 84 71 ?42: sty <decidx4098 FE9AAC A4 73 ldy <decidx24099 FE9AAE 8A txa4100 FE9AAF 49 FF eor #$FF4101 FE9AB1 29 80 and #$804102 FE9AB3 AA tax4103 FE9AB4 C0 24 cpy #$244104 FE9AB6 F0 04 beq ?444105 FE9AB8 C0 3C cpy #$3C4106 FE9ABA D0 A2 bne ?324107 FE9ABC A4 71 ?44: ldy <decidx4108 FE9ABE B9 59 FF ?46: lda !TTNUMSTR-1,y4109 FE9AC1 88 dey4110 FE9AC2 C9 30 cmp #$304111 FE9AC4 F0 F8 beq ?464112 FE9AC6 C9 2E cmp #$2ETue Jul 17 11:00:22 2018 Page 284113 FE9AC8 F0 01 beq ?484114 FE9ACA C8 iny4115 FE9ACB A9 2B ?48: lda #$2B4116 FE9ACD A6 5E ldx <expval4117 FE9ACF F0 2E beq ?604118 FE9AD1 10 08 bpl ?504119 FE9AD3 A9 00 lda #$004120 FE9AD5 38 sec4121 FE9AD6 E5 5E sbc <expval4122 FE9AD8 AA tax4123 FE9AD9 A9 2D lda #$2D4124 FE9ADB 99 5B FF ?50: sta !TTNUMSTR+1,y4125 FE9ADE A9 45 lda #$454126 FE9AE0 99 5A FF sta !TTNUMSTR,y4127 FE9AE3 8A txa4128 FE9AE4 A2 2F ldx #$2F4129 FE9AE6 38 sec4130 FE9AE7 E8 ?52: inx4131 FE9AE8 E9 0A sbc #$0A4132 FE9AEA B0 FB bcs ?524133 FE9AEC 69 3A adc #$3A4134 FE9AEE 99 5D FF sta !TTNUMSTR+3,y4135 FE9AF1 8A txa4136 FE9AF2 99 5C FF sta !TTNUMSTR+2,y4137 FE9AF5 A9 00 lda #$004138 FE9AF7 99 5E FF sta !TTNUMSTR+4,y4139 FE9AFA F0 08 beq ?624140 FE9AFC 99 59 FF ?58: sta !TTNUMSTR-1,y4141 FE9AFF A9 00 ?60: lda #$004142 FE9B01 99 5A FF sta !TTNUMSTR,y4143 FE9B04 AB ?62: plb ; restore data bank4144 FE9B05 A9 5A lda #<TTNUMSTR4145 FE9B07 A0 FF ldy #>TTNUMSTR4146 FE9B09 60 rts41474148 ;---------------------------------------------------------------------------4149 ; sqrt() function implementation - sqrt(x) = exp(ln(x) * 0.5)4150 ;---------------------------------------------------------------------------41514152 ; compute SQRT([FAC #1])4153 FE9B0A XBF71:4154 FE9B0A 20 CE 97 jsr XBC0C ; round FAC #1 and move to FAC #24155 FE9B0D A9 AC lda #<FCON054156 FE9B0F A0 9C ldy #>FCON054157 FE9B11 20 49 97 jsr XBBA2K ; move AYK to FAC #141584159 ;---------------------------------------------------------------------------4160 ; power() function implementation - power(x,y) = x^y = exp(ln(x) * y)4161 ;---------------------------------------------------------------------------41624163 ; compute [FAC #2] ^ [FAC #1] - A = FAC #1 exponent4164 FE9B14 XBF7B:4165 FE9B14 F0 3F beq XBFED ; if null go to exp(0)4166 FE9B16 A5 69 lda <argexp4167 FE9B18 D0 03 bne ?024168 FE9B1A 4C A1 94 jmp XB8F94169 FE9B1D 20 82 97 ?02: jsr XBBX3 ; round FAC #1 and move to FAC3TMPTue Jul 17 11:00:22 2018 Page 294170 FE9B20 A5 6E lda <argsgn4171 FE9B22 10 0F bpl ?044172 FE9B24 20 9F 98 jsr XBCCC ; to int4173 FE9B27 A9 F2 lda #<FAC3TMP4174 FE9B29 A0 FF ldy #>FAC3TMP4175 FE9B2B 20 66 98 jsr XBC5BT ; compare in temp. fac bank4176 FE9B2E D0 03 bne ?044177 FE9B30 98 tya4178 FE9B31 A4 5B ldy <faclsb4179 FE9B33 20 C0 97 ?04: jsr XBBFE ; copy FAC #1 to #FAC #24180 FE9B36 5A phy4181 FE9B37 20 8E 95 jsr XB9EA ; compute ln()4182 FE9B3A A9 F2 lda #<FAC3TMP4183 FE9B3C A0 FF ldy #>FAC3TMP4184 FE9B3E A6 A2 ldx <vbnk ; bank that old FAC3TMP4185 FE9B40 20 DF 95 jsr XBA28 ; mult with FAC3TMP4186 FE9B43 20 55 9B jsr XBFED ; compute exp()4187 FE9B46 68 pla4188 FE9B47 4A lsr a4189 FE9B48 90 0A bcc XBFBE ; if even41904191 XBFB4: ; change sign4192 FE9B4A A5 61 lda <facexp4193 FE9B4C F0 06 beq XBFBE4194 FE9B4E A5 66 lda <facsgn4195 FE9B50 49 FF eor #$FF4196 FE9B52 85 66 sta <facsgn4197 FE9B54 60 XBFBE: rts41984199 ;---------------------------------------------------------------------------4200 ; exp() function implementation4201 ;---------------------------------------------------------------------------42024203 ; compute exp([FAC #1])4204 FE9B55 XBFED:4205 FE9B55 A9 DA lda #<XBFBF ; 1/LN(2)4206 FE9B57 A0 9C ldy #>XBFBF4207 FE9B59 A2 FE ldx #^XBFBF ; current program bank4208 FE9B5B 20 DF 95 jsr XBA28 ; mult 1/LN(2)4209 FE9B5E A5 70 lda <facbits4210 FE9B60 69 50 adc #$504211 FE9B62 90 03 bcc ?024212 FE9B64 20 E5 97 jsr XBC23 ; increment mantissa4213 FE9B67 85 59 ?02: sta <argbits4214 FE9B69 20 D1 97 jsr XBC0F ; copy FAC # 1 to FAC #2 (no rounding before)4215 FE9B6C A5 61 lda <facexp4216 FE9B6E C9 88 cmp #$884217 FE9B70 90 03 bcc ?064218 FE9B72 20 5C 96 ?04: jsr XBAD44219 FE9B75 20 9F 98 ?06: jsr XBCCC ; compute integer part4220 FE9B78 A5 5B lda <faclsb4221 FE9B7A 18 clc4222 FE9B7B 69 81 adc #$814223 FE9B7D F0 F3 beq ?044224 FE9B7F 38 sec4225 FE9B80 E9 01 sbc #$014226 FE9B82 48 phaTue Jul 17 11:00:22 2018 Page 304227 FE9B83 A2 05 ldx #$05 ; exchange FAC #1 and FAC #24228 FE9B85 B5 69 ?08: lda <argexp,x4229 FE9B87 B4 61 ldy <facexp,x4230 FE9B89 95 61 sta <facexp,x4231 FE9B8B 94 69 sty <argexp,x4232 FE9B8D CA dex4233 FE9B8E 10 F5 bpl ?084234 FE9B90 A5 59 lda <argbits4235 FE9B92 85 70 sta <facbits4236 FE9B94 20 F6 93 jsr XB853 ; [FAC #2] - [FAC #1]4237 FE9B97 20 4A 9B jsr XBFB4 ; negate FAC #14238 FE9B9A A9 17 lda #<XBFC4 ; compute exp() polynomial4239 FE9B9C A0 9D ldy #>XBFC44240 FE9B9E 20 C3 9B jsr XE0594241 FE9BA1 64 6F stz <arisgn4242 FE9BA3 68 pla4243 FE9BA4 4C 40 96 jmp XBAB9 ; add exponent4244 FE9BA74245 ;---------------------------------------------------------------------------4246 ; series evaluation implementation4247 ;---------------------------------------------------------------------------42484249 FE9BA7 XE043:4250 FE9BA7 85 71 sta <fbufpt4251 FE9BA9 84 72 sty <fbufpt+14252 FE9BAB 20 7A 97 jsr XBBCA ; round FAC #1 and move it to FAC2TMP4253 FE9BAE A9 ED lda #<FAC2TMP4254 FE9BB0 A0 FF ldy #>FAC2TMP4255 FE9BB2 A6 A2 ldx <vbnk ; bank that old FAC2TMP4256 FE9BB4 20 DF 95 jsr XBA28 ; compute X^24257 FE9BB7 20 C7 9B jsr XE05D ; compute P(X^2)4258 FE9BBA A9 ED lda #<FAC2TMP4259 FE9BBC A0 FF ldy #>FAC2TMP4260 FE9BBE A6 A2 ldx <vbnk ; bank that old FAC2TMP4261 FE9BC0 4C DF 95 jmp XBA28 ; compute X*P(X^2)4262 FE9BC34263 FE9BC3 XE059:4264 FE9BC3 85 71 sta <fbufpt4265 FE9BC5 84 72 sty <fbufpt+14266 FE9BC7 XE05D:4267 FE9BC7 20 72 97 jsr XBBC7 ; round FAC #1 and move it to FAC1TMP4268 FE9BCA A9 FE lda #^XE0434269 FE9BCC 85 73 sta <fbufpt+2 ; polynomial coefficient in program bank4270 FE9BCE B7 71 lda [fbufpt],y4271 FE9BD0 85 67 sta <sercnt ; number of coefficients4272 FE9BD2 A5 71 lda <fbufpt4273 FE9BD4 1A inc a4274 FE9BD5 D0 02 bne ?024275 FE9BD7 E6 72 inc <fbufpt+14276 FE9BD9 85 71 ?02: sta <fbufpt4277 FE9BDB A4 72 ldy <fbufpt+14278 FE9BDD A6 73 ldx <fbufpt+24279 FE9BDF 20 DF 95 ?04: jsr XBA28 ; [FAC #1] * coeff[n]4280 FE9BE2 A5 71 lda <fbufpt4281 FE9BE4 A4 72 ldy <fbufpt+14282 FE9BE6 18 clc4283 FE9BE7 69 05 adc #$05Tue Jul 17 11:00:22 2018 Page 314284 FE9BE9 90 01 bcc ?064285 FE9BEB C8 iny4286 FE9BEC 85 71 ?06: sta <fbufpt4287 FE9BEE 84 72 sty <fbufpt+14288 FE9BF0 20 09 94 jsr XB867 ; [FAC #1] + coeff[n+1]4289 FE9BF3 A9 E8 lda #<FAC1TMP4290 FE9BF5 A0 FF ldy #>FAC1TMP4291 FE9BF7 A6 A2 ldx <vbnk ; bank that old FAC1TMP4292 FE9BF9 C6 67 dec <sercnt4293 FE9BFB D0 E2 bne ?044294 FE9BFD 60 XE0X: rts42954296 ;---------------------------------------------------------------------------4297 ; arctan() function implementation4298 ;---------------------------------------------------------------------------42994300 FE9BFE XE30E:4301 FE9BFE A5 66 lda <facsgn4302 FE9C00 48 pha4303 FE9C01 10 03 bpl ?024304 FE9C03 20 4A 9B jsr XBFB4 ; change sign4305 FE9C06 A5 61 ?02: lda <facexp4306 FE9C08 48 pha4307 FE9C09 C9 81 cmp #$814308 FE9C0B 90 09 bcc ?04 ; x <= 14309 FE9C0D A9 A2 lda #<XB9BC ; 1.04310 FE9C0F A0 9C ldy #>XB9BC4311 FE9C11 A2 FE ldx #^XE30E4312 FE9C13 20 91 96 jsr XBB0F ; x = 1/x4313 FE9C16 A9 6E ?04: lda #<XE33E ; polynomial table4314 FE9C18 A0 9D ldy #>XE33E4315 FE9C1A 20 A7 9B jsr XE0434316 FE9C1D 68 pla4317 FE9C1E C9 81 cmp #$814318 FE9C20 90 07 bcc ?064319 FE9C22 A9 40 lda #<XE2E0 ; pi / 24320 FE9C24 A0 9D ldy #>XE2E04321 FE9C26 20 F3 93 jsr XB850 ; pi/2 - [FAC 31]4322 FE9C29 68 ?06: pla4323 FE9C2A 10 D1 bpl XE0X ; exit4324 FE9C2C 4C 4A 9B jmp XBFB4 ; change sign43254326 ;---------------------------------------------------------------------------4327 ; cos() function implementation -- cos(x) = sin(x + pi/2)4328 ;---------------------------------------------------------------------------43294330 ; compute cos([FAC #1])4331 FE9C2F XE264:4332 FE9C2F A9 40 lda #<XE2E0 ; PI/24333 FE9C31 A0 9D ldy #>XE2E04334 FE9C33 20 09 94 jsr XB86743354336 ;---------------------------------------------------------------------------4337 ; sin() function implementation4338 ;---------------------------------------------------------------------------43394340 ; compute sin([FAC #1])Tue Jul 17 11:00:22 2018 Page 324341 FE9C36 XE26B:4342 FE9C36 20 CE 97 jsr XBC0C ; round FAC #1 and move to FAC #24343 FE9C39 A9 45 lda #<XE2E5 ; 2 * PI4344 FE9C3B A0 9D ldy #>XE2E54345 FE9C3D A6 6E ldx <argsgn4346 FE9C3F 20 8A 96 jsr XBB07 ; divide by 2*pi (reduction angle)4347 FE9C42 20 CE 97 jsr XBC0C ; round FAC #1 and move to FAC #24348 FE9C45 20 9F 98 jsr XBCCC ; to int4349 FE9C48 64 6F stz <arisgn4350 FE9C4A 20 F6 93 jsr XB853 ; now argument is reduced to 2 * pi max.4351 FE9C4D A9 4A lda #<XE2EA ; 0.254352 FE9C4F A0 9D ldy #>XE2EA4353 FE9C51 20 F3 93 jsr XB850 ; 0.25 - [FAC #1]4354 FE9C54 A5 66 lda <facsgn4355 FE9C56 48 pha4356 FE9C57 10 0D bpl XE29D4357 FE9C59 20 ED 93 jsr XB849 ; add 0.54358 FE9C5C A5 66 lda <facsgn4359 FE9C5E 30 09 bmi XE2A04360 FE9C60 A5 5A lda <cprmsk4361 FE9C62 49 FF eor #$FF4362 FE9C64 85 5A sta <cprmsk4363 FE9C66 XE29D:4364 FE9C66 20 4A 9B jsr XBFB4 ; change sign4365 FE9C69 XE2A0:4366 FE9C69 A9 4A lda #<XE2EA ; 0.254367 FE9C6B A0 9D ldy #>XE2EA4368 FE9C6D 20 09 94 jsr XB867 ; add 0.254369 FE9C70 68 pla4370 FE9C71 10 03 bpl ?024371 FE9C73 20 4A 9B jsr XBFB4 ; change sign4372 FE9C76 A9 4F ?02: lda #<XE2EF4373 FE9C78 A0 9D ldy #>XE2EF4374 FE9C7A 4C A7 9B jmp XE043 ; compute sin() polynomial43754376 ;---------------------------------------------------------------------------4377 ; tan() function implementation -- tan(x) = sin(x) / cos(x)4378 ;---------------------------------------------------------------------------43794380 ; compute tan([FAC #1])4381 FE9C7D XE2B4:4382 FE9C7D 20 7A 97 jsr XBBCA ; move FAC #1 rounded to FAC1TMP4383 FE9C80 64 5A stz <cprmsk4384 FE9C82 20 36 9C jsr XE26B ; compute sin()4385 FE9C85 20 82 97 jsr XBBX3 ; move FAC #1 rounded to FAC3TMP4386 FE9C88 A9 ED lda #<FAC2TMP4387 FE9C8A A0 FF ldy #>FAC2TMP4388 FE9C8C 20 41 97 jsr XBBA2T4389 FE9C8F 64 66 stz <facsgn4390 FE9C91 A5 5A lda <cprmsk4391 FE9C93 20 9F 9C jsr ?024392 FE9C96 A9 F2 lda #<FAC3TMP4393 FE9C98 A0 FF ldy #>FAC3TMP4394 FE9C9A A6 A2 ldx <vbnk4395 FE9C9C 4C 91 96 jmp XBB0F ; div4396 FE9C9F 48 ?02 pha4397 FE9CA0 80 C4 bra XE29DTue Jul 17 11:00:22 2018 Page 3343984399 ;---------------------------------------------------------------------------4400 ; some floating point costants4401 ;---------------------------------------------------------------------------44024403 FE9CA2 81 00 00 00 00 XB9BC: .DB $81,$00,$00,$00,$00 ; 1.04404 FE9CA7 84 20 00 00 00 FCON10: .DB $84,$20,$00,$00,$00 ; 10.04405 FE9CAC 80 00 00 00 00 FCON05: .DB $80,$00,$00,$00,$00 ; 0.544064407 FE9CB1 03 XB9C1: .DB $03 ; LN polynomial table degree 44408 FE9CB2 7F 5E 56 CB 79 .DB $7F,$5E,$56,$CB,$79 ; B9C24409 FE9CB7 80 13 9B 0B 64 .DB $80,$13,$9B,$0B,$64 ; B9C74410 FE9CBC 80 76 38 93 16 .DB $80,$76,$38,$93,$16 ; B9CC4411 FE9CC1 82 38 AA 3B 20 .DB $82,$38,$AA,$3B,$20 ; B9D144124413 FE9CC6 80 35 04 F3 34 XB9D6: .DB $80,$35,$04,$F3,$34 ; 0.5 * SQRT(2)4414 FE9CCB 81 35 04 F3 34 XB9DB: .DB $81,$35,$04,$F3,$34 ; SQRT(2)4415 FE9CD0 80 80 00 00 00 XB9E0: .DB $80,$80,$00,$00,$00 ; -0.54416 FE9CD5 80 31 72 17 F8 XB9E5: .DB $80,$31,$72,$17,$F8 ; LN(2)4417 FE9CDA 81 38 AA 3B 29 XBFBF: .DB $81,$38,$AA,$3B,$29 ; 1/LN(2)4418 FE9CDF 82 13 5D 8D DE XLN10 .DB $82,$13,$5D,$8D,$DE ; LN(10)44194420 FE9CE4 9B 3E BC 1F FD XBDB3: .DB $9B,$3E,$BC,$1F,$FD ; used by float to string conversion4421 FE9CE9 9E 6E 6B 27 FD XBDB8: .DB $9E,$6E,$6B,$27,$FD ; 1E84422 FE9CEE 9E 6E 6B 28 00 XBDBD: .DB $9E,$6E,$6B,$28,$00 ; 1E944234424 FE9CF3 FA 0A 1F 00 XBF16: .DB $FA,$0A,$1F,$00 ; -1000000004425 FE9CF7 00 98 96 80 .DB $00,$98,$96,$80 ; 100000004426 FE9CFB FF F0 BD C0 .DB $FF,$F0,$BD,$C0 ; -10000004427 FE9CFF 00 01 86 A0 .DB $00,$01,$86,$A0 ; 1000004428 FE9D03 FF FF D8 F0 .DB $FF,$FF,$D8,$F0 ; -100004429 FE9D07 00 00 03 E8 .DB $00,$00,$03,$E8 ; 10004430 FE9D0B FF FF FF 9C .DB $FF,$FF,$FF,$9C ; -1004431 FE9D0F 00 00 00 0A .DB $00,$00,$00,$0A ; 104432 FE9D13 FF FF FF FF .DB $FF,$FF,$FF,$FF ; -144334434 FE9D17 07 XBFC4: .DB $07 ; EXP polynomial table degree 84435 FE9D18 71 34 58 3E 56 .DB $71,$34,$58,$3E,$564436 FE9D1D 74 16 7E B3 1B .DB $74,$16,$7E,$B3,$1B4437 FE9D22 77 2F EE E3 85 .DB $77,$2F,$EE,$E3,$854438 FE9D27 7A 1D 84 1C 2A .DB $7A,$1D,$84,$1C,$2A4439 FE9D2C 7C 63 59 58 0A .DB $7C,$63,$59,$58,$0A4440 FE9D31 7E 75 FD E7 C6 .DB $7E,$75,$FD,$E7,$C64441 FE9D36 80 31 72 18 10 .DB $80,$31,$72,$18,$104442 FE9D3B 81 00 00 00 00 .DB $81,$00,$00,$00,$0044434444 FE9D40 81 49 0F DA A2 XE2E0 .DB $81,$49,$0F,$DA,$A2 ; pi * 0.54445 FE9D45 83 49 0F DA A2 XE2E5 .DB $83,$49,$0F,$DA,$A2 ; pi * 24446 FE9D4A 7F 00 00 00 00 XE2EA .DB $7F,$00,$00,$00,$00 ; 0.2544474448 FE9D4F 05 XE2EF .DB $05 ; SIN polynomial table degree 64449 FE9D50 84 E6 1A 2D 1B .DB $84,$E6,$1A,$2D,$1B4450 FE9D55 86 28 07 FB F8 .DB $86,$28,$07,$FB,$F84451 FE9D5A 87 99 68 89 01 .DB $87,$99,$68,$89,$014452 FE9D5F 87 23 35 DF E1 .DB $87,$23,$35,$DF,$E14453 FE9D64 86 A5 5D E7 28 .DB $86,$A5,$5D,$E7,$284454 FE9D69 83 49 0F DA A2 .DB $83,$49,$0F,$DA,$A2Tue Jul 17 11:00:22 2018 Page 3444554456 FE9D6E 0B XE33E: .DB $0B ; ARCTAN polynomial table degree 124457 FE9D6F 76 B3 83 BD D3 .DB $76,$B3,$83,$BD,$D34458 FE9D74 79 1E F4 A6 F5 .DB $79,$1E,$F4,$A6,$F54459 FE9D79 7B 83 FC B0 10 .DB $7B,$83,$FC,$B0,$104460 FE9D7E 7C 0C 1F 67 CA .DB $7C,$0C,$1F,$67,$CA4461 FE9D83 7C DE 53 CB C1 .DB $7C,$DE,$53,$CB,$C14462 FE9D88 7D 14 64 70 4C .DB $7D,$14,$64,$70,$4C4463 FE9D8D 7D B7 EA 51 7A .DB $7D,$B7,$EA,$51,$7A4464 FE9D92 7D 63 30 88 7E .DB $7D,$63,$30,$88,$7E4465 FE9D97 7E 92 44 99 3A .DB $7E,$92,$44,$99,$3A4466 FE9D9C 7E 4C CC 91 C7 .DB $7E,$4C,$CC,$91,$C74467 FE9DA1 7F AA AA AA 13 .DB $7F,$AA,$AA,$AA,$134468 FE9DA6 81 00 00 00 00 .DB $81,$00,$00,$00,$00Lines Assembled : 4446 Errors : 0