Find an Integer Square Root on Amiga | Pattern Library

Sometimes the useful answer is the whole-number distance, with no floating point in sight.

Amiga shows the working diagnostic row 1F 27 00 E8 03 00 00 00. — A completed native caller displays 1F 27 00 E8 03 00 00 00. The public field order below identifies these bytes; the capture is a state diagnostic, not a recording of a finished game.

Use it in a project

NUMBER_LO/NUMBER_HI hold an unsigned 16-bit number. Return ROOT=floor(sqrt(number)) and the remainder number−ROOT² as a low-byte-first word. Inputs are preserved. The remainder can reach 510, so it is not a single byte.

Include component.inc once (MOS callers include the shared 6502 implementation through it). Reserve the public fields in this exact order: ROOT, REMAINDER_LO, REMAINDER_HI, NUMBER_LO, NUMBER_HI. They occupy 5 contiguous bytes. Also reserve the private workspace: T0, T1, T2, T3 (4 bytes). Array names in the explanation refer to the numbered consecutive bindings, such as C0, C1 and so on.

Call integer_square_root_init once at the intended reset boundary. It clears inputs as well as outputs and workspace. Set your starting state and configuration afterwards, then call integer_square_root_step whenever the described event occurs. The helper reads no device, waits for no frame and installs no interrupt handler.

Ownership	Fields
Retain between calls	none
Supply before service	NUMBER_LO and NUMBER_HI
Read after service	ROOT, REMAINDER_LO and REMAINDER_HI

The first n odd numbers sum to n². Subtract 1, 3, 5 and so on until the next subtraction would go negative. The number of accepted subtractions is the root.

Reusable initialisation and service

component.inc

; SPDX-License-Identifier: MIT
; Find an Integer Square Root. See README.md for the byte interface and ownership.
; init clears the public fields and private workspace; set inputs afterwards.
integer_square_root_init:
    move.b #0,d0
    move.b d0,ROOT
    move.b d0,REMAINDER_LO
    move.b d0,REMAINDER_HI
    move.b d0,NUMBER_LO
    move.b d0,NUMBER_HI
    move.b d0,T0
    move.b d0,T1
    move.b d0,T2
    move.b d0,T3
    rts
integer_square_root_step:
; ROOT = 0
    move.b #0,d0
    move.b d0,ROOT
; REMAINDER_LO = NUMBER_LO
    move.b NUMBER_LO,d0
    move.b d0,REMAINDER_LO
; REMAINDER_HI = NUMBER_HI
    move.b NUMBER_HI,d0
    move.b d0,REMAINDER_HI
; T0 = 1
    move.b #1,d0
    move.b d0,T0
; T1 = 0
    move.b #0,d0
    move.b d0,T1
; T2 = 1
    move.b #1,d0
    move.b d0,T2
; while T2 == 1:
integer_square_root_L1:
    move.b T2,d0
    cmp.b #1,d0
    beq integer_square_root_L2
    jmp integer_square_root_L3
integer_square_root_L2:
; if REMAINDER_HI < T1:
    move.b REMAINDER_HI,d0
    move.b T1,d1
    cmp.b d1,d0
    bcs integer_square_root_L4
    jmp integer_square_root_L5
integer_square_root_L4:
; return
    rts
    jmp integer_square_root_L6
integer_square_root_L5:
integer_square_root_L6:
; if REMAINDER_HI == T1:
    move.b REMAINDER_HI,d0
    move.b T1,d1
    cmp.b d1,d0
    beq integer_square_root_L7
    jmp integer_square_root_L8
integer_square_root_L7:
; if REMAINDER_LO < T0:
    move.b REMAINDER_LO,d0
    move.b T0,d1
    cmp.b d1,d0
    bcs integer_square_root_L10
    jmp integer_square_root_L11
integer_square_root_L10:
; return
    rts
    jmp integer_square_root_L12
integer_square_root_L11:
integer_square_root_L12:
    jmp integer_square_root_L9
integer_square_root_L8:
integer_square_root_L9:
; T3 = 0
    move.b #0,d0
    move.b d0,T3
; if REMAINDER_LO < T0:
    move.b REMAINDER_LO,d0
    move.b T0,d1
    cmp.b d1,d0
    bcs integer_square_root_L13
    jmp integer_square_root_L14
integer_square_root_L13:
; T3 = 1
    move.b #1,d0
    move.b d0,T3
    jmp integer_square_root_L15
integer_square_root_L14:
integer_square_root_L15:
; REMAINDER_LO = REMAINDER_LO - T0
    move.b REMAINDER_LO,d0
    move.b T0,d1
    sub.b d1,d0
    move.b d0,REMAINDER_LO
; REMAINDER_HI = REMAINDER_HI - T1
    move.b REMAINDER_HI,d0
    move.b T1,d1
    sub.b d1,d0
    move.b d0,REMAINDER_HI
