WIP FINT/FDIV for BK0010

Author: nzeemin

BasicCompiler-VS2022.vcxproj

@@ -107,6 +107,7 @@
       <PreprocessorDefinitions>WIN32;NDEBUG;_CONSOLE;%(PreprocessorDefinitions)</PreprocessorDefinitions>
       <ConformanceMode>true</ConformanceMode>
       <RuntimeLibrary>MultiThreaded</RuntimeLibrary>
+      <LanguageStandard>stdcpp17</LanguageStandard>
     </ClCompile>
     <Link>
       <SubSystem>Console</SubSystem>

runtime-BK0010.tmac

@@ -225,7 +225,7 @@ IMULOV:	MOV	(SP), R5	; 
 	MOV	#1706., R0	; Ошибка: переполнение
 	JMP	ERRR
 ;
-; DAUG ? умножение двух 16-разрядных положительных чисел (R1 * R0).
+; DAUG - умножение двух 16-разрядных положительных чисел (R1 * R0).
 ; Вход: R1, R0 ? положительные числа (0..32767).
 ; Выход: результат в R5 (старшее слово) и R1 (младшее).
 ; Если результат не помещается в 32 бита (переполнение), возвращает R5 = -1.
@@ -441,6 +441,41 @@ FFIX:
 ;
 	JMP	@R5		; возврат
 
+;#####################################################################
+;## FINT
+;## Need FFIX
+;## Need FCMP
+; Для плавающего числа (float) получает целое, меньшее чем аргумент.
+; Для положительных чисел идёт к нулю, для отрицательных чисел идёт от нуля.
+; На основе кода процедуры INT из реализации Бейсик Вильнюс для УКНЦ.
+; На входе: (SP) - адрес возврата
+;           (SP+2)(SP+4) - плвающее число (старшее слово, младшее слово)
+; Результат: (SP)(SP+2) - плвающее число (старшее слово, младшее слово)
+FINT:
+	TST	2(SP)		; проверить знак числа (старшее слово)
+	BPL	FFIX		; положительное или 0 => так же как FIX
+; Отрицательное число: нужно округлить вниз, от нуля
+	CMP	-(SP),-(SP)	; резервируем место на стеке
+	MOV	10(SP),-(SP)	; дублируем исходное число на стек
+	MOV	10(SP),-(SP)
+	MOV	2(SP),-(SP)	; дублируем исходное число ещё раз
+	MOV	2(SP),-(SP)
+; На стеке: исх.число, исх.число, (SP+10)(SP+12) резерв, адрес возврата, исх.число
+	CALL	FFIX		; отбрасываем дробную часть
+	MOV	(SP),10(SP)	; копируем результат FIX в зарезерв. место на стеке
+	MOV	2(SP),12(SP)
+	CALL	FCMP		; сравниваем два floats на стеке, снимает 4 слова
+; На стеке: результат FIX, адрес возврата, исходное число
+	BEQ	1$		; числа равны => нет дробной части
+	CLR	-(SP)		; младшее слово константы 1.0
+	MOV	#40200,-(SP)	; старшее слово константы 1.0
+	CALL	FSUB		; вычитаем 1.0 из результата FIX
+1$:
+; На стеке: результат INT, (SP+4) адрес возврата, (SP+6)(SP+10) исходное число
+	MOV	(SP)+,4(SP)	; копируем результат поверх исходного числа
+	MOV	(SP)+,4(SP)	; и сдвигаем указатель стека
+	RETURN			; возврат
+
 ;#####################################################################
 ;## FINT
 ;## Need FFIX
@@ -712,7 +747,7 @@ VISKAS:	MOVB	#40,(R2)+	; 
 	MOV	R5, (SP)	; адрес возврата
 	JMP	WRSTR1		; выводим полученную строку R0, R2 = длина
 ;
-; DG4DL5 ? умножение 24-битного числа (R1:R0) на 4/5 (0.8)
+; DG4DL5 - умножение 24-битного числа (R1:R0) на 4/5 (0.8)
 ; путём многократного деления на 2 и сложения с исходным значением.
 ; Используется для извлечения десятичных цифр при выводе чисел.
 ; Вход: R0 ? старшее слово, R1 ? младшее слово (24-битное число).
@@ -1035,7 +1070,7 @@ FSUB:
 
 ;#####################################################################
 ;## FMUL
-; Умножение для плавающих чисел, параметры B и A на стеке
+; Умножение для плавающих чисел, параметры B и A на стеке.
 ; На основе кода процедуры MLR из реализации Бейсик Вильнюс для БК-0010.
 ; На входе: (SP) - адрес возврата
 ;           (SP+2)(SP+4) - плавающее число B (старшее, младшее слово)
@@ -1163,10 +1198,140 @@ FMUL$P:	ASR	R4
 
 ;#####################################################################
 ;## FDIV
-; Деление для float, эмуляция FIS, параметры B и A на стеке
+; Деление для плавающих чисел, параметры B и A на стеке.
+; На основе кода процедуры DVR из реализации Бейсик Вильнюс для БК-0010.
+; На входе: (SP) - адрес возврата
+;           (SP+2)(SP+4) - плавающее число B (старшее, младшее слово)
+;           (SP+6)(SP+10) - плавающее число A (старшее, младшее слово)
+; Результат: (SP)(SP+2) - плавающее число (старшее слово, младшее слово)
 FDIV:
