Аттракцион повышенной точности
Когда в 1967 году советским вычислительным математикам досталась очень большая (аж до 192 Кб адресуемой памяти, на более ранних, которые они могли видеть, было меньше 50 Кб) очень быстрая (аж до миллиона операций в секунду, на более ранних - не больше 50 тысяч) вычислительная машина, они тут же захотели сделать себе библиотеку программ, полезных для повседневной работы.
Программой №1 была "прокрутка" - выполнение программы по одной команде и печать адреса и кода каждой команды и значений используемых в ней регистров на АЦПУ. При программировании в машинных кодах - незаменимая вещь.
Программами №2 и №3 были "автокод" (ассемблер) и "программа стыковки". Тоже дело полезное для повседневной работы.
А вот программой №4 была "печать с форматом". Программа печатает в формате, утвержденном нормоконтролем, программы с адресами команд и названиями кодов операций. Кстати о приоритетах.
Дальше шли:
А сразу за ними в наборе была стандартная программа "увеличение точности" (хотя, казалось бы, 40 бит мантиссы должно быть достаточно для большинства применений). Насколько я понимаю, если у кого результаты вычислений оказывались бессмысленными из-за потери точности при вычислениях с плавающей точкой,
то рекомендовалось не сразу заморачиваться выяснениями, как исправить алгоритм, чтобы ликвидировать потерю точности, а сначала выяснить, становится ли хоть сколь-нибудь лучше при увеличении точности.
Но как? Тут тебе не Фортран, где можно объявить переменные DOUBLE PRECISION вместо REAL, и не C++, родной.
А вот как: программа увеличения точности интерпретировала машинный код, выполняя все команды, кроме операций с плавающей точкой, как есть, а всю плавающую точку - с заданной увеличенной точностью. Нужно ли было при этом программисту изменять работу с адресами переменных и индексами массивов, или программа сама выбирала в буферной памяти, где хранить "хвосты" переменных, требующих большей точности, ещё предстоит разобраться.
Работало это, понятное дело, гораздо медленнее, чем могло бы при использовании библиотеки повышенной точности с отдельными вызовами, но, по-видимому, народ устраивало.
Программой №1 была "прокрутка" - выполнение программы по одной команде и печать адреса и кода каждой команды и значений используемых в ней регистров на АЦПУ. При программировании в машинных кодах - незаменимая вещь.
Программами №2 и №3 были "автокод" (ассемблер) и "программа стыковки". Тоже дело полезное для повседневной работы.
А вот программой №4 была "печать с форматом". Программа печатает в формате, утвержденном нормоконтролем, программы с адресами команд и названиями кодов операций. Кстати о приоритетах.
Дальше шли:
- перемножение матриц
- решение системы линейных алгебраических уравнений
- вычисление определенных интегралов по формуле Симпсона
- "упорядочение массива" (то ли слово "сортировка" тогда ещё не придумали, то ли конкретно "упорядочение" значило сохранение исходного массива нетронутым, а запись результата в другой массив - каким алгоритмом это делалось, ещё предстоит разобраться; чует моё сердце, что пузырьком - машина-то быстрая, а массивы маленькие)
- ассоциативный (по маске и значению) поиск слова в массиве
- поиск максимального числа в массиве (поиска минимального не нажили)
- решение системы дифференциальных уравнений 1-го порядка методом Рунге-Кутта
- слияние двух массивов ([a, b, c, ...], [d, e, f, ...] -> [a, d, b, e, c, f, ...] - зачем бы это так уж часто было нужно?)
- обращение матрицы
- квадратичная интерполяция
- буферизация ввода-вывода на внешние носители
- пакетный отладчик с весьма развитым скриптовым языком
А сразу за ними в наборе была стандартная программа "увеличение точности" (хотя, казалось бы, 40 бит мантиссы должно быть достаточно для большинства применений). Насколько я понимаю, если у кого результаты вычислений оказывались бессмысленными из-за потери точности при вычислениях с плавающей точкой,
то рекомендовалось не сразу заморачиваться выяснениями, как исправить алгоритм, чтобы ликвидировать потерю точности, а сначала выяснить, становится ли хоть сколь-нибудь лучше при увеличении точности.
Но как? Тут тебе не Фортран, где можно объявить переменные DOUBLE PRECISION вместо REAL, и не C++, родной.
А вот как: программа увеличения точности интерпретировала машинный код, выполняя все команды, кроме операций с плавающей точкой, как есть, а всю плавающую точку - с заданной увеличенной точностью. Нужно ли было при этом программисту изменять работу с адресами переменных и индексами массивов, или программа сама выбирала в буферной памяти, где хранить "хвосты" переменных, требующих большей точности, ещё предстоит разобраться.
Работало это, понятное дело, гораздо медленнее, чем могло бы при использовании библиотеки повышенной точности с отдельными вызовами, но, по-видимому, народ устраивало.