; REMAINDER_HI = REMAINDER_HI - T3
    move.b REMAINDER_HI,d0
    move.b T3,d1
    sub.b d1,d0
    move.b d0,REMAINDER_HI
; ROOT = ROOT + 1
    move.b ROOT,d0
    add.b #1,d0
    move.b d0,ROOT
; T3 = T0
    move.b T0,d0
    move.b d0,T3
; T0 = T0 + 2
    move.b T0,d0
    add.b #2,d0
    move.b d0,T0
; if T0 < T3:
    move.b T0,d0
    move.b T3,d1
    cmp.b d1,d0
    bcs integer_square_root_L16
    jmp integer_square_root_L17
integer_square_root_L16:
; T1 = T1 + 1
    move.b T1,d0
    add.b #1,d0
    move.b d0,T1
    jmp integer_square_root_L18
integer_square_root_L17:
integer_square_root_L18:
    jmp integer_square_root_L1
integer_square_root_L3:
    rts
component_end:

Machine contract

D0–D3/A0–A1/CCR may change; D4–D7/A2–A6 are preserved. A direct BSR/JSR needs four return-address stack bytes. The helper makes no nested calls and adds no other stack use. Its named scratch is private workspace, not a pointer, a ROM variable or shared zero-page space that the game may borrow concurrently. The samples use fixed symbol bindings. Scope or rename every label and binding, including the component_end measurement marker, and allocate a separate block for each independent instance. Different helpers can also use common field names, so combining them needs separate namespaces or explicitly renamed bindings. Give each block one mutation owner.

The A500 OCS diagnostic cold-boots an OFS disk, owns the Copper and bitplane and disables interrupt/DMA sources. The helper does not take over hardware or allocate Chip RAM.

Adapt the caller

demo.asm supplies the bindings and include order. exercise.inc first runs table cases, then a separate client retaining state between calls. A game can call the public step entry directly and remove the diagnostic packet and display. Do not copy the diagnostic’s hardware takeover into an existing game loop.

Build with make ASM198X=/path/to/asm198x and BUILD198X=/path/to/build198x. The Makefile produces normal, custom, deliberately wrong and repaired programs. Tools must already be installed; the sample downloads nothing.

Bindings and include order

demo.asm

    section code,code_c
CUSTOM = 0
BROKEN = 0
COYOTE_STEPS = 4
BUFFER_STEPS = 4
TEXT_CAP = 8
PROPERTY_COUNT = 16
STREAM_BYTES = 13
CRC_POLY = 7
    include "caller.inc"
    include "exercise.inc"
component_start:
    include "component.inc"
    even
    include "data.inc"
    include "cases.inc"
    even
    include "display.inc"
    include "../../framework/turn-deadline/font.inc"
code_end:
STATE:
ROOT: dc.b 0
REMAINDER_LO: dc.b 0
REMAINDER_HI: dc.b 0
NUMBER_LO: dc.b 0
NUMBER_HI: dc.b 0
T0: dc.b 0
T1: dc.b 0
T2: dc.b 0
T3: dc.b 0
    even
ready: dc.w 0
preview: dcb.b 8,0
OUTPUT: dcb.b 1300,0
    even
signature: dc.b "198XWORLDKIT000",0
    even
copperlist:
    dc.w $08e,$2c81,$090,$2cc1,$092,$0038,$094,$00d0
    dc.w $100,$1200,$102,0,$104,0,$108,0,$10a,0
    dc.w $180,$012,$182,$fff
copbpl:
    dc.w $0e0,0,$0e2,0
    dc.w $ffff,$fffe
plane: dcb.b 10240,0

A client retaining state between calls

exercise.inc, excerpt

preview_rules:
    ifne CUSTOM
; NUMBER_LO = 255
    move.b #255,d0
    move.b d0,NUMBER_LO
; NUMBER_HI = 255
    move.b #255,d0
    move.b d0,NUMBER_HI
    ifne BROKEN
; NUMBER_HI = 0
    move.b #0,d0
    move.b d0,NUMBER_HI
    endif
    bsr checked_step
    bsr emit_state
    rts
    else
; NUMBER_LO = 232
    move.b #232,d0
    move.b d0,NUMBER_LO
; NUMBER_HI = 3
    move.b #3,d0
    move.b d0,NUMBER_HI
    ifne BROKEN
; NUMBER_HI = 0
    move.b #0,d0
    move.b d0,NUMBER_HI
    endif
    bsr checked_step
    bsr emit_state
    endif
    rts

Customise it

Exercise the maximum input and its two-byte remainder. The maintained custom caller changes its inputs or client sequence rather than hiding a different implementation behind the same label. Compare every state of the two clients, not only their final display.

Custom client selection

custom.asm

    section code,code_c
