[Users-fr] question à propos de la fonction inv(.)

[Paul BIGNIER]

Bonjour Sylvain,

La "tolérance de conditionnement" est l'inverse de la racine carrée de 
la précision machine.

  Soit en termes de programmation : assert_checktrue( cond(A) < 
1/sqrt(%eps) )

Le conditionnement est calculé en norme 1 ou infinie par la fonction 
dgecon() que vous pouvez consulter ici : 
http://www.netlib.no/netlib/lapack/double/dgecon.f

Bonne journée,
Cordialement,
Paul.


On 10/09/2013 05:22 PM, sylvain benito wrote:
> Bonjour,
>
> Je me permet de vous écrire car j'aurais besoin d'un certain nombre de 
> précisions concernant l'utilisation de la fonction d'inversion 
> matricielle inv(.).
> Pour une matrice carrée A de taille n*n, lorsque la matrice est mal 
> conditionnée, l'appel de inv(A) fait apparaître un message d'erreur 
> suivant: "Attention: la matrice est presque singulière ou mal 
> conditionnée. rcond=...".
> J'imagine que ce message présuppose que le conditionnement de A est 
> comparé à une certaine tolerance fixée (que je noterai tol, avec 
> 0<=tol<=1): rcond(A) est évalué et si rcond(A)<=tol alors le message 
> s'affiche.
>
> Ma question est la suivante:
> 1) je voudrais savoir quelle est la valeur de tolerance en question;
> 2) disposez-vous d'une valeur de tolérance pour le conditionnement en 
> norme l2 ()?
> 3) sur quelles bases avez-vous décidé de fixer la tolérance à une 
> valeur plutôt qu'une autre? En particulier, vous appuyez-vous sur des 
> références scientifiques ou des travaux particuliers?
>
> D'avance merci pour votre éclairage,
>
> Bien scientifiquement,
> Sylvain Bénito
>
>
> _______________________________________________
> users-fr mailing list
> users-fr at lists.scilab.org
> http://lists.scilab.org/mailman/listinfo/users-fr


-- 
Paul BIGNIER
Scilab Engineer & Xcos Developer
-----------------------------------------------------------
Scilab Enterprises
143bis rue Yves Le Coz - 78000 Versailles, France
Phone: +33.1.80.77.04.69
http://www.scilab-enterprises.com

-------------- section suivante --------------
Une pièce jointe HTML a été nettoyée...
URL: <https://lists.scilab.org/pipermail/users-fr/attachments/20131011/102ad4cb/attachment.htm>