Где ошибка в моем коде полиномиального вычисления?

Question

Моя цель состоит в том, чтобы сравнить два разных метода для вычисления полинома ... но если вы запустите его на Mac (мой компьютер - MacBook Air), вы обнаружите, что два результата различны ... .но ....если вы удаляете часть "/ * ... * /" ИЛИ удаляете два "//" перед "для ...", он работает нормально ...

плюс ... он отлично работает на Linux...

Может кто-нибудь сказать мне, почему? ..

вот моя программа:

#define MAX_NUM 10
#define TEST_TIME 1
#define test_type double

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

test_type normal_method(test_type *a, test_type x, int degree)
{
    test_type s = 0.0;
    test_type e = 1.0;

    for (int i = 0; i <= degree; i++) {
        s += a[i] * e;
        e *= x;
    }

    printf("%lf\n", s);

    return s;
}

test_type horner_method(test_type *a, test_type x, int degree)
{
    test_type s = a[degree];

    for (int i = degree - 1; i >= 0; i--) {
        s *= x;
        s += a[i];
    }

    printf("%lf\n", s);

    return s;
}

void generate_data(test_type *a, test_type *x, int degree)
{
    srand( time(NULL) );
    for (int i = 0; i <= degree; i++) {
        a[i] = (test_type)(rand() % MAX_NUM + 1) / (test_type)(rand() % MAX_NUM + 1);
        *x = (test_type)(rand() % MAX_NUM + 1) / (test_type)(rand() % MAX_NUM + 1);
    }

}

int main()
{
    const int degree = 10;
    test_type a[degree];
    test_type x;

    generate_data(a, &x, degree);

//Just by clear the /**/ below, the program will work fine....
/*    
    printf("x = %lf\n", x);
    for (int i = 0; i <= degree; i++) {
        printf("%lf\n", a[i]);
    }
*/
    clock_t begin, end;

//  Or clear the two // below, it will work fine too....

    begin = clock();
//    for (int i = 0; i < TEST_TIME; i++) 
        normal_method(a, x, degree);
    end = clock();
    printf("The normal method used %d clock times\n", end - begin);

    begin = clock();
//    for (int i = 0; i < TEST_TIME; i++)
        horner_method(a, x, degree);
    end = clock();
    printf("The horner method used %d clock times\n", end - begin);

    return 0;
}

Answer 1

Вы обращаетесь к памяти, которая находится за пределами массива, который вы создаете в main, и используете его в других функциях. Это происходит как минимум в трех местах.

• Границы цикла в normal_method переходят от нуля к десяти:

  for (int i = 0; i <= degree; i++) {  // a[degree] is out of bounds
  

• В первой строке horner_method вы обращаетесь к памяти, которая находится за пределами вашего массива:

  test_type s = a[degree];  // a[degree] is out of bounds
  

• Границы цикла в generate_data неверны так же, как и в normal_method:

  for (int i = 0; i <= degree; i++) {  // a[degree] is out of bounds
  

Answer 2

Не следует использовать %lf для печати double. Просто %f в порядке. Вы, вероятно, перепутали это с scanf, который нуждается в l.

С опциями предупреждения компилятора -Wall -Wextra -Wformat=2 gcc должен сообщить вам, что не так в вашем коде.