Печать страницы - Функция от одного аргумента возвращает разные значения ?

Название: Функция от одного аргумента возвращает разные значения ?
Отправлено: Белый пони от Март 05, 2011, 10:56

Здравствуйте!
Совсем не могу понять как так.

вот функция fX(r):

Код:

double MainWidget::fX( double r)
{
double f = sqrt( r*r + l13*l13 - 2 * r * l13 * cos( beta23 - falpha(r)) ) - r*sqI13;
return f;
}

double MainWidget::falpha( double r)
{
double a = acos ( (r*r + l12*l12 - r*r*I1/I2) / (2 * r * l12) );
if( isnan(a)) a = 0;
return a;
}

При этом такой кусок кода:

Код:

void MainWidget::toFile()
{
double hi;
...

std::cerr << " hi = " << hi << std::endl;
std::cerr << "fX ( "<< hi << " ) = " << fX( hi ) << std::endl;
hi = 233.623;
std::cerr << "fX ( "<< hi << " ) = " << fX( hi ) << std::endl;
}

выдаёт такой результат:

Код:

hi = 233.623
fX ( 233.623 ) = 34.367
fX ( 233.623 ) = -71.3663

Как так получается ??? Я не понимаю каких-нибудь основ С++ ?

(остальные переменные - l12, l13, beta23, I1, I2, sqI13 - члены класса MainWidget, вычисляемые один раз в конструкторе)

Название: Re: Функция от одного аргумента возвращает разные значения ?
Отправлено: SimpleSunny от Март 05, 2011, 11:13

Зайдите под отладчиком в эти функции и посмотрите что не так (:

К тому же аргументы могут быть не совсем одинаковыми, при выводе могли потеряться разряды у double.

Название: Re: Функция от одного аргумента возвращает разные значения ?
Отправлено: Белый пони от Март 05, 2011, 14:03

Цитата: SimpleSunny от Март 05, 2011, 11:13

К тому же аргументы могут быть не совсем одинаковыми, при выводе могли потеряться разряды у double.

Да, так и было :)
Хотя всего четвёртый знак после запятой, но у арккосинуса по краям бесконечная производная, и итоговая ошибка получалась большой.

Ошибка была в том, что выражение

Код:

if( isnan(a)) a = 0;

учитывает только случай аргумент > +1.

отдельно расписал условия для -1 и +1 и всё заработало.

Спасибо :)

Russian Qt Forum

Программирование => С/C++ => Тема начата: Белый пони от Март 05, 2011, 10:56