Когда это правда?

Да, я думаю, что я действительно живу мечтой. У меня есть следующий фрагмент кода, который я компилирую и запускаю на машине AIX:

AIX 3 5
PowerPC_POWER5 processor type
IBM XL C/C++ for AIX, V10.1
Version: 10.01.0000.0003


#include <stdio.h>
#include <math.h>

#define RADIAN(x) ((x) * acos(0.0) / 90.0)

double nearest_distance(double radius,double lon1, double lat1, double lon2, double lat2){
    double rlat1=RADIAN(lat1);
    double rlat2=RADIAN(lat2);
    double rlon1=lon1;
    double rlon2=lon2;
    double a=0,b=0,c=0;

    a = sin(rlat1)*sin(rlat2)+ cos(rlat1)*cos(rlat2)*cos(rlon2-rlon1);
    printf("%lf\n",a);
    if (a > 1) {
      printf("aaaaaaaaaaaaaaaa\n");
    }
    b = acos(a);
    c = radius * b;

    return radius*(acos(sin(rlat1)*sin(rlat2)+
        cos(rlat1)*cos(rlat2)*cos(rlon2-rlon1)));

}

int main(int argc, char** argv) {
  nearest_distance(6367.47,10,64,10,64);
  return 0;
}

Теперь значение «а» после расчета сообщается как «1». И на этой машине AIX, похоже, 1> 1 истинно, так как введено мое «если» !!! И мой акос о том, что я считаю «1», возвращает NanQ, поскольку 1 больше 1. Могу я спросить, как это вообще возможно? Я не знаю, что думать больше!

Код прекрасно работает на других архитектурах, где «a» действительно принимает значение того, что я считаю равным 1, а acos (a) - 0.

Ответы на вопрос(4)

Ваш ответ на вопрос