Форум программистов, компьютерный форум CyberForum.ru

Floating point overflow - C++

Восстановить пароль Регистрация
 
Alligieri
CEO
Эксперт C++
 Аватар для Alligieri
2247 / 1236 / 24
Регистрация: 16.03.2009
Сообщений: 3,581
18.05.2009, 21:56     Floating point overflow #1
в последней строке кода функции double uzag_coef(double inp_1,double inp_2) возникает ошибка переполнения, при чем 1 раз скомпилилось нормально, а остальные - появилась ошибка. Если есть у кого идеи или даже решение - помогите

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
#include <stdio.h>
#include <conio.h>
#include <math.h>
 
double uzag_coef(double inp_1,double inp_2)
{
int index_c;
double result;
if (inp_2==0) return 1;
if (inp_2<1) { result=result*(inp_1-inp_2+1); return result;}
for (index_c=1;index_c<inp_2;index_c++) result=result*(inp_1-index_c);
result*=(inp_1-inp_2+1);
return result;
}
 
int main()
{
static double x[257];
static double Y[257];
static double poly[257];
static double an[257];
static double J[257];
int index;
int x_counter;
 
int N=256;
double h;
 
double p=0.5;
double q=0.5;
 
FILE* outFile;
 
x[0]=-1;
x[256]=2;
h=fabs((x[256])-(x[0]))/(double) N;
 
for (index=1;index<256;index++)
x[index]=x[index-1]+h;
 
index=0;
x_counter=0;
 
while (x_counter!=1)
{Y[index]=1;
poly[index]=(1/6)*(x[index]*x[index]*x[index]-3*p*(q+(double)N)*x[index]*x[index]+(3*p*p*((double)(N-1))*((double)(N-2))+3*p*((double)(N-2))+2)*x[index]-p*p*p*((double)N)*((double)(N-1))*((double)(N-2)));
J[index]=uzag_coef((double)N,x[index])*pow(p,x[index])*pow(q,((double)N)-x[index]);
an[index]=pow((uzag_coef((double)N,(double)index)*pow(p*q,index)),-1)*J[index]*Y[index]*poly[index];
index++;
if (x[index]>0) x_counter=1;
}
 
x_counter=0;
 
while (x_counter!=1)
{Y[index]=-(x[index]*x[index])+1;
poly[index]=(1/6)*(x[index]*x[index]*x[index]-3*p*(q+(double)N)*x[index]*x[index]+(3*p*p*((double)(N-1))*((double)(N-2))+3*p*((double)(N-2))+2)*x[index]-p*p*p*((double)N)*((double)(N-1))*((double)(N-2)));
J[index]=uzag_coef((double)N,x[index])*pow(p,x[index])*pow(q,((double)N)-x[index]);
an[index]=pow((uzag_coef((double)N,(double)index)*pow(p*q,index)),-1)*J[index]*Y[index]*poly[index];
index++;
if (x[index]>1) x_counter=1;
}
 
x_counter=0;
 
while (x_counter!=1)
{Y[index]=sqrt(x[index]-1);
poly[index]=(1/6)*(x[index]*x[index]*x[index]-3*p*(q+(double)N)*x[index]*x[index]+(3*p*p*((double)(N-1))*((double)(N-2))+3*p*((double)(N-2))+2)*x[index]-p*p*p*((double)N)*((double)(N-1))*((double)(N-2)));
J[index]=uzag_coef((double)N,x[index])*pow(p,x[index])*pow(q,((double)N)-x[index]);
an[index]=pow((uzag_coef((double)N,(double)index)*pow(p*q,index)),-1)*J[index]*Y[index]*poly[index];
index++;
if (x[index]>=2) x_counter=1;
}
 
outFile=fopen("full_res.out","w");
for (index=0;index<257;index++)
fprintf(outFile,"%f \n",an[index]);
 
fclose(outFile);
 
printf("All OK!");
 
getch();
return 0;
}
Добавлено через 6 минут 22 секунды
Сам же и нашел - не определил переменную result и поэтому программа использовала мусор который находился в ячейке памяти.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
18.05.2009, 21:56     Floating point overflow
Посмотрите здесь:

C++ Floating point error:Overfull
C++ Floating point overflow
Illegal use of floating point C++
C++ floating point stack check at 0x00402f75
floating point error C++
C++ Ошибка: Illegal use of floating point
C++ Illegal use of floating point
Floating point overflow C++
Много ошибок illegal use of floating point C++
Exception class overflow with message floating point overflow C++
C++ Floating point exception (core dumped)
Ошибка Illegal use of floating point C++

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 01:41. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru