Probleme sur un programme en C
OverdOzed
Inscrit:
31/12/2011 23:13
Post(s): 319
Salut à tous blenderien, je sais que ce sujet n'a pas trop sa place sur un forum en rapport avec la 3d mais les autres forums sur lesquels j'ai posté et également des gens à qui j'ai demandé ne savent pas me répondre, et comme j'ai vu que sur le clan il y a des codeurs je me suis dit pourquoi pas !
Si personne ne sait me répondre pourriez vous me renvoyer vers des forums plus expert.
Bon alors voilà le problème pou m’entraîner au langage C j'ai décidé il n'y a pas longtemps de coder un programme de calcul de racine carré par dichotomie, vous allez me dire " hé mais ça sert à rien tu peux faire sqrt bla bla " Le truc c'est que c'est pour m’entraîner pas pour calculer des racines carrées, Le programme fonctionne a deux détail près:
1. La racine carré de 36 est 5.997799 donc il tombe presque juste mais pas totalement
2. Quand on demande une précision trop élevé le programme s'arrète
Voila le code

Citation :

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

int main()
{
int n = 0;
float a = 0, b = 0, r = 0, eps = 0, delta = 0, m = 0;
printf("Entrez le nombre dont vous voulez la racine carré:\n");
scanf("%d", &n);
printf("Entrez la première borne:\n");
scanf("%f", &a);
printf("Entrez la deuxiéme borne:\n");
scanf("%f", &b);
printf("Entrez le degré de précision sous la forme 0.000 ... 1:\n");
scanf("%f", &eps);
while (b - a > eps)
{
delta = b - a;
m = delta / 2;
r = a + m;
if (r * r > n)
{
b = b - m;
}
else
{
a = a + m;
}
}
printf("La racine carré de %d est %f\n", n, r);
printf("Info complémentaire: La valeur de la première borne est maintenant %f, celle de la deuxième borne est de %f, celle de delta de %f, celle de m %f et pour epsilon %f\n", a, b, delta, m, eps);

return 0;
}

je n'ai pas pris la peine de le commenter car je le trouve relativement simple, cependant si vous voulez des commentaires dites-le.
Merci d'avance.

edit;
: je ne suis pas sur d'avoir accès à Internet les jours qui vont suivre donc ne vous étonné pas si je ne répond pas .

Contribution le : 08/08/2012 15:59
_________________
Mon site Internet
Créer un fichier PDF de la contribution Imprimer


Re: Probleme sur un programme en C ( Calcul de racine carré )
OverdOzed
Inscrit:
29/04/2007 19:11
De Nîmes...300 jours de soleil par an, inondations le reste du temps
Post(s): 3774
Pour 1/, c'est normal : tu as un résultat approché qui converge vers le résultat réel. Donc, si tu n'itères pas assez, tu arrêtes le calcul avant d'avoir atteint une précision acceptable. Le calcul par dichotomie ne tombe jamais juste, sauf si par hasard tu as choisi des bornes qui tombent bien (avec 4 et 12, ça tomberait pile, par exemple).

Pour 2/, que veut dire "une précision trop élevée" ? Combien ? Et "le programme s'arrête" => comment ? Un message d'erreur quelconque ? Je ne vois pas ce qui peut faire planter un truc aussi simple...

Par ailleurs, le code peut se simplifier un chouilla en se passant de delta et m, qui n'apportent rien (ces variables ne servent que comme résultats intermédiaires, écrits et lus une seule fois) :
while (b - a > eps)
{
   r = (a + b)/2;
   if (r * r > n)
   {
      b = r;
   }
   else
   {
      a = r;
   }
}

Contribution le : 08/08/2012 16:52
_________________
|C'est en forgeant qu'on devient forgeron, c'est en mouchant qu'on devient moucheron et c'est en sciant que Léonard devint scie.
Créer un fichier PDF de la contribution Imprimer


Re: Probleme sur un programme en C ( Calcul de racine carré )
Touriste
Inscrit:
08/08/2012 17:26
Post(s): 1
Bonjour Aspirine62400.
Je plussoie DaWaaaaghBabal le calcul de dichotomique ne peut pas tomber juste. Tu ne peux pas vraiment régler ton problème 1/ la seule chose que tu peux faire c'est vérifier si, par un coup de bol extraordinaire tu ne tomerait pas pile sur le bon nombre en rajoutant dans ton while quelque chose comme
Citation :

else if (r*r==n)
{
a=r;
b=r;
m=0;
delta=0;
break;
}

mais ça ne règle pas beaucoup de cas...

