Как рассчитать остаток от деления в сборке SPARC? - PullRequest
2 голосов
/ 10 октября 2011

Вот псевдокод, который вычисляет деление двух натуральных чисел.
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 ]

3 голосов
/ 10 октября 2011

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

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

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

1 голос
/ 10 октября 2011

Я не говорю на языке 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;
}
...