Discussion :

[souviron34]

Bonjour à tous

J'avais déposé il y a longtemps dans la rubrique Contribuez du forum Algorithme un code Fortran (IV !!!) de fit d'ellipses.

J'ai eu besoin de traduire en C les fonctions de moindre-carrés.

Pourriez-vous svp vérifier que mon portage est correct ? (seule modif réelle : profiter de la souplesse du C pour généraliser à N dimensions)

Merci

Code Fortran ici (première partie du code)

Code C ici

[dva2tlse]

Bonjour ou bonsoir Souviron,
je suis plus physicien que matheux (ingénieur) donc je n'ai pas plongé dans le code; je connais un peu le fortran et "Hello Word !" en C sans vraiment plus, mais il existe probablement des bibliothèques de fonctions mathématiques qui font ça bien. (Il est vrai qu'il faut les trouver et apprendre à les utiliser, mais moi c'est ce que je ferais)
Bon courage
À+,
David

[souviron34]

Problème résolu.. voir lien fin de la discussion citée pour le code C (ou rubrique Contribuez de la section Algorithmes)

[pvincent]

Comme dva2tlse, j'aurai plutôt recours à des bibliothèques existantes, sauf si le but est d'avoir un code simple et indépendant.

Je pense par exemple que CLapack (http://www.netlib.org/clapack/) ou les variantes optimisées libres (Atlas) ou propriétaires (e.g. MKL pour Intel) ont l'avantage d'utiliser des algorithmes plus complexes (e.g. décomposition LU ou QR au lieu de l'élimination de Gauss pour résoudre un système linéaire) mais qui sont plus performants dès que le volume des données devient important.