CUSTOM = 1
BROKEN = 0
COYOTE_STEPS = 2
BUFFER_STEPS = 2
TEXT_CAP = 4
PROPERTY_COUNT = 16
STREAM_BYTES = 13
CRC_POLY = 49
    include "caller.inc"
    include "exercise.inc"
component_start:
    include "component.inc"
    even
    include "data.inc"
    include "cases.inc"
    even
    include "display.inc"
    include "../../framework/turn-deadline/font.inc"
code_end:
STATE:
ROOT: dc.b 0
REMAINDER_LO: dc.b 0
REMAINDER_HI: dc.b 0
NUMBER_LO: dc.b 0
NUMBER_HI: dc.b 0
T0: dc.b 0
T1: dc.b 0
T2: dc.b 0
T3: dc.b 0
    even
ready: dc.w 0
preview: dcb.b 8,0
OUTPUT: dcb.b 1300,0
    even
signature: dc.b "198XWORLDKIT000",0
    even
copperlist:
    dc.w $08e,$2c81,$090,$2cc1,$092,$0038,$094,$00d0
    dc.w $100,$1200,$102,0,$104,0,$108,0,$10a,0
    dc.w $180,$012,$182,$fff
copbpl:
    dc.w $0e0,0,$0e2,0
    dc.w $ffff,$fffe
plane: dcb.b 10240,0

The custom diagnostic row is FF FE 01 FF FF 00 00 00. The full client trace in the source explains how it gets there.

Mistake and repair

The wrong caller keeps only the low input byte. Large distances collapse into small unrelated values. Supply the complete squared value, including its high byte.

The BROKEN branch lives in the sequential client, not the reusable helper. wrong.asm selects it; repaired.asm restores the normal call sequence. The repaired executable matches the normal program (for Spectrum, compare loaded CODE rather than the tape filename).

Select the deliberately wrong client

wrong.asm

    section code,code_c
CUSTOM = 0
BROKEN = 1
COYOTE_STEPS = 4
BUFFER_STEPS = 4
TEXT_CAP = 8
PROPERTY_COUNT = 16
STREAM_BYTES = 13
CRC_POLY = 7
    include "caller.inc"
    include "exercise.inc"
component_start:
    include "component.inc"
    even
    include "data.inc"
    include "cases.inc"
    even
    include "display.inc"
    include "../../framework/turn-deadline/font.inc"
code_end:
STATE:
ROOT: dc.b 0
REMAINDER_LO: dc.b 0
REMAINDER_HI: dc.b 0
NUMBER_LO: dc.b 0
NUMBER_HI: dc.b 0
T0: dc.b 0
T1: dc.b 0
T2: dc.b 0
T3: dc.b 0
    even
ready: dc.w 0
preview: dcb.b 8,0
OUTPUT: dcb.b 1300,0
    even
signature: dc.b "198XWORLDKIT000",0
    even
copperlist:
    dc.w $08e,$2c81,$090,$2cc1,$092,$0038,$094,$00d0
    dc.w $100,$1200,$102,0,$104,0,$108,0,$10a,0
    dc.w $180,$012,$182,$fff
copbpl:
    dc.w $0e0,0,$0e2,0
    dc.w $ffff,$fffe
plane: dcb.b 10240,0

Amiga shows the deliberate error row 0F 07 00 E8 00 00 00 00. — The wrong client completes with 0F 07 00 E8 00 00 00 00; its changed inputs and state explain the result.

Cost and alternatives

The measured helper extent is 358 code bytes; public state is 5 bytes and private scratch is 4 bytes. There are no immutable helper tables. Caller data, the test cases, return-address stack, alignment and display buffers cost additional space. These are assembled sizes, not cycle measurements.

Zero to 255 subtraction stages. A bitwise restoring square root has a fixed number of bit-pair stages but is less direct to read. Use a table only when the input range and memory cost justify it.

This readable version subtracts consecutive odd numbers and can take 255 iterations. Collision tests usually compare squared distances instead, avoiding the square root entirely.

Execution evidence

Four variants run on A500 OCS PAL, Kickstart1.3 OFS cold boot, using Asm198x 0.0.72, Emu198x 0.44.2 and Build198x 0.2.9. Each executes 64 boundary/varied/error table cases, dirty initialisation and the separate retained-state client. An independent mathematical/container model checks every public byte after every call. This is bounded example coverage, not exhaustive correctness for every possible state.

The records retain loaded immutable code, completion, register sentinels, full packets, source/tool/firmware hashes and captures. The public gate validates these retained observations and rejects altered evidence; it does not run restricted firmware. No firmware is distributed. The eight displayed bytes are the first eight public fields, zero-padded when shorter; the packet includes every field.

Where it fits

Byte products: Squaring a byte is a direct way to produce a full-width input for this routine.

The Craft · Maths for Games · Source and native evidence.