Abstract
We demonstrated that supernodal techniques were more efficient than traditional methods for factorization and inversion of a coefficient matrix of mixed model equations (MME), which are often required in residual maximum likelihood (REML). Supernodal left-looking and inverse multifrontal algorithms were employed for sparse factorization and inversion, respectively. The approximate minimum degree or multilevel nested dissection was used for ordering. A new computer package, Yet Another MME Solver (YAMS), was developed and compared with FSPAK with respect to computing time and size of temporary memory for 13 test matrices. The matrices were produced by fitting animal models to dairy data and by using simulations from sire, sire-maternal grand sire, maternal and dominance models for phenotypic data and animal model for genomic data. The order of matrices ranged from 32,840 to 1,048,872. The YAMS software factorized and inverted the matrices up to 13 and 10 times faster than FSPAK, respectively, when an appropriate ordering strategy was applied. The YAMS package required at most 282 MB and 512 MB of temporary memory for factorization and inversion, respectively. Processing time per iteration in average information REML was reduced, using YAMS. The YAMS package is freely available on request by contacting the corresponding author.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
More From: Journal of animal breeding and genetics = Zeitschrift fur Tierzuchtung und Zuchtungsbiologie
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.