evidemment c le contraire que je voulais ecrire :D
qd tu passe de double à float les naN sont plus proche. Ce qui parait normal dans l'hypothese de la precision qui fout la merde!
j'edit cette phrase desuite.
avoue que sinon t'aurai rien eu a dire ! et puis ca a fait un petit moment de suspens ! Bref j'essaye de justifier le fait que je me suis trompé ;)
n'as-t'on pas du tout de NaN ? ou alors se produit'il trop loin pour qu'on l'observe ?
Moi pas comprendre.
Ce que je peux dire c'est
1) J'ai lancé une myriade fois le code avec matrice aléatoire avec Math.random() mais en forçant, la non nullité de coordonnées et je n'ai pas vu seul Nan
2) j'ai fait pareil avec Math.random() * 0.0001 et je n'ai pas vu un NaN
3) j'ai fait pareil avec Math.random() * 0.0000001 et j'ai vu des Nan
4) j'ai donc raison depuis longtemps: le code et les fonctions mathématique sont convenables avec des coordonnées non nulles. Mais au voisinage des zéros la division par zéro est problèmatique ce qui est quand même normal, Java ou pas.
5) je voudrais bien savoir où facturer pour le million d'Euro.
par contre j'ai un code tout pret avec une matrice tout jolie qui donne du naN à 26. Le code est adapté pour lire les fichiers txt ou il y a la matrice d'ecrite.
Donc y'a le code Geotest, et la matrice qu'il faut mettre dans un fichier test_matrix_30.txt, et puis creer une methode qui lit cette matrice (moi je l'ai dcette methode, mais secret defense je la garde, et en plus ca marche pas chez certain)
desolé pour le post immense, mais j'ai pas reussi à mettre sur mon serveur web, il plante ce con, pil poil quand j'en ai besoin !
Le code de Geotest.java
import java.io.File;
import java.util.ArrayList;
class Geotest {
final int NB_POINTS = 30;
// type Matrice de départ : 0 = c'est nous qui la créons ; 1 = aléatoire ; 2=on la lit dans un fichier
final int typeMat = 2;
/* typeAffichage :
* 0 = on affiche toutes les étapes de calcul
* 1 = on affiche juste résultats
* 2 = on affiche juste quand y'a du NaN */
public int typeAffichage = 2;
// n = nombre de points
// m = dimension de l'espace
private int n = NB_POINTS;
private int m = NB_POINTS - 1;
// Matrice des positions initiales
public double[][] coordonneesInitiales = new double[n][m];;
// Matrice des positions reconstruites par l'algorithme
public double[][] coordonnees = new double[n][m];
// Matrice des distances
public double[][] distances = new double[n][n];
// vecteur pour controle de division par zero
// zeroPos[0]=position en i ; zeroPos[1]=position en j ; zeroPos[2] = valeur de coord
public double[] zeroPos = { -1,-1,-1};
// petite methode pour calculer un carre
double carre(double x) {
return x * x;
}
// Calcul de la distance euclidienne pour la matrice de position
double dist(double[][] coords, int a, int b) {
double sum = 0;
for (int i = 0; i < m; i++) {
sum = sum + carre(coords[a][i] - coords[b][i]);
}
return Math.sqrt(sum);
}
public void creationMatrice(int type){
if(typeAffichage==0) System.out.println("Matrice Initiale");
switch(type) {
case 0 :
// ON cree nous même la matrice
double[][] temp =
{
{ 0.0, 0.0, 0.0 } ,
{ 2.0, 0.0, 0.0 },
{ 7.0, 0.0, 0.0 },
{ 0.0, 1.0, 2.0 }
};
for(int i=0;i<n;i++){
if(typeAffichage==0) System.out.print(i + "=> ");
for(int j=0;j<m;j++){
coordonneesInitiales[i][j]= temp[i][j];
if(typeAffichage==0) System.out.print(coordonneesInitiales[i][j] + " ");
}
if(typeAffichage==0) System.out.println();
}
break;
case 1 :
// On crée une matrice triangulaire aléatoire
for (int i = 0; i < n; i++) {
if(typeAffichage==0) System.out.print(i + "=> ");
for (int j = 0; j < i; j++) {
coordonneesInitiales[i][j] = (Math.random()*10);
if(typeAffichage==0) System.out.print(coordonneesInitiales[i][j] + " ");
}
if(typeAffichage==0) System.out.println();
}
break;
case 2 :
coordonneesInitiales = readMatrixFromFile(new File("test_matrix_30.txt"));
}
}
public void executeTest() {
creationMatrice(typeMat);
// On calcul les distances euclidiennes entre les points
if(typeAffichage==0) System.out.println("Matrice distance");
for (int i = 0; i < n; i++) {
//System.out.print(i + "=> ");
for (int j = 0; j < m; j++) {
distances[i][j] = dist(coordonneesInitiales, i, j);
distances[j][i] = dist(coordonneesInitiales, i, j);
if(typeAffichage==0) System.out.print(distances[i][j] + " ");
}
if(typeAffichage==0) System.out.println();
}
if(typeAffichage==0) System.out.println("Matrice Résultat (intermediaire)");
if(typeAffichage==0) System.out.print("1=> ");
// Le premier point a tout ces coordonnées nulles
// Le deuxième n'a qu'une seule coordonnée nulle
coordonnees[1][0] = distances[1][0];
if(typeAffichage==0) System.out.println(coordonneesInitiales[1][0] + " ");
/* On lance notre algorithme pour trouver les coordonnées des points 3 à n */
double zero = 0;
for (int WichPoint = 2; WichPoint < n; WichPoint++) {
int i = WichPoint; // je met i juste pour simplifier le code, et Wichpoint pour comprendre de
// quoi je parle
if(typeAffichage==0) System.out.print(i + "=> ");
// La norme du vecteur sera sa distance au point 1
double norme = carre(distances[0][i]);
// On calcul d'abord les n-2 premières coordonnées du point courant
for (int p = 0; p < i - 1; p++) {
coordonnees[i][p] = norme - carre(distances[p + 1][i]) + carre(distances[p + 1][0]);
if(typeAffichage==0) System.out.print("p=" + p + " norme-d²=" + coordonnees[i][p]);
double sumCoord = 0;
for (int k = 0; k <= p - 1; k++) {
sumCoord += coordonnees[i][k] * coordonnees[p + 1][k];
}
if(typeAffichage==0) System.out.print(" sumCoord=" + sumCoord + " ");
coordonnees[i][p] = coordonnees[i][p] - 2 * sumCoord;
coordonnees[i][p] = coordonnees[i][p] / (2 * coordonnees[p + 1][p]);
if(typeAffichage==0) System.out.print(" x=" + coordonnees[i][p] + " ");
if(typeAffichage==0) System.out.println();
}
if(typeAffichage==0) System.out.println();
// On calcul enfin la dernière coord du point courant en utilisant les coord calculés à
// l'instant
double sumSQR = 0;
for (int j = 0; j < i - 1; j++) {
sumSQR += carre(coordonnees[i][j]);
}
// Le dernier calcul devrait tjs etre vrai car norme > sumSQR (mathematiquement j'entend)
coordonnees[i][i - 1] = Math.sqrt(norme - sumSQR);
// on test si un des résultats s'approche de zero = controle division par zero
if(Math.abs(coordonnees[i][i-1])<0.0000001){
System.out.println(" petit resultat= " + coordonnees[i][i-1]);
zeroPos[0]=i;
zeroPos[1]=i-1;
zeroPos[2]=coordonnees[i][i-1];
}
// Or y'a du NaN ???? donc norme < sumSQR et c pas normal !!!
// On test le dernier calcul voir quand apparait ce NaN
if (Double.isNaN(coordonnees[i][i - 1]) == true && zero == 0) {
zero++;
System.out.println("Nan a " + i);
System.out.println(" Position du zero : i=" + zeroPos[0] + " j=" + zeroPos[1] + " coord="+ zeroPos[2]);
if(typeAffichage==0) typeAffichage=2; // on chg la valeur pour arreter l'affichage au NaN
}
}
if(typeAffichage==1) System.out.println("Matrice resultat");
for (int i = 0; i < n; i++) {
if(typeAffichage==1) System.out.print(i + "=> ");
for (int j = 0; j < i; j++) {
if(typeAffichage==1) System.out.print(coordonnees[i][j] + " ");
}
if(typeAffichage==1) System.out.println();
}
}
public static void main(String[] args) {
Geotest test = new Geotest();
test.executeTest();
}
}
La matrice à mettre dans un fichier : test_matrix_30.txt
n'as-t'on pas du tout de NaN ? ou alors se produit'il trop loin pour qu'on l'observe ?
Moi pas comprendre.
ben plus le nombre de points est grand plus on des chances de touver un NaN, et dans ton cas si tu fait avec une matrice de 100 tu va rien trouver, alors qu'avec n=1000 tu va en trouver du NaN
Quote:
1) J'ai lancé une myriade fois le code avec matrice aléatoire avec Math.random() mais en forçant, la non nullité de coordonnées et je n'ai pas vu seul Nan
tu as lancé avec quel nombre de points ?
Quote:
3) j'ai fait pareil avec Math.random() * 0.0000001 et j'ai vu des Nan
si tu force la division par zero, tu va en trouver.
Quote:
4) j'ai donc raison depuis longtemps: le code et les fonctions mathématique sont convenables avec des coordonnées non nulles. Mais au voisinage des zéros la division par zéro est problèmatique ce qui est quand même normal, Java ou pas.
certe! mais dans une execution normal, on est jamais proche de zero. Mon avant dernier code avec le test de non nullite t'en convaincra!
Quote:
5) je voudrais bien savoir où facturer pour le million d'Euro.
alors le nom c'est Jean Vouarnac, 13 place de la sainte redemption, 11111 desole cedex
Sans vouloir ton peiner, il est bien compliqué ton code BBReader :shock:
Au fait pendant que j'y pense: mets tes classes dans des PACKAGES screugneugneu (arf ces scientifiques) :evil: :evil:
Donc bien compliqué disais je ce code, avec des buffers psychédéliques ... final byte[] r_buffer = new byte[500000];
Rien que ça.... :shock:
Pourquoi faire simple quand on peut faire compliqué... (arf ces scientifiques)
Bref quand j'ai vu ce code, mauvaise langue comme je suis, je me suis dit:" va pas marcher ce tromblon". Bref, téméraire, je fais quand même un copier/coller et quoi ?
Alors toi, comanktafé pour avoir du NaN avec un code qui marche pas.. Sasréti pas des fois que tu as lancé une autre classe que celle que tu croyais, avec ta manie de pas les mettre dans un package...
Et puis je ne vois pas ce que tu cherches, obstiné que tu es, de plus que ce que je t'ai dis dans mon post où je dis "J'ai lancé une myriade fois le code". Les choses ne peuvent pas être plus claires...
Tes formules marchent sauf au voisinage de 0, parce que là tu t'étais planté. Parce qu'on à pas le droit de diviser par zéro quoi... C'est tout. (arf...)
Mes excuses, j'ai pas voulu me faire chié, en postant le code de lecture de fichier que j'utilise, mais il marche chez moi et pas chez les autres, bref je pense que vous s'aurez faire un code pour lire cette matrice
Quote:
Et puis je ne vois pas ce que tu cherches, obstiné que tu es, de plus que ce que je t'ai dis dans mon post où je dis "J'ai lancé une myriade fois le code". Les choses ne peuvent pas être plus claires...
Lance mon code en mode matrice aléatoire avec 20 points t'aura pas de naN , avec 500 t'en aura,c'est tout. donc je demandais avec combien de point t'avais testé ?
Et c'est toi qu'es borné, les valeurs au voisinage de zero, y'en a pratiquelment aucune quand on lance l'algo avec des matrices aléatoires et pourtant on a du Nan!
j'ai pas voulu me faire chié, en postant le code de lecture de fichier que j'utilise, mais il marche chez moi et pas chez les autres
cool 8)
Quote:
bref je pense que vous s'aurez faire un code pour lire cette matrice
Ben tiens, on a que ça à faire :lol:
Que tu ne veuilles pas te faire chier est quelque chose que je conçois volontiers. Par contre quand d'autres t'aident, (bénévolement, n'est-ce pas M Vouarnac ;) ) il faut essayer de savoir jusqu'où tu peux user les patiences :lol:
Quote:
Lance mon code en mode matrice aléatoire avec 20 points, t'aura pas de naN avec 500 t'en aura
Effectivement.
Quote:
Et c'est toi qu'es borné, les valeurs au voisinage de zero, y'en a pratiquelment aucune quand on lance l'algo avec des matrices aléatoires et pourtant on a du Nan!
:lol:
Borné et obtu, tout le monde le sait ici.
Le borné dit quand même que plus les valeurs sont petites plus les NaN apparaissent vite.
Ce que je remarque aussi, c'est que les Nan apparaissent (chez moi du moins) en bas et à droite du triangle de la matrice triangulaire résultat, où ils y forment à leur tour un triangle.
On dirait que tu as inventés le NaN de Sierpinski :lol:
Sans doute le matheux averti que tu es, et que n'étant pas allé à l'école je ne suis pas, saura comprendre le pourquoi du comment.
Mais à ce niveau, je ne pense pas que le problème soit dans Java ou dans le code stricto sensu.
bref je pense que vous s'aurez faire un code pour lire cette matrice
Je voulez flatter vos ego en fait ;) mais effectivement lu comme tu la lu ca fait vraiment : demerdez vous a debbuger mon code ! et plus vite que ca!
Mille excuses encore une fois !
Les valeurs petites qui donnent du NaN c'est ok, mathématiquement je sais pourquoi, mais c loin d'etre La raison du NaN.
après : NaN de Sierpinski , aucune idée de ce que c'est, je vais voir mes amis wikipedia et google, ils auront la réponse.
Pas de problème :)
Et puis sur ce forum j'en ai vu d'autres.
Quote:
mais c loin d'etre La raison du NaN.
Mais là je dirais que c'est dans tes formules que ça se passe. Et non dans leur codagre. Alors je crains de ne pas t'être d'une grande aide. Même si sur le fond le problème m'intéresse.
D'ailleurs question naïve... quand tu cherches des coordonnées, à partir de distances, les solutions ne sont pas uniques. Tes formules ne donnent qu'une solution unique. Mais je pense que ceci est chose évidente pour toi et que tu y as pensé.
Quote:
NaN de Sierpinski , aucune idée de ce que c'est
Un triangle de Sierpinski construit avec des NaN :lol:
Comme je ne sais pas faire des maths, et comme je ne fais qu'une chose de bien, c'est de tracer l'exécution du programme, je vous propose de modifier légèrement la boucle du calcul de sumSQR avec la code que voici :
double sumSQRant = 0;
for (int j = 0; j < i - 1; j++) {
sumSQR += carre(coordonnees[i][j]);
if (Double.isNaN(sumSQR)) {
System.out.println("NaN dans le carré sumSQRant, coordonnees[i][j], carre(coordonnees[i][j]):"+sumSQRant+"|"+coordonnees[i][j]+"|"+(carre(coordonnees[i][j])));
}
sumSQRant = sumSQR;
}
Bref, je cherche simplement ou se déclenche le NaN. (je vais finir par aller manger Indien à force de parler de NaN .... j'ai pas pu m'empêcher de la faire celle-là).
J'ai déclenché plusieurs éxécutions aléatoires, car j'ai la flemme d'écrire la fonction de lecture du fichier, mais bon, je ferais peut-être ça plus tard...
Et voilà ce que j'obtiens comme résultat :
NaN dans le carré i, j, sumSQRant, coordonnees[i][j], carre(coordonnees[i][j]):37|35|15.706202875865124|NaN|NaN
NaN dans le carré i, j, sumSQRant, coordonnees[i][j], carre(coordonnees[i][j]):38|35|11.689362043699964|NaN|NaN
NaN dans le carré i, j, sumSQRant, coordonnees[i][j], carre(coordonnees[i][j]):38|36|NaN|NaN|NaN
NaN dans le carré i, j, sumSQRant, coordonnees[i][j], carre(coordonnees[i][j]):39|35|10.739514060053297|NaN|NaN
NaN dans le carré i, j, sumSQRant, coordonnees[i][j], carre(coordonnees[i][j]):39|36|NaN|NaN|NaN
NaN dans le carré i, j, sumSQRant, coordonnees[i][j], carre(coordonnees[i][j]):39|37|NaN|NaN|NaN
Le Nan se propage car il s'appuit sur un Nan précédement calculé. Le Nan précédement calculé provient de la racine carré d'un nombre négatif.
J'ai fait afficher norme et sumSQR côte à côte et l'on voit bien le phénomène se produire (sumSQR>norme)
Même en utilisant des BigDécimaux, le problème se pose également.
Donc si tout est juste (formule math, algo etc...), soit java ne peut pas t'apporte la précision que nécessite tes calculs, donc change de language pour faire ton problème, soit, puisque le problème temps a disparaître avec des nombres plus grands, transpose tes unités de mesure si c'est possible (*100 ou *1000) et ton problème sera résolue. En tout cas, je lirai la suite de cette aventure avec plaisir, mais là moi, je jette l'éponge
transposer mes unités de mesure c'est une idée auquel je vais réfléchir !
Hier dans le ReR de retour j'ai eu une idée sur comment le naN pourrai apparaitre qui me semble bonne. Faut que je fasse un peu de math pour le prouver donc je vais voir ca aujourd'hui, je vous tiens au courant si y'a des news!
transposer mes unités de mesure c'est une idée auquel je vais réfléchir !
Déjà fait. Mais je me demande si je l'ai dit hier.
J'ai lancé un calcul sur mille points aléatoires, mais sans valuer nulle tirés par Math.random() *1000 et j'ai vu du NaN. Sans doute comme dit willbback un Nan qui se propage en effet, d'où une apparition en remplissage d'un triangle en bas à droite de la matrice résultat.
Quote:
Hier dans le ReR de retour
Allo ?
Quote:
que je fasse un peu de math pour le prouver donc je vais voir ca aujourd'hui, je vous tiens au courant si y'a des news!
Oui tiens nous au courant. Ca m'intéresse, bien que, je dois l'avouer, je n'ai toujours pas regardé tes formules de près faute de temps :oops: En admettant par ailleurs que j'en sois capable ce qui n'est pas sûr du tout :)
Bon j'annonce solenellement que je tiens quelque chose de carrèment bien, et qui expliquerai tout !
Mais d'abord faisons encore une fois un résumé de ce que je veux faire :
J'ai une matrice de distance et je veux trouver une matrice de position.
J'ai donc un programme de test, qui lui tire une matrice au hasard (triangulaire) calcul les distances et a partir de celle ci trouve la matrice position. Ce programme plante au bout d'un moment en ressortant du NaN.
Une facon de controler la justesse de ce programme est de verifier que les distances entre les points reconstruit sont bien les mêmes que celle de la matrice de distance de départ. Et ca je le verifie parfaitement à 10^-15 près en moyenne.
Il y a un mois je vérifiais si mon programme était correcte d'une autre manière : Je regardais la différence entre la matrice de positionInitiale et la matrice de positionFinale. Et cette étude nous montrait une divergence (desole je me rend compte que je suis con je vous l'ai jamais dit). C'est a dire que sur la diagonale des matrices position plus on avance et plus la difference est grande (10^-17 puis 10-16 puis 10^-15 etc ...).
Alors pourquoi je ne vous l'avais pas dit, et que maintenant je ressors cette histoire de divergence des positions du chapeau ?
C'est parce que les distances, elles, ne montrent aucune divergence. Cette matrice de distance étant la seule accessible dans mon programme futur, et donc la seule matrice interessante, et sur lequel je dois me baser, dans mon futur programme la matrice positionInitiales ne sera pas accessible.
J'avais donc éludé la question des divergences de position en me disant que y'avait pas de raison pour que la matrice positionInitiales et la Finale soit les mêmes . alors que si !!!
J'ai montré ce matin avec les math que c obligatoirement les mêmes! Donc la divergence qu'on observe reflete totalement une divergence dans l'algorithme!!!!
La blague qui fait qu'on le voit pas est que : Les positions divergent, mais la structure de mes équations font que les distances ne connaissent aucune divergence! et cela est du au fait que mes calculs utilisent les distances pour trouver les positions.
En résumé : Moi y'en a avoir une matrice diagonale de position au départ -> moi y'en a retrouver la même à la fin, mais ca y'en pas marcher car ca diverge. Mais moi y'en a pas le voir avec ma facon de verifier les calculs, car les distances entre les positions y'en a rester constante à cause de la structure de mes équations!
Bref moi y'en a tres tres content d'avoir trouvé ca !!!
C la fête, merci à tous pour votre aide. Ca m'a permis de bien defricher ce qui pouvait etre et ce qui ne pouvait pas être.
Maintenant si le sujet vous interesse vraiment la question nouvelle est de savoir, comment controler la divergence ???
Car dans mon programme finale je n'aurai aucun acces à la matrice positionInitiales, il faut que je me demerde juste avec la matrice de distance pour controler cette divergence! or cette matrice ne subit aucune divergence .... :roll:
Ah la recherche, j'adore ! maintenant je sais ou je dois aller, j'ai plus l'impression d'etre emprisonné par l'ordi, ca soulage.
On the road again .... :P
J'ai absolument rien compris à ton explication.... Ce qui veut dire que ce n'est plus un problème de précision de Java, c'est ça ? :?
J'ai refait des essais et c'est bien sumSQR qui devient supérieur à norme dans certain cas. Un NaN apparaît alors dans la matrice coordonnee, puis est utilisé dans les calculs de la boucle suivante, provancant alors un autre Nan et ainsi de suite. C'est pour celà que les Nan sont bien rangé dans la matrice.
Le fichier qui fait des Nan à tous les coups sans 0 null part :
"norme - sumSQR"= -0.11681206539908295 (en mode debug avec "i"= 38
)
Pas de division par zero, rien un gros problème d'algorithmie, c'est tout.
Bref ce fut passionnant.
Bon courrage pour la suite de tes aventures.
Ah si, juste une remarque : il est inutile de faire des calculs avec des précisions astronomiques mais avec une précision légèrement inférieure à celles des mesures (1/100) en principe le résultat reste juste. Je crois que c'est un reste de mes cours de physique
lance ce joli code, et le secret sur la divergence tu comprendra! Pour le point situé juste avant le NaN je calcule les différences entre ce qu'on trouve et ce qu'on aurait du avoir, et tu vois que les erreur sont de plus en plus grande:
class Geotest {
final int NB_POINTS = 50;
// type Matrice de départ : 0 = c'est nous qui la créons ; 1 = aléatoire ; 2=on la lit dans un fichier
final int typeMat = 1;
/* typeAffichage :
* 0 = on affiche toutes les étapes de calcul
* 1 = on affiche juste résultats
* 2 = on affiche juste quand y'a du NaN
* 3 = difference entre somme des distance init et finale
* 4 = on affiche difference entre matrice position depart et matrice position finale */
public int typeAffichage = 4;
int quelLigne=9;
// n = nombre de points
// m = dimension de l'espace
private int n = NB_POINTS;
private int m = NB_POINTS - 1;
// Matrice des positions initiales
public double[][] coordonneesInitiales = new double[n][m];
// Matrice des positions reconstruites par l'algorithme
public double[][] coordonnees = new double[n][m];
// Matrice des distances
public double[][] distances = new double[n][n];
// vecteur pour controle de division par zero
// zeroPos[0]=position en i ; zeroPos[1]=position en j ; zeroPos[2] = valeur de coord
public double[] zeroPos = { -1,-1,-1};
// petite methode pour calculer un carre
double carre(double x) {
return x * x;
}
// Calcul de la distance euclidienne pour la matrice de position
double dist(double[][] coords, int a, int b) {
double sum = 0;
for (int i = 0; i < m; i++) {
sum = sum + carre(coords[a][i] - coords[b][i]);
}
return Math.sqrt(sum);
}
public void creationMatrice(int type){
if(typeAffichage==0) System.out.println("Matrice Initiale");
switch(type) {
case 0 :
// ON cree nous même la matrice
double[][] temp =
{
{ 0.0, 0.0, 0.0 } ,
{ 2.0, 0.0, 0.0 },
{ 7.0, 0.0, 0.0 },
{ 0.0, 1.0, 2.0 }
};
for(int i=0;i<n;i++){
if(typeAffichage==0) System.out.print(i + "=> ");
for(int j=0;j<m;j++){
coordonneesInitiales[i][j]= temp[i][j];
if(typeAffichage==0) System.out.print(coordonneesInitiales[i][j] + " ");
}
if(typeAffichage==0) System.out.println();
}
break;
case 1 :
// On crée une matrice triangulaire aléatoire
for (int i = 0; i < n; i++) {
if(typeAffichage==0) System.out.print(i + "=> ");
for (int j = 0; j < i; j++) {
coordonneesInitiales[i][j] = (Math.random()*0.1);
if(typeAffichage==0) System.out.print(coordonneesInitiales[i][j] + " ");
}
if(typeAffichage==0) System.out.println();
}
break;
}
}
public void executeTest() {
creationMatrice(typeMat);
// On calcul les distances euclidiennes entre les points
if(typeAffichage==0) System.out.println("Matrice distance");
for (int i = 0; i < n; i++) {
//System.out.print(i + "=> ");
for (int j = 0; j < m; j++) {
distances[i][j] = dist(coordonneesInitiales, i, j);
distances[j][i] = dist(coordonneesInitiales, i, j);
if(typeAffichage==0) System.out.print(distances[i][j] + " ");
}
if(typeAffichage==0) System.out.println();
}
if(typeAffichage==0) System.out.println("Matrice Résultat (intermediaire)");
if(typeAffichage==0) System.out.print("1=> ");
// Le premier point a tout ces coordonnées nulles
// Le deuxième n'a qu'une seule coordonnée nulle
coordonnees[1][0] = distances[1][0];
if(typeAffichage==0) System.out.println(coordonneesInitiales[1][0] + " ");
/* On lance notre algorithme pour trouver les coordonnées des points 3 à n */
double zero = 0;
for (int WichPoint = 2; WichPoint < n; WichPoint++) {
int i = WichPoint; // je met i juste pour simplifier le code, et Wichpoint pour comprendre de
// quoi je parle
if(typeAffichage==0) System.out.print(i + "=> ");
// La norme du vecteur sera sa distance au point 1
double norme = carre(distances[0][i]);
// On calcul d'abord les n-2 premières coordonnées du point courant
for (int p = 0; p < i - 1; p++) {
coordonnees[i][p] = norme - carre(distances[p + 1][i]) + carre(distances[p + 1][0]);
if(typeAffichage==0) System.out.print("p=" + p + " norme-d²=" + coordonnees[i][p]);
double sumCoord = 0;
for (int k = 0; k <= p - 1; k++) {
sumCoord += coordonnees[i][k] * coordonnees[p + 1][k];
}
if(typeAffichage==0) System.out.print(" sumCoord=" + sumCoord + " ");
coordonnees[i][p] = coordonnees[i][p] - 2 * sumCoord;
coordonnees[i][p] = coordonnees[i][p] / (2 * coordonnees[p + 1][p]);
// on controle la divergence en soustrayant ou additionnant une erreur qui pourrait apparaitre
if(typeAffichage==0) System.out.print(" x=" + coordonnees[i][p] + " ");
if(typeAffichage==0) System.out.println();
}
if(typeAffichage==0) System.out.println();
// On calcul enfin la dernière coord du point courant en utilisant les coord calculés à
// l'instant
double sumSQR = 0;
for (int j = 0; j < i - 1; j++) {
sumSQR += carre(coordonnees[i][j]);
}
// Le dernier calcul devrait tjs etre vrai car norme > sumSQR (mathematiquement j'entend)
coordonnees[i][i - 1] = Math.sqrt(norme - sumSQR);
if(i%2==0) coordonnees[i][i - 1]=coordonnees[i][i - 1]-Math.pow(10,-17);
else coordonnees[i][i - 1]=coordonnees[i][i - 1]+Math.pow(10,-17);
// on test si un des résultats s'approche de zero = controle division par zero
if(Math.abs(coordonnees[i][i-1])<0.0000001){
System.out.println(" petit resultat= " + coordonnees[i][i-1]);
zeroPos[0]=i;
zeroPos[1]=i-1;
zeroPos[2]=coordonnees[i][i-1];
}
// Or y'a du NaN ???? donc norme < sumSQR et c pas normal !!!
// On test le dernier calcul voir quand apparait ce NaN
if (Double.isNaN(coordonnees[i][i - 1]) == true && zero == 0) {
zero++;
quelLigne=i-1;
System.out.println("Nan a " + i);
System.out.println(" Position du zero : i=" + zeroPos[0] + " j=" + zeroPos[1] + " coord="+ zeroPos[2]);
if(typeAffichage==0) typeAffichage=2; // on chg la valeur pour arreter l'affichage au NaN
}
}
if(typeAffichage==1) System.out.println("Matrice resultat");
for (int i = 0; i < n; i++) {
if(typeAffichage==1) System.out.print(i + "=> ");
for (int j = 0; j < i; j++) {
if(typeAffichage==1) System.out.print(coordonnees[i][j] + " ");
}
if(typeAffichage==1) System.out.println();
}
// sommme des distances
if(typeAffichage==3){
for(int i=0;i<quelLigne;i++){
double sumInit=0;
double sumFinale=0;
System.out.print(i + "=> ");
for(int j=0;j<i;j++){
sumInit+=distances[i][j];
sumFinale+=dist(coordonnees,i,j);
}
if(typeAffichage==3) System.out.println(" erreur somme des distances= " + (sumInit-sumFinale));
}
}
// difference de coordonnées
if(typeAffichage==4){
for(int j=0;j<quelLigne-1;j++){
System.out.print("j="+j+" erreur="+ (coordonneesInitiales[quelLigne-1][j]-coordonnees[quelLigne-1][j]));
System.out.println();
}
}
}
public static void main(String[] args) {
Geotest test = new Geotest();
test.executeTest();
}
}
change maintenant typeAffichage=4 en typeAffichage=3 et tu verra que les distances ne divergent pas.
Ah :!: :lol:
Y a qu'un parisien qui se croit au centre du monde pour penser que c'est évident pour tout monde.
Moi quand j'ai vu ReR (et non RER), dans le contexte de la conversation, je me suis demandé quel coin obscur de l'informatique j'avais l'infortune d'ignorer. :lol:
Ah... la région parisienne... quand je pense que j'ai habité cet asile de dingues pendant 12 longues années. Mais je me suis évadé! Pour me réfugier dans un trou du sud de la France on on sait à peine ce que c'est qu'une voiture :) Alors le RER...
Le RER... c'était dans une autre vie pour moi ça :)
Pour ton problème je crois que j'ai compris au moins l'idée dans les grandes lignes. Mais ce que je ne comprends pas c'est pourquoi le calcul diverge. Si à l'occasion tu a sle temps d'écrire deux ou trois lignes bien vulgarisées à l'attention du primate que je suis :)
evidemment c le contraire que je voulais ecrire :D
qd tu passe de double à float les naN sont plus proche. Ce qui parait normal dans l'hypothese de la precision qui fout la merde!
j'edit cette phrase desuite.
avoue que sinon t'aurai rien eu a dire ! et puis ca a fait un petit moment de suspens ! Bref j'essaye de justifier le fait que je me suis trompé ;)
Moi pas comprendre.
Ce que je peux dire c'est
1) J'ai lancé une myriade fois le code avec matrice aléatoire avec Math.random() mais en forçant, la non nullité de coordonnées et je n'ai pas vu seul Nan
2) j'ai fait pareil avec Math.random() * 0.0001 et je n'ai pas vu un NaN
3) j'ai fait pareil avec Math.random() * 0.0000001 et j'ai vu des Nan
4) j'ai donc raison depuis longtemps: le code et les fonctions mathématique sont convenables avec des coordonnées non nulles. Mais au voisinage des zéros la division par zéro est problèmatique ce qui est quand même normal, Java ou pas.
5) je voudrais bien savoir où facturer pour le million d'Euro.
par contre j'ai un code tout pret avec une matrice tout jolie qui donne du naN à 26. Le code est adapté pour lire les fichiers txt ou il y a la matrice d'ecrite.
Donc y'a le code Geotest, et la matrice qu'il faut mettre dans un fichier test_matrix_30.txt, et puis creer une methode qui lit cette matrice (moi je l'ai dcette methode, mais secret defense je la garde, et en plus ca marche pas chez certain)
desolé pour le post immense, mais j'ai pas reussi à mettre sur mon serveur web, il plante ce con, pil poil quand j'en ai besoin !
Le code de Geotest.java
import java.io.File; import java.util.ArrayList; class Geotest { final int NB_POINTS = 30; // type Matrice de départ : 0 = c'est nous qui la créons ; 1 = aléatoire ; 2=on la lit dans un fichier final int typeMat = 2; /* typeAffichage : * 0 = on affiche toutes les étapes de calcul * 1 = on affiche juste résultats * 2 = on affiche juste quand y'a du NaN */ public int typeAffichage = 2; // n = nombre de points // m = dimension de l'espace private int n = NB_POINTS; private int m = NB_POINTS - 1; // Matrice des positions initiales public double[][] coordonneesInitiales = new double[n][m];; // Matrice des positions reconstruites par l'algorithme public double[][] coordonnees = new double[n][m]; // Matrice des distances public double[][] distances = new double[n][n]; // vecteur pour controle de division par zero // zeroPos[0]=position en i ; zeroPos[1]=position en j ; zeroPos[2] = valeur de coord public double[] zeroPos = { -1,-1,-1}; // petite methode pour calculer un carre double carre(double x) { return x * x; } // Calcul de la distance euclidienne pour la matrice de position double dist(double[][] coords, int a, int b) { double sum = 0; for (int i = 0; i < m; i++) { sum = sum + carre(coords[a][i] - coords[b][i]); } return Math.sqrt(sum); } public void creationMatrice(int type){ if(typeAffichage==0) System.out.println("Matrice Initiale"); switch(type) { case 0 : // ON cree nous même la matrice double[][] temp = { { 0.0, 0.0, 0.0 } , { 2.0, 0.0, 0.0 }, { 7.0, 0.0, 0.0 }, { 0.0, 1.0, 2.0 } }; for(int i=0;i<n;i++){ if(typeAffichage==0) System.out.print(i + "=> "); for(int j=0;j<m;j++){ coordonneesInitiales[i][j]= temp[i][j]; if(typeAffichage==0) System.out.print(coordonneesInitiales[i][j] + " "); } if(typeAffichage==0) System.out.println(); } break; case 1 : // On crée une matrice triangulaire aléatoire for (int i = 0; i < n; i++) { if(typeAffichage==0) System.out.print(i + "=> "); for (int j = 0; j < i; j++) { coordonneesInitiales[i][j] = (Math.random()*10); if(typeAffichage==0) System.out.print(coordonneesInitiales[i][j] + " "); } if(typeAffichage==0) System.out.println(); } break; case 2 : coordonneesInitiales = readMatrixFromFile(new File("test_matrix_30.txt")); } } public void executeTest() { creationMatrice(typeMat); // On calcul les distances euclidiennes entre les points if(typeAffichage==0) System.out.println("Matrice distance"); for (int i = 0; i < n; i++) { //System.out.print(i + "=> "); for (int j = 0; j < m; j++) { distances[i][j] = dist(coordonneesInitiales, i, j); distances[j][i] = dist(coordonneesInitiales, i, j); if(typeAffichage==0) System.out.print(distances[i][j] + " "); } if(typeAffichage==0) System.out.println(); } if(typeAffichage==0) System.out.println("Matrice Résultat (intermediaire)"); if(typeAffichage==0) System.out.print("1=> "); // Le premier point a tout ces coordonnées nulles // Le deuxième n'a qu'une seule coordonnée nulle coordonnees[1][0] = distances[1][0]; if(typeAffichage==0) System.out.println(coordonneesInitiales[1][0] + " "); /* On lance notre algorithme pour trouver les coordonnées des points 3 à n */ double zero = 0; for (int WichPoint = 2; WichPoint < n; WichPoint++) { int i = WichPoint; // je met i juste pour simplifier le code, et Wichpoint pour comprendre de // quoi je parle if(typeAffichage==0) System.out.print(i + "=> "); // La norme du vecteur sera sa distance au point 1 double norme = carre(distances[0][i]); // On calcul d'abord les n-2 premières coordonnées du point courant for (int p = 0; p < i - 1; p++) { coordonnees[i][p] = norme - carre(distances[p + 1][i]) + carre(distances[p + 1][0]); if(typeAffichage==0) System.out.print("p=" + p + " norme-d²=" + coordonnees[i][p]); double sumCoord = 0; for (int k = 0; k <= p - 1; k++) { sumCoord += coordonnees[i][k] * coordonnees[p + 1][k]; } if(typeAffichage==0) System.out.print(" sumCoord=" + sumCoord + " "); coordonnees[i][p] = coordonnees[i][p] - 2 * sumCoord; coordonnees[i][p] = coordonnees[i][p] / (2 * coordonnees[p + 1][p]); if(typeAffichage==0) System.out.print(" x=" + coordonnees[i][p] + " "); if(typeAffichage==0) System.out.println(); } if(typeAffichage==0) System.out.println(); // On calcul enfin la dernière coord du point courant en utilisant les coord calculés à // l'instant double sumSQR = 0; for (int j = 0; j < i - 1; j++) { sumSQR += carre(coordonnees[i][j]); } // Le dernier calcul devrait tjs etre vrai car norme > sumSQR (mathematiquement j'entend) coordonnees[i][i - 1] = Math.sqrt(norme - sumSQR); // on test si un des résultats s'approche de zero = controle division par zero if(Math.abs(coordonnees[i][i-1])<0.0000001){ System.out.println(" petit resultat= " + coordonnees[i][i-1]); zeroPos[0]=i; zeroPos[1]=i-1; zeroPos[2]=coordonnees[i][i-1]; } // Or y'a du NaN ???? donc norme < sumSQR et c pas normal !!! // On test le dernier calcul voir quand apparait ce NaN if (Double.isNaN(coordonnees[i][i - 1]) == true && zero == 0) { zero++; System.out.println("Nan a " + i); System.out.println(" Position du zero : i=" + zeroPos[0] + " j=" + zeroPos[1] + " coord="+ zeroPos[2]); if(typeAffichage==0) typeAffichage=2; // on chg la valeur pour arreter l'affichage au NaN } } if(typeAffichage==1) System.out.println("Matrice resultat"); for (int i = 0; i < n; i++) { if(typeAffichage==1) System.out.print(i + "=> "); for (int j = 0; j < i; j++) { if(typeAffichage==1) System.out.print(coordonnees[i][j] + " "); } if(typeAffichage==1) System.out.println(); } } public static void main(String[] args) { Geotest test = new Geotest(); test.executeTest(); } }La matrice à mettre dans un fichier : test_matrix_30.txt
[/code]ben plus le nombre de points est grand plus on des chances de touver un NaN, et dans ton cas si tu fait avec une matrice de 100 tu va rien trouver, alors qu'avec n=1000 tu va en trouver du NaN
tu as lancé avec quel nombre de points ?
si tu force la division par zero, tu va en trouver.
certe! mais dans une execution normal, on est jamais proche de zero. Mon avant dernier code avec le test de non nullite t'en convaincra!
Sans vouloir ton peiner, il est bien compliqué ton code BBReader :shock:
Au fait pendant que j'y pense: mets tes classes dans des PACKAGES screugneugneu (arf ces scientifiques) :evil: :evil:
Donc bien compliqué disais je ce code, avec des buffers psychédéliques ...
final byte[] r_buffer = new byte[500000];Rien que ça.... :shock:
Pourquoi faire simple quand on peut faire compliqué... (arf ces scientifiques)
Bref quand j'ai vu ce code, mauvaise langue comme je suis, je me suis dit:" va pas marcher ce tromblon". Bref, téméraire, je fais quand même un copier/coller et quoi ?
Exception in thread "main" java.lang.NumberFormatException: For input string: "0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0"
at sun.misc.FloatingDecimal.readJavaFormatString(FloatingDecimal.java:1224)
at java.lang.Double.parseDouble(Double.java:510)
at readtest.Geotest.readMatrixFromFile(Geotest.java:202)
at readtest.Geotest.creationMatrice(Geotest.java:84)
at readtest.Geotest.executeTest(Geotest.java:90)
at readtest.Geotest.main(Geotest.java:177)
Java Result: 1
:twisted: :lol:
Alors toi, comanktafé pour avoir du NaN avec un code qui marche pas.. Sasréti pas des fois que tu as lancé une autre classe que celle que tu croyais, avec ta manie de pas les mettre dans un package...
Et puis je ne vois pas ce que tu cherches, obstiné que tu es, de plus que ce que je t'ai dis dans mon post où je dis "J'ai lancé une myriade fois le code". Les choses ne peuvent pas être plus claires...
Tes formules marchent sauf au voisinage de 0, parce que là tu t'étais planté. Parce qu'on à pas le droit de diviser par zéro quoi... C'est tout. (arf...)
MON CHEEEQQUUUEEEEEEE !!!!!!!!!!!!
:twisted: :twisted: :twisted:
Mes excuses, j'ai pas voulu me faire chié, en postant le code de lecture de fichier que j'utilise, mais il marche chez moi et pas chez les autres, bref je pense que vous s'aurez faire un code pour lire cette matrice
Lance mon code en mode matrice aléatoire avec 20 points t'aura pas de naN , avec 500 t'en aura,c'est tout. donc je demandais avec combien de point t'avais testé ?
Et c'est toi qu'es borné, les valeurs au voisinage de zero, y'en a pratiquelment aucune quand on lance l'algo avec des matrices aléatoires et pourtant on a du Nan!
cool 8)
Ben tiens, on a que ça à faire :lol:
Que tu ne veuilles pas te faire chier est quelque chose que je conçois volontiers. Par contre quand d'autres t'aident, (bénévolement, n'est-ce pas M Vouarnac ;) ) il faut essayer de savoir jusqu'où tu peux user les patiences :lol:
Effectivement.
:lol:
Borné et obtu, tout le monde le sait ici.
Le borné dit quand même que plus les valeurs sont petites plus les NaN apparaissent vite.
Ce que je remarque aussi, c'est que les Nan apparaissent (chez moi du moins) en bas et à droite du triangle de la matrice triangulaire résultat, où ils y forment à leur tour un triangle.
On dirait que tu as inventés le NaN de Sierpinski :lol:
Sans doute le matheux averti que tu es, et que n'étant pas allé à l'école je ne suis pas, saura comprendre le pourquoi du comment.
Mais à ce niveau, je ne pense pas que le problème soit dans Java ou dans le code stricto sensu.
desole pour le :
Je voulez flatter vos ego en fait ;) mais effectivement lu comme tu la lu ca fait vraiment : demerdez vous a debbuger mon code ! et plus vite que ca!
Mille excuses encore une fois !
Les valeurs petites qui donnent du NaN c'est ok, mathématiquement je sais pourquoi, mais c loin d'etre La raison du NaN.
après : NaN de Sierpinski , aucune idée de ce que c'est, je vais voir mes amis wikipedia et google, ils auront la réponse.
Pas de problème :)
Et puis sur ce forum j'en ai vu d'autres.
Mais là je dirais que c'est dans tes formules que ça se passe. Et non dans leur codagre. Alors je crains de ne pas t'être d'une grande aide. Même si sur le fond le problème m'intéresse.
D'ailleurs question naïve... quand tu cherches des coordonnées, à partir de distances, les solutions ne sont pas uniques. Tes formules ne donnent qu'une solution unique. Mais je pense que ceci est chose évidente pour toi et que tu y as pensé.
Un triangle de Sierpinski construit avec des NaN :lol:
Comme je ne sais pas faire des maths, et comme je ne fais qu'une chose de bien, c'est de tracer l'exécution du programme, je vous propose de modifier légèrement la boucle du calcul de sumSQR avec la code que voici :
double sumSQRant = 0; for (int j = 0; j < i - 1; j++) { sumSQR += carre(coordonnees[i][j]); if (Double.isNaN(sumSQR)) { System.out.println("NaN dans le carré sumSQRant, coordonnees[i][j], carre(coordonnees[i][j]):"+sumSQRant+"|"+coordonnees[i][j]+"|"+(carre(coordonnees[i][j]))); } sumSQRant = sumSQR; }Bref, je cherche simplement ou se déclenche le NaN. (je vais finir par aller manger Indien à force de parler de NaN .... j'ai pas pu m'empêcher de la faire celle-là).
J'ai déclenché plusieurs éxécutions aléatoires, car j'ai la flemme d'écrire la fonction de lecture du fichier, mais bon, je ferais peut-être ça plus tard...
Et voilà ce que j'obtiens comme résultat :
Le Nan se propage car il s'appuit sur un Nan précédement calculé. Le Nan précédement calculé provient de la racine carré d'un nombre négatif.
J'ai fait afficher norme et sumSQR côte à côte et l'on voit bien le phénomène se produire (sumSQR>norme)
Même en utilisant des BigDécimaux, le problème se pose également.
Donc si tout est juste (formule math, algo etc...), soit java ne peut pas t'apporte la précision que nécessite tes calculs, donc change de language pour faire ton problème, soit, puisque le problème temps a disparaître avec des nombres plus grands, transpose tes unités de mesure si c'est possible (*100 ou *1000) et ton problème sera résolue. En tout cas, je lirai la suite de cette aventure avec plaisir, mais là moi, je jette l'éponge
transposer mes unités de mesure c'est une idée auquel je vais réfléchir !
Hier dans le ReR de retour j'ai eu une idée sur comment le naN pourrai apparaitre qui me semble bonne. Faut que je fasse un peu de math pour le prouver donc je vais voir ca aujourd'hui, je vous tiens au courant si y'a des news!
Déjà fait. Mais je me demande si je l'ai dit hier.
J'ai lancé un calcul sur mille points aléatoires, mais sans valuer nulle tirés par Math.random() *1000 et j'ai vu du NaN. Sans doute comme dit willbback un Nan qui se propage en effet, d'où une apparition en remplissage d'un triangle en bas à droite de la matrice résultat.
Allo ?
Oui tiens nous au courant. Ca m'intéresse, bien que, je dois l'avouer, je n'ai toujours pas regardé tes formules de près faute de temps :oops: En admettant par ailleurs que j'en sois capable ce qui n'est pas sûr du tout :)
ben le RER quoi ?!
Bon j'annonce solenellement que je tiens quelque chose de carrèment bien, et qui expliquerai tout !
Mais d'abord faisons encore une fois un résumé de ce que je veux faire :
J'ai une matrice de distance et je veux trouver une matrice de position.
J'ai donc un programme de test, qui lui tire une matrice au hasard (triangulaire) calcul les distances et a partir de celle ci trouve la matrice position. Ce programme plante au bout d'un moment en ressortant du NaN.
Une facon de controler la justesse de ce programme est de verifier que les distances entre les points reconstruit sont bien les mêmes que celle de la matrice de distance de départ. Et ca je le verifie parfaitement à 10^-15 près en moyenne.
Il y a un mois je vérifiais si mon programme était correcte d'une autre manière : Je regardais la différence entre la matrice de positionInitiale et la matrice de positionFinale. Et cette étude nous montrait une divergence (desole je me rend compte que je suis con je vous l'ai jamais dit). C'est a dire que sur la diagonale des matrices position plus on avance et plus la difference est grande (10^-17 puis 10-16 puis 10^-15 etc ...).
Alors pourquoi je ne vous l'avais pas dit, et que maintenant je ressors cette histoire de divergence des positions du chapeau ?
C'est parce que les distances, elles, ne montrent aucune divergence. Cette matrice de distance étant la seule accessible dans mon programme futur, et donc la seule matrice interessante, et sur lequel je dois me baser, dans mon futur programme la matrice positionInitiales ne sera pas accessible.
J'avais donc éludé la question des divergences de position en me disant que y'avait pas de raison pour que la matrice positionInitiales et la Finale soit les mêmes . alors que si !!!
J'ai montré ce matin avec les math que c obligatoirement les mêmes! Donc la divergence qu'on observe reflete totalement une divergence dans l'algorithme!!!!
La blague qui fait qu'on le voit pas est que : Les positions divergent, mais la structure de mes équations font que les distances ne connaissent aucune divergence! et cela est du au fait que mes calculs utilisent les distances pour trouver les positions.
En résumé : Moi y'en a avoir une matrice diagonale de position au départ -> moi y'en a retrouver la même à la fin, mais ca y'en pas marcher car ca diverge. Mais moi y'en a pas le voir avec ma facon de verifier les calculs, car les distances entre les positions y'en a rester constante à cause de la structure de mes équations!
Bref moi y'en a tres tres content d'avoir trouvé ca !!!
C la fête, merci à tous pour votre aide. Ca m'a permis de bien defricher ce qui pouvait etre et ce qui ne pouvait pas être.
Maintenant si le sujet vous interesse vraiment la question nouvelle est de savoir, comment controler la divergence ???
Car dans mon programme finale je n'aurai aucun acces à la matrice positionInitiales, il faut que je me demerde juste avec la matrice de distance pour controler cette divergence! or cette matrice ne subit aucune divergence .... :roll:
Ah la recherche, j'adore ! maintenant je sais ou je dois aller, j'ai plus l'impression d'etre emprisonné par l'ordi, ca soulage.
On the road again .... :P
J'ai absolument rien compris à ton explication.... Ce qui veut dire que ce n'est plus un problème de précision de Java, c'est ça ? :?
J'ai refait des essais et c'est bien sumSQR qui devient supérieur à norme dans certain cas. Un NaN apparaît alors dans la matrice coordonnee, puis est utilisé dans les calculs de la boucle suivante, provancant alors un autre Nan et ainsi de suite. C'est pour celà que les Nan sont bien rangé dans la matrice.
Le fichier qui fait des Nan à tous les coups sans 0 null part :
"norme - sumSQR"= -0.11681206539908295 (en mode debug avec "i"= 38
)
Pas de division par zero, rien un gros problème d'algorithmie, c'est tout.
Bref ce fut passionnant.
Bon courrage pour la suite de tes aventures.
Ah si, juste une remarque : il est inutile de faire des calculs avec des précisions astronomiques mais avec une précision légèrement inférieure à celles des mesures (1/100) en principe le résultat reste juste. Je crois que c'est un reste de mes cours de physique
lance ce joli code, et le secret sur la divergence tu comprendra! Pour le point situé juste avant le NaN je calcule les différences entre ce qu'on trouve et ce qu'on aurait du avoir, et tu vois que les erreur sont de plus en plus grande:
class Geotest { final int NB_POINTS = 50; // type Matrice de départ : 0 = c'est nous qui la créons ; 1 = aléatoire ; 2=on la lit dans un fichier final int typeMat = 1; /* typeAffichage : * 0 = on affiche toutes les étapes de calcul * 1 = on affiche juste résultats * 2 = on affiche juste quand y'a du NaN * 3 = difference entre somme des distance init et finale * 4 = on affiche difference entre matrice position depart et matrice position finale */ public int typeAffichage = 4; int quelLigne=9; // n = nombre de points // m = dimension de l'espace private int n = NB_POINTS; private int m = NB_POINTS - 1; // Matrice des positions initiales public double[][] coordonneesInitiales = new double[n][m]; // Matrice des positions reconstruites par l'algorithme public double[][] coordonnees = new double[n][m]; // Matrice des distances public double[][] distances = new double[n][n]; // vecteur pour controle de division par zero // zeroPos[0]=position en i ; zeroPos[1]=position en j ; zeroPos[2] = valeur de coord public double[] zeroPos = { -1,-1,-1}; // petite methode pour calculer un carre double carre(double x) { return x * x; } // Calcul de la distance euclidienne pour la matrice de position double dist(double[][] coords, int a, int b) { double sum = 0; for (int i = 0; i < m; i++) { sum = sum + carre(coords[a][i] - coords[b][i]); } return Math.sqrt(sum); } public void creationMatrice(int type){ if(typeAffichage==0) System.out.println("Matrice Initiale"); switch(type) { case 0 : // ON cree nous même la matrice double[][] temp = { { 0.0, 0.0, 0.0 } , { 2.0, 0.0, 0.0 }, { 7.0, 0.0, 0.0 }, { 0.0, 1.0, 2.0 } }; for(int i=0;i<n;i++){ if(typeAffichage==0) System.out.print(i + "=> "); for(int j=0;j<m;j++){ coordonneesInitiales[i][j]= temp[i][j]; if(typeAffichage==0) System.out.print(coordonneesInitiales[i][j] + " "); } if(typeAffichage==0) System.out.println(); } break; case 1 : // On crée une matrice triangulaire aléatoire for (int i = 0; i < n; i++) { if(typeAffichage==0) System.out.print(i + "=> "); for (int j = 0; j < i; j++) { coordonneesInitiales[i][j] = (Math.random()*0.1); if(typeAffichage==0) System.out.print(coordonneesInitiales[i][j] + " "); } if(typeAffichage==0) System.out.println(); } break; } } public void executeTest() { creationMatrice(typeMat); // On calcul les distances euclidiennes entre les points if(typeAffichage==0) System.out.println("Matrice distance"); for (int i = 0; i < n; i++) { //System.out.print(i + "=> "); for (int j = 0; j < m; j++) { distances[i][j] = dist(coordonneesInitiales, i, j); distances[j][i] = dist(coordonneesInitiales, i, j); if(typeAffichage==0) System.out.print(distances[i][j] + " "); } if(typeAffichage==0) System.out.println(); } if(typeAffichage==0) System.out.println("Matrice Résultat (intermediaire)"); if(typeAffichage==0) System.out.print("1=> "); // Le premier point a tout ces coordonnées nulles // Le deuxième n'a qu'une seule coordonnée nulle coordonnees[1][0] = distances[1][0]; if(typeAffichage==0) System.out.println(coordonneesInitiales[1][0] + " "); /* On lance notre algorithme pour trouver les coordonnées des points 3 à n */ double zero = 0; for (int WichPoint = 2; WichPoint < n; WichPoint++) { int i = WichPoint; // je met i juste pour simplifier le code, et Wichpoint pour comprendre de // quoi je parle if(typeAffichage==0) System.out.print(i + "=> "); // La norme du vecteur sera sa distance au point 1 double norme = carre(distances[0][i]); // On calcul d'abord les n-2 premières coordonnées du point courant for (int p = 0; p < i - 1; p++) { coordonnees[i][p] = norme - carre(distances[p + 1][i]) + carre(distances[p + 1][0]); if(typeAffichage==0) System.out.print("p=" + p + " norme-d²=" + coordonnees[i][p]); double sumCoord = 0; for (int k = 0; k <= p - 1; k++) { sumCoord += coordonnees[i][k] * coordonnees[p + 1][k]; } if(typeAffichage==0) System.out.print(" sumCoord=" + sumCoord + " "); coordonnees[i][p] = coordonnees[i][p] - 2 * sumCoord; coordonnees[i][p] = coordonnees[i][p] / (2 * coordonnees[p + 1][p]); // on controle la divergence en soustrayant ou additionnant une erreur qui pourrait apparaitre if(typeAffichage==0) System.out.print(" x=" + coordonnees[i][p] + " "); if(typeAffichage==0) System.out.println(); } if(typeAffichage==0) System.out.println(); // On calcul enfin la dernière coord du point courant en utilisant les coord calculés à // l'instant double sumSQR = 0; for (int j = 0; j < i - 1; j++) { sumSQR += carre(coordonnees[i][j]); } // Le dernier calcul devrait tjs etre vrai car norme > sumSQR (mathematiquement j'entend) coordonnees[i][i - 1] = Math.sqrt(norme - sumSQR); if(i%2==0) coordonnees[i][i - 1]=coordonnees[i][i - 1]-Math.pow(10,-17); else coordonnees[i][i - 1]=coordonnees[i][i - 1]+Math.pow(10,-17); // on test si un des résultats s'approche de zero = controle division par zero if(Math.abs(coordonnees[i][i-1])<0.0000001){ System.out.println(" petit resultat= " + coordonnees[i][i-1]); zeroPos[0]=i; zeroPos[1]=i-1; zeroPos[2]=coordonnees[i][i-1]; } // Or y'a du NaN ???? donc norme < sumSQR et c pas normal !!! // On test le dernier calcul voir quand apparait ce NaN if (Double.isNaN(coordonnees[i][i - 1]) == true && zero == 0) { zero++; quelLigne=i-1; System.out.println("Nan a " + i); System.out.println(" Position du zero : i=" + zeroPos[0] + " j=" + zeroPos[1] + " coord="+ zeroPos[2]); if(typeAffichage==0) typeAffichage=2; // on chg la valeur pour arreter l'affichage au NaN } } if(typeAffichage==1) System.out.println("Matrice resultat"); for (int i = 0; i < n; i++) { if(typeAffichage==1) System.out.print(i + "=> "); for (int j = 0; j < i; j++) { if(typeAffichage==1) System.out.print(coordonnees[i][j] + " "); } if(typeAffichage==1) System.out.println(); } // sommme des distances if(typeAffichage==3){ for(int i=0;i<quelLigne;i++){ double sumInit=0; double sumFinale=0; System.out.print(i + "=> "); for(int j=0;j<i;j++){ sumInit+=distances[i][j]; sumFinale+=dist(coordonnees,i,j); } if(typeAffichage==3) System.out.println(" erreur somme des distances= " + (sumInit-sumFinale)); } } // difference de coordonnées if(typeAffichage==4){ for(int j=0;j<quelLigne-1;j++){ System.out.print("j="+j+" erreur="+ (coordonneesInitiales[quelLigne-1][j]-coordonnees[quelLigne-1][j])); System.out.println(); } } } public static void main(String[] args) { Geotest test = new Geotest(); test.executeTest(); } }change maintenant typeAffichage=4 en typeAffichage=3 et tu verra que les distances ne divergent pas.
Ah :!: :lol:
Y a qu'un parisien qui se croit au centre du monde pour penser que c'est évident pour tout monde.
Moi quand j'ai vu ReR (et non RER), dans le contexte de la conversation, je me suis demandé quel coin obscur de l'informatique j'avais l'infortune d'ignorer. :lol:
Ah... la région parisienne... quand je pense que j'ai habité cet asile de dingues pendant 12 longues années. Mais je me suis évadé! Pour me réfugier dans un trou du sud de la France on on sait à peine ce que c'est qu'une voiture :) Alors le RER...
Le RER... c'était dans une autre vie pour moi ça :)
Pour ton problème je crois que j'ai compris au moins l'idée dans les grandes lignes. Mais ce que je ne comprends pas c'est pourquoi le calcul diverge. Si à l'occasion tu a sle temps d'écrire deux ou trois lignes bien vulgarisées à l'attention du primate que je suis :)
j'ai fait expres de repondre : Ben le RER quoi ! je sentais que ca allait te faire grommeler !
sinon pour l'explication, je te fait des que j'ai un peu de temps, la je suis en pleine bourre post decouverte :D
Je vois :lol:
Pages