Pour 2/ après quelque tests je pense avoir compris le problème.Mon hypothèse est qu'il provient de la précision du format float en C. Quand tu demandes une trop grande précision ta variable eps vaut en fait 0 et ta boucle while devient infinie donc le programme s'arrête.

Solutions :
Tu peux partiellement régler ce problème en définissant tes variables en tant que long double ,qui est un type plus précis que float. Cependant, si tu demande une précision encore plus grande, tu auras le même problème.

En bref tu peux atténuer le problème mais pas le résoudre complètement à moins, peut être, d'utiliser une bibliothèque spécialisée dans le calcul en précision arbitraire de ce style (à vérifier je n'ai pas testé). Mais ça parait excessif vu la simplicité de ton programme.

Enfin si tu veux t'amuser tu peux essayer de coder la même chose avec d'autre algorithmes.Il me semble que la Méthode de Newton est bien plus rapide que la méthode dichotomique pour ce genre de calculs.

Voilà j'espère t'avoir un peu aidé (et ne pas avoir dis trop de bêtises)

Contribution le : 08/08/2012 18:30
Créer un fichier PDF de la contribution Imprimer


Re: Probleme sur un programme en C ( Calcul de racine carré )
OverdOzed
Inscrit:
31/12/2011 23:13
Post(s): 319
Tout d'abord je vous remercie tous les deux pour vos réponses qui me sont très utiles, @ DaWaaaaghBabal et Vaydre, c'est bien ce que je craignais en effet, impossible de tomber juste avec le calcul par dichotomie, merci de me le confirmer, @ DaWaaaaghBabal, c'est vrai que j'aurai pu simplifier le code de cette façon, merci de me le faire remarquer, @ Vaydre, oui apparemment c'est ce qui se passe la boucle devient infinie, d'ailleurs sa en a tous les symptômes ! J'aurai du penser à utiliser un double mais comme tu le dis toi même ce n'est que repousser le problème, pour la bibliothèque ça devient un peu compliqué compte tenu de mon niveau mais c'est vrai que c'est intéressant, j'y reviendrai, quant à la méthode de Newton je ne la connaissait pas et merci de m'en avoir parlé, je vais m'y intéresser.
Encore merci à tous les deux quant a Vaydre je vois que c'est ton premier post donc bienvenue sur le
@+

Contribution le : 09/08/2012 14:58
_________________
Mon site Internet
Créer un fichier PDF de la contribution Imprimer


Re: Probleme sur un programme en C ( Calcul de racine carré )
OverdOzed
Inscrit:
23/02/2006 18:10
De Alpes-Maritimes
Post(s): 3080
Salut
C'est un problème d'algortihmique (numérique).

Comme dit, la dichotomie ne permet pas d'avoir un résultat précis dans certains cas mais c'est une méthode rapide. Pour les problèmes classiques il existe plusieurs méthodes et il faut choisir la meilleure selon ce qu'on recherche (performances vs résultat).
Ici, la méthode de Newton est efficace pour trouver une racine (carrée). Tu peux y jeter un oeil.

Il n'est pas utile d'utiliser les double parfois. La mémoire utilisée est doublée et il n'y a pas forcément de meilleurs résultats (l'expérience le montre). AMHA, GMP n'apporte rien dans ce cas, à part pour apprendre à s'en servir ce qui n'est pas inintéressant ! Si la méthode est foireuse, elle le reste. Ce n'est pas comme cela qu'il faut penser.

Ceci dit, pour le 2/ quand tu dis précision élevée, pense qu'un float, de par sa structure, ne peut aller au delà de 10-7 (il vaut mieux restreindre à 10-6 pour une meilleure précision). Donc si tu choisis quelque chose de plus petit il le considère comme étant 0 ou va savoir quoi...

Avant de lancer la boucle, fais donc un test sur la valeur rentrée par l'utilisateur... De toute façon, il faut toujours vérifier les données rentrées. C'est bien pour cela que l'interaction avec un utilisateur se fait uniquement à la fin du développement quand on est sûr que son moteur fonctionne sans bug.

Dans ce genre de programme, il faut généralement utiliser un garde fou (compteur et un nombre de tours max) qui permet de ne pas faire planter le programme en cas de forever. Tu conclues ton programme avec un test sur ce compteur : soit il est égal au nombre max et la racine n'est pas trouvée, sinon elle est approchée à la précision souhaitée.
Eventuellement l'utilisateur peut choisir le nombre d'iterations maximales (toujours en traitant cette donnée par la suite dans ton programme).

