Ошибка конверсии (unsigned int) double ? Опции

[ZSN1]

Вычисляю sin табличным методом
(LPC2478 в среде Keil 4.23)

кусок кода

Код

....
#define Tabl_acos 0x80000000
...
...
double tsin (double arg)
{
double argout;
unsigned int *adr;

    adr = (unsigned int *)((unsigned int)Tabl_sin + (unsigned int)(fabs(arg)*400));
    argout = *(double *)adr;
    return (argout);
}

Указатель adr при вычислении имеет совершенно фантастическое значение !

А в коде ниже (unsigned int)fabs(arg)*400 работает правильно, но теряется точность, поскольку приведение к uint происходит до умножения на 400.

Код

....
#define Tabl_acos 0x80000000
...
...
double tsin (double arg)
{
double argout;
unsigned int *adr;

    adr = (unsigned int *)((unsigned int)Tabl_sin + (unsigned int)fabs(arg)*400);
    argout = *(double *)adr;
    return (argout);
}

Пробовал разные ключи компилятора - не помогает
Заметил, что ошибка возникает, когда значение arg*400 превышает 255.

Глюк компилятора?

[ZSN1]

Сделал втупую

Код

  value = arccos (0);     // result = 90
  value = arccos (0.465); // result = 62.289
  value = arccos (0.5);   // result = 60
  value = arccos (0.563); // result = -4967.265 (incorrect!)
  value = arccos (0.564); // result = -4967.265 (incorrect!)
  value = arccos (0.565); // result = -4967.265 (incorrect!)
  value = arccos (0.566); // result = 0 (incorrect!)
  value = arccos (0.567); // result = -4967.265 (incorrect!)
  value = arccos (0.568); // result = -4967.265 (incorrect!)
  value = arccos (0.569); // result = 0 (incorrect!)
  value = arccos (0.57);  // result = 0 (incorrect!)
  value = arccos (0.58);  // result = 54.549
  value = arccos (0.59);  // result = 53.842
  value = arccos (1);     // result = 0

результат поставил вообще в тупик ....

таблицу проверял 100 раз
выводил на экран значения от 0.5 до 0.6
Возможно в других диапазонах тоже самое

Только это всё равно ничего не изменит.
У меня есть версия. Возможно, стек не выравнен по границе 8 байт, и из-за этого не работает умножение double*double.

А как стек выровнять? Добавить __align(32)?

Сообщение отредактировал ZSN1 - Aug 20 2015, 10:51

[Obam]

А как стек выровнять? Добавить __align(32)?

Чёй-то 32? С какой радости?