ASREML does not use Moore-Penrose inverses but simply zeroes
singular rows. I do not have any Moore Penrose inverse code
Arthur
On Wed, 5 Dec 2001, crepieux wrote:
> Dear ASREMLers,
> sorry to ask a question not really about ASREML (but the step before).
> I have troubles to find a good (an efficient) C program to compute
> moore-penrose inverses for (about) 1000*1000 square matrices.
> Does someone have the code of one of them, or an internet link to find one?
> Thanks for your help
> Best regards,
>
> Sebastien Crepieux
>
>
> Sébastien CREPIEUX, Ph.D. student
> INRA domaine de Crouelle
> Unité d'amélioration et santé des plantes
> 234, Avenue du Brezet
> 63 039 Clermont-Ferrand Cedex
> France
> Tel: (0033) 4 73 62 43 90
>
> --
> Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml
>
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
Arthur Gilmour PhD mailto:Arthur.Gilmour@agric.nsw.gov.au
Principal Research Scientist (Biometrics) fax: <61> 2 6391 3899
NSW Agriculture <61> 2 6391 3922
Orange Agricultural Institute telephone work: <61> 2 6391 3815
Forest Rd, ORANGE, 2800, AUSTRALIA home: <61> 2 6364 3288
mobile: <61> 4 1964 7675
ASREML is available by anonymous ftp from pub/aar on ftp.res.bbsrc.ac.uk
or point your web browser at ftp://ftp.res.bbsrc.ac.uk/pub/aar/
To join the asreml discussion list, send the message
subscribe
mailto:asreml-request@chiswick.anprod.CSIRO.au
To send messages to the list, mailto:asreml@chiswick.anprod.CSIRO.au
Asreml list archive: http://www.chiswick.anprod.csiro.au/lists/asreml
<> <> <> <> <> <> <>
"Blessed are the peacemakers,
for they shall be called sons of God" Jesus; Matthew 5:9
<><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><><>
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml