problèmes de précision

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : problèmes de précision

bibi218

Bonjour à tous

Voilà, comme d'habitude, j'ai un problème.
J'ai programmé quelques fonctions, on faisant toujours attention à déclarer mes variables en double, et non en float (j'ai besoin de grands nombres, et surtout d'au moins 10 décimales).
Seulement, il apparait que lors d'une bête addition, on me renvoie un truc du genre 10E20 + 10E20 = -Inf (je suis pas sûr de la puissance, mais bon, ça devrait pas arriver avec des exposants aussi petits en double).
D'où ma question : quand on déclare une variable en double, peut-elle être implicitement castée en float (avec la perte de précision que cela entraîne), ou alors est-ce que le double de la déclaration "l'emporte" ???
Par exemple :

double x,y;
y=(un double quoi);
x=1.0/y; ----> x est double ou float maintenant ???

Publicité

bibi218

Juste une petite chose .... pour les calculs, j'utilise pas mal de #define pour déclarer des constantes, dans un fichier .h
Peut-on spécifier le format de la constante par la même occasion ?

jagstang

Pa Capona ಠ_ಠ

non elle est en double là y'a pas de doute.

---------------
What if I were smiling and running into your arms? Would you see then what I see now?

bibi218

La variable est stockée en double, mais le résultat est-il obtenu par exemple en faisant (double) (1.0/y) ou ((double) 1.0) /y ??? ce qui ne me semble pas tout à fait équivalent

bibi218

Pour etre un peu plus clair, voila ou ca coince ...

Code :

void odeint(double ystart[], unsigned long nvar, double t1, double t2, double eps, double h1, double hmin, int *nok, int *nbad, void (*derivs)(double *,
double *, unsigned long, double, double **, double **, unsigned long *), void (*rkqs)(double *, double *, unsigned long , double *, double , double , double *, double *, double
*, void (*)(double *, double *, unsigned long, double, double **, double **, unsigned long *), double , double **, double
**, unsigned long *), double dx, double **a, double **a1, unsigned long *indx)
{
unsigned long nstp,i,n;
double tsav,t,hnext,hdid,h;
double *yscal,*y,*dydt;
n=nvar/4;
yscal=dvector(1,nvar);
y=dvector(1,nvar);
dydt=dvector(1,nvar);
t=t1;
h=SIGN(h1,t2-t1);
*nok = (*nbad) = kount = 0;
for (i=1;i<=nvar;i++) y[i]=ystart[i];
if (kmax > 0) tsav=t-dtsav*2.0; // stockage du premier pas
for (nstp=1;nstp<=MAXSTP;nstp++) {
(*derivs)(y,dydt,n,dx,a,a1,indx); ---> ca foire ici, la fonction marche quand je l'utilise a part, mais ici, j'ai du Nan qui apparait, avec les memes arguments en entree
for (i=1;i<=nvar;i++) yscal[i]=fabs(y[i])+fabs(dydt[i]*h)+TINY2;
if (kmax > 0 && kount < kmax-1 && fabs(t-tsav) > fabs(dtsav)) {
tp[++kount]=t; // stockage des valeurs intermediaires
for (i=1;i<=nvar;i++) yp[i][kount]=y[i]; tsav=t;
}
if ((t+h-t2)*(t+h-t1) > 0.0) h=t2-t; // si l'on sort du domaine, on reduit le pas
(*rkqs)(y,dydt,n,&t,h,eps,yscal,&hdid,&hnext,derivs,dx,a,a1,indx);
if (hdid == h) ++(*nok);
else ++(*nbad);
if ((t-t2)*(t2-t1) >= 0.0) {
for (i=1;i<=nvar;i++) ystart[i]=y[i];
if (kmax) {
tp[++kount]=t; // sauvegarde des dernieres valeurs
for (i=1;i<=nvar;i++) yp[i][kount]=y[i];
}
free_dvector(dydt,1,nvar);
free_dvector(y,1,nvar);
free_dvector(yscal,1,nvar);
return;
}
if (fabs(hnext) <= hmin) nrerror("Step size too small in odeint" );
h=hnext;
}
nrerror("Too many steps in routine odeint" );
}

bibi218

desole pour la declaration de fonction, c'est vraiment bourrin, mais j'ai rien trouve de mieux que ca .....

Taz

bisounours-codeur

les typedef abruti

FORUM HardWare.fr

Programmation

problèmes de précision

Sujets relatifs
[IE & Mozilla] A peine commencé déja des problèmes	problèmes de Threads .
Shell sh/tcsh : problèmes de caractères spéciaux	Problemes d erreurs en ADA pour un devoir à rendre Help
Problemes d erreurs en ADA pour un devoir à rendre Help	Problèmes de cookies avec IE mais pas Firefox
[CSS] Mac : Problèmes divers - PART II : IE Mac et les marges	[PHP] Une petite précision sur date()
Problèmes de création de noeud avec JavaCC	[Visual C++]Problèmes d'affichage d'images
Plus de sujets relatifs à : problèmes de précision

Page générée en 0.056 secondes