[Caml-list] inversion de matrice en caml

[Caml-list] inversion de matrice en caml

[thomas Icart]

Je suis en math-spe .Pour mon tipe sur la geometrie des nombres ,j'aurai 
besoin d'un programme qui puisse  m'inverser des matrices de taille au moins 
15*15 dans un temps assez raisonnable.
Je vous remercie d'avance.
Thomas ICART.

_________________________________________________________________
Discutez en ligne avec vos amis, essayez MSN Messenger : 
http://messenger.msn.com
-------------------
To unsubscribe, mail caml-list-request@inria.fr Archives: http://caml.inria.fr
Bug reports: http://caml.inria.fr/bin/caml-bugs FAQ: http://caml.inria.fr/FAQ/
Beginner's list: http://groups.yahoo.com/group/ocaml_beginners

Re: [Caml-list] inversion de matrice en caml

[Laurent Chéno]

Le lundi 27 mai 2002, à 04:06 , thomas Icart a écrit :

> Je suis en math-spe .Pour mon tipe sur la geometrie des nombres 
> ,j'aurai besoin d'un programme qui puisse  m'inverser des matrices 
> de taille au moins 15*15 dans un temps assez raisonnable.
> Je vous remercie d'avance.
> Thomas ICART.
>

;-)   tu t'y prends un petit peu tard, non ?   ;-)

Bon, cela dit, tout dépend, comme tu le sais sans doute un peu, du 
type de tes matrices ? ont-elles quelque propriété particulière que 
ce soit ? (matrices symétriques, complexes/réelles, matrices 
creuses, ...)

Quelle précision attends-tu dans les résultats ? ordre de grandeur 
des coefficients ? précision ?

Sait-on a priori quelque chose sur le déterminant des matrices 
considérées ? sur leur conditionnement ?


Laurent


--
Laurent Chéno <mailto:laurent.cheno@ac-paris.fr>
Prof. de mathématiques en PSI* et de l'option informatique en 
MP/MP* au lycée Louis-le-Grand, Paris
http://pauillac.inria.fr/~cheno/

-------------------
To unsubscribe, mail caml-list-request@inria.fr Archives: http://caml.inria.fr
Bug reports: http://caml.inria.fr/bin/caml-bugs FAQ: http://caml.inria.fr/FAQ/
Beginner's list: http://groups.yahoo.com/group/ocaml_beginners

Re: [Caml-list] inversion de matrice en caml

[Oleg]

On Monday 27 May 2002 10:06 am, thomas Icart wrote:
> Je suis en math-spe .Pour mon tipe sur la geometrie des nombres ,j'aurai
> besoin d'un programme qui puisse  m'inverser des matrices de taille au
> moins 15*15 dans un temps assez raisonnable.
> Je vous remercie d'avance.
> Thomas ICART.
>

Typically, responsible programmers will use LAPACK ( 
http://www.netlib.org/lapack/lug/lapack_lug.html ) for all their linear 
algebra needs. LAPACK in turn uses BLAS (http://www.netlib.org/blas/).
If speed is important to you, and your hardware vendor does not supply 
hardware-specific BLAS, consider a self-optimizing version of BLAS called 
ATLAS ( http://math-atlas.sourceforge.net/ ). Up to 10x speed-ups are 
possible with ATLAS compared to generic BLAS. 

If you don't want to "roll your own" interface to LAPACK (as it is in 
FORTRAN), consider LACAML (found at 
http://www.ai.univie.ac.at/~markus/home/ocaml_sources.html ). I haven't tried 
it yet.

BTW, in most cases matrix inversion is not required. To solve A*X=B (where X 
is either a vector or a matrix) use xGESV subroutine (x = C,Z,S,D - denotes 
complex, double complex, float, double respectively). If you insist on 
finding the inverse, use xGETRI.

HTH

Oleg
-------------------
To unsubscribe, mail caml-list-request@inria.fr Archives: http://caml.inria.fr
Bug reports: http://caml.inria.fr/bin/caml-bugs FAQ: http://caml.inria.fr/FAQ/
Beginner's list: http://groups.yahoo.com/group/ocaml_beginners

Re: [Caml-list] inversion de matrice en caml

[Markus Mottl]

On Mon, 27 May 2002, Oleg wrote:
> If you don't want to "roll your own" interface to
> LAPACK (as it is in FORTRAN), consider LACAML (found at
> http://www.ai.univie.ac.at/~markus/home/ocaml_sources.html ). I haven't
> tried it yet.

The library currently only contains functions that I have really needed
so far (mostly functions for linear regression + vector operations),
but I have been very careful to write it in such a way that people can
easily add further ones by following the existing design.  If anybody
wants certain functions added, just tell me or even better: send me the
glue code :-)

Regards,
Markus Mottl

-- 
Markus Mottl                                             markus@oefai.at
Austrian Research Institute
for Artificial Intelligence                  http://www.oefai.at/~markus
-------------------
To unsubscribe, mail caml-list-request@inria.fr Archives: http://caml.inria.fr
Bug reports: http://caml.inria.fr/bin/caml-bugs FAQ: http://caml.inria.fr/FAQ/
Beginner's list: http://groups.yahoo.com/group/ocaml_beginners