Квадратный корень

[Destr]

В эл.журнале Bugs #01 наткнулся на сабжевую статью: http://www.zxpress.ru/article.php?id=1819&skip=1

Вот одна из процедур которая меня заинтересовала:

4. Неявное табличное вычисление. ( Осторожно - взрыв мозгов ! )

LD A,H ; 4 : 1
CP 16 ; 7 : 2
JR NC,M1 ; 12/7: 2
OR 128 ; 7 : 2
LD H,A ; 4 : 1
LD A,(HL) ; 7 : 1
RET ; 10 : 1
M1: LD L,H ; 4 : 1
LD H,127 ; 7 : 2
LD H,(HL) ; 7 : 1
LD L,A ; 4 : 1
ADD A,(HL); 7 : 1
RET ; 10 : 1
----------------------
17 байт.

Пpи HL = от 0 до 4095, Вpемя = 46 тактов.
Пpи HL = от 4096 до 65535, Вpемя = 62 такта.

Pазмеp таблицы 5120 . ( адрес в памяти 32768 )
Таблица состоит:
0000-4095 - ответы для HL от 0-4095
4096-4351 - интерполирование диапазона 4096-19199
4352-4607 - интерполирование диапазона 19200-33791
4608-4863 - интерполирование диапазона 33792-34815
4864-5119 - интерполирование диапазона 34816-65535

Почему именно так - можете и не спрашивать.
Дополнительная 256-байтная таблица с адреса 32512.
Состав:
DS 75 ,144
DS 57 ,145
DS 4 ,146
DS 120,147

===== НУ, ЧТО - БАШНЮ ОТОРВАЛО ??? =====

Сейчас отпустит, когда расскажу о точности вычисления:

Абсолютно точно = 50046 чисел. 77%
Ошибка на -1 = 10321 число. 15%
Ошибка на +1 = 5199 чисел. 8%

Ну, что сказать в свое оправдание - надо лучше интерполировать !
( Но так впадлу ... да и времени нет совсем ... )

Вопрос: Кто-нибудь подскажет как сформировать таблицу?
(первые 4096 байт понятно, а вот остальное?)

[Advertiser]

[jerri]

у меня вторая половина сомнения вызывает

; ld a,h
M1: LD L,H ; 4 : 1
LD H,127 ; 7 : 2
LD H,(HL) ; 7 : 1
;здесь понятно - берем адрес таблицы для диапазона
LD L,A ; 4 : 1
; и здесь понятно - отбрасываем младший байт

ADD A,(HL); 7 : 1
; а вот это непонятно - зачем прибавлять? почему просто не взять?

RET ; 10 : 1

для примера - число #1000
ld a,#10

M1: LD L,#10 ; 4 : 1
LD H,127 ; 7 : 2
LD H,(HL) ; 7 : 1
LD l,#10 ; 4 : 1
ADD #10,(HL); 7 : 1
;по идее в таблице должно быть число #30
;что по факту - хз
RET ; 10 : 1[/QUOTE]
;а вычислять - наверное умножая числа

[Destr]

а вот это непонятно - зачем прибавлять? почему просто не взять?

Вот и я о том-же думаю.
Видимо эта часть таблицы так хитро сформирована, что корень получается для большего диапазона чисел чем 256, т.е. значения для целого ряда, но для получения конкретного результата нужно прибавить старшую часть...
Как это сформировать - ума не приложу %(
Экперементирую...

[jerri]

он же пишет что есть погрешность

[Destr]

он же пишет что есть погрешность

+-1 = ерунда.
К тому-же грит что если хорошо интерполировать то всё будет гут

[jerri]

а таблицы кстати есть в комплекте?

[Destr]

а таблицы кстати есть в комплекте?

В том то и вопрос как их сформировать, готовых я лично не нашел...

В самом образе журнала нашлись два файла, оба являются архивами Hrust. Но распаковыватся не хотят.

[jerri]

что пишут?

[Destr]

что пишут?

"Неизвестный формат или ошибка в архиве".

P.S. Зайди в асю что-ли, а то опять флейм какой-то получается...

[jerri]

да тут наверное можно и самому рассчитать
перемножить да посмотреть