Ошибка конверсии (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.

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

[zltigo]

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

Нет.
400ULL

[scifi]

400ULL

Ага, щас. Вообще-то 400 перед умножением автоматически приводится к double. Так что ULL ничего не изменит.
Скорее уж 400.0, чтобы убрать неявное приведение типа.

[zltigo]

Ага, щас. Вообще-то 400 перед умножением автоматически приводится к double. Так что ULL ничего не изменит.
Скорее уж 400.0, чтобы убрать неявное приведение типа.

Не обратил внимание, что там double. Да, 400.0

[scifi]

Не обратил внимание, что там double. Да, 400.0

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