Good day.
Changes I'd like to propose have been uploaded to
http://www.jsoftware.com/svn/addons/trunk/math/lapack/
I suppose they are not optimal. But they work and extend LAPACK addon
functionality. Anyway, criticism and improvement suggestions are welcome.
Summary of changes:
- some mistypings are fixed,
- some present interfaces are enhanced,
- some new interfaces are added.
Details:
- lapack.ijs
- Verbs extended to dyadic version:
idmat, ltmat, utmat, ltri, utri, sltri, sutri.
Now they accept optional diagonal number as x argument, default x
value is 0 (backward compatibility).
Exploited by: dgerqf, zgerqf.
- Verbs added to verify vector-or-matrix arguments:
isvector, ismatrixorvector, vmatrixorvector.
Exploited by: ztrtrs.ijs, dtrtrs.ijs.
- dgeqrf.ijs, zgeqrf.ijs
- cosmetic changes
- "non-square input matrix" case implemented
- dgeqrf, zgerqf
- "error" verb call fixed
- dgeqrfQ
- Early datatype branching. It should prevent wrong Q calculation which
presumes non-complex output from dgeqrf. Indeed, dgeqrf may return
either non-complex or complex datatype depending on datatype of input.
- dgesv.ijs
- cosmetic changes
- dgesv
- "error" verb call fixed
- dgetrf.ijs
- mistyping fixed
- cosmetic changes
- dgetrf
- "info" exit code checking added
- dgesvd.ijs
- cosmetic changes
- dgesvd
- "error" verb call fixed
- dgeev.ijs, zgeev.ijs
- cosmetic changes
- dgeev
- "info" exit code checking added
- dpotrf.ijs
- cosmetic changes
- dpotrf
- Early datatype check to avoid vsymposdef for complex input matrix.
- "error" verb call fixed
- new interfaces
- dtrtrs.ijs, ztrtrs.ijs
- subject to proposal: ability to accept vector RHS. Some kind of J-ism,
which supports both cases:
(? 5 1 $ 10) %. (? 5 5 $ 10)
(? 5 $ 10) %. (? 5 5 $ 10)
- dgerqf.ijs, zgerqf.ijs
- zpotrf.ijs
- dgees.ijs, zgees.ijs
- zgesv.ijs
Tested under base library 6.2.9:
9!:14 ''
j602/beta/2007-09-19/23:00
P.S. Have a Happy New Year!
--
WBR
Igor
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