Contribution le : 10/08/2012 02:06
Créer un fichier PDF de la contribution Imprimer


Re: Probleme sur un programme en C ( Calcul de racine carré )
OverdOzed
Inscrit:
31/12/2011 23:13
Post(s): 319
Ce que tu dit Bibi09 est également intéressant, je ne connaissais pas le principe du garde-fou et c'est intéressant, je vais également mettre un test sur la valeur en entrée, bonne idée.
Encore Merci de votre aide.

Contribution le : 10/08/2012 11:40
_________________
Mon site Internet
Créer un fichier PDF de la contribution Imprimer



 Haut   Précédent   Suivant




Enregistrer votre réponse
Compte*
Nom   Mot de passe   Authentification
Message:*



[Recherche avancée]



Sujets récemment répondus
Forums Sujets Réponses Lus Dernières contributions
Questions & Réponses [résolu] Découper un objet selon un autre 4 109 Aujourd'hui 18:51:31
mamain83 
Questions & Réponses X-Ray uniquement en mode Solid 2 46 Aujourd'hui 12:15:45
Horemheb 
Hors Sujet !! les tutos de Moonboots    [1][2][3]...[24] 239 24882 Hier 21:57:38
moonboots 
Questions & Réponses Solution rendu saccade    [1][2] 14 224 Hier 20:41:12
Elise1002 
Questions & Réponses Texture baké devient noir    [1][2] 14 126 Hier 07:58:33
moonboots 
Questions & Réponses Raccourcis clavier qui ne marchent plus v 2.93.4 0 51 26/11 19:16:31
GFC 
The Blender Clan 'tchat Conseil pour débutante 4 368 25/11 16:24:59
Ksuhma 
Moteur de jeu GameBlender et alternatives [WIP] DeadSigns FPS Unity - Version alpha disponible + discord    [1][2][3]...[68] 672 152226 25/11 11:59:06
Hook 
Questions & Réponses Récupérer la couleur en sortie de shader 8 175 25/11 09:33:29
bibi 
Questions & Réponses Peut-on entrer les coordonnées du point visé de la caméra ? 4 144 24/11 19:16:16
CBY 
Questions & Réponses [résolu] snap fixe sur les cm ou mm 2 112 24/11 17:47:40
mamain83 
Questions & Réponses debutant- engrenage en pointe 8 259 24/11 17:30:55
CBY 
The Blender Clan 'tchat ANNONCE IMPORTANTE : LE BLENDER CLAN REOUVRE !! Etat des lieux sur le présent et le futur :)    [1][2][3] 20 47083 24/11 01:55:17
jor760 
Questions & Réponses [non résolu] Fusion 360 - recherche d'un connaisseur 0 106 22/11 18:24:48
Redstar 
The Blender Clan 'tchat [abandonné] Clavier Corsair K55 RGB Pro / Raccourcis clavier    [1][2] 14 459 21/11 19:46:01
CBY 
Questions & Réponses Paramétrage des Particle System 0 116 21/11 10:05:15
Horemheb 
Questions & Réponses Gestion des Hook 4 179 21/11 09:26:17
Horemheb 
Questions & Réponses Tissu élastique    [1][2] 14 490 21/11 09:21:27
Horemheb 
Questions & Réponses csv>python>blender 8 1253 18/11 18:46:06
Cdric 
Questions & Réponses Editer un nom avec la police 'DroidSans.ttf' en hindi 0 127 18/11 09:55:16
Melodicpinpon 

Qui est en ligne
183 utilisateur(s) en ligne (dont 70 sur Forums)

Membre(s): 6
Invité(s): 177


LucileNann, WendyFrale, JuliannSae, SZKEzekiel, MeredithCa, RandalLayt, plus...
Nouveaux membres
Waylon0033 28/11/2021
ElishaBaud 28/11/2021
TerrieJelk 28/11/2021
KittyMokar 28/11/2021
JoseBess10 28/11/2021
DannyWette 28/11/2021
SwenWhitem 28/11/2021
ReinaReinh 28/11/2021
BookerGrah 28/11/2021
AllisonMcK 28/11/2021
Dernier Ajout
2020-09-24.jpg

Evènements à venir
Dec 29
Anniversaire d'ebrain
Jan 8
BUG de Lyon
Fev 15
Anniversaire de Dany
plus 246 plus d'élément(s)
 Par Mickaël Guédon [ebrain] © 2003-2021 The Blender Clan - hébergé par TuxFamily - Site déclaré à la CNIL sous le numéro 1155445