Maths:Racine Nième [Résolu]
Taur33
Membre Dernière intervention
20 mai 2011
- - Dernière réponse :
BruNews
Modérateur Dernière intervention
7 novembre 2014
- 12 déc. 2010 à 10:56
- Messages postés
- 85
- Date d'inscription
- vendredi 24 septembre 2010
- Statut
- Messages postés
- 21054
- Date d'inscription
- jeudi 23 janvier 2003
- Statut
A voir également:
- Racine quatrième
- Limite de la racine nième - Meilleures réponses
- Limite racine nième de n - Meilleures réponses
- Racine carrée par la méthode de newton - Codes sources - Visual Basic / VB.NET (Maths)
- Trouver une racine treeees precise - Codes sources - Visual Basic / VB.NET (Maths)
- Visual Basic / VB.NET : Racine carré sans sqrt - CodeS SourceS - Guide
- Calcul des racines d ' une équation du second degré - Codes sources - Visual Basic / VB.NET (Maths)
- Graphique de fonction du quatrieme° - Codes sources - Visual Basic / VB.NET (Maths)
5 réponses
- Messages postés
- 417
- Date d'inscription
- vendredi 22 février 2008
- Statut
Bonsoir Taur33,
Je n'ai pas de réponse mais des liens:
http://f-cpu.seul.org/new/manuel_arith.pdf
c'est assez théorique et compliqué mais si tu cherche une méthode sans division elle est peut être sur ce tuto.(methode CMOS ... proche de celle de ton PC)
Mais il existe d'autres méthodes:
http://membres.multimania.fr/ericmer/Racines/racines.htm
Sinon bonne prog CGSI3
Je n'ai pas de réponse mais des liens:
http://f-cpu.seul.org/new/manuel_arith.pdf
c'est assez théorique et compliqué mais si tu cherche une méthode sans division elle est peut être sur ce tuto.(methode CMOS ... proche de celle de ton PC)
Mais il existe d'autres méthodes:
http://membres.multimania.fr/ericmer/Racines/racines.htm
Sinon bonne prog CGSI3
- Messages postés
- 21054
- Date d'inscription
- jeudi 23 janvier 2003
- Statut
Tout ceci est de l'ordre de la masturbation matheuse, ce qui est une pure nuisance en informatique.
Exemple:
"Si cela ne suffit toujours pas à rendre R(N - 1) supérieur à T, on supprime la ligne du +101, on remultiplie de nouveau les R(N) par pow(10, N), etc..."
Si c'est pour faire des puissances à répétition, je n'en fais qu'une et 1 SEULE:
R = pow(N, 1/degre);
ciao...
BruNews, MVP VC++
Exemple:
"Si cela ne suffit toujours pas à rendre R(N - 1) supérieur à T, on supprime la ligne du +101, on remultiplie de nouveau les R(N) par pow(10, N), etc..."
Si c'est pour faire des puissances à répétition, je n'en fais qu'une et 1 SEULE:
R = pow(N, 1/degre);
ciao...
BruNews, MVP VC++
- Messages postés
- 85
- Date d'inscription
- vendredi 24 septembre 2010
- Statut
Merci de votre aide CGSI3 ET BRUNEWS
le premier est en effet interessant quand au second il utilise des divisions ou
alors ce n'est que pour la racine carrée.
La méthode fonctionne avec des tranches de Racine nième donc n (le nombre de chiffres max par tranche) ne dépasse pas un Integer(sauf si on cherche Racine 3 000 000 000 de 157894 par exemple)
Peut etre mais dans ce cas pourquoi trouve-t-on une méthode de calcul validée par Wikipédia qui pourrait ne pas etre valide ?
Certes racine nième de x revient à faire x^1/n comme tu le dis BruNews mais là faut
gérer une division ou créer une fonction pow(x,n) qui puisse gérer des nombres de + de 100 chiffres et des exposants décimaux.(si tu as un bout de code je suis preneur,j'ai un algo rapide mais avec x et n entiers)
Cordialement.
Je n'ai pas de réponse mais des liens:
http://f-cpu.seul.org/new/manuel_arith.pdf
Mais il existe d'autres méthodes:
http://membres.multimania.fr/ericmer/Racines/racines.htm
le premier est en effet interessant quand au second il utilise des divisions ou
alors ce n'est que pour la racine carrée.
BigInteger: Entier signé très grand (sans limite supérieure ou inférieure) (VB2010)
Je crois que tu as un pb de limite de calcul et j'ai cette idée.
La méthode fonctionne avec des tranches de Racine nième donc n (le nombre de chiffres max par tranche) ne dépasse pas un Integer(sauf si on cherche Racine 3 000 000 000 de 157894 par exemple)
Tout ceci est de l'ordre de la masturbation matheuse, ce qui est une pure nuisance en informatique.
Si c'est pour faire des puissances à répétition, je n'en fais qu'une et 1 SEULE:
R = pow(N, 1/degre);
Peut etre mais dans ce cas pourquoi trouve-t-on une méthode de calcul validée par Wikipédia qui pourrait ne pas etre valide ?
Certes racine nième de x revient à faire x^1/n comme tu le dis BruNews mais là faut
gérer une division ou créer une fonction pow(x,n) qui puisse gérer des nombres de + de 100 chiffres et des exposants décimaux.(si tu as un bout de code je suis preneur,j'ai un algo rapide mais avec x et n entiers)
Cordialement.
- Messages postés
- 21054
- Date d'inscription
- jeudi 23 janvier 2003
- Statut
Les traitements lourds de type pow() et/ou nombre de taille indéterminée ne se font qu'en langage compilé (C et ASM).
Exemple d'une fonction pow() qui prend float32 en base et exponentielle en mode x64.
Cette pow() est très rapide mais APPROCHEE à quelques millièmes prés.
Pour du float64 ce sera nettement plus compliqué, si taille indéterminée alors...
bnPowF PROC ; XMM0 float fVal, XMM1 float fexp
lea rdx, [rsp - 24]
xorps xmm5, xmm5
movd xmm2, dword ptr BNPACKinv_mant_mask
movaps [rdx], xmm7
cmpss xmm5, xmm0, 1
maxss xmm0, dword ptr BNPACKmin_norm_pos
movd xmm7, dword ptr BNPACK1PZ
movaps xmm3, xmm0
movaps [rdx - 16], xmm6
andps xmm0, xmm2
orps xmm0, xmm7
movd eax, xmm5
movaps xmm4, xmm0
movd xmm2, dword ptr BNPACKDW7F
subss xmm0, xmm7
addss xmm4, xmm7
psrld xmm3, 23
rcpss xmm4, xmm4
mulss xmm0, xmm4
psubd xmm3, xmm2
addss xmm0, xmm0
movaps xmm2, xmm0
mulss xmm0, xmm0
movd xmm4, dword ptr BNPACKlog_p0
movd xmm6, dword ptr BNPACKlog_q0
mulss xmm4, xmm0
movss xmm5, dword ptr BNPACKlog_p1
mulss xmm6, xmm0
movd xmm7, dword ptr BNPACKlog_q1
addss xmm4, xmm5
addss xmm6, xmm7
movd xmm5, dword ptr BNPACKlog_p2
mulss xmm4, xmm0
movd xmm7, dword ptr BNPACKlog_q2
mulss xmm6, xmm0
addss xmm4, xmm5
movd xmm5, dword ptr BNPACKlog2_c0
addss xmm6, xmm7
cvtdq2ps xmm7, xmm3
mulss xmm0, xmm4
rcpss xmm6, xmm6
mulss xmm0, xmm6
movd xmm4, dword ptr BNPACKexp2_hi
mulss xmm0, xmm2
movd xmm6, dword ptr BNPACKexp2_lo
mulss xmm2, xmm5
mulss xmm0, xmm5
addss xmm2, xmm7
movd xmm3, dword ptr BNPACKZP5
addss xmm0, xmm2
xorps xmm2, xmm2
movd xmm5, dword ptr BNPACKDW1
mulss xmm0, xmm1
minss xmm0, xmm4
movd xmm4, dword ptr BNPACKexp2_p0
maxss xmm0, xmm6
movd xmm6, dword ptr BNPACKexp2_q0
addss xmm3, xmm0
cmpnltss xmm2, xmm3
pand xmm2, xmm5
cvttps2dq xmm3, xmm3
psubd xmm3, xmm2
cvtdq2ps xmm2, xmm3
subss xmm0, xmm2
movaps xmm2, xmm0
mulss xmm0, xmm0
paddd xmm3, BNPACKDW7F
mulss xmm4, xmm0
mulss xmm6, xmm0
addss xmm4, dword ptr BNPACKexp2_p1
addss xmm6, dword ptr BNPACKexp2_q1
mulss xmm4, xmm0
addss xmm4, dword ptr BNPACKexp2_p2
mulss xmm2, xmm4
movd xmm0, dword ptr BNPACK1PZ
subss xmm6, xmm2
pslld xmm3, 23
rcpss xmm6, xmm6
movd xmm5, eax
mulss xmm2, xmm6
andps xmm3, xmm5
movaps xmm7, [rdx]
addss xmm2, xmm2
movaps xmm6, [rdx - 16]
addss xmm0, xmm2
mulss xmm0, xmm3
ret 0
bnPowF ENDP
ciao...
BruNews, MVP VC++
Exemple d'une fonction pow() qui prend float32 en base et exponentielle en mode x64.
Cette pow() est très rapide mais APPROCHEE à quelques millièmes prés.
Pour du float64 ce sera nettement plus compliqué, si taille indéterminée alors...
bnPowF PROC ; XMM0 float fVal, XMM1 float fexp
lea rdx, [rsp - 24]
xorps xmm5, xmm5
movd xmm2, dword ptr BNPACKinv_mant_mask
movaps [rdx], xmm7
cmpss xmm5, xmm0, 1
maxss xmm0, dword ptr BNPACKmin_norm_pos
movd xmm7, dword ptr BNPACK1PZ
movaps xmm3, xmm0
movaps [rdx - 16], xmm6
andps xmm0, xmm2
orps xmm0, xmm7
movd eax, xmm5
movaps xmm4, xmm0
movd xmm2, dword ptr BNPACKDW7F
subss xmm0, xmm7
addss xmm4, xmm7
psrld xmm3, 23
rcpss xmm4, xmm4
mulss xmm0, xmm4
psubd xmm3, xmm2
addss xmm0, xmm0
movaps xmm2, xmm0
mulss xmm0, xmm0
movd xmm4, dword ptr BNPACKlog_p0
movd xmm6, dword ptr BNPACKlog_q0
mulss xmm4, xmm0
movss xmm5, dword ptr BNPACKlog_p1
mulss xmm6, xmm0
movd xmm7, dword ptr BNPACKlog_q1
addss xmm4, xmm5
addss xmm6, xmm7
movd xmm5, dword ptr BNPACKlog_p2
mulss xmm4, xmm0
movd xmm7, dword ptr BNPACKlog_q2
mulss xmm6, xmm0
addss xmm4, xmm5
movd xmm5, dword ptr BNPACKlog2_c0
addss xmm6, xmm7
cvtdq2ps xmm7, xmm3
mulss xmm0, xmm4
rcpss xmm6, xmm6
mulss xmm0, xmm6
movd xmm4, dword ptr BNPACKexp2_hi
mulss xmm0, xmm2
movd xmm6, dword ptr BNPACKexp2_lo
mulss xmm2, xmm5
mulss xmm0, xmm5
addss xmm2, xmm7
movd xmm3, dword ptr BNPACKZP5
addss xmm0, xmm2
xorps xmm2, xmm2
movd xmm5, dword ptr BNPACKDW1
mulss xmm0, xmm1
minss xmm0, xmm4
movd xmm4, dword ptr BNPACKexp2_p0
maxss xmm0, xmm6
movd xmm6, dword ptr BNPACKexp2_q0
addss xmm3, xmm0
cmpnltss xmm2, xmm3
pand xmm2, xmm5
cvttps2dq xmm3, xmm3
psubd xmm3, xmm2
cvtdq2ps xmm2, xmm3
subss xmm0, xmm2
movaps xmm2, xmm0
mulss xmm0, xmm0
paddd xmm3, BNPACKDW7F
mulss xmm4, xmm0
mulss xmm6, xmm0
addss xmm4, dword ptr BNPACKexp2_p1
addss xmm6, dword ptr BNPACKexp2_q1
mulss xmm4, xmm0
addss xmm4, dword ptr BNPACKexp2_p2
mulss xmm2, xmm4
movd xmm0, dword ptr BNPACK1PZ
subss xmm6, xmm2
pslld xmm3, 23
rcpss xmm6, xmm6
movd xmm5, eax
mulss xmm2, xmm6
andps xmm3, xmm5
movaps xmm7, [rdx]
addss xmm2, xmm2
movaps xmm6, [rdx - 16]
addss xmm0, xmm2
mulss xmm0, xmm3
ret 0
bnPowF ENDP
ciao...
BruNews, MVP VC++