-;TODO
-	RETURN ;STUB
+	MOV	(SP)+, FDIV$X+2	; снимаем и подставляем адрес возврата
+	MOV	R5,-(SP)	; сохранение регистров
+	MOV	R4,-(SP)
+; Подготовительные действия
+	CLR	R0		; Подготовить место
+	CLR	R1		; для мантиссы
+	CLR	-(SP)		; и знака результата
+	ASL	12(SP)		; Запомнить в стеке
+	ROL	(SP)		; знак A
+	CLR	-(SP)		; Выделить место для порядка результата
+	TST	10(SP)		; делитель 0 - ошибка
+	BNE	2$
+	TRAP	11.
+2$:	BISB	15(SP),(SP)	; запомнить знак A
+	BEQ	FDIV$0		; A = 0 => результат 0
+	BISB	14(SP),R0
+	SWAB	R0		; заслать в R0-R1
+	SEC			; мантиссу A,
+	ROR	R0		; присоединив спереди
+	BISB	17(SP),R0	; значащий разряд
+	BISB	16(SP),R1
+	SWAB	R1
+	CLR	R2		; подготовить место
+	CLR	R3		; для мантиссы делителя
+	ASL	10(SP)		; сформировать знак
+	ADC	2(SP)		; результата
+	BISB	11(SP),R2	; подсчитать порядок
+	SUB	R2,(SP)		; результата
+	CLR	R2
+	BISB	10(SP),R2
+	SWAB	R2		; заслать в R2-R3
+	SEC			; мантиссу B
+	ROR	R2
+	BISB	13(SP),R2
+	BISB	12(SP),R3
+	SWAB	R3
+; Переход на деление мантисс
+	CMP	R0,R2		; если A > B,
+	BCS	FDIV$Z		; начинать вычитывать B
+	BHI	FDIV$M		; с первого разряда A
+	CMP	R1,R3		; если A < B,
+	BHI	FDIV$M		; начинать вычитывать
+	BNE	FDIV$Z		; cо второго разряда A
+	CLR	14(SP)		; если A = B,
+	INC	(SP)		; подготовить
+	CLR	R5		; формирование результата
+	BR	FDIV$D		; c мантиссой 1
+;
+;REZ0: Подготовка формирования результата 0
+FDIV$0:	CMP	(SP)+,(SP)+
+	BR	FDIV$N
+;
+;EPERP: Обработка переполнения
+FDIV$O:	TRAP	6
+;
+; Обработка исчезновения порядка
+EISNYK:	TST	(SP)+
+;NUL: Формирование результата 0
+FDIV$N:	CLR	10(SP)
+	CLR	12(SP)
+	BR	FDIV$E
+;
+;DAUG: Деление мантисс
+FDIV$M:	ROR	R0		; погашение
+	ROR	R1		; начального сдвига
+	INC	(SP)		; подпрограммы DALYT
+;MAZIAU:
+FDIV$Z:	MOV	#11,R4		; сформирование
+	CALL	FDIV$T		; старшей части
+	MOVB	R5,14(SP)	; мантиссы результата
+	TST	R4		; если мантиссы не сравнялись,
+	BEQ	20$		; продолжать деление
+	CLR	R5		; младшая часть мантиссы
+	BR	FDIV$D		; результата = 0
+20$:	MOV	#20,R4		; продолжение деления
+	CALL	FDIV$T
+;
+;END: Формирование результата
+FDIV$D:	MOV	(SP)+,R4	; формирование порядка
+	ADD	#200,R4
+	BLE	EISNYK		; исчезновение порядка?
+	CMP	#377,R4
+	BLT	FDIV$O		; переполнение?
+	MOVB	R4,13(SP)	; Присоединение порядка к мантиссе результата
+	ROR	(SP)+		; присоединение знака
+	ROR	10(SP)
+	ROR	R5
+	ADC	R5		; округление
+	MOV	R5,12(SP)
+	ADC	10(SP)
+	BCS	FDIV$O		; переполнение?
+	BVS	FDIV$O
+;ENDD:
+;ENDM: Завершение
+FDIV$E:	MOV	(SP)+,R4	; восстановление регистров
+	MOV	(SP)+,R5
+	CMP	(SP)+,(SP)+	; удаление числа B из стека
+FDIV$X:	JMP	@#000000	; возврат
+;
+;DALYT: Подпрограмма деления мантисс
+FDIV$T:	ASL	R5		; сдвиг результата
+	ASL	R1		; и мантиссы A
+	ROL	R0
+	BCS	30$		; C=1 => A > B => вычитать
+	CMP	R2,R0		; сравнить A и B
+	BHI	32$		; B > A => сдвинуть A
+	BCS	30$		; B < A => вычитать B
+	CMP	R3,R1		; сравнить младшие
+	BHI	32$		; части
+	BEQ	36$		; A = B => окончить
+;ATIMT:
+30$:	SUB	R3,R1		; вычитание A-B
+	SBC	R0		; сравнить A и B
+	SUB	R2,R0
+	INC	R5		; Формирование единицы в результате
+;SLINKT:
+32$:	DEC	R4		; результат сформирован?
+	BGT	FDIV$T		; нет - продолжить
+	RETURN			; нет - окончить
+;LYGU:
+36$:	INC	R5		; делимое и делитель
+	BR	41$		; сравнились
+40$:	ASL	R5		; заполнить
+41$:	DEC	R4		; конец результата
+	BGT	40$		; нолями
+	INC	R4		; признак окончания деления
+	RETURN
 
 ;#####################################################################
 ;## STCP