Введение: Легендарный сопропроцессор и его математическое наследие

Когда компания Intel в 1980 году выпустила математический сопропроцессор 8087, это стало настоящей революцией для мира персональных компьютеров. Впервые x86-совместимые системы получили аппаратную поддержку чисел с плавающей запятой, что освободило центральный процессор от тяжелейших вычислений. Архитектура x87 задала стандарт для целой эпохи, определив формат IEEE 754 и набор тригонометрических инструкций, таких как FPTAN, FSIN и FCOS.

Среди всех математических функций вычисление тангенса (FPTAN) всегда стояло особняком. Инженеры Intel столкнулись с жесткими ограничениями по площади кремния, энергопотреблению и тактовой частоте. Возникает закономерный вопрос: как инженеры 70-х годов заставили кремниевый чип вычислять сложные трансцендентные функции с высокой точностью? Долгое время в сообществе бытовало мнение, что процессоры Intel использовали стандартный алгоритм CORDIC (Coordinate Rotation Digital Computer). Однако детальный анализ и реверс-инжиниринг микрокода 8087 показывают гораздо более сложную, элегантную и гибридную картину.

Миф о чистом CORDIC: Почему традиционный подход не подошел Intel

Алгоритм CORDIC, разработанный Джеком Вотом в 1959 году, кажется идеальным кандидатом для аппаратной реализации тригонометрии. Он требует только операций сдвига и сложения/вычитания (shift-add), избегая дорогих аппаратных умножителей. Однако у классического CORDIC есть существенные недостатки, которые сделали его непригодным для архитектуры Intel 8087 в чистом виде:

  • Ограниченный диапазон сходимости: Базовый CORDIC сходится только для углов в диапазоне примерно от -99 до +99 градусов. Для обработки полных периодов требуется предварительное приведение аргумента.
  • Низкая скорость сходимости на бит: Для достижения 64-битной точности (формат Extended Precision в x87) требуется больше шагов итерации, что увеличивает задержку (latency) выполнения инструкции.
  • Накладные расходы на масштабирование: После серии вращений CORDIC масштабирует результат на фиксированный коэффициент $K$, который требует дополнительной коррекции.

Интелу требовалось нечто гораздо более мощное и универсальное, способное обрабатывать числа двойной точности и расширенного формата (80 бит) с минимальной погрешностью (ulp — unit in the last place). Именно поэтому реальный алгоритм тангенса в Intel 8087 представляет собой сложный гибрид табличных методов, полиномиальных аппроксимаций и специализированных итеративных шагов.

Архитектура микрокода: Как устроен конвейер FPTAN

Инструкция FPTAN в Intel 8087 выполняет сразу две задачи: она вычисляет и тангенс аргумента $X$, и помещает в стек сопропроцессора число 1.0 (в качестве коэффициента для тригонометрического тождества $\tan(x) = y/x$), если операция прошла успешно. Весь процесс делится на несколько последовательных микроархитектурных фаз:

  1. Анализ и редукция диапазона (Range Reduction): Проверка входного значения на специальные случаи (NaN, бесконечность, нуль). Уменьшение аргумента по модулю $\pi/2$ с использованием констант высокой точности, хранящихся в ПЗУ микрокода.
  2. Выбор ветвления: В зависимости от попадания редуцированного аргумента в определенные поддиапазоны, процессор выбирает оптимальный математический метод — будь то усеченный ряд Тейлора, рациональная аппроксимация (метод Паде) или модифицированные итерации.
  3. Финальное масштабирование и округление: Приведение результата к требуемому формату IEEE 754 с учетом текущего режима округления (rounding mode), заданного в управляющем слове (Control Word).
    4. ; Концептуальная схема работы FPTAN на уровне микрокода
1. ST0 = Input_Argument
2. Call Check_Special_Cases(ST0)
3. ST0, flags = Range_Reduce_Pi_Over_2(ST0)
4. If |ST0| < Threshold:
     Jump Polynomial_Approximation
   Else:
     Jump CORDIC_Like_Iteration
5. Apply_Rounding_And_Flags()

    Гибридный математический аппарат: Полиномы и модифицированный CORDIC

    Вместо слепого следования одной математической модели, микрокод 8087 применяет адаптивный подход. Для малых углов разложение в ряд или рациональные дроби работают значительно быстрее и точнее, чем итеративные методы. Лишь для определенных интервалов задействуются алгоритмы, родственные CORDIC, но оптимизированные под аппаратную базу сопропроцессора.

    Инженеры Intel использовали ПЗУ микрокода не просто как хранилище инструкций, а как многомерный массив предвычисленных коэффициентов, которые минимизировали количество арифметических тактов процессора.

    Рассмотрим, как выглядит этот баланс на практике. Аппаратный умножитель 8087 работал по принципу побитового сдвига и сложения, что делало умножение относительно быстрым (хотя и медленнее сложения). Поэтому комбинация нескольких умножений с заранее зашитыми константами часто оказывалась эффективнее длинных циклов CORDIC.

    Для вычисления $\tan(x)$ при малых $x$ используется аппроксимация полиномами:

    tan(x) ≈ x + c3*x^3 + c5*x^5 + c7*x^7

    Где константы $c_3, c_5, c_7$ зашиты в ПЗУ с точностью, превышающей 80 бит, чтобы избежать накопления ошибок округления на промежуточных этапах.

    Борьба за точность: Обработка исключительных ситуаций и аргументов

    Одним из самых сложных аспектов реверс-инжиниринга 8087 является понимание того, как процессор обрабатывает краевые случаи. Математические функции с плавающей запятой уязвимы к потере точности при редукции аргумента (catastrophic cancellation), особенно когда входной угол $X$ представляет собой огромное число, близкое к целому кратному $\pi/2$.

    В Intel 8087 эта проблема решалась с помощью расширенного внутреннего регистра шириной более 80 бит для хранения константы $\pi/2$. Когда пользователь передавал аргумент, превышающий по модулю $2^{63}$, стандартные методы редукции давали сбой. Микрокод сопроцессора включал специальные подпрограммы точного деления по модулю, предотвращающие выдачу неверных результатов для периодических функций.

    Если аргумент превышал допустимый диапазон (около $10^{13}$), флаг C1 в статусной строке (Status Word) устанавливался в единицу, а генерация исключения зависела от масок в управляющем регистре. Это демонстрирует, что алгоритм тангенса в 8087 — это не просто формула, а целая система управления состоянием вычислений.

    Заключение: Наследие архитектуры x87 в современных процессорах

    Исследование и реверс-инжиниринг алгоритма вычисления тангенса в Intel 8087 показывают, насколько глубоко продуманной была эта кремниевая архитектура. Вопреки мифам о чистом CORDIC, инженеры Intel создали гибкий гибридный движок, объединивший достоинства табличных методов, быстрых полиномиальных разложений и итеративных уточнений.

    Удивительно, но базовые принципы, заложенные в микрокод Intel 8087 в конце 1970-х годов, прошли сквозь десятилетия. Современные процессоры архитектур x86 и x86-64 продолжают поддерживать инструкции x87 (включая FPTAN) для обратной совместимости, хотя внутренне они часто транслируют их в микрооперации или обрабатывают с помощью векторных модулей SSE/AVX. Понимание того, как инженеры прошлого справлялись с жесткими аппаратными ограничениями ради достижения максимальной математической точности, остается бесценным опытом для современных разработчиков компиляторов и системных архитекторов.