Integer Math for 6502#

00010          .LI OFF
00030 * RECHENROUTINEN FUER GANZZAHLEN
00040 *
00050 *
00060 *
00070 *
00080 *
00090 *
00100 MLPCND   .EQ $D0
00110 MLPLER   .EQ $D1
00120 PROD     .EQ $D2
00130 *
00140 DVDND    .EQ $D4
00150 DVSOR    .EQ $D6
00160 RMNDR    .EQ $D8
00170 *
00180 RADCND   .EQ $DA
00190 ROOT     .EQ $DC
00200 TEMP     .EQ $DE
00210 SIGN     .EQ $DF
00220 *
00230 *
00240 *
00250 *
00260 *
00270 *
00280 SQUARE   LDA MLPCND    Multiplikant
00290          BPL POSITIV   =|MLPCND|
00300          SEC           Multiplikator
00310          LDA #$0
00320          SBC MLPCND
00330          STA MLPCND
00340 POSITIV  STA MLPLER
00350 *
00360 MULT     LDA #0
00370          LDX #8        8 Bits
00380 MLOOP    LSR MLPLER    Multiplikator/2
00390          BCC NOADD     Bit=0
00400          CLC           addiere
00410          ADC MLPCND    Multiplikanten
00420 NOADD    ROR           schiebe Bits
00430          ROR PROD      ins Produktregister
00440          DEX           Naechstes
00450          BNE MLOOP     Bit
00460          STA PROD+1    MSB Produkt
00470          RTS
00480 ------------------------------
00490 MULTV    LDA MLPCND    Vorzeichen
00500          EOR MLPLER    des Produktes
00510          AND #$80      bestimmen
00520          STA SIGN
00530          LDX #MLPCND   Berechnung
00540          JSR ABS8      von
00550          LDX #MLPLER   |MLPLER|*|MLPCND|
00560          JSR ABS8
00570          JSR MULT
00580          LDA SIGN      Ergebnis
00590          BPL .1        positiv?
00600          LDX #PROD     Vorzeichen des
00610          JSR CHS16     Produktes aendern
00620 .1       RTS
00630 ------------------------------
00640 ABS8     LDA $0,X      Absolutwert einer 8 Bit Zahl
00650          BPL AO        Zahl Positiv
00660 CHS8     LDA #0        Vorzeichen aendern
00670          SEC
00680          SBC $0,X
00690          STA $0,X
00700 AO       RTS
00710 *
00720 ABS16    LDA $1,X      Wie ABS8 Jedoch fuer 16 Bit Zahl
00730          BPL AO
00740 CHS16    JSR CHS8      Wie CHS8
00750          LDA #0
00760          SBC $1,X
00770          STA $1,X
00780          RTS
00790 ------------------------------
00800 DIVIDEV  LDA DVDND+1   Vorzeichen des
00810          EOR DVSOR+1   Quotienten
00820          AND #$80      bestimmen
00830          STA SIGN
00840          LDX #DVDND    Berechnung
00850          JSR ABS16     des Quotienten
00860          LDX #DVSOR    aus |DVDND|/|DVSOR|
00870          JSR ABS16
00880          JSR DIVIDE
00890          BCS .1        Division durch Null?
00900          LDA SIGN      Vorzeichen
00910          BPL .1        aendern wenn
00920          LDX #DVDND    Ergebnis
00930          JSR CHS16     negativ sein muss
00940          CLC
00950 .1       RTS
00960 ------------------------------
00970 DIVIDE   LDA DVSOR     Divisor=0
00980          ORA DVSOR+1   dann Error
00990          BNE DIVOK
01000          SEC
01010          RTS
01020 *
01030 DIVOK    LDA #$00      Werte vorbesetzen
01040          STA RMNDR
01050          STA RMNDR+1
01060 *
01070          LDX #$10      16 Bit
01080 DLOOP    ROL DVDND     Bit von DVDND nach
01090          ROL DVDND+1
01100          ROL RMNDR     RMNDR
01110          ROL RMNDR+1   schieben
01120          SEC
01130          LDA RMNDR     Wenn RMNDR-DVSOR>=0
01140          SBC DVSOR     dann
01150          TAY           RMNDR=
01160          LDA RMNDR+1   RMNDR-DVSOR
01170          SBC DVSOR+1
01180          BCC DECCNT
01190          STY RMNDR
01200          STA RMNDR+1
01210 DECCNT   DEX           Naechstes Bit
01220          BNE DLOOP
01221 *
01230          ROL DVDND     Anschaetzung
01240          ROL DVDND+1   fuer
01250          ASL RMNDR     Rundung
01260          ROL RMNDR+1
01270          BCS ROUND
01280          SEC
01290          LDA DVSOR
01300          SBC RMNDR
01310          LDA DVSOR+1
01320          SBC RMNDR+1
01330          BCS NOCHNG
01340 ROUND    INC DVDND     Ergebnis
01350          BNE NOCHNG    steht in
01360          INC DVDND+1   DVDND
01370 NOCHNG   CLC           No Error
01380          RTS
01390 *
01400 ------------------------------
01410 *
01420 SQRT     LDX #$8       Berechnung
01430          LDA #0        der Qudrat-
01440          STA ROOT      wurzel einer
01450          STA ROOT+1    16 Bit Zahl
01460          STA TEMP      Nach dem
01470          STA TEMP+1    Newton-
01480 SQRT1    ASL ROOT      Verfahren
01490          ROL ROOT+1
01500          INC ROOT
01510          BNE NEXT1
01520          INC ROOT+1
01530 NEXT1    ASL RADCND
01540          ROL RADCND+1
01550          ROL TEMP
01560          ROL TEMP+1
01570          ASL RADCND
01580          ROL RADCND+1
01590          ROL TEMP
01600          ROL TEMP+1
01610          SEC
01620          LDA TEMP
01630          SBC ROOT
01640          TAY
01650          LDA TEMP+1
01660          SBC ROOT+1
01670          BCC RESTOR
01680          STA TEMP+1
01690          STY TEMP
01700          INC ROOT
01710          BNE NEXT2
01720          INC ROOT+1
01730 NEXT2    DEX
01740          BNE SQRT1
01750          JMP FINIS
01760 RESTOR   SEC
01770          LDA ROOT
01780          SBC #1
01790          STA ROOT
01800          BCS NEXT3
01810          DEC ROOT+1
01820 NEXT3    DEX
01830          BNE SQRT1
01840 FINIS    ROR ROOT+1
01850          ROR ROOT
01860          RTS
01870 *
01880 WURZEL   LDA RADCND+1  Einsprung
01890          BMI WERR      in Wurzel-
01891          JSR SQRT      Routine
01900          JSR SQRT
01910          CLC           wenn
01920          RTS           Radikant positiv
01930 WERR     SEC           sonst Error
01940          RTS