Programowanie: Dzielenie 16-bitowe z wynikiem 32-bitowym
From Atariki
(Różnice między wersjami)
Wersja z dnia 14:16, 21 gru 2005 KMK (Dyskusja | wkład) (→6502 - ditto) ← Previous diff |
Aktualna wersja KMK (Dyskusja | wkład) (w sumie po co dwa razy odejmowanie ...) |
||
Linia 9: | Linia 9: | ||
beq ?e | beq ?e | ||
+ | ; <-- czesc pomijalna | ||
lda fr0 ;dzielnik wiekszy od dzielnej | lda fr0 ;dzielnik wiekszy od dzielnej | ||
cmp fr1 ;daje w wyniku zero, a dzielna | cmp fr1 ;daje w wyniku zero, a dzielna | ||
Linia 23: | Linia 24: | ||
sta fr0+1 | sta fr0+1 | ||
rts | rts | ||
+ | ; --> | ||
?div lda #$00 ;kasujemy reszte | ?div lda #$00 ;kasujemy reszte | ||
Linia 35: | Linia 37: | ||
rol fr0+3 | rol fr0+3 | ||
+ | sec | ||
lda fr0+2 | lda fr0+2 | ||
- | cmp fr1 | + | sbc fr1 |
+ | tay | ||
lda fr0+3 | lda fr0+3 | ||
sbc fr1+1 | sbc fr1+1 | ||
bcc ?s | bcc ?s | ||
- | lda fr0+2 | + | sty fr0+2 |
- | sbc fr1 | + | |
- | sta fr0+2 | + | |
- | lda fr0+3 | + | |
- | sbc fr1+1 | + | |
sta fr0+3 | sta fr0+3 | ||
Linia 77: | Linia 77: | ||
beq ?e | beq ?e | ||
+ | ; <-- czesc pomijalna | ||
lda fr0 ;dzielnik wiekszy od dzielnej | lda fr0 ;dzielnik wiekszy od dzielnej | ||
cmp fr1 ;daje w wyniku zero, a dzielna | cmp fr1 ;daje w wyniku zero, a dzielna | ||
Linia 85: | Linia 86: | ||
stz fr0 | stz fr0 | ||
rts | rts | ||
+ | ; --> | ||
?div stz fr0+2 ;kasujemy reszte | ?div stz fr0+2 ;kasujemy reszte | ||
Linia 93: | Linia 95: | ||
rol fr0+2 | rol fr0+2 | ||
+ | sec | ||
lda fr0+2 | lda fr0+2 | ||
- | cmp fr1 | + | sbc fr1 |
bcc ?s | bcc ?s | ||
- | lda fr0+2 | ||
- | sbc fr1 | ||
sta fr0+2 | sta fr0+2 | ||
Aktualna wersja
Procedura dzieli 16-bitową liczbę bez znaku umieszczoną we FR0 ($D4) przez takąż liczbę umieszczoną w FR1 ($E0). 16-bitowy wynik całkowity umieszczany jest w FR0, a 16-bitowa reszta z dzielenia w FR0+2. Opuszczenie procedury z ustawionym znacznikiem C sygnalizuje przepełnienie.
[Edytuj]
6502
int16div lda fr1 ;dzielnik musi byc wiekszy od zera ora fr1+1 beq ?e ; <-- czesc pomijalna lda fr0 ;dzielnik wiekszy od dzielnej cmp fr1 ;daje w wyniku zero, a dzielna lda fr0+1 ;stanowi reszte sbc fr1+1 bcs ?div lda fr0 sta fr0+2 lda fr0+1 sta fr0+3 lda #$00 sta fr0 sta fr0+1 rts ; --> ?div lda #$00 ;kasujemy reszte sta fr0+2 sta fr0+3 ldx #16 ?l asl fr0 rol fr0+1 rol fr0+2 rol fr0+3 sec lda fr0+2 sbc fr1 tay lda fr0+3 sbc fr1+1 bcc ?s sty fr0+2 sta fr0+3 inc fr0 bne ?s inc fr0+1 bne ?s ?e sec rts ?s dex bne ?l clc rts
Główna pętla zajmuje - w zależności od składników - od 624 do 1152 cykli maszynowych.
[Edytuj]
65C816
int16div .aw .ib rep #$20 ;CPU musi byc w trybie natywnym sep #$10 lda fr1 ;dzielnik musi byc wiekszy od zera beq ?e ; <-- czesc pomijalna lda fr0 ;dzielnik wiekszy od dzielnej cmp fr1 ;daje w wyniku zero, a dzielna bcs ?div ;stanowi reszte lda fr0 sta fr0+2 stz fr0 rts ; --> ?div stz fr0+2 ;kasujemy reszte ldx #16 ?l asl fr0 rol fr0+2 sec lda fr0+2 sbc fr1 bcc ?s sta fr0+2 inc fr0 bne ?s ?e sec rts ?s dex bne ?l clc rts
Elegancka procedura biblioteczna powinna jeszcze przechowywać na stosie i przed powrotem do miejsca wywołania przywracać rozmiar rejestrów. W powyższym przykładzie pominięto to, żeby nie zaciemniać zagadnienia.
[Edytuj]
Uwagi
Procedura jest "kompatybilna" z procedurą mnożenia (quod vide).