Как рассчитать напоминание о делении в сборке SPARC?

Вот псевдокод, который вычисляет деление двух натуральных чисел.
HR-регистр сохраняет напоминание, а LR сохраняет дивиденды. (и в итоге сохраняет root)
Однако я думаю, что этот алгоритм имеет некоторые проблемы.
Поскольку этот алгоритм иногда не восстанавливает вычитание (деление является продолжением вычитания).
Например 6 / 3 (0110 / 011)
Этот алгоритм вычитает -3 еще раз. (Такая ситуация никогда не возникает, когда мы вычисляем это деление вручную)
Поэтому я думаю, что этот алгоритм имеет некоторые проблемы.
Ты не согласен со мной? Как рассчитать деление напоминания в сборке?

for i = 1 to num_of_bits do
(HR LR) << 1
if (HR >= 0) then
   HR = HR - DIVISOR
else
   HR = HR + DIVISOR
endif
if (HR > 0) then LR(lsb) = 1 endif
endfor
Источник

2 ответа

Решение

Я не говорю на языке SPARC, но я говорю на C. Вот пример реализации алгоритма для деления 16/8=8,8:

#include <stdio.h>

typedef unsigned char uint8;
typedef unsigned int uint;

int u8div(uint8* dividendh, uint8* dividendl, uint8 divisor)
{
  int i;

  if (*dividendh >= divisor)
    return 0; // overflow

  for (i = 0; i < 8; i++)
  {
    if (*dividendh >= 0x80)
    {
      *dividendh = (*dividendh << 1) | (*dividendl >> (8 - 1));
      *dividendl <<= 1;

      *dividendh -= divisor;
      *dividendl |= 1;
    }
    else
    {
      *dividendh = (*dividendh << 1) | (*dividendl >> (8 - 1));
      *dividendl <<= 1;

      if (*dividendh >= divisor)
      {
        *dividendh -= divisor;
        *dividendl |= 1;
      }
    }
  }

  return 1;
}

int u8div2(uint8* dividendh, uint8* dividendl, uint8 divisor)
{
  uint dividend = (*dividendh << 8) | *dividendl;

  if (*dividendh >= divisor)
    return 0; // overflow

  *dividendl = dividend / divisor;
  *dividendh = dividend % divisor;

  return 1;
}

int main(void)
{
  uint dividendh, dividendl, divisor;

  for (dividendh = 0; dividendh <= 0xFF; dividendh++)
    for (dividendl = 0; dividendl <= 0xFF; dividendl++)
      for (divisor = 0; divisor <= 0xFF; divisor++)
      {
        uint8 divh = dividendh, divl = dividendl, divr = divisor;
        uint8 divh2 = dividendh, divl2 = dividendl;

        printf("0x%04X/0x%02X=", (divh << 8) | divl, divr);

        if (u8div(&divh, &divl, divr))
          printf("0x%02X.0x%02X", divl, divh);
        else
          printf("ovf");

        printf(" ");

        if (u8div2(&divh2, &divl2, divr))
          printf("0x%02X.0x%02X", divl2, divh2);
        else
          printf("ovf");

        if ((divl != divl2) || (divh != divh2))
          printf(" err"); // "err" will be printed if u8div() computes incorrect result

        printf("\n");
      }

  return 0;
}
Источник

Несколько реализаций алгоритма деления (который также вычисляет остаток) можно найти в приложении E руководства по архитектуре SPARC.

Более новые версии архитектуры SPARC включают в себя операторы разделения UDIV и SDIV.

Дальнейшее внедрение можно найти здесь.

Источник
Другие вопросы по тегам