ruby-lapack 1.3
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- data/COPYING +56 -0
- data/GPL +340 -0
- data/README.rdoc +61 -0
- data/Rakefile +115 -0
- data/dev/common.rb +9 -0
- data/dev/defs/cbbcsd +297 -0
- data/dev/defs/cbdsqr +196 -0
- data/dev/defs/cgbbrd +174 -0
- data/dev/defs/cgbcon +114 -0
- data/dev/defs/cgbequ +121 -0
- data/dev/defs/cgbequb +128 -0
- data/dev/defs/cgbrfs +182 -0
- data/dev/defs/cgbrfsx +418 -0
- data/dev/defs/cgbsv +134 -0
- data/dev/defs/cgbsvx +356 -0
- data/dev/defs/cgbsvxx +539 -0
- data/dev/defs/cgbtf2 +110 -0
- data/dev/defs/cgbtrf +109 -0
- data/dev/defs/cgbtrs +106 -0
- data/dev/defs/cgebak +95 -0
- data/dev/defs/cgebal +125 -0
- data/dev/defs/cgebd2 +162 -0
- data/dev/defs/cgebrd +178 -0
- data/dev/defs/cgecon +90 -0
- data/dev/defs/cgeequ +107 -0
- data/dev/defs/cgeequb +114 -0
- data/dev/defs/cgees +169 -0
- data/dev/defs/cgeesx +209 -0
- data/dev/defs/cgeev +154 -0
- data/dev/defs/cgeevx +263 -0
- data/dev/defs/cgegs +204 -0
- data/dev/defs/cgegv +261 -0
- data/dev/defs/cgehd2 +115 -0
- data/dev/defs/cgehrd +136 -0
- data/dev/defs/cgelq2 +89 -0
- data/dev/defs/cgelqf +113 -0
- data/dev/defs/cgels +157 -0
- data/dev/defs/cgelsd +211 -0
- data/dev/defs/cgelss +156 -0
- data/dev/defs/cgelsx +155 -0
- data/dev/defs/cgelsy +188 -0
- data/dev/defs/cgeql2 +91 -0
- data/dev/defs/cgeqlf +116 -0
- data/dev/defs/cgeqp3 +127 -0
- data/dev/defs/cgeqpf +118 -0
- data/dev/defs/cgeqr2 +89 -0
- data/dev/defs/cgeqr2p +89 -0
- data/dev/defs/cgeqrf +114 -0
- data/dev/defs/cgeqrfp +114 -0
- data/dev/defs/cgerfs +164 -0
- data/dev/defs/cgerfsx +394 -0
- data/dev/defs/cgerq2 +91 -0
- data/dev/defs/cgerqf +116 -0
- data/dev/defs/cgesc2 +87 -0
- data/dev/defs/cgesdd +203 -0
- data/dev/defs/cgesv +97 -0
- data/dev/defs/cgesvd +195 -0
- data/dev/defs/cgesvx +333 -0
- data/dev/defs/cgesvxx +519 -0
- data/dev/defs/cgetc2 +82 -0
- data/dev/defs/cgetf2 +76 -0
- data/dev/defs/cgetrf +76 -0
- data/dev/defs/cgetri +86 -0
- data/dev/defs/cgetrs +91 -0
- data/dev/defs/cggbak +119 -0
- data/dev/defs/cggbal +144 -0
- data/dev/defs/cgges +247 -0
- data/dev/defs/cggesx +311 -0
- data/dev/defs/cggev +197 -0
- data/dev/defs/cggevx +353 -0
- data/dev/defs/cggglm +160 -0
- data/dev/defs/cgghrd +171 -0
- data/dev/defs/cgglse +157 -0
- data/dev/defs/cggqrf +195 -0
- data/dev/defs/cggrqf +194 -0
- data/dev/defs/cggsvd +324 -0
- data/dev/defs/cggsvp +240 -0
- data/dev/defs/cgtcon +112 -0
- data/dev/defs/cgtrfs +193 -0
- data/dev/defs/cgtsv +97 -0
- data/dev/defs/cgtsvx +275 -0
- data/dev/defs/cgttrf +97 -0
- data/dev/defs/cgttrs +119 -0
- data/dev/defs/cgtts2 +108 -0
- data/dev/defs/chbev +124 -0
- data/dev/defs/chbevd +194 -0
- data/dev/defs/chbevx +234 -0
- data/dev/defs/chbgst +137 -0
- data/dev/defs/chbgv +163 -0
- data/dev/defs/chbgvd +231 -0
- data/dev/defs/chbgvx +265 -0
- data/dev/defs/chbtrd +132 -0
- data/dev/defs/checon +91 -0
- data/dev/defs/cheequb +87 -0
- data/dev/defs/cheev +111 -0
- data/dev/defs/cheevd +178 -0
- data/dev/defs/cheevr +325 -0
- data/dev/defs/cheevx +225 -0
- data/dev/defs/chegs2 +96 -0
- data/dev/defs/chegst +96 -0
- data/dev/defs/chegv +155 -0
- data/dev/defs/chegvd +222 -0
- data/dev/defs/chegvx +270 -0
- data/dev/defs/cherfs +170 -0
- data/dev/defs/cherfsx +380 -0
- data/dev/defs/chesv +147 -0
- data/dev/defs/chesvx +261 -0
- data/dev/defs/chesvxx +481 -0
- data/dev/defs/chetd2 +143 -0
- data/dev/defs/chetf2 +141 -0
- data/dev/defs/chetrd +162 -0
- data/dev/defs/chetrf +151 -0
- data/dev/defs/chetri +83 -0
- data/dev/defs/chetrs +90 -0
- data/dev/defs/chetrs2 +97 -0
- data/dev/defs/chfrk +138 -0
- data/dev/defs/chgeqz +258 -0
- data/dev/defs/chla_transtype +32 -0
- data/dev/defs/chpcon +85 -0
- data/dev/defs/chpev +110 -0
- data/dev/defs/chpevd +179 -0
- data/dev/defs/chpevx +204 -0
- data/dev/defs/chpgst +80 -0
- data/dev/defs/chpgv +145 -0
- data/dev/defs/chpgvd +221 -0
- data/dev/defs/chpgvx +247 -0
- data/dev/defs/chprfs +156 -0
- data/dev/defs/chpsv +134 -0
- data/dev/defs/chpsvx +248 -0
- data/dev/defs/chptrd +118 -0
- data/dev/defs/chptrf +120 -0
- data/dev/defs/chptri +77 -0
- data/dev/defs/chptrs +84 -0
- data/dev/defs/chsein +217 -0
- data/dev/defs/chseqr +263 -0
- data/dev/defs/cla_gbamv +154 -0
- data/dev/defs/cla_gbrcond_c +146 -0
- data/dev/defs/cla_gbrcond_x +140 -0
- data/dev/defs/cla_gbrfsx_extended +391 -0
- data/dev/defs/cla_gbrpvgrw +94 -0
- data/dev/defs/cla_geamv +142 -0
- data/dev/defs/cla_gercond_c +128 -0
- data/dev/defs/cla_gercond_x +122 -0
- data/dev/defs/cla_gerfsx_extended +376 -0
- data/dev/defs/cla_heamv +142 -0
- data/dev/defs/cla_hercond_c +125 -0
- data/dev/defs/cla_hercond_x +119 -0
- data/dev/defs/cla_herfsx_extended +375 -0
- data/dev/defs/cla_herpvgrw +105 -0
- data/dev/defs/cla_lin_berr +79 -0
- data/dev/defs/cla_porcond_c +116 -0
- data/dev/defs/cla_porcond_x +110 -0
- data/dev/defs/cla_porfsx_extended +366 -0
- data/dev/defs/cla_porpvgrw +85 -0
- data/dev/defs/cla_rpvgrw +76 -0
- data/dev/defs/cla_syamv +142 -0
- data/dev/defs/cla_syrcond_c +125 -0
- data/dev/defs/cla_syrcond_x +119 -0
- data/dev/defs/cla_syrfsx_extended +375 -0
- data/dev/defs/cla_syrpvgrw +105 -0
- data/dev/defs/cla_wwaddw +53 -0
- data/dev/defs/clabrd +187 -0
- data/dev/defs/clacgv +44 -0
- data/dev/defs/clacn2 +91 -0
- data/dev/defs/clacon +75 -0
- data/dev/defs/clacp2 +78 -0
- data/dev/defs/clacpy +78 -0
- data/dev/defs/clacrm +88 -0
- data/dev/defs/clacrt +73 -0
- data/dev/defs/cladiv +37 -0
- data/dev/defs/claed0 +121 -0
- data/dev/defs/claed7 +237 -0
- data/dev/defs/claed8 +211 -0
- data/dev/defs/claein +117 -0
- data/dev/defs/claesy +81 -0
- data/dev/defs/claev2 +83 -0
- data/dev/defs/clag2z +73 -0
- data/dev/defs/clags2 +116 -0
- data/dev/defs/clagtm +115 -0
- data/dev/defs/clahef +127 -0
- data/dev/defs/clahqr +159 -0
- data/dev/defs/clahr2 +153 -0
- data/dev/defs/clahrd +141 -0
- data/dev/defs/claic1 +101 -0
- data/dev/defs/clals0 +251 -0
- data/dev/defs/clalsa +267 -0
- data/dev/defs/clalsd +160 -0
- data/dev/defs/clangb +92 -0
- data/dev/defs/clange +84 -0
- data/dev/defs/clangt +77 -0
- data/dev/defs/clanhb +99 -0
- data/dev/defs/clanhe +93 -0
- data/dev/defs/clanhf +211 -0
- data/dev/defs/clanhp +85 -0
- data/dev/defs/clanhs +78 -0
- data/dev/defs/clanht +69 -0
- data/dev/defs/clansb +97 -0
- data/dev/defs/clansp +83 -0
- data/dev/defs/clansy +92 -0
- data/dev/defs/clantb +108 -0
- data/dev/defs/clantp +93 -0
- data/dev/defs/clantr +109 -0
- data/dev/defs/clapll +67 -0
- data/dev/defs/clapmr +73 -0
- data/dev/defs/clapmt +73 -0
- data/dev/defs/claqgb +124 -0
- data/dev/defs/claqge +107 -0
- data/dev/defs/claqhb +105 -0
- data/dev/defs/claqhe +98 -0
- data/dev/defs/claqhp +89 -0
- data/dev/defs/claqp2 +115 -0
- data/dev/defs/claqps +146 -0
- data/dev/defs/claqr0 +203 -0
- data/dev/defs/claqr1 +66 -0
- data/dev/defs/claqr2 +243 -0
- data/dev/defs/claqr3 +238 -0
- data/dev/defs/claqr4 +203 -0
- data/dev/defs/claqr5 +221 -0
- data/dev/defs/claqsb +105 -0
- data/dev/defs/claqsp +89 -0
- data/dev/defs/claqsy +98 -0
- data/dev/defs/clar1v +199 -0
- data/dev/defs/clar2v +88 -0
- data/dev/defs/clarcm +88 -0
- data/dev/defs/clarf +97 -0
- data/dev/defs/clarfb +139 -0
- data/dev/defs/clarfg +71 -0
- data/dev/defs/clarfgp +69 -0
- data/dev/defs/clarft +133 -0
- data/dev/defs/clarfx +88 -0
- data/dev/defs/clargv +88 -0
- data/dev/defs/clarnv +63 -0
- data/dev/defs/clarrv +259 -0
- data/dev/defs/clarscl2 +61 -0
- data/dev/defs/clartg +66 -0
- data/dev/defs/clartv +82 -0
- data/dev/defs/clarz +111 -0
- data/dev/defs/clarzb +150 -0
- data/dev/defs/clarzt +151 -0
- data/dev/defs/clascl +102 -0
- data/dev/defs/clascl2 +61 -0
- data/dev/defs/claset +77 -0
- data/dev/defs/clasr +169 -0
- data/dev/defs/classq +71 -0
- data/dev/defs/claswp +81 -0
- data/dev/defs/clasyf +127 -0
- data/dev/defs/clatbs +206 -0
- data/dev/defs/clatdf +133 -0
- data/dev/defs/clatps +193 -0
- data/dev/defs/clatrd +168 -0
- data/dev/defs/clatrs +202 -0
- data/dev/defs/clatrz +106 -0
- data/dev/defs/clatzm +124 -0
- data/dev/defs/clauu2 +68 -0
- data/dev/defs/clauum +68 -0
- data/dev/defs/cpbcon +99 -0
- data/dev/defs/cpbequ +96 -0
- data/dev/defs/cpbrfs +166 -0
- data/dev/defs/cpbstf +116 -0
- data/dev/defs/cpbsv +136 -0
- data/dev/defs/cpbsvx +314 -0
- data/dev/defs/cpbtf2 +105 -0
- data/dev/defs/cpbtrf +103 -0
- data/dev/defs/cpbtrs +98 -0
- data/dev/defs/cpftrf +183 -0
- data/dev/defs/cpftri +175 -0
- data/dev/defs/cpftrs +185 -0
- data/dev/defs/cpocon +87 -0
- data/dev/defs/cpoequ +80 -0
- data/dev/defs/cpoequb +80 -0
- data/dev/defs/cporfs +160 -0
- data/dev/defs/cporfsx +371 -0
- data/dev/defs/cposv +105 -0
- data/dev/defs/cposvx +281 -0
- data/dev/defs/cposvxx +471 -0
- data/dev/defs/cpotf2 +75 -0
- data/dev/defs/cpotrf +73 -0
- data/dev/defs/cpotri +67 -0
- data/dev/defs/cpotrs +79 -0
- data/dev/defs/cppcon +85 -0
- data/dev/defs/cppequ +82 -0
- data/dev/defs/cpprfs +146 -0
- data/dev/defs/cppsv +115 -0
- data/dev/defs/cppsvx +283 -0
- data/dev/defs/cpptrf +81 -0
- data/dev/defs/cpptri +58 -0
- data/dev/defs/cpptrs +84 -0
- data/dev/defs/cpstf2 +108 -0
- data/dev/defs/cpstrf +108 -0
- data/dev/defs/cptcon +84 -0
- data/dev/defs/cpteqr +116 -0
- data/dev/defs/cptrfs +161 -0
- data/dev/defs/cptsv +89 -0
- data/dev/defs/cptsvx +210 -0
- data/dev/defs/cpttrf +59 -0
- data/dev/defs/cpttrs +98 -0
- data/dev/defs/cptts2 +89 -0
- data/dev/defs/crot +72 -0
- data/dev/defs/cspcon +85 -0
- data/dev/defs/cspmv +121 -0
- data/dev/defs/cspr +98 -0
- data/dev/defs/csprfs +156 -0
- data/dev/defs/cspsv +134 -0
- data/dev/defs/cspsvx +248 -0
- data/dev/defs/csptrf +121 -0
- data/dev/defs/csptri +77 -0
- data/dev/defs/csptrs +84 -0
- data/dev/defs/csrscl +49 -0
- data/dev/defs/cstedc +191 -0
- data/dev/defs/cstegr +233 -0
- data/dev/defs/cstein +159 -0
- data/dev/defs/cstemr +302 -0
- data/dev/defs/csteqr +103 -0
- data/dev/defs/csycon +91 -0
- data/dev/defs/csyconv +90 -0
- data/dev/defs/csyequb +104 -0
- data/dev/defs/csymv +126 -0
- data/dev/defs/csyr +102 -0
- data/dev/defs/csyrfs +170 -0
- data/dev/defs/csyrfsx +380 -0
- data/dev/defs/csysv +149 -0
- data/dev/defs/csysvx +261 -0
- data/dev/defs/csysvxx +485 -0
- data/dev/defs/csyswapr +70 -0
- data/dev/defs/csytf2 +140 -0
- data/dev/defs/csytrf +156 -0
- data/dev/defs/csytri +83 -0
- data/dev/defs/csytri2 +110 -0
- data/dev/defs/csytri2x +90 -0
- data/dev/defs/csytrs +90 -0
- data/dev/defs/csytrs2 +97 -0
- data/dev/defs/ctbcon +109 -0
- data/dev/defs/ctbrfs +164 -0
- data/dev/defs/ctbtrs +113 -0
- data/dev/defs/ctfsm +259 -0
- data/dev/defs/ctftri +183 -0
- data/dev/defs/ctfttp +172 -0
- data/dev/defs/ctfttr +182 -0
- data/dev/defs/ctgevc +194 -0
- data/dev/defs/ctgex2 +158 -0
- data/dev/defs/ctgexc +176 -0
- data/dev/defs/ctgsen +406 -0
- data/dev/defs/ctgsja +344 -0
- data/dev/defs/ctgsna +282 -0
- data/dev/defs/ctgsy2 +235 -0
- data/dev/defs/ctgsyl +273 -0
- data/dev/defs/ctpcon +95 -0
- data/dev/defs/ctprfs +150 -0
- data/dev/defs/ctptri +79 -0
- data/dev/defs/ctptrs +98 -0
- data/dev/defs/ctpttf +172 -0
- data/dev/defs/ctpttr +73 -0
- data/dev/defs/ctrcon +103 -0
- data/dev/defs/ctrevc +188 -0
- data/dev/defs/ctrexc +99 -0
- data/dev/defs/ctrrfs +158 -0
- data/dev/defs/ctrsen +234 -0
- data/dev/defs/ctrsna +223 -0
- data/dev/defs/ctrsyl +126 -0
- data/dev/defs/ctrti2 +76 -0
- data/dev/defs/ctrtri +75 -0
- data/dev/defs/ctrtrs +107 -0
- data/dev/defs/ctrttf +181 -0
- data/dev/defs/ctrttp +72 -0
- data/dev/defs/ctzrqf +104 -0
- data/dev/defs/ctzrzf +128 -0
- data/dev/defs/cunbdb +270 -0
- data/dev/defs/cuncsd +283 -0
- data/dev/defs/cung2l +85 -0
- data/dev/defs/cung2r +85 -0
- data/dev/defs/cungbr +129 -0
- data/dev/defs/cunghr +97 -0
- data/dev/defs/cungl2 +84 -0
- data/dev/defs/cunglq +100 -0
- data/dev/defs/cungql +101 -0
- data/dev/defs/cungqr +101 -0
- data/dev/defs/cungr2 +85 -0
- data/dev/defs/cungrq +101 -0
- data/dev/defs/cungtr +95 -0
- data/dev/defs/cunm2l +130 -0
- data/dev/defs/cunm2r +130 -0
- data/dev/defs/cunmbr +179 -0
- data/dev/defs/cunmhr +157 -0
- data/dev/defs/cunml2 +130 -0
- data/dev/defs/cunmlq +143 -0
- data/dev/defs/cunmql +143 -0
- data/dev/defs/cunmqr +143 -0
- data/dev/defs/cunmr2 +130 -0
- data/dev/defs/cunmr3 +152 -0
- data/dev/defs/cunmrq +143 -0
- data/dev/defs/cunmrz +157 -0
- data/dev/defs/cunmtr +152 -0
- data/dev/defs/cupgtr +87 -0
- data/dev/defs/cupmtr +120 -0
- data/dev/defs/dbbcsd +297 -0
- data/dev/defs/dbdsdc +194 -0
- data/dev/defs/dbdsqr +203 -0
- data/dev/defs/ddisna +84 -0
- data/dev/defs/dgbbrd +167 -0
- data/dev/defs/dgbcon +114 -0
- data/dev/defs/dgbequ +121 -0
- data/dev/defs/dgbequb +128 -0
- data/dev/defs/dgbrfs +182 -0
- data/dev/defs/dgbrfsx +418 -0
- data/dev/defs/dgbsv +134 -0
- data/dev/defs/dgbsvx +353 -0
- data/dev/defs/dgbsvxx +536 -0
- data/dev/defs/dgbtf2 +110 -0
- data/dev/defs/dgbtrf +109 -0
- data/dev/defs/dgbtrs +106 -0
- data/dev/defs/dgebak +95 -0
- data/dev/defs/dgebal +125 -0
- data/dev/defs/dgebd2 +162 -0
- data/dev/defs/dgebrd +178 -0
- data/dev/defs/dgecon +90 -0
- data/dev/defs/dgeequ +107 -0
- data/dev/defs/dgeequb +114 -0
- data/dev/defs/dgees +188 -0
- data/dev/defs/dgeesx +251 -0
- data/dev/defs/dgeev +166 -0
- data/dev/defs/dgeevx +282 -0
- data/dev/defs/dgegs +207 -0
- data/dev/defs/dgegv +286 -0
- data/dev/defs/dgehd2 +115 -0
- data/dev/defs/dgehrd +136 -0
- data/dev/defs/dgejsv +862 -0
- data/dev/defs/dgelq2 +89 -0
- data/dev/defs/dgelqf +113 -0
- data/dev/defs/dgels +158 -0
- data/dev/defs/dgelsd +193 -0
- data/dev/defs/dgelss +149 -0
- data/dev/defs/dgelsx +148 -0
- data/dev/defs/dgelsy +181 -0
- data/dev/defs/dgeql2 +91 -0
- data/dev/defs/dgeqlf +116 -0
- data/dev/defs/dgeqp3 +120 -0
- data/dev/defs/dgeqpf +111 -0
- data/dev/defs/dgeqr2 +89 -0
- data/dev/defs/dgeqr2p +89 -0
- data/dev/defs/dgeqrf +114 -0
- data/dev/defs/dgeqrfp +114 -0
- data/dev/defs/dgerfs +164 -0
- data/dev/defs/dgerfsx +394 -0
- data/dev/defs/dgerq2 +91 -0
- data/dev/defs/dgerqf +116 -0
- data/dev/defs/dgesc2 +86 -0
- data/dev/defs/dgesdd +197 -0
- data/dev/defs/dgesv +97 -0
- data/dev/defs/dgesvd +188 -0
- data/dev/defs/dgesvj +308 -0
- data/dev/defs/dgesvx +333 -0
- data/dev/defs/dgesvxx +516 -0
- data/dev/defs/dgetc2 +82 -0
- data/dev/defs/dgetf2 +76 -0
- data/dev/defs/dgetrf +76 -0
- data/dev/defs/dgetri +86 -0
- data/dev/defs/dgetrs +91 -0
- data/dev/defs/dggbak +119 -0
- data/dev/defs/dggbal +145 -0
- data/dev/defs/dgges +261 -0
- data/dev/defs/dggesx +342 -0
- data/dev/defs/dggev +206 -0
- data/dev/defs/dggevx +371 -0
- data/dev/defs/dggglm +160 -0
- data/dev/defs/dgghrd +174 -0
- data/dev/defs/dgglse +157 -0
- data/dev/defs/dggqrf +195 -0
- data/dev/defs/dggrqf +194 -0
- data/dev/defs/dggsvd +319 -0
- data/dev/defs/dggsvp +233 -0
- data/dev/defs/dgsvj0 +200 -0
- data/dev/defs/dgsvj1 +220 -0
- data/dev/defs/dgtcon +119 -0
- data/dev/defs/dgtrfs +193 -0
- data/dev/defs/dgtsv +100 -0
- data/dev/defs/dgtsvx +275 -0
- data/dev/defs/dgttrf +97 -0
- data/dev/defs/dgttrs +119 -0
- data/dev/defs/dgtts2 +106 -0
- data/dev/defs/dhgeqz +282 -0
- data/dev/defs/dhsein +234 -0
- data/dev/defs/dhseqr +279 -0
- data/dev/defs/disnan +29 -0
- data/dev/defs/dla_gbamv +152 -0
- data/dev/defs/dla_gbrcond +152 -0
- data/dev/defs/dla_gbrfsx_extended +393 -0
- data/dev/defs/dla_gbrpvgrw +90 -0
- data/dev/defs/dla_geamv +142 -0
- data/dev/defs/dla_gercond +134 -0
- data/dev/defs/dla_gerfsx_extended +379 -0
- data/dev/defs/dla_lin_berr +78 -0
- data/dev/defs/dla_porcond +122 -0
- data/dev/defs/dla_porfsx_extended +366 -0
- data/dev/defs/dla_porpvgrw +81 -0
- data/dev/defs/dla_rpvgrw +72 -0
- data/dev/defs/dla_syamv +142 -0
- data/dev/defs/dla_syrcond +131 -0
- data/dev/defs/dla_syrfsx_extended +375 -0
- data/dev/defs/dla_syrpvgrw +101 -0
- data/dev/defs/dla_wwaddw +53 -0
- data/dev/defs/dlabad +43 -0
- data/dev/defs/dlabrd +186 -0
- data/dev/defs/dlacn2 +96 -0
- data/dev/defs/dlacon +80 -0
- data/dev/defs/dlacpy +78 -0
- data/dev/defs/dladiv +57 -0
- data/dev/defs/dlae2 +66 -0
- data/dev/defs/dlaebz +292 -0
- data/dev/defs/dlaed0 +147 -0
- data/dev/defs/dlaed1 +141 -0
- data/dev/defs/dlaed2 +194 -0
- data/dev/defs/dlaed3 +167 -0
- data/dev/defs/dlaed4 +113 -0
- data/dev/defs/dlaed5 +78 -0
- data/dev/defs/dlaed6 +101 -0
- data/dev/defs/dlaed7 +241 -0
- data/dev/defs/dlaed8 +227 -0
- data/dev/defs/dlaed9 +137 -0
- data/dev/defs/dlaeda +147 -0
- data/dev/defs/dlaein +141 -0
- data/dev/defs/dlaev2 +83 -0
- data/dev/defs/dlaexc +107 -0
- data/dev/defs/dlag2 +123 -0
- data/dev/defs/dlag2s +81 -0
- data/dev/defs/dlags2 +108 -0
- data/dev/defs/dlagtf +128 -0
- data/dev/defs/dlagtm +115 -0
- data/dev/defs/dlagts +135 -0
- data/dev/defs/dlagv2 +128 -0
- data/dev/defs/dlahqr +172 -0
- data/dev/defs/dlahr2 +153 -0
- data/dev/defs/dlahrd +141 -0
- data/dev/defs/dlaic1 +101 -0
- data/dev/defs/dlaln2 +187 -0
- data/dev/defs/dlals0 +250 -0
- data/dev/defs/dlalsa +268 -0
- data/dev/defs/dlalsd +152 -0
- data/dev/defs/dlamrg +66 -0
- data/dev/defs/dlaneg +87 -0
- data/dev/defs/dlangb +93 -0
- data/dev/defs/dlange +84 -0
- data/dev/defs/dlangt +77 -0
- data/dev/defs/dlanhs +78 -0
- data/dev/defs/dlansb +97 -0
- data/dev/defs/dlansf +178 -0
- data/dev/defs/dlansp +83 -0
- data/dev/defs/dlanst +69 -0
- data/dev/defs/dlansy +92 -0
- data/dev/defs/dlantb +108 -0
- data/dev/defs/dlantp +93 -0
- data/dev/defs/dlantr +109 -0
- data/dev/defs/dlanv2 +84 -0
- data/dev/defs/dlapll +69 -0
- data/dev/defs/dlapmr +73 -0
- data/dev/defs/dlapmt +73 -0
- data/dev/defs/dlapy2 +30 -0
- data/dev/defs/dlapy3 +34 -0
- data/dev/defs/dlaqgb +124 -0
- data/dev/defs/dlaqge +107 -0
- data/dev/defs/dlaqp2 +115 -0
- data/dev/defs/dlaqps +146 -0
- data/dev/defs/dlaqr0 +225 -0
- data/dev/defs/dlaqr1 +79 -0
- data/dev/defs/dlaqr2 +252 -0
- data/dev/defs/dlaqr3 +247 -0
- data/dev/defs/dlaqr4 +225 -0
- data/dev/defs/dlaqr5 +230 -0
- data/dev/defs/dlaqsb +105 -0
- data/dev/defs/dlaqsp +89 -0
- data/dev/defs/dlaqsy +98 -0
- data/dev/defs/dlaqtr +134 -0
- data/dev/defs/dlar1v +199 -0
- data/dev/defs/dlar2v +86 -0
- data/dev/defs/dlarf +93 -0
- data/dev/defs/dlarfb +139 -0
- data/dev/defs/dlarfg +71 -0
- data/dev/defs/dlarfgp +69 -0
- data/dev/defs/dlarft +133 -0
- data/dev/defs/dlarfx +89 -0
- data/dev/defs/dlargv +74 -0
- data/dev/defs/dlarnv +61 -0
- data/dev/defs/dlarra +106 -0
- data/dev/defs/dlarrb +167 -0
- data/dev/defs/dlarrc +99 -0
- data/dev/defs/dlarrd +290 -0
- data/dev/defs/dlarre +265 -0
- data/dev/defs/dlarrf +168 -0
- data/dev/defs/dlarrj +138 -0
- data/dev/defs/dlarrk +108 -0
- data/dev/defs/dlarrr +62 -0
- data/dev/defs/dlarrv +259 -0
- data/dev/defs/dlarscl2 +61 -0
- data/dev/defs/dlartg +64 -0
- data/dev/defs/dlartgp +62 -0
- data/dev/defs/dlartgs +57 -0
- data/dev/defs/dlartv +82 -0
- data/dev/defs/dlaruv +59 -0
- data/dev/defs/dlarz +109 -0
- data/dev/defs/dlarzb +150 -0
- data/dev/defs/dlarzt +151 -0
- data/dev/defs/dlas2 +71 -0
- data/dev/defs/dlascl +102 -0
- data/dev/defs/dlascl2 +61 -0
- data/dev/defs/dlasd0 +133 -0
- data/dev/defs/dlasd1 +180 -0
- data/dev/defs/dlasd2 +254 -0
- data/dev/defs/dlasd3 +213 -0
- data/dev/defs/dlasd4 +123 -0
- data/dev/defs/dlasd5 +88 -0
- data/dev/defs/dlasd6 +300 -0
- data/dev/defs/dlasd7 +264 -0
- data/dev/defs/dlasd8 +144 -0
- data/dev/defs/dlasda +276 -0
- data/dev/defs/dlasdq +187 -0
- data/dev/defs/dlasdt +76 -0
- data/dev/defs/dlaset +81 -0
- data/dev/defs/dlasq1 +75 -0
- data/dev/defs/dlasq2 +73 -0
- data/dev/defs/dlasq3 +141 -0
- data/dev/defs/dlasq4 +107 -0
- data/dev/defs/dlasq5 +100 -0
- data/dev/defs/dlasq6 +90 -0
- data/dev/defs/dlasr +169 -0
- data/dev/defs/dlasrt +53 -0
- data/dev/defs/dlassq +68 -0
- data/dev/defs/dlasv2 +100 -0
- data/dev/defs/dlaswp +81 -0
- data/dev/defs/dlasy2 +145 -0
- data/dev/defs/dlasyf +126 -0
- data/dev/defs/dlat2s +83 -0
- data/dev/defs/dlatbs +206 -0
- data/dev/defs/dlatdf +131 -0
- data/dev/defs/dlatps +192 -0
- data/dev/defs/dlatrd +168 -0
- data/dev/defs/dlatrs +202 -0
- data/dev/defs/dlatrz +106 -0
- data/dev/defs/dlatzm +123 -0
- data/dev/defs/dlauu2 +68 -0
- data/dev/defs/dlauum +68 -0
- data/dev/defs/dopgtr +87 -0
- data/dev/defs/dopmtr +120 -0
- data/dev/defs/dorbdb +270 -0
- data/dev/defs/dorcsd +264 -0
- data/dev/defs/dorg2l +85 -0
- data/dev/defs/dorg2r +85 -0
- data/dev/defs/dorgbr +129 -0
- data/dev/defs/dorghr +97 -0
- data/dev/defs/dorgl2 +84 -0
- data/dev/defs/dorglq +100 -0
- data/dev/defs/dorgql +101 -0
- data/dev/defs/dorgqr +101 -0
- data/dev/defs/dorgr2 +85 -0
- data/dev/defs/dorgrq +101 -0
- data/dev/defs/dorgtr +95 -0
- data/dev/defs/dorm2l +130 -0
- data/dev/defs/dorm2r +130 -0
- data/dev/defs/dormbr +176 -0
- data/dev/defs/dormhr +157 -0
- data/dev/defs/dorml2 +130 -0
- data/dev/defs/dormlq +143 -0
- data/dev/defs/dormql +143 -0
- data/dev/defs/dormqr +143 -0
- data/dev/defs/dormr2 +130 -0
- data/dev/defs/dormr3 +152 -0
- data/dev/defs/dormrq +143 -0
- data/dev/defs/dormrz +157 -0
- data/dev/defs/dormtr +152 -0
- data/dev/defs/dpbcon +98 -0
- data/dev/defs/dpbequ +96 -0
- data/dev/defs/dpbrfs +166 -0
- data/dev/defs/dpbstf +115 -0
- data/dev/defs/dpbsv +136 -0
- data/dev/defs/dpbsvx +314 -0
- data/dev/defs/dpbtf2 +105 -0
- data/dev/defs/dpbtrf +103 -0
- data/dev/defs/dpbtrs +98 -0
- data/dev/defs/dpftrf +161 -0
- data/dev/defs/dpftri +154 -0
- data/dev/defs/dpftrs +164 -0
- data/dev/defs/dpocon +87 -0
- data/dev/defs/dpoequ +80 -0
- data/dev/defs/dpoequb +80 -0
- data/dev/defs/dporfs +160 -0
- data/dev/defs/dporfsx +370 -0
- data/dev/defs/dposv +105 -0
- data/dev/defs/dposvx +281 -0
- data/dev/defs/dposvxx +468 -0
- data/dev/defs/dpotf2 +75 -0
- data/dev/defs/dpotrf +73 -0
- data/dev/defs/dpotri +67 -0
- data/dev/defs/dpotrs +79 -0
- data/dev/defs/dppcon +85 -0
- data/dev/defs/dppequ +82 -0
- data/dev/defs/dpprfs +146 -0
- data/dev/defs/dppsv +115 -0
- data/dev/defs/dppsvx +283 -0
- data/dev/defs/dpptrf +81 -0
- data/dev/defs/dpptri +58 -0
- data/dev/defs/dpptrs +84 -0
- data/dev/defs/dpstf2 +108 -0
- data/dev/defs/dpstrf +108 -0
- data/dev/defs/dptcon +84 -0
- data/dev/defs/dpteqr +117 -0
- data/dev/defs/dptrfs +141 -0
- data/dev/defs/dptsv +89 -0
- data/dev/defs/dptsvx +203 -0
- data/dev/defs/dpttrf +59 -0
- data/dev/defs/dpttrs +88 -0
- data/dev/defs/dptts2 +77 -0
- data/dev/defs/drscl +49 -0
- data/dev/defs/dsbev +117 -0
- data/dev/defs/dsbevd +169 -0
- data/dev/defs/dsbevx +231 -0
- data/dev/defs/dsbgst +130 -0
- data/dev/defs/dsbgv +156 -0
- data/dev/defs/dsbgvd +206 -0
- data/dev/defs/dsbgvx +259 -0
- data/dev/defs/dsbtrd +132 -0
- data/dev/defs/dsfrk +134 -0
- data/dev/defs/dsgesv +170 -0
- data/dev/defs/dspcon +92 -0
- data/dev/defs/dspev +103 -0
- data/dev/defs/dspevd +155 -0
- data/dev/defs/dspevx +197 -0
- data/dev/defs/dspgst +80 -0
- data/dev/defs/dspgv +139 -0
- data/dev/defs/dspgvd +188 -0
- data/dev/defs/dspgvx +241 -0
- data/dev/defs/dsposv +173 -0
- data/dev/defs/dsprfs +156 -0
- data/dev/defs/dspsv +134 -0
- data/dev/defs/dspsvx +249 -0
- data/dev/defs/dsptrd +118 -0
- data/dev/defs/dsptrf +120 -0
- data/dev/defs/dsptri +77 -0
- data/dev/defs/dsptrs +84 -0
- data/dev/defs/dstebz +229 -0
- data/dev/defs/dstedc +165 -0
- data/dev/defs/dstegr +233 -0
- data/dev/defs/dstein +152 -0
- data/dev/defs/dstemr +285 -0
- data/dev/defs/dsteqr +103 -0
- data/dev/defs/dsterf +54 -0
- data/dev/defs/dstev +88 -0
- data/dev/defs/dstevd +138 -0
- data/dev/defs/dstevr +273 -0
- data/dev/defs/dstevx +193 -0
- data/dev/defs/dsycon +98 -0
- data/dev/defs/dsyconv +90 -0
- data/dev/defs/dsyequb +104 -0
- data/dev/defs/dsyev +104 -0
- data/dev/defs/dsyevd +157 -0
- data/dev/defs/dsyevr +300 -0
- data/dev/defs/dsyevx +218 -0
- data/dev/defs/dsygs2 +96 -0
- data/dev/defs/dsygst +96 -0
- data/dev/defs/dsygv +148 -0
- data/dev/defs/dsygvd +197 -0
- data/dev/defs/dsygvx +263 -0
- data/dev/defs/dsyrfs +170 -0
- data/dev/defs/dsyrfsx +380 -0
- data/dev/defs/dsysv +149 -0
- data/dev/defs/dsysvx +261 -0
- data/dev/defs/dsysvxx +481 -0
- data/dev/defs/dsyswapr +70 -0
- data/dev/defs/dsytd2 +142 -0
- data/dev/defs/dsytf2 +143 -0
- data/dev/defs/dsytrd +162 -0
- data/dev/defs/dsytrf +156 -0
- data/dev/defs/dsytri +83 -0
- data/dev/defs/dsytri2 +110 -0
- data/dev/defs/dsytri2x +90 -0
- data/dev/defs/dsytrs +90 -0
- data/dev/defs/dsytrs2 +97 -0
- data/dev/defs/dtbcon +109 -0
- data/dev/defs/dtbrfs +164 -0
- data/dev/defs/dtbtrs +113 -0
- data/dev/defs/dtfsm +243 -0
- data/dev/defs/dtftri +163 -0
- data/dev/defs/dtfttp +151 -0
- data/dev/defs/dtfttr +176 -0
- data/dev/defs/dtgevc +265 -0
- data/dev/defs/dtgex2 +191 -0
- data/dev/defs/dtgexc +190 -0
- data/dev/defs/dtgsen +422 -0
- data/dev/defs/dtgsja +343 -0
- data/dev/defs/dtgsna +354 -0
- data/dev/defs/dtgsy2 +252 -0
- data/dev/defs/dtgsyl +274 -0
- data/dev/defs/dtpcon +95 -0
- data/dev/defs/dtprfs +150 -0
- data/dev/defs/dtptri +79 -0
- data/dev/defs/dtptrs +98 -0
- data/dev/defs/dtpttf +151 -0
- data/dev/defs/dtpttr +73 -0
- data/dev/defs/dtrcon +103 -0
- data/dev/defs/dtrevc +191 -0
- data/dev/defs/dtrexc +113 -0
- data/dev/defs/dtrrfs +158 -0
- data/dev/defs/dtrsen +284 -0
- data/dev/defs/dtrsna +239 -0
- data/dev/defs/dtrsyl +133 -0
- data/dev/defs/dtrti2 +76 -0
- data/dev/defs/dtrtri +75 -0
- data/dev/defs/dtrtrs +107 -0
- data/dev/defs/dtrttf +171 -0
- data/dev/defs/dtrttp +72 -0
- data/dev/defs/dtzrqf +104 -0
- data/dev/defs/dtzrzf +128 -0
- data/dev/defs/dzsum1 +49 -0
- data/dev/defs/icmax1 +55 -0
- data/dev/defs/ieeeck +49 -0
- data/dev/defs/ilaclc +46 -0
- data/dev/defs/ilaclr +46 -0
- data/dev/defs/iladiag +30 -0
- data/dev/defs/iladlc +46 -0
- data/dev/defs/iladlr +46 -0
- data/dev/defs/ilaenv +133 -0
- data/dev/defs/ilaprec +33 -0
- data/dev/defs/ilaslc +46 -0
- data/dev/defs/ilaslr +46 -0
- data/dev/defs/ilatrans +32 -0
- data/dev/defs/ilauplo +30 -0
- data/dev/defs/ilaver +35 -0
- data/dev/defs/ilazlc +46 -0
- data/dev/defs/ilazlr +46 -0
- data/dev/defs/iparmq +177 -0
- data/dev/defs/izmax1 +55 -0
- data/dev/defs/lsamen +50 -0
- data/dev/defs/sbbcsd +297 -0
- data/dev/defs/sbdsdc +192 -0
- data/dev/defs/sbdsqr +203 -0
- data/dev/defs/scsum1 +49 -0
- data/dev/defs/sdisna +84 -0
- data/dev/defs/sgbbrd +167 -0
- data/dev/defs/sgbcon +114 -0
- data/dev/defs/sgbequ +121 -0
- data/dev/defs/sgbequb +128 -0
- data/dev/defs/sgbrfs +182 -0
- data/dev/defs/sgbrfsx +419 -0
- data/dev/defs/sgbsv +134 -0
- data/dev/defs/sgbsvx +356 -0
- data/dev/defs/sgbsvxx +539 -0
- data/dev/defs/sgbtf2 +110 -0
- data/dev/defs/sgbtrf +109 -0
- data/dev/defs/sgbtrs +106 -0
- data/dev/defs/sgebak +95 -0
- data/dev/defs/sgebal +125 -0
- data/dev/defs/sgebd2 +162 -0
- data/dev/defs/sgebrd +178 -0
- data/dev/defs/sgecon +90 -0
- data/dev/defs/sgeequ +107 -0
- data/dev/defs/sgeequb +114 -0
- data/dev/defs/sgees +188 -0
- data/dev/defs/sgeesx +251 -0
- data/dev/defs/sgeev +166 -0
- data/dev/defs/sgeevx +282 -0
- data/dev/defs/sgegs +207 -0
- data/dev/defs/sgegv +286 -0
- data/dev/defs/sgehd2 +115 -0
- data/dev/defs/sgehrd +136 -0
- data/dev/defs/sgejsv +860 -0
- data/dev/defs/sgelq2 +89 -0
- data/dev/defs/sgelqf +113 -0
- data/dev/defs/sgels +158 -0
- data/dev/defs/sgelsd +194 -0
- data/dev/defs/sgelss +149 -0
- data/dev/defs/sgelsx +148 -0
- data/dev/defs/sgelsy +181 -0
- data/dev/defs/sgeql2 +91 -0
- data/dev/defs/sgeqlf +116 -0
- data/dev/defs/sgeqp3 +120 -0
- data/dev/defs/sgeqpf +111 -0
- data/dev/defs/sgeqr2 +89 -0
- data/dev/defs/sgeqr2p +89 -0
- data/dev/defs/sgeqrf +114 -0
- data/dev/defs/sgeqrfp +114 -0
- data/dev/defs/sgerfs +164 -0
- data/dev/defs/sgerfsx +394 -0
- data/dev/defs/sgerq2 +91 -0
- data/dev/defs/sgerqf +116 -0
- data/dev/defs/sgesc2 +86 -0
- data/dev/defs/sgesdd +197 -0
- data/dev/defs/sgesv +97 -0
- data/dev/defs/sgesvd +188 -0
- data/dev/defs/sgesvj +304 -0
- data/dev/defs/sgesvx +333 -0
- data/dev/defs/sgesvxx +519 -0
- data/dev/defs/sgetc2 +82 -0
- data/dev/defs/sgetf2 +76 -0
- data/dev/defs/sgetrf +76 -0
- data/dev/defs/sgetri +86 -0
- data/dev/defs/sgetrs +91 -0
- data/dev/defs/sggbak +119 -0
- data/dev/defs/sggbal +144 -0
- data/dev/defs/sgges +261 -0
- data/dev/defs/sggesx +342 -0
- data/dev/defs/sggev +206 -0
- data/dev/defs/sggevx +371 -0
- data/dev/defs/sggglm +160 -0
- data/dev/defs/sgghrd +174 -0
- data/dev/defs/sgglse +157 -0
- data/dev/defs/sggqrf +195 -0
- data/dev/defs/sggrqf +194 -0
- data/dev/defs/sggsvd +319 -0
- data/dev/defs/sggsvp +233 -0
- data/dev/defs/sgsvj0 +200 -0
- data/dev/defs/sgsvj1 +220 -0
- data/dev/defs/sgtcon +119 -0
- data/dev/defs/sgtrfs +193 -0
- data/dev/defs/sgtsv +100 -0
- data/dev/defs/sgtsvx +275 -0
- data/dev/defs/sgttrf +97 -0
- data/dev/defs/sgttrs +119 -0
- data/dev/defs/sgtts2 +106 -0
- data/dev/defs/shgeqz +282 -0
- data/dev/defs/shsein +234 -0
- data/dev/defs/shseqr +279 -0
- data/dev/defs/sisnan +29 -0
- data/dev/defs/sla_gbamv +152 -0
- data/dev/defs/sla_gbrcond +152 -0
- data/dev/defs/sla_gbrfsx_extended +391 -0
- data/dev/defs/sla_gbrpvgrw +90 -0
- data/dev/defs/sla_geamv +142 -0
- data/dev/defs/sla_gercond +134 -0
- data/dev/defs/sla_gerfsx_extended +376 -0
- data/dev/defs/sla_lin_berr +78 -0
- data/dev/defs/sla_porcond +122 -0
- data/dev/defs/sla_porfsx_extended +363 -0
- data/dev/defs/sla_porpvgrw +81 -0
- data/dev/defs/sla_rpvgrw +72 -0
- data/dev/defs/sla_syamv +142 -0
- data/dev/defs/sla_syrcond +131 -0
- data/dev/defs/sla_syrfsx_extended +372 -0
- data/dev/defs/sla_syrpvgrw +101 -0
- data/dev/defs/sla_wwaddw +53 -0
- data/dev/defs/slabad +43 -0
- data/dev/defs/slabrd +186 -0
- data/dev/defs/slacn2 +96 -0
- data/dev/defs/slacon +80 -0
- data/dev/defs/slacpy +78 -0
- data/dev/defs/sladiv +57 -0
- data/dev/defs/slae2 +66 -0
- data/dev/defs/slaebz +292 -0
- data/dev/defs/slaed0 +147 -0
- data/dev/defs/slaed1 +141 -0
- data/dev/defs/slaed2 +194 -0
- data/dev/defs/slaed3 +167 -0
- data/dev/defs/slaed4 +113 -0
- data/dev/defs/slaed5 +78 -0
- data/dev/defs/slaed6 +101 -0
- data/dev/defs/slaed7 +241 -0
- data/dev/defs/slaed8 +227 -0
- data/dev/defs/slaed9 +137 -0
- data/dev/defs/slaeda +147 -0
- data/dev/defs/slaein +141 -0
- data/dev/defs/slaev2 +83 -0
- data/dev/defs/slaexc +107 -0
- data/dev/defs/slag2 +123 -0
- data/dev/defs/slag2d +74 -0
- data/dev/defs/slags2 +108 -0
- data/dev/defs/slagtf +128 -0
- data/dev/defs/slagtm +115 -0
- data/dev/defs/slagts +135 -0
- data/dev/defs/slagv2 +128 -0
- data/dev/defs/slahqr +172 -0
- data/dev/defs/slahr2 +153 -0
- data/dev/defs/slahrd +142 -0
- data/dev/defs/slaic1 +101 -0
- data/dev/defs/slaln2 +187 -0
- data/dev/defs/slals0 +250 -0
- data/dev/defs/slalsa +268 -0
- data/dev/defs/slalsd +150 -0
- data/dev/defs/slamrg +66 -0
- data/dev/defs/slaneg +87 -0
- data/dev/defs/slangb +93 -0
- data/dev/defs/slange +84 -0
- data/dev/defs/slangt +77 -0
- data/dev/defs/slanhs +78 -0
- data/dev/defs/slansb +97 -0
- data/dev/defs/slansf +179 -0
- data/dev/defs/slansp +83 -0
- data/dev/defs/slanst +69 -0
- data/dev/defs/slansy +92 -0
- data/dev/defs/slantb +108 -0
- data/dev/defs/slantp +93 -0
- data/dev/defs/slantr +109 -0
- data/dev/defs/slanv2 +84 -0
- data/dev/defs/slapll +69 -0
- data/dev/defs/slapmr +73 -0
- data/dev/defs/slapmt +73 -0
- data/dev/defs/slapy2 +30 -0
- data/dev/defs/slapy3 +34 -0
- data/dev/defs/slaqgb +124 -0
- data/dev/defs/slaqge +107 -0
- data/dev/defs/slaqp2 +115 -0
- data/dev/defs/slaqps +146 -0
- data/dev/defs/slaqr0 +225 -0
- data/dev/defs/slaqr1 +79 -0
- data/dev/defs/slaqr2 +252 -0
- data/dev/defs/slaqr3 +247 -0
- data/dev/defs/slaqr4 +225 -0
- data/dev/defs/slaqr5 +230 -0
- data/dev/defs/slaqsb +105 -0
- data/dev/defs/slaqsp +89 -0
- data/dev/defs/slaqsy +98 -0
- data/dev/defs/slaqtr +134 -0
- data/dev/defs/slar1v +199 -0
- data/dev/defs/slar2v +86 -0
- data/dev/defs/slarf +93 -0
- data/dev/defs/slarfb +139 -0
- data/dev/defs/slarfg +71 -0
- data/dev/defs/slarfgp +69 -0
- data/dev/defs/slarft +133 -0
- data/dev/defs/slarfx +89 -0
- data/dev/defs/slargv +74 -0
- data/dev/defs/slarnv +61 -0
- data/dev/defs/slarra +106 -0
- data/dev/defs/slarrb +167 -0
- data/dev/defs/slarrc +99 -0
- data/dev/defs/slarrd +290 -0
- data/dev/defs/slarre +265 -0
- data/dev/defs/slarrf +165 -0
- data/dev/defs/slarrj +138 -0
- data/dev/defs/slarrk +108 -0
- data/dev/defs/slarrr +62 -0
- data/dev/defs/slarrv +259 -0
- data/dev/defs/slarscl2 +61 -0
- data/dev/defs/slartg +64 -0
- data/dev/defs/slartgp +62 -0
- data/dev/defs/slartgs +57 -0
- data/dev/defs/slartv +82 -0
- data/dev/defs/slaruv +59 -0
- data/dev/defs/slarz +109 -0
- data/dev/defs/slarzb +150 -0
- data/dev/defs/slarzt +151 -0
- data/dev/defs/slas2 +71 -0
- data/dev/defs/slascl +102 -0
- data/dev/defs/slascl2 +61 -0
- data/dev/defs/slasd0 +131 -0
- data/dev/defs/slasd1 +179 -0
- data/dev/defs/slasd2 +254 -0
- data/dev/defs/slasd3 +213 -0
- data/dev/defs/slasd4 +123 -0
- data/dev/defs/slasd5 +88 -0
- data/dev/defs/slasd6 +300 -0
- data/dev/defs/slasd7 +264 -0
- data/dev/defs/slasd8 +144 -0
- data/dev/defs/slasda +275 -0
- data/dev/defs/slasdq +187 -0
- data/dev/defs/slasdt +76 -0
- data/dev/defs/slaset +81 -0
- data/dev/defs/slasq1 +75 -0
- data/dev/defs/slasq2 +73 -0
- data/dev/defs/slasq3 +141 -0
- data/dev/defs/slasq4 +107 -0
- data/dev/defs/slasq5 +100 -0
- data/dev/defs/slasq6 +90 -0
- data/dev/defs/slasr +169 -0
- data/dev/defs/slasrt +53 -0
- data/dev/defs/slassq +68 -0
- data/dev/defs/slasv2 +100 -0
- data/dev/defs/slaswp +81 -0
- data/dev/defs/slasy2 +145 -0
- data/dev/defs/slasyf +126 -0
- data/dev/defs/slatbs +206 -0
- data/dev/defs/slatdf +131 -0
- data/dev/defs/slatps +192 -0
- data/dev/defs/slatrd +168 -0
- data/dev/defs/slatrs +202 -0
- data/dev/defs/slatrz +106 -0
- data/dev/defs/slatzm +123 -0
- data/dev/defs/slauu2 +68 -0
- data/dev/defs/slauum +68 -0
- data/dev/defs/sopgtr +87 -0
- data/dev/defs/sopmtr +120 -0
- data/dev/defs/sorbdb +270 -0
- data/dev/defs/sorcsd +264 -0
- data/dev/defs/sorg2l +85 -0
- data/dev/defs/sorg2r +85 -0
- data/dev/defs/sorgbr +129 -0
- data/dev/defs/sorghr +97 -0
- data/dev/defs/sorgl2 +84 -0
- data/dev/defs/sorglq +100 -0
- data/dev/defs/sorgql +101 -0
- data/dev/defs/sorgqr +101 -0
- data/dev/defs/sorgr2 +85 -0
- data/dev/defs/sorgrq +101 -0
- data/dev/defs/sorgtr +95 -0
- data/dev/defs/sorm2l +130 -0
- data/dev/defs/sorm2r +130 -0
- data/dev/defs/sormbr +176 -0
- data/dev/defs/sormhr +157 -0
- data/dev/defs/sorml2 +130 -0
- data/dev/defs/sormlq +143 -0
- data/dev/defs/sormql +143 -0
- data/dev/defs/sormqr +143 -0
- data/dev/defs/sormr2 +130 -0
- data/dev/defs/sormr3 +152 -0
- data/dev/defs/sormrq +143 -0
- data/dev/defs/sormrz +157 -0
- data/dev/defs/sormtr +152 -0
- data/dev/defs/spbcon +98 -0
- data/dev/defs/spbequ +96 -0
- data/dev/defs/spbrfs +166 -0
- data/dev/defs/spbstf +115 -0
- data/dev/defs/spbsv +136 -0
- data/dev/defs/spbsvx +314 -0
- data/dev/defs/spbtf2 +105 -0
- data/dev/defs/spbtrf +103 -0
- data/dev/defs/spbtrs +98 -0
- data/dev/defs/spftrf +161 -0
- data/dev/defs/spftri +154 -0
- data/dev/defs/spftrs +164 -0
- data/dev/defs/spocon +87 -0
- data/dev/defs/spoequ +80 -0
- data/dev/defs/spoequb +80 -0
- data/dev/defs/sporfs +160 -0
- data/dev/defs/sporfsx +370 -0
- data/dev/defs/sposv +105 -0
- data/dev/defs/sposvx +281 -0
- data/dev/defs/sposvxx +471 -0
- data/dev/defs/spotf2 +75 -0
- data/dev/defs/spotrf +73 -0
- data/dev/defs/spotri +67 -0
- data/dev/defs/spotrs +79 -0
- data/dev/defs/sppcon +85 -0
- data/dev/defs/sppequ +82 -0
- data/dev/defs/spprfs +146 -0
- data/dev/defs/sppsv +115 -0
- data/dev/defs/sppsvx +283 -0
- data/dev/defs/spptrf +81 -0
- data/dev/defs/spptri +58 -0
- data/dev/defs/spptrs +84 -0
- data/dev/defs/spstf2 +108 -0
- data/dev/defs/spstrf +108 -0
- data/dev/defs/sptcon +84 -0
- data/dev/defs/spteqr +117 -0
- data/dev/defs/sptrfs +141 -0
- data/dev/defs/sptsv +89 -0
- data/dev/defs/sptsvx +203 -0
- data/dev/defs/spttrf +59 -0
- data/dev/defs/spttrs +88 -0
- data/dev/defs/sptts2 +77 -0
- data/dev/defs/srscl +49 -0
- data/dev/defs/ssbev +117 -0
- data/dev/defs/ssbevd +169 -0
- data/dev/defs/ssbevx +231 -0
- data/dev/defs/ssbgst +130 -0
- data/dev/defs/ssbgv +156 -0
- data/dev/defs/ssbgvd +206 -0
- data/dev/defs/ssbgvx +259 -0
- data/dev/defs/ssbtrd +132 -0
- data/dev/defs/ssfrk +134 -0
- data/dev/defs/sspcon +92 -0
- data/dev/defs/sspev +103 -0
- data/dev/defs/sspevd +154 -0
- data/dev/defs/sspevx +197 -0
- data/dev/defs/sspgst +80 -0
- data/dev/defs/sspgv +139 -0
- data/dev/defs/sspgvd +188 -0
- data/dev/defs/sspgvx +241 -0
- data/dev/defs/ssprfs +156 -0
- data/dev/defs/sspsv +134 -0
- data/dev/defs/sspsvx +249 -0
- data/dev/defs/ssptrd +118 -0
- data/dev/defs/ssptrf +120 -0
- data/dev/defs/ssptri +77 -0
- data/dev/defs/ssptrs +84 -0
- data/dev/defs/sstebz +229 -0
- data/dev/defs/sstedc +164 -0
- data/dev/defs/sstegr +233 -0
- data/dev/defs/sstein +152 -0
- data/dev/defs/sstemr +285 -0
- data/dev/defs/ssteqr +103 -0
- data/dev/defs/ssterf +54 -0
- data/dev/defs/sstev +88 -0
- data/dev/defs/sstevd +138 -0
- data/dev/defs/sstevr +275 -0
- data/dev/defs/sstevx +193 -0
- data/dev/defs/ssycon +98 -0
- data/dev/defs/ssyconv +90 -0
- data/dev/defs/ssyequb +104 -0
- data/dev/defs/ssyev +104 -0
- data/dev/defs/ssyevd +157 -0
- data/dev/defs/ssyevr +302 -0
- data/dev/defs/ssyevx +218 -0
- data/dev/defs/ssygs2 +96 -0
- data/dev/defs/ssygst +96 -0
- data/dev/defs/ssygv +148 -0
- data/dev/defs/ssygvd +197 -0
- data/dev/defs/ssygvx +263 -0
- data/dev/defs/ssyrfs +170 -0
- data/dev/defs/ssyrfsx +380 -0
- data/dev/defs/ssysv +149 -0
- data/dev/defs/ssysvx +261 -0
- data/dev/defs/ssysvxx +484 -0
- data/dev/defs/ssyswapr +70 -0
- data/dev/defs/ssytd2 +142 -0
- data/dev/defs/ssytf2 +143 -0
- data/dev/defs/ssytrd +162 -0
- data/dev/defs/ssytrf +156 -0
- data/dev/defs/ssytri +83 -0
- data/dev/defs/ssytri2 +110 -0
- data/dev/defs/ssytri2x +90 -0
- data/dev/defs/ssytrs +90 -0
- data/dev/defs/ssytrs2 +97 -0
- data/dev/defs/stbcon +109 -0
- data/dev/defs/stbrfs +164 -0
- data/dev/defs/stbtrs +113 -0
- data/dev/defs/stfsm +243 -0
- data/dev/defs/stftri +163 -0
- data/dev/defs/stfttp +151 -0
- data/dev/defs/stfttr +174 -0
- data/dev/defs/stgevc +265 -0
- data/dev/defs/stgex2 +191 -0
- data/dev/defs/stgexc +190 -0
- data/dev/defs/stgsen +421 -0
- data/dev/defs/stgsja +343 -0
- data/dev/defs/stgsna +354 -0
- data/dev/defs/stgsy2 +252 -0
- data/dev/defs/stgsyl +274 -0
- data/dev/defs/stpcon +95 -0
- data/dev/defs/stprfs +150 -0
- data/dev/defs/stptri +79 -0
- data/dev/defs/stptrs +98 -0
- data/dev/defs/stpttf +151 -0
- data/dev/defs/stpttr +73 -0
- data/dev/defs/strcon +103 -0
- data/dev/defs/strevc +191 -0
- data/dev/defs/strexc +113 -0
- data/dev/defs/strrfs +158 -0
- data/dev/defs/strsen +284 -0
- data/dev/defs/strsna +239 -0
- data/dev/defs/strsyl +133 -0
- data/dev/defs/strti2 +76 -0
- data/dev/defs/strtri +75 -0
- data/dev/defs/strtrs +107 -0
- data/dev/defs/strttf +171 -0
- data/dev/defs/strttp +72 -0
- data/dev/defs/stzrqf +104 -0
- data/dev/defs/stzrzf +128 -0
- data/dev/defs/xerbla +40 -0
- data/dev/defs/xerbla_array +65 -0
- data/dev/defs/zbbcsd +297 -0
- data/dev/defs/zbdsqr +196 -0
- data/dev/defs/zcgesv +176 -0
- data/dev/defs/zcposv +182 -0
- data/dev/defs/zdrscl +49 -0
- data/dev/defs/zgbbrd +174 -0
- data/dev/defs/zgbcon +114 -0
- data/dev/defs/zgbequ +121 -0
- data/dev/defs/zgbequb +128 -0
- data/dev/defs/zgbrfs +182 -0
- data/dev/defs/zgbrfsx +418 -0
- data/dev/defs/zgbsv +134 -0
- data/dev/defs/zgbsvx +356 -0
- data/dev/defs/zgbsvxx +536 -0
- data/dev/defs/zgbtf2 +110 -0
- data/dev/defs/zgbtrf +109 -0
- data/dev/defs/zgbtrs +106 -0
- data/dev/defs/zgebak +95 -0
- data/dev/defs/zgebal +125 -0
- data/dev/defs/zgebd2 +162 -0
- data/dev/defs/zgebrd +178 -0
- data/dev/defs/zgecon +90 -0
- data/dev/defs/zgeequ +107 -0
- data/dev/defs/zgeequb +114 -0
- data/dev/defs/zgees +170 -0
- data/dev/defs/zgeesx +209 -0
- data/dev/defs/zgeev +154 -0
- data/dev/defs/zgeevx +263 -0
- data/dev/defs/zgegs +205 -0
- data/dev/defs/zgegv +261 -0
- data/dev/defs/zgehd2 +115 -0
- data/dev/defs/zgehrd +136 -0
- data/dev/defs/zgelq2 +89 -0
- data/dev/defs/zgelqf +113 -0
- data/dev/defs/zgels +157 -0
- data/dev/defs/zgelsd +211 -0
- data/dev/defs/zgelss +156 -0
- data/dev/defs/zgelsx +155 -0
- data/dev/defs/zgelsy +188 -0
- data/dev/defs/zgeql2 +91 -0
- data/dev/defs/zgeqlf +116 -0
- data/dev/defs/zgeqp3 +127 -0
- data/dev/defs/zgeqpf +118 -0
- data/dev/defs/zgeqr2 +89 -0
- data/dev/defs/zgeqr2p +89 -0
- data/dev/defs/zgeqrf +114 -0
- data/dev/defs/zgeqrfp +114 -0
- data/dev/defs/zgerfs +164 -0
- data/dev/defs/zgerfsx +394 -0
- data/dev/defs/zgerq2 +91 -0
- data/dev/defs/zgerqf +116 -0
- data/dev/defs/zgesc2 +87 -0
- data/dev/defs/zgesdd +203 -0
- data/dev/defs/zgesv +97 -0
- data/dev/defs/zgesvd +195 -0
- data/dev/defs/zgesvx +333 -0
- data/dev/defs/zgesvxx +516 -0
- data/dev/defs/zgetc2 +82 -0
- data/dev/defs/zgetf2 +76 -0
- data/dev/defs/zgetrf +76 -0
- data/dev/defs/zgetri +86 -0
- data/dev/defs/zgetrs +91 -0
- data/dev/defs/zggbak +119 -0
- data/dev/defs/zggbal +144 -0
- data/dev/defs/zgges +247 -0
- data/dev/defs/zggesx +311 -0
- data/dev/defs/zggev +197 -0
- data/dev/defs/zggevx +353 -0
- data/dev/defs/zggglm +160 -0
- data/dev/defs/zgghrd +171 -0
- data/dev/defs/zgglse +157 -0
- data/dev/defs/zggqrf +195 -0
- data/dev/defs/zggrqf +194 -0
- data/dev/defs/zggsvd +324 -0
- data/dev/defs/zggsvp +240 -0
- data/dev/defs/zgtcon +112 -0
- data/dev/defs/zgtrfs +193 -0
- data/dev/defs/zgtsv +97 -0
- data/dev/defs/zgtsvx +275 -0
- data/dev/defs/zgttrf +97 -0
- data/dev/defs/zgttrs +119 -0
- data/dev/defs/zgtts2 +108 -0
- data/dev/defs/zhbev +124 -0
- data/dev/defs/zhbevd +194 -0
- data/dev/defs/zhbevx +234 -0
- data/dev/defs/zhbgst +137 -0
- data/dev/defs/zhbgv +163 -0
- data/dev/defs/zhbgvd +231 -0
- data/dev/defs/zhbgvx +265 -0
- data/dev/defs/zhbtrd +132 -0
- data/dev/defs/zhecon +91 -0
- data/dev/defs/zheequb +87 -0
- data/dev/defs/zheev +111 -0
- data/dev/defs/zheevd +178 -0
- data/dev/defs/zheevr +325 -0
- data/dev/defs/zheevx +225 -0
- data/dev/defs/zhegs2 +96 -0
- data/dev/defs/zhegst +96 -0
- data/dev/defs/zhegv +155 -0
- data/dev/defs/zhegvd +222 -0
- data/dev/defs/zhegvx +270 -0
- data/dev/defs/zherfs +170 -0
- data/dev/defs/zherfsx +380 -0
- data/dev/defs/zhesv +147 -0
- data/dev/defs/zhesvx +261 -0
- data/dev/defs/zhesvxx +478 -0
- data/dev/defs/zhetd2 +143 -0
- data/dev/defs/zhetf2 +141 -0
- data/dev/defs/zhetrd +162 -0
- data/dev/defs/zhetrf +151 -0
- data/dev/defs/zhetri +83 -0
- data/dev/defs/zhetrs +90 -0
- data/dev/defs/zhetrs2 +97 -0
- data/dev/defs/zhfrk +138 -0
- data/dev/defs/zhgeqz +258 -0
- data/dev/defs/zhpcon +85 -0
- data/dev/defs/zhpev +110 -0
- data/dev/defs/zhpevd +180 -0
- data/dev/defs/zhpevx +204 -0
- data/dev/defs/zhpgst +80 -0
- data/dev/defs/zhpgv +145 -0
- data/dev/defs/zhpgvd +221 -0
- data/dev/defs/zhpgvx +247 -0
- data/dev/defs/zhprfs +156 -0
- data/dev/defs/zhpsv +134 -0
- data/dev/defs/zhpsvx +248 -0
- data/dev/defs/zhptrd +118 -0
- data/dev/defs/zhptrf +120 -0
- data/dev/defs/zhptri +77 -0
- data/dev/defs/zhptrs +84 -0
- data/dev/defs/zhsein +217 -0
- data/dev/defs/zhseqr +263 -0
- data/dev/defs/zla_gbamv +154 -0
- data/dev/defs/zla_gbrcond_c +146 -0
- data/dev/defs/zla_gbrcond_x +140 -0
- data/dev/defs/zla_gbrfsx_extended +391 -0
- data/dev/defs/zla_gbrpvgrw +94 -0
- data/dev/defs/zla_geamv +142 -0
- data/dev/defs/zla_gercond_c +128 -0
- data/dev/defs/zla_gercond_x +122 -0
- data/dev/defs/zla_gerfsx_extended +378 -0
- data/dev/defs/zla_heamv +142 -0
- data/dev/defs/zla_hercond_c +125 -0
- data/dev/defs/zla_hercond_x +119 -0
- data/dev/defs/zla_herfsx_extended +375 -0
- data/dev/defs/zla_herpvgrw +105 -0
- data/dev/defs/zla_lin_berr +79 -0
- data/dev/defs/zla_porcond_c +116 -0
- data/dev/defs/zla_porcond_x +110 -0
- data/dev/defs/zla_porfsx_extended +366 -0
- data/dev/defs/zla_porpvgrw +85 -0
- data/dev/defs/zla_rpvgrw +76 -0
- data/dev/defs/zla_syamv +142 -0
- data/dev/defs/zla_syrcond_c +125 -0
- data/dev/defs/zla_syrcond_x +119 -0
- data/dev/defs/zla_syrfsx_extended +375 -0
- data/dev/defs/zla_syrpvgrw +105 -0
- data/dev/defs/zla_wwaddw +53 -0
- data/dev/defs/zlabrd +187 -0
- data/dev/defs/zlacgv +44 -0
- data/dev/defs/zlacn2 +91 -0
- data/dev/defs/zlacon +75 -0
- data/dev/defs/zlacp2 +78 -0
- data/dev/defs/zlacpy +78 -0
- data/dev/defs/zlacrm +88 -0
- data/dev/defs/zlacrt +73 -0
- data/dev/defs/zladiv +37 -0
- data/dev/defs/zlaed0 +121 -0
- data/dev/defs/zlaed7 +237 -0
- data/dev/defs/zlaed8 +211 -0
- data/dev/defs/zlaein +117 -0
- data/dev/defs/zlaesy +81 -0
- data/dev/defs/zlaev2 +83 -0
- data/dev/defs/zlag2c +82 -0
- data/dev/defs/zlags2 +116 -0
- data/dev/defs/zlagtm +115 -0
- data/dev/defs/zlahef +127 -0
- data/dev/defs/zlahqr +159 -0
- data/dev/defs/zlahr2 +153 -0
- data/dev/defs/zlahrd +141 -0
- data/dev/defs/zlaic1 +101 -0
- data/dev/defs/zlals0 +251 -0
- data/dev/defs/zlalsa +267 -0
- data/dev/defs/zlalsd +160 -0
- data/dev/defs/zlangb +92 -0
- data/dev/defs/zlange +84 -0
- data/dev/defs/zlangt +77 -0
- data/dev/defs/zlanhb +99 -0
- data/dev/defs/zlanhe +93 -0
- data/dev/defs/zlanhf +211 -0
- data/dev/defs/zlanhp +85 -0
- data/dev/defs/zlanhs +78 -0
- data/dev/defs/zlanht +69 -0
- data/dev/defs/zlansb +97 -0
- data/dev/defs/zlansp +83 -0
- data/dev/defs/zlansy +92 -0
- data/dev/defs/zlantb +108 -0
- data/dev/defs/zlantp +93 -0
- data/dev/defs/zlantr +109 -0
- data/dev/defs/zlapll +67 -0
- data/dev/defs/zlapmr +73 -0
- data/dev/defs/zlapmt +73 -0
- data/dev/defs/zlaqgb +124 -0
- data/dev/defs/zlaqge +107 -0
- data/dev/defs/zlaqhb +105 -0
- data/dev/defs/zlaqhe +98 -0
- data/dev/defs/zlaqhp +89 -0
- data/dev/defs/zlaqp2 +115 -0
- data/dev/defs/zlaqps +141 -0
- data/dev/defs/zlaqr0 +203 -0
- data/dev/defs/zlaqr1 +66 -0
- data/dev/defs/zlaqr2 +243 -0
- data/dev/defs/zlaqr3 +238 -0
- data/dev/defs/zlaqr4 +203 -0
- data/dev/defs/zlaqr5 +221 -0
- data/dev/defs/zlaqsb +105 -0
- data/dev/defs/zlaqsp +89 -0
- data/dev/defs/zlaqsy +98 -0
- data/dev/defs/zlar1v +199 -0
- data/dev/defs/zlar2v +88 -0
- data/dev/defs/zlarcm +88 -0
- data/dev/defs/zlarf +97 -0
- data/dev/defs/zlarfb +139 -0
- data/dev/defs/zlarfg +71 -0
- data/dev/defs/zlarfgp +69 -0
- data/dev/defs/zlarft +133 -0
- data/dev/defs/zlarfx +88 -0
- data/dev/defs/zlargv +88 -0
- data/dev/defs/zlarnv +63 -0
- data/dev/defs/zlarrv +259 -0
- data/dev/defs/zlarscl2 +61 -0
- data/dev/defs/zlartg +66 -0
- data/dev/defs/zlartv +82 -0
- data/dev/defs/zlarz +111 -0
- data/dev/defs/zlarzb +150 -0
- data/dev/defs/zlarzt +151 -0
- data/dev/defs/zlascl +102 -0
- data/dev/defs/zlascl2 +61 -0
- data/dev/defs/zlaset +77 -0
- data/dev/defs/zlasr +169 -0
- data/dev/defs/zlassq +71 -0
- data/dev/defs/zlaswp +81 -0
- data/dev/defs/zlasyf +127 -0
- data/dev/defs/zlat2c +85 -0
- data/dev/defs/zlatbs +206 -0
- data/dev/defs/zlatdf +133 -0
- data/dev/defs/zlatps +193 -0
- data/dev/defs/zlatrd +168 -0
- data/dev/defs/zlatrs +202 -0
- data/dev/defs/zlatrz +106 -0
- data/dev/defs/zlatzm +124 -0
- data/dev/defs/zlauu2 +68 -0
- data/dev/defs/zlauum +68 -0
- data/dev/defs/zpbcon +99 -0
- data/dev/defs/zpbequ +96 -0
- data/dev/defs/zpbrfs +166 -0
- data/dev/defs/zpbstf +116 -0
- data/dev/defs/zpbsv +136 -0
- data/dev/defs/zpbsvx +314 -0
- data/dev/defs/zpbtf2 +105 -0
- data/dev/defs/zpbtrf +103 -0
- data/dev/defs/zpbtrs +98 -0
- data/dev/defs/zpftrf +182 -0
- data/dev/defs/zpftri +175 -0
- data/dev/defs/zpftrs +185 -0
- data/dev/defs/zpocon +87 -0
- data/dev/defs/zpoequ +80 -0
- data/dev/defs/zpoequb +80 -0
- data/dev/defs/zporfs +160 -0
- data/dev/defs/zporfsx +370 -0
- data/dev/defs/zposv +105 -0
- data/dev/defs/zposvx +281 -0
- data/dev/defs/zposvxx +468 -0
- data/dev/defs/zpotf2 +75 -0
- data/dev/defs/zpotrf +73 -0
- data/dev/defs/zpotri +67 -0
- data/dev/defs/zpotrs +79 -0
- data/dev/defs/zppcon +85 -0
- data/dev/defs/zppequ +82 -0
- data/dev/defs/zpprfs +146 -0
- data/dev/defs/zppsv +115 -0
- data/dev/defs/zppsvx +283 -0
- data/dev/defs/zpptrf +81 -0
- data/dev/defs/zpptri +58 -0
- data/dev/defs/zpptrs +84 -0
- data/dev/defs/zpstf2 +108 -0
- data/dev/defs/zpstrf +108 -0
- data/dev/defs/zptcon +84 -0
- data/dev/defs/zpteqr +116 -0
- data/dev/defs/zptrfs +161 -0
- data/dev/defs/zptsv +89 -0
- data/dev/defs/zptsvx +210 -0
- data/dev/defs/zpttrf +59 -0
- data/dev/defs/zpttrs +98 -0
- data/dev/defs/zptts2 +89 -0
- data/dev/defs/zrot +72 -0
- data/dev/defs/zspcon +85 -0
- data/dev/defs/zspmv +119 -0
- data/dev/defs/zspr +98 -0
- data/dev/defs/zsprfs +156 -0
- data/dev/defs/zspsv +134 -0
- data/dev/defs/zspsvx +248 -0
- data/dev/defs/zsptrf +121 -0
- data/dev/defs/zsptri +77 -0
- data/dev/defs/zsptrs +84 -0
- data/dev/defs/zstedc +192 -0
- data/dev/defs/zstegr +233 -0
- data/dev/defs/zstein +159 -0
- data/dev/defs/zstemr +302 -0
- data/dev/defs/zsteqr +103 -0
- data/dev/defs/zsycon +91 -0
- data/dev/defs/zsyconv +90 -0
- data/dev/defs/zsyequb +104 -0
- data/dev/defs/zsymv +126 -0
- data/dev/defs/zsyr +102 -0
- data/dev/defs/zsyrfs +170 -0
- data/dev/defs/zsyrfsx +380 -0
- data/dev/defs/zsysv +149 -0
- data/dev/defs/zsysvx +261 -0
- data/dev/defs/zsysvxx +482 -0
- data/dev/defs/zsyswapr +70 -0
- data/dev/defs/zsytf2 +140 -0
- data/dev/defs/zsytrf +156 -0
- data/dev/defs/zsytri +83 -0
- data/dev/defs/zsytri2 +111 -0
- data/dev/defs/zsytri2x +90 -0
- data/dev/defs/zsytrs +90 -0
- data/dev/defs/zsytrs2 +97 -0
- data/dev/defs/ztbcon +109 -0
- data/dev/defs/ztbrfs +164 -0
- data/dev/defs/ztbtrs +113 -0
- data/dev/defs/ztfsm +259 -0
- data/dev/defs/ztftri +183 -0
- data/dev/defs/ztfttp +172 -0
- data/dev/defs/ztfttr +182 -0
- data/dev/defs/ztgevc +194 -0
- data/dev/defs/ztgex2 +158 -0
- data/dev/defs/ztgexc +176 -0
- data/dev/defs/ztgsen +406 -0
- data/dev/defs/ztgsja +344 -0
- data/dev/defs/ztgsna +282 -0
- data/dev/defs/ztgsy2 +235 -0
- data/dev/defs/ztgsyl +273 -0
- data/dev/defs/ztpcon +95 -0
- data/dev/defs/ztprfs +150 -0
- data/dev/defs/ztptri +79 -0
- data/dev/defs/ztptrs +98 -0
- data/dev/defs/ztpttf +172 -0
- data/dev/defs/ztpttr +73 -0
- data/dev/defs/ztrcon +103 -0
- data/dev/defs/ztrevc +188 -0
- data/dev/defs/ztrexc +99 -0
- data/dev/defs/ztrrfs +158 -0
- data/dev/defs/ztrsen +234 -0
- data/dev/defs/ztrsna +223 -0
- data/dev/defs/ztrsyl +126 -0
- data/dev/defs/ztrti2 +76 -0
- data/dev/defs/ztrtri +75 -0
- data/dev/defs/ztrtrs +107 -0
- data/dev/defs/ztrttf +181 -0
- data/dev/defs/ztrttp +72 -0
- data/dev/defs/ztzrqf +104 -0
- data/dev/defs/ztzrzf +128 -0
- data/dev/defs/zunbdb +270 -0
- data/dev/defs/zuncsd +283 -0
- data/dev/defs/zung2l +85 -0
- data/dev/defs/zung2r +85 -0
- data/dev/defs/zungbr +129 -0
- data/dev/defs/zunghr +97 -0
- data/dev/defs/zungl2 +84 -0
- data/dev/defs/zunglq +100 -0
- data/dev/defs/zungql +101 -0
- data/dev/defs/zungqr +101 -0
- data/dev/defs/zungr2 +85 -0
- data/dev/defs/zungrq +101 -0
- data/dev/defs/zungtr +95 -0
- data/dev/defs/zunm2l +130 -0
- data/dev/defs/zunm2r +130 -0
- data/dev/defs/zunmbr +177 -0
- data/dev/defs/zunmhr +157 -0
- data/dev/defs/zunml2 +130 -0
- data/dev/defs/zunmlq +143 -0
- data/dev/defs/zunmql +143 -0
- data/dev/defs/zunmqr +143 -0
- data/dev/defs/zunmr2 +130 -0
- data/dev/defs/zunmr3 +152 -0
- data/dev/defs/zunmrq +143 -0
- data/dev/defs/zunmrz +157 -0
- data/dev/defs/zunmtr +152 -0
- data/dev/defs/zupgtr +87 -0
- data/dev/defs/zupmtr +120 -0
- data/dev/make_csrc.rb +857 -0
- data/dev/mkdoc.rb +265 -0
- data/dev/parse.rb +1972 -0
- data/doc/bd.html +16 -0
- data/doc/c.html +36 -0
- data/doc/cbd.html +161 -0
- data/doc/cgb.html +1865 -0
- data/doc/cge.html +5261 -0
- data/doc/cgg.html +2027 -0
- data/doc/cgt.html +711 -0
- data/doc/chb.html +1031 -0
- data/doc/che.html +3165 -0
- data/doc/chg.html +201 -0
- data/doc/chp.html +1696 -0
- data/doc/chs.html +386 -0
- data/doc/cpb.html +994 -0
- data/doc/cpo.html +1520 -0
- data/doc/cpp.html +770 -0
- data/doc/cpt.html +706 -0
- data/doc/csp.html +905 -0
- data/doc/cst.html +742 -0
- data/doc/csy.html +2194 -0
- data/doc/ctb.html +284 -0
- data/doc/ctg.html +1544 -0
- data/doc/ctp.html +553 -0
- data/doc/ctr.html +1281 -0
- data/doc/ctz.html +211 -0
- data/doc/cun.html +2553 -0
- data/doc/cup.html +166 -0
- data/doc/d.html +35 -0
- data/doc/dbd.html +304 -0
- data/doc/ddi.html +87 -0
- data/doc/dgb.html +1857 -0
- data/doc/dge.html +7267 -0
- data/doc/dgg.html +2102 -0
- data/doc/dgt.html +713 -0
- data/doc/dhg.html +225 -0
- data/doc/dhs.html +414 -0
- data/doc/di.html +14 -0
- data/doc/dop.html +166 -0
- data/doc/dor.html +2540 -0
- data/doc/dpb.html +992 -0
- data/doc/dpo.html +1517 -0
- data/doc/dpp.html +770 -0
- data/doc/dpt.html +675 -0
- data/doc/dsb.html +995 -0
- data/doc/dsp.html +1777 -0
- data/doc/dst.html +1422 -0
- data/doc/dsy.html +3433 -0
- data/doc/dtb.html +284 -0
- data/doc/dtg.html +1730 -0
- data/doc/dtp.html +532 -0
- data/doc/dtr.html +1346 -0
- data/doc/dtz.html +211 -0
- data/doc/gb.html +16 -0
- data/doc/ge.html +16 -0
- data/doc/gg.html +16 -0
- data/doc/gt.html +16 -0
- data/doc/hb.html +14 -0
- data/doc/he.html +14 -0
- data/doc/hg.html +16 -0
- data/doc/hp.html +14 -0
- data/doc/hs.html +16 -0
- data/doc/index.html +53 -0
- data/doc/op.html +14 -0
- data/doc/or.html +14 -0
- data/doc/others.html +1142 -0
- data/doc/pb.html +16 -0
- data/doc/po.html +16 -0
- data/doc/pp.html +16 -0
- data/doc/pt.html +16 -0
- data/doc/s.html +35 -0
- data/doc/sb.html +14 -0
- data/doc/sbd.html +303 -0
- data/doc/sdi.html +87 -0
- data/doc/sgb.html +1863 -0
- data/doc/sge.html +7263 -0
- data/doc/sgg.html +2102 -0
- data/doc/sgt.html +713 -0
- data/doc/shg.html +225 -0
- data/doc/shs.html +414 -0
- data/doc/sop.html +166 -0
- data/doc/sor.html +2540 -0
- data/doc/sp.html +16 -0
- data/doc/spb.html +992 -0
- data/doc/spo.html +1520 -0
- data/doc/spp.html +770 -0
- data/doc/spt.html +675 -0
- data/doc/ssb.html +995 -0
- data/doc/ssp.html +1647 -0
- data/doc/sst.html +1423 -0
- data/doc/ssy.html +3438 -0
- data/doc/st.html +16 -0
- data/doc/stb.html +284 -0
- data/doc/stg.html +1729 -0
- data/doc/stp.html +532 -0
- data/doc/str.html +1346 -0
- data/doc/stz.html +211 -0
- data/doc/sy.html +16 -0
- data/doc/tb.html +16 -0
- data/doc/tg.html +16 -0
- data/doc/tp.html +16 -0
- data/doc/tr.html +16 -0
- data/doc/tz.html +16 -0
- data/doc/un.html +14 -0
- data/doc/up.html +14 -0
- data/doc/z.html +36 -0
- data/doc/zbd.html +161 -0
- data/doc/zgb.html +1862 -0
- data/doc/zge.html +5258 -0
- data/doc/zgg.html +2027 -0
- data/doc/zgt.html +711 -0
- data/doc/zhb.html +1031 -0
- data/doc/zhe.html +3162 -0
- data/doc/zhg.html +201 -0
- data/doc/zhp.html +1697 -0
- data/doc/zhs.html +386 -0
- data/doc/zpb.html +994 -0
- data/doc/zpo.html +1517 -0
- data/doc/zpp.html +770 -0
- data/doc/zpt.html +706 -0
- data/doc/zsp.html +905 -0
- data/doc/zst.html +743 -0
- data/doc/zsy.html +2191 -0
- data/doc/ztb.html +284 -0
- data/doc/ztg.html +1544 -0
- data/doc/ztp.html +553 -0
- data/doc/ztr.html +1281 -0
- data/doc/ztz.html +211 -0
- data/doc/zun.html +2553 -0
- data/doc/zup.html +166 -0
- data/ext/cbbcsd.c +279 -0
- data/ext/cbdsqr.c +178 -0
- data/ext/cgbbrd.c +153 -0
- data/ext/cgbcon.c +94 -0
- data/ext/cgbequ.c +94 -0
- data/ext/cgbequb.c +92 -0
- data/ext/cgbrfs.c +157 -0
- data/ext/cgbrfsx.c +245 -0
- data/ext/cgbsv.c +111 -0
- data/ext/cgbsvx.c +282 -0
- data/ext/cgbsvxx.c +285 -0
- data/ext/cgbtf2.c +89 -0
- data/ext/cgbtrf.c +89 -0
- data/ext/cgbtrs.c +107 -0
- data/ext/cgebak.c +97 -0
- data/ext/cgebal.c +87 -0
- data/ext/cgebd2.c +108 -0
- data/ext/cgebrd.c +123 -0
- data/ext/cgecon.c +74 -0
- data/ext/cgeequ.c +84 -0
- data/ext/cgeequb.c +84 -0
- data/ext/cgees.c +138 -0
- data/ext/cgeesx.c +148 -0
- data/ext/cgeev.c +128 -0
- data/ext/cgeevx.c +169 -0
- data/ext/cgegs.c +162 -0
- data/ext/cgegv.c +167 -0
- data/ext/cgehd2.c +88 -0
- data/ext/cgehrd.c +103 -0
- data/ext/cgelq2.c +82 -0
- data/ext/cgelqf.c +99 -0
- data/ext/cgels.c +133 -0
- data/ext/cgelsd.c +150 -0
- data/ext/cgelss.c +147 -0
- data/ext/cgelsx.c +135 -0
- data/ext/cgelsy.c +162 -0
- data/ext/cgeql2.c +84 -0
- data/ext/cgeqlf.c +99 -0
- data/ext/cgeqp3.c +125 -0
- data/ext/cgeqpf.c +110 -0
- data/ext/cgeqr2.c +84 -0
- data/ext/cgeqr2p.c +84 -0
- data/ext/cgeqrf.c +99 -0
- data/ext/cgeqrfp.c +99 -0
- data/ext/cgerfs.c +149 -0
- data/ext/cgerfsx.c +215 -0
- data/ext/cgerq2.c +82 -0
- data/ext/cgerqf.c +99 -0
- data/ext/cgesc2.c +104 -0
- data/ext/cgesdd.c +131 -0
- data/ext/cgesv.c +103 -0
- data/ext/cgesvd.c +142 -0
- data/ext/cgesvx.c +274 -0
- data/ext/cgesvxx.c +277 -0
- data/ext/cgetc2.c +85 -0
- data/ext/cgetf2.c +81 -0
- data/ext/cgetrf.c +81 -0
- data/ext/cgetri.c +99 -0
- data/ext/cgetrs.c +99 -0
- data/ext/cggbak.c +109 -0
- data/ext/cggbal.c +124 -0
- data/ext/cgges.c +188 -0
- data/ext/cggesx.c +226 -0
- data/ext/cggev.c +167 -0
- data/ext/cggevx.c +222 -0
- data/ext/cggglm.c +152 -0
- data/ext/cgghrd.c +163 -0
- data/ext/cgglse.c +167 -0
- data/ext/cggqrf.c +133 -0
- data/ext/cggrqf.c +137 -0
- data/ext/cggsvd.c +180 -0
- data/ext/cggsvp.c +170 -0
- data/ext/cgtcon.c +117 -0
- data/ext/cgtrfs.c +205 -0
- data/ext/cgtsv.c +138 -0
- data/ext/cgtsvx.c +252 -0
- data/ext/cgttrf.c +128 -0
- data/ext/cgttrs.c +133 -0
- data/ext/cgtts2.c +130 -0
- data/ext/chbev.c +106 -0
- data/ext/chbevd.c +154 -0
- data/ext/chbevx.c +156 -0
- data/ext/chbgst.c +116 -0
- data/ext/chbgv.c +136 -0
- data/ext/chbgvd.c +184 -0
- data/ext/chbgvx.c +185 -0
- data/ext/chbtrd.c +126 -0
- data/ext/checon.c +83 -0
- data/ext/cheequb.c +78 -0
- data/ext/cheev.c +106 -0
- data/ext/cheevd.c +139 -0
- data/ext/cheevr.c +186 -0
- data/ext/cheevx.c +156 -0
- data/ext/chegs2.c +91 -0
- data/ext/chegst.c +91 -0
- data/ext/chegv.c +136 -0
- data/ext/chegvd.c +169 -0
- data/ext/chegvx.c +186 -0
- data/ext/cherfs.c +149 -0
- data/ext/cherfsx.c +214 -0
- data/ext/chesv.c +119 -0
- data/ext/chesvx.c +179 -0
- data/ext/chesvxx.c +254 -0
- data/ext/chetd2.c +97 -0
- data/ext/chetf2.c +81 -0
- data/ext/chetrd.c +109 -0
- data/ext/chetrf.c +93 -0
- data/ext/chetri.c +88 -0
- data/ext/chetrs.c +99 -0
- data/ext/chetrs2.c +102 -0
- data/ext/chfrk.c +105 -0
- data/ext/chgeqz.c +204 -0
- data/ext/chla_transtype.c +47 -0
- data/ext/chpcon.c +81 -0
- data/ext/chpev.c +101 -0
- data/ext/chpevd.c +149 -0
- data/ext/chpevx.c +140 -0
- data/ext/chpgst.c +90 -0
- data/ext/chpgv.c +128 -0
- data/ext/chpgvd.c +166 -0
- data/ext/chpgvx.c +166 -0
- data/ext/chprfs.c +145 -0
- data/ext/chpsv.c +106 -0
- data/ext/chpsvx.c +159 -0
- data/ext/chptrd.c +96 -0
- data/ext/chptrf.c +80 -0
- data/ext/chptri.c +85 -0
- data/ext/chptrs.c +97 -0
- data/ext/chsein.c +181 -0
- data/ext/chseqr.c +141 -0
- data/ext/cla_gbamv.c +123 -0
- data/ext/cla_gbrcond_c.c +138 -0
- data/ext/cla_gbrcond_x.c +134 -0
- data/ext/cla_gbrfsx_extended.c +291 -0
- data/ext/cla_gbrpvgrw.c +83 -0
- data/ext/cla_geamv.c +113 -0
- data/ext/cla_gercond_c.c +130 -0
- data/ext/cla_gercond_x.c +126 -0
- data/ext/cla_gerfsx_extended.c +277 -0
- data/ext/cla_heamv.c +112 -0
- data/ext/cla_hercond_c.c +130 -0
- data/ext/cla_hercond_x.c +126 -0
- data/ext/cla_herfsx_extended.c +279 -0
- data/ext/cla_herpvgrw.c +103 -0
- data/ext/cla_lin_berr.c +80 -0
- data/ext/cla_porcond_c.c +118 -0
- data/ext/cla_porcond_x.c +114 -0
- data/ext/cla_porfsx_extended.c +267 -0
- data/ext/cla_porpvgrw.c +91 -0
- data/ext/cla_rpvgrw.c +75 -0
- data/ext/cla_syamv.c +111 -0
- data/ext/cla_syrcond_c.c +130 -0
- data/ext/cla_syrcond_x.c +126 -0
- data/ext/cla_syrfsx_extended.c +279 -0
- data/ext/cla_syrpvgrw.c +103 -0
- data/ext/cla_wwaddw.c +98 -0
- data/ext/clabrd.c +128 -0
- data/ext/clacgv.c +71 -0
- data/ext/clacn2.c +99 -0
- data/ext/clacon.c +76 -0
- data/ext/clacp2.c +73 -0
- data/ext/clacpy.c +73 -0
- data/ext/clacrm.c +86 -0
- data/ext/clacrt.c +104 -0
- data/ext/cladiv.c +53 -0
- data/ext/claed0.c +130 -0
- data/ext/claed7.c +243 -0
- data/ext/claed8.c +194 -0
- data/ext/claein.c +109 -0
- data/ext/claesy.c +70 -0
- data/ext/claev2.c +67 -0
- data/ext/clag2z.c +72 -0
- data/ext/clags2.c +88 -0
- data/ext/clagtm.c +128 -0
- data/ext/clahef.c +93 -0
- data/ext/clahqr.c +131 -0
- data/ext/clahr2.c +108 -0
- data/ext/clahrd.c +108 -0
- data/ext/claic1.c +86 -0
- data/ext/clals0.c +197 -0
- data/ext/clalsa.c +266 -0
- data/ext/clalsd.c +141 -0
- data/ext/clangb.c +72 -0
- data/ext/clange.c +70 -0
- data/ext/clangt.c +83 -0
- data/ext/clanhb.c +74 -0
- data/ext/clanhe.c +68 -0
- data/ext/clanhf.c +76 -0
- data/ext/clanhp.c +70 -0
- data/ext/clanhs.c +66 -0
- data/ext/clanht.c +71 -0
- data/ext/clansb.c +74 -0
- data/ext/clansp.c +70 -0
- data/ext/clansy.c +68 -0
- data/ext/clantb.c +76 -0
- data/ext/clantp.c +76 -0
- data/ext/clantr.c +78 -0
- data/ext/clapll.c +101 -0
- data/ext/clapmr.c +93 -0
- data/ext/clapmt.c +97 -0
- data/ext/claqgb.c +113 -0
- data/ext/claqge.c +105 -0
- data/ext/claqhb.c +93 -0
- data/ext/claqhe.c +93 -0
- data/ext/claqhp.c +90 -0
- data/ext/claqp2.c +154 -0
- data/ext/claqps.c +204 -0
- data/ext/claqr0.c +141 -0
- data/ext/claqr1.c +72 -0
- data/ext/claqr2.c +170 -0
- data/ext/claqr3.c +170 -0
- data/ext/claqr4.c +141 -0
- data/ext/claqr5.c +175 -0
- data/ext/claqsb.c +97 -0
- data/ext/claqsp.c +90 -0
- data/ext/claqsy.c +93 -0
- data/ext/clar1v.c +169 -0
- data/ext/clar2v.c +145 -0
- data/ext/clarcm.c +82 -0
- data/ext/clarf.c +98 -0
- data/ext/clarfb.c +119 -0
- data/ext/clarfg.c +80 -0
- data/ext/clarfgp.c +80 -0
- data/ext/clarft.c +101 -0
- data/ext/clarfx.c +90 -0
- data/ext/clargv.c +110 -0
- data/ext/clarnv.c +79 -0
- data/ext/clarrv.c +267 -0
- data/ext/clarscl2.c +78 -0
- data/ext/clartg.c +59 -0
- data/ext/clartv.c +126 -0
- data/ext/clarz.c +102 -0
- data/ext/clarzb.c +123 -0
- data/ext/clarzt.c +101 -0
- data/ext/clascl.c +93 -0
- data/ext/clascl2.c +78 -0
- data/ext/claset.c +84 -0
- data/ext/clasr.c +106 -0
- data/ext/classq.c +66 -0
- data/ext/claswp.c +90 -0
- data/ext/clasyf.c +93 -0
- data/ext/clatbs.c +126 -0
- data/ext/clatdf.c +115 -0
- data/ext/clatps.c +120 -0
- data/ext/clatrd.c +101 -0
- data/ext/clatrs.c +122 -0
- data/ext/clatrz.c +83 -0
- data/ext/clatzm.c +128 -0
- data/ext/clauu2.c +73 -0
- data/ext/clauum.c +73 -0
- data/ext/cpbcon.c +78 -0
- data/ext/cpbequ.c +79 -0
- data/ext/cpbrfs.c +141 -0
- data/ext/cpbstf.c +77 -0
- data/ext/cpbsv.c +103 -0
- data/ext/cpbsvx.c +197 -0
- data/ext/cpbtf2.c +77 -0
- data/ext/cpbtrf.c +77 -0
- data/ext/cpbtrs.c +91 -0
- data/ext/cpftrf.c +78 -0
- data/ext/cpftri.c +78 -0
- data/ext/cpftrs.c +93 -0
- data/ext/cpocon.c +74 -0
- data/ext/cpoequ.c +71 -0
- data/ext/cpoequb.c +71 -0
- data/ext/cporfs.c +137 -0
- data/ext/cporfsx.c +202 -0
- data/ext/cposv.c +99 -0
- data/ext/cposvx.c +193 -0
- data/ext/cposvxx.c +231 -0
- data/ext/cpotf2.c +73 -0
- data/ext/cpotrf.c +73 -0
- data/ext/cpotri.c +73 -0
- data/ext/cpotrs.c +87 -0
- data/ext/cppcon.c +74 -0
- data/ext/cppequ.c +75 -0
- data/ext/cpprfs.c +135 -0
- data/ext/cppsv.c +100 -0
- data/ext/cppsvx.c +187 -0
- data/ext/cpptrf.c +74 -0
- data/ext/cpptri.c +74 -0
- data/ext/cpptrs.c +89 -0
- data/ext/cpstf2.c +91 -0
- data/ext/cpstrf.c +91 -0
- data/ext/cptcon.c +77 -0
- data/ext/cpteqr.c +122 -0
- data/ext/cptrfs.c +157 -0
- data/ext/cptsv.c +115 -0
- data/ext/cptsvx.c +167 -0
- data/ext/cpttrf.c +89 -0
- data/ext/cpttrs.c +97 -0
- data/ext/cptts2.c +94 -0
- data/ext/crot.c +103 -0
- data/ext/cspcon.c +81 -0
- data/ext/cspmv.c +111 -0
- data/ext/cspr.c +92 -0
- data/ext/csprfs.c +145 -0
- data/ext/cspsv.c +106 -0
- data/ext/cspsvx.c +159 -0
- data/ext/csptrf.c +80 -0
- data/ext/csptri.c +85 -0
- data/ext/csptrs.c +97 -0
- data/ext/csrscl.c +75 -0
- data/ext/cstedc.c +173 -0
- data/ext/cstegr.c +184 -0
- data/ext/cstein.c +130 -0
- data/ext/cstemr.c +189 -0
- data/ext/csteqr.c +122 -0
- data/ext/csycon.c +83 -0
- data/ext/csyconv.c +80 -0
- data/ext/csyequb.c +78 -0
- data/ext/csymv.c +111 -0
- data/ext/csyr.c +91 -0
- data/ext/csyrfs.c +149 -0
- data/ext/csyrfsx.c +214 -0
- data/ext/csysv.c +125 -0
- data/ext/csysvx.c +179 -0
- data/ext/csysvxx.c +254 -0
- data/ext/csyswapr.c +78 -0
- data/ext/csytf2.c +81 -0
- data/ext/csytrf.c +93 -0
- data/ext/csytri.c +88 -0
- data/ext/csytri2.c +104 -0
- data/ext/csytri2x.c +92 -0
- data/ext/csytrs.c +99 -0
- data/ext/csytrs2.c +102 -0
- data/ext/ctbcon.c +82 -0
- data/ext/ctbrfs.c +123 -0
- data/ext/ctbtrs.c +99 -0
- data/ext/ctfsm.c +107 -0
- data/ext/ctftri.c +82 -0
- data/ext/ctfttp.c +75 -0
- data/ext/ctfttr.c +76 -0
- data/ext/ctgevc.c +152 -0
- data/ext/ctgex2.c +167 -0
- data/ext/ctgexc.c +168 -0
- data/ext/ctgsen.c +240 -0
- data/ext/ctgsja.c +223 -0
- data/ext/ctgsna.c +160 -0
- data/ext/ctgsy2.c +172 -0
- data/ext/ctgsyl.c +186 -0
- data/ext/ctpcon.c +78 -0
- data/ext/ctprfs.c +119 -0
- data/ext/ctptri.c +78 -0
- data/ext/ctptrs.c +97 -0
- data/ext/ctpttf.c +75 -0
- data/ext/ctpttr.c +72 -0
- data/ext/ctrcon.c +78 -0
- data/ext/ctrevc.c +150 -0
- data/ext/ctrexc.c +107 -0
- data/ext/ctrrfs.c +119 -0
- data/ext/ctrsen.c +150 -0
- data/ext/ctrsna.c +133 -0
- data/ext/ctrsyl.c +112 -0
- data/ext/ctrti2.c +77 -0
- data/ext/ctrtri.c +77 -0
- data/ext/ctrtrs.c +95 -0
- data/ext/ctrttf.c +73 -0
- data/ext/ctrttp.c +69 -0
- data/ext/ctzrqf.c +79 -0
- data/ext/ctzrzf.c +97 -0
- data/ext/cunbdb.c +228 -0
- data/ext/cuncsd.c +200 -0
- data/ext/cung2l.c +88 -0
- data/ext/cung2r.c +88 -0
- data/ext/cungbr.c +111 -0
- data/ext/cunghr.c +107 -0
- data/ext/cungl2.c +86 -0
- data/ext/cunglq.c +103 -0
- data/ext/cungql.c +103 -0
- data/ext/cungqr.c +103 -0
- data/ext/cungr2.c +86 -0
- data/ext/cungrq.c +103 -0
- data/ext/cungtr.c +103 -0
- data/ext/cunm2l.c +110 -0
- data/ext/cunm2r.c +110 -0
- data/ext/cunmbr.c +135 -0
- data/ext/cunmhr.c +129 -0
- data/ext/cunml2.c +106 -0
- data/ext/cunmlq.c +121 -0
- data/ext/cunmql.c +125 -0
- data/ext/cunmqr.c +125 -0
- data/ext/cunmr2.c +106 -0
- data/ext/cunmr3.c +110 -0
- data/ext/cunmrq.c +121 -0
- data/ext/cunmrz.c +125 -0
- data/ext/cunmtr.c +125 -0
- data/ext/cupgtr.c +87 -0
- data/ext/cupmtr.c +112 -0
- data/ext/dbbcsd.c +283 -0
- data/ext/dbdsdc.c +147 -0
- data/ext/dbdsqr.c +178 -0
- data/ext/ddisna.c +71 -0
- data/ext/dgbbrd.c +150 -0
- data/ext/dgbcon.c +94 -0
- data/ext/dgbequ.c +94 -0
- data/ext/dgbequb.c +92 -0
- data/ext/dgbrfs.c +157 -0
- data/ext/dgbrfsx.c +245 -0
- data/ext/dgbsv.c +111 -0
- data/ext/dgbsvx.c +282 -0
- data/ext/dgbsvxx.c +285 -0
- data/ext/dgbtf2.c +89 -0
- data/ext/dgbtrf.c +89 -0
- data/ext/dgbtrs.c +107 -0
- data/ext/dgebak.c +97 -0
- data/ext/dgebal.c +87 -0
- data/ext/dgebd2.c +108 -0
- data/ext/dgebrd.c +123 -0
- data/ext/dgecon.c +74 -0
- data/ext/dgeequ.c +84 -0
- data/ext/dgeequb.c +84 -0
- data/ext/dgees.c +144 -0
- data/ext/dgeesx.c +166 -0
- data/ext/dgeev.c +133 -0
- data/ext/dgeevx.c +177 -0
- data/ext/dgegs.c +167 -0
- data/ext/dgegv.c +167 -0
- data/ext/dgehd2.c +88 -0
- data/ext/dgehrd.c +103 -0
- data/ext/dgejsv.c +155 -0
- data/ext/dgelq2.c +82 -0
- data/ext/dgelqf.c +99 -0
- data/ext/dgels.c +133 -0
- data/ext/dgelsd.c +145 -0
- data/ext/dgelss.c +144 -0
- data/ext/dgelsx.c +132 -0
- data/ext/dgelsy.c +159 -0
- data/ext/dgeql2.c +84 -0
- data/ext/dgeqlf.c +99 -0
- data/ext/dgeqp3.c +122 -0
- data/ext/dgeqpf.c +107 -0
- data/ext/dgeqr2.c +84 -0
- data/ext/dgeqr2p.c +84 -0
- data/ext/dgeqrf.c +99 -0
- data/ext/dgeqrfp.c +99 -0
- data/ext/dgerfs.c +149 -0
- data/ext/dgerfsx.c +215 -0
- data/ext/dgerq2.c +82 -0
- data/ext/dgerqf.c +99 -0
- data/ext/dgesc2.c +104 -0
- data/ext/dgesdd.c +128 -0
- data/ext/dgesv.c +103 -0
- data/ext/dgesvd.c +139 -0
- data/ext/dgesvj.c +152 -0
- data/ext/dgesvx.c +274 -0
- data/ext/dgesvxx.c +277 -0
- data/ext/dgetc2.c +85 -0
- data/ext/dgetf2.c +81 -0
- data/ext/dgetrf.c +81 -0
- data/ext/dgetri.c +99 -0
- data/ext/dgetrs.c +99 -0
- data/ext/dggbak.c +109 -0
- data/ext/dggbal.c +124 -0
- data/ext/dgges.c +194 -0
- data/ext/dggesx.c +227 -0
- data/ext/dggev.c +167 -0
- data/ext/dggevx.c +225 -0
- data/ext/dggglm.c +152 -0
- data/ext/dgghrd.c +163 -0
- data/ext/dgglse.c +167 -0
- data/ext/dggqrf.c +133 -0
- data/ext/dggrqf.c +137 -0
- data/ext/dggsvd.c +177 -0
- data/ext/dggsvp.c +167 -0
- data/ext/dgsvj0.c +178 -0
- data/ext/dgsvj1.c +182 -0
- data/ext/dgtcon.c +120 -0
- data/ext/dgtrfs.c +205 -0
- data/ext/dgtsv.c +138 -0
- data/ext/dgtsvx.c +252 -0
- data/ext/dgttrf.c +128 -0
- data/ext/dgttrs.c +133 -0
- data/ext/dgtts2.c +130 -0
- data/ext/dhgeqz.c +209 -0
- data/ext/dhsein.c +201 -0
- data/ext/dhseqr.c +149 -0
- data/ext/disnan.c +47 -0
- data/ext/dla_gbamv.c +125 -0
- data/ext/dla_gbrcond.c +138 -0
- data/ext/dla_gbrfsx_extended.c +289 -0
- data/ext/dla_gbrpvgrw.c +83 -0
- data/ext/dla_geamv.c +115 -0
- data/ext/dla_gercond.c +130 -0
- data/ext/dla_gerfsx_extended.c +277 -0
- data/ext/dla_lin_berr.c +80 -0
- data/ext/dla_porcond.c +118 -0
- data/ext/dla_porfsx_extended.c +267 -0
- data/ext/dla_porpvgrw.c +91 -0
- data/ext/dla_rpvgrw.c +75 -0
- data/ext/dla_syamv.c +109 -0
- data/ext/dla_syrcond.c +130 -0
- data/ext/dla_syrfsx_extended.c +279 -0
- data/ext/dla_syrpvgrw.c +103 -0
- data/ext/dla_wwaddw.c +98 -0
- data/ext/dlabad.c +50 -0
- data/ext/dlabrd.c +128 -0
- data/ext/dlacn2.c +102 -0
- data/ext/dlacon.c +79 -0
- data/ext/dlacpy.c +73 -0
- data/ext/dladiv.c +62 -0
- data/ext/dlae2.c +58 -0
- data/ext/dlaebz.c +214 -0
- data/ext/dlaed0.c +123 -0
- data/ext/dlaed1.c +129 -0
- data/ext/dlaed2.c +185 -0
- data/ext/dlaed3.c +157 -0
- data/ext/dlaed4.c +86 -0
- data/ext/dlaed5.c +83 -0
- data/ext/dlaed6.c +86 -0
- data/ext/dlaed7.c +244 -0
- data/ext/dlaed8.c +202 -0
- data/ext/dlaed9.c +107 -0
- data/ext/dlaeda.c +156 -0
- data/ext/dlaein.c +139 -0
- data/ext/dlaev2.c +64 -0
- data/ext/dlaexc.c +114 -0
- data/ext/dlag2.c +87 -0
- data/ext/dlag2s.c +72 -0
- data/ext/dlags2.c +86 -0
- data/ext/dlagtf.c +136 -0
- data/ext/dlagtm.c +128 -0
- data/ext/dlagts.c +135 -0
- data/ext/dlagv2.c +128 -0
- data/ext/dlahqr.c +139 -0
- data/ext/dlahr2.c +108 -0
- data/ext/dlahrd.c +108 -0
- data/ext/dlaic1.c +85 -0
- data/ext/dlaln2.c +116 -0
- data/ext/dlals0.c +197 -0
- data/ext/dlalsa.c +266 -0
- data/ext/dlalsd.c +138 -0
- data/ext/dlamrg.c +76 -0
- data/ext/dlaneg.c +79 -0
- data/ext/dlangb.c +74 -0
- data/ext/dlange.c +70 -0
- data/ext/dlangt.c +83 -0
- data/ext/dlanhs.c +66 -0
- data/ext/dlansb.c +74 -0
- data/ext/dlansf.c +76 -0
- data/ext/dlansp.c +72 -0
- data/ext/dlanst.c +71 -0
- data/ext/dlansy.c +70 -0
- data/ext/dlantb.c +78 -0
- data/ext/dlantp.c +76 -0
- data/ext/dlantr.c +78 -0
- data/ext/dlanv2.c +78 -0
- data/ext/dlapll.c +101 -0
- data/ext/dlapmr.c +93 -0
- data/ext/dlapmt.c +97 -0
- data/ext/dlapy2.c +51 -0
- data/ext/dlapy3.c +55 -0
- data/ext/dlaqgb.c +113 -0
- data/ext/dlaqge.c +105 -0
- data/ext/dlaqp2.c +154 -0
- data/ext/dlaqps.c +204 -0
- data/ext/dlaqr0.c +149 -0
- data/ext/dlaqr1.c +78 -0
- data/ext/dlaqr2.c +178 -0
- data/ext/dlaqr3.c +178 -0
- data/ext/dlaqr4.c +149 -0
- data/ext/dlaqr5.c +196 -0
- data/ext/dlaqsb.c +97 -0
- data/ext/dlaqsp.c +90 -0
- data/ext/dlaqsy.c +93 -0
- data/ext/dlaqtr.c +110 -0
- data/ext/dlar1v.c +169 -0
- data/ext/dlar2v.c +145 -0
- data/ext/dlarf.c +97 -0
- data/ext/dlarfb.c +119 -0
- data/ext/dlarfg.c +79 -0
- data/ext/dlarfgp.c +79 -0
- data/ext/dlarft.c +101 -0
- data/ext/dlarfx.c +89 -0
- data/ext/dlargv.c +110 -0
- data/ext/dlarnv.c +79 -0
- data/ext/dlarra.c +120 -0
- data/ext/dlarrb.c +174 -0
- data/ext/dlarrc.c +92 -0
- data/ext/dlarrd.c +186 -0
- data/ext/dlarre.c +217 -0
- data/ext/dlarrf.c +172 -0
- data/ext/dlarrj.c +143 -0
- data/ext/dlarrk.c +93 -0
- data/ext/dlarrr.c +78 -0
- data/ext/dlarrv.c +267 -0
- data/ext/dlarscl2.c +78 -0
- data/ext/dlartg.c +57 -0
- data/ext/dlartgp.c +57 -0
- data/ext/dlartgs.c +58 -0
- data/ext/dlartv.c +126 -0
- data/ext/dlaruv.c +75 -0
- data/ext/dlarz.c +101 -0
- data/ext/dlarzb.c +123 -0
- data/ext/dlarzt.c +101 -0
- data/ext/dlas2.c +58 -0
- data/ext/dlascl.c +93 -0
- data/ext/dlascl2.c +78 -0
- data/ext/dlasd0.c +116 -0
- data/ext/dlasd1.c +158 -0
- data/ext/dlasd2.c +224 -0
- data/ext/dlasd3.c +198 -0
- data/ext/dlasd4.c +89 -0
- data/ext/dlasd5.c +86 -0
- data/ext/dlasd6.c +232 -0
- data/ext/dlasd7.c +221 -0
- data/ext/dlasd8.c +169 -0
- data/ext/dlasda.c +217 -0
- data/ext/dlasdq.c +182 -0
- data/ext/dlasdt.c +78 -0
- data/ext/dlaset.c +82 -0
- data/ext/dlasq1.c +92 -0
- data/ext/dlasq2.c +70 -0
- data/ext/dlasq3.c +134 -0
- data/ext/dlasq4.c +103 -0
- data/ext/dlasq5.c +90 -0
- data/ext/dlasq6.c +82 -0
- data/ext/dlasr.c +106 -0
- data/ext/dlasrt.c +70 -0
- data/ext/dlassq.c +66 -0
- data/ext/dlasv2.c +70 -0
- data/ext/dlaswp.c +90 -0
- data/ext/dlasy2.c +122 -0
- data/ext/dlasyf.c +93 -0
- data/ext/dlat2s.c +72 -0
- data/ext/dlatbs.c +126 -0
- data/ext/dlatdf.c +115 -0
- data/ext/dlatps.c +120 -0
- data/ext/dlatrd.c +101 -0
- data/ext/dlatrs.c +122 -0
- data/ext/dlatrz.c +83 -0
- data/ext/dlatzm.c +127 -0
- data/ext/dlauu2.c +73 -0
- data/ext/dlauum.c +73 -0
- data/ext/dopgtr.c +87 -0
- data/ext/dopmtr.c +112 -0
- data/ext/dorbdb.c +228 -0
- data/ext/dorcsd.c +193 -0
- data/ext/dorg2l.c +88 -0
- data/ext/dorg2r.c +88 -0
- data/ext/dorgbr.c +111 -0
- data/ext/dorghr.c +107 -0
- data/ext/dorgl2.c +86 -0
- data/ext/dorglq.c +103 -0
- data/ext/dorgql.c +103 -0
- data/ext/dorgqr.c +103 -0
- data/ext/dorgr2.c +86 -0
- data/ext/dorgrq.c +103 -0
- data/ext/dorgtr.c +103 -0
- data/ext/dorm2l.c +110 -0
- data/ext/dorm2r.c +110 -0
- data/ext/dormbr.c +135 -0
- data/ext/dormhr.c +129 -0
- data/ext/dorml2.c +106 -0
- data/ext/dormlq.c +121 -0
- data/ext/dormql.c +125 -0
- data/ext/dormqr.c +125 -0
- data/ext/dormr2.c +106 -0
- data/ext/dormr3.c +110 -0
- data/ext/dormrq.c +121 -0
- data/ext/dormrz.c +125 -0
- data/ext/dormtr.c +125 -0
- data/ext/dpbcon.c +78 -0
- data/ext/dpbequ.c +79 -0
- data/ext/dpbrfs.c +141 -0
- data/ext/dpbstf.c +77 -0
- data/ext/dpbsv.c +103 -0
- data/ext/dpbsvx.c +197 -0
- data/ext/dpbtf2.c +77 -0
- data/ext/dpbtrf.c +77 -0
- data/ext/dpbtrs.c +91 -0
- data/ext/dpftrf.c +78 -0
- data/ext/dpftri.c +78 -0
- data/ext/dpftrs.c +93 -0
- data/ext/dpocon.c +74 -0
- data/ext/dpoequ.c +71 -0
- data/ext/dpoequb.c +71 -0
- data/ext/dporfs.c +137 -0
- data/ext/dporfsx.c +202 -0
- data/ext/dposv.c +99 -0
- data/ext/dposvx.c +193 -0
- data/ext/dposvxx.c +231 -0
- data/ext/dpotf2.c +73 -0
- data/ext/dpotrf.c +73 -0
- data/ext/dpotri.c +73 -0
- data/ext/dpotrs.c +87 -0
- data/ext/dppcon.c +74 -0
- data/ext/dppequ.c +75 -0
- data/ext/dpprfs.c +135 -0
- data/ext/dppsv.c +100 -0
- data/ext/dppsvx.c +187 -0
- data/ext/dpptrf.c +74 -0
- data/ext/dpptri.c +74 -0
- data/ext/dpptrs.c +89 -0
- data/ext/dpstf2.c +91 -0
- data/ext/dpstrf.c +91 -0
- data/ext/dptcon.c +77 -0
- data/ext/dpteqr.c +122 -0
- data/ext/dptrfs.c +150 -0
- data/ext/dptsv.c +115 -0
- data/ext/dptsvx.c +164 -0
- data/ext/dpttrf.c +89 -0
- data/ext/dpttrs.c +93 -0
- data/ext/dptts2.c +90 -0
- data/ext/drscl.c +75 -0
- data/ext/dsbev.c +103 -0
- data/ext/dsbevd.c +136 -0
- data/ext/dsbevx.c +153 -0
- data/ext/dsbgst.c +113 -0
- data/ext/dsbgv.c +133 -0
- data/ext/dsbgvd.c +166 -0
- data/ext/dsbgvx.c +193 -0
- data/ext/dsbtrd.c +126 -0
- data/ext/dsfrk.c +105 -0
- data/ext/dsgesv.c +111 -0
- data/ext/dspcon.c +84 -0
- data/ext/dspev.c +98 -0
- data/ext/dspevd.c +131 -0
- data/ext/dspevx.c +137 -0
- data/ext/dspgst.c +90 -0
- data/ext/dspgv.c +125 -0
- data/ext/dspgvd.c +158 -0
- data/ext/dspgvx.c +164 -0
- data/ext/dsposv.c +107 -0
- data/ext/dsprfs.c +145 -0
- data/ext/dspsv.c +106 -0
- data/ext/dspsvx.c +159 -0
- data/ext/dsptrd.c +96 -0
- data/ext/dsptrf.c +80 -0
- data/ext/dsptri.c +85 -0
- data/ext/dsptrs.c +97 -0
- data/ext/dstebz.c +131 -0
- data/ext/dstedc.c +155 -0
- data/ext/dstegr.c +184 -0
- data/ext/dstein.c +130 -0
- data/ext/dstemr.c +189 -0
- data/ext/dsteqr.c +122 -0
- data/ext/dsterf.c +89 -0
- data/ext/dstev.c +107 -0
- data/ext/dstevd.c +140 -0
- data/ext/dstevr.c +184 -0
- data/ext/dstevx.c +154 -0
- data/ext/dsycon.c +86 -0
- data/ext/dsyconv.c +80 -0
- data/ext/dsyequb.c +78 -0
- data/ext/dsyev.c +103 -0
- data/ext/dsyevd.c +121 -0
- data/ext/dsyevr.c +168 -0
- data/ext/dsyevx.c +153 -0
- data/ext/dsygs2.c +91 -0
- data/ext/dsygst.c +91 -0
- data/ext/dsygv.c +133 -0
- data/ext/dsygvd.c +151 -0
- data/ext/dsygvx.c +183 -0
- data/ext/dsyrfs.c +149 -0
- data/ext/dsyrfsx.c +214 -0
- data/ext/dsysv.c +125 -0
- data/ext/dsysvx.c +179 -0
- data/ext/dsysvxx.c +254 -0
- data/ext/dsyswapr.c +78 -0
- data/ext/dsytd2.c +97 -0
- data/ext/dsytf2.c +81 -0
- data/ext/dsytrd.c +109 -0
- data/ext/dsytrf.c +93 -0
- data/ext/dsytri.c +88 -0
- data/ext/dsytri2.c +104 -0
- data/ext/dsytri2x.c +92 -0
- data/ext/dsytrs.c +99 -0
- data/ext/dsytrs2.c +102 -0
- data/ext/dtbcon.c +82 -0
- data/ext/dtbrfs.c +123 -0
- data/ext/dtbtrs.c +99 -0
- data/ext/dtfsm.c +106 -0
- data/ext/dtftri.c +82 -0
- data/ext/dtfttp.c +75 -0
- data/ext/dtfttr.c +76 -0
- data/ext/dtgevc.c +149 -0
- data/ext/dtgex2.c +176 -0
- data/ext/dtgexc.c +183 -0
- data/ext/dtgsen.c +248 -0
- data/ext/dtgsja.c +223 -0
- data/ext/dtgsna.c +160 -0
- data/ext/dtgsy2.c +178 -0
- data/ext/dtgsyl.c +186 -0
- data/ext/dtpcon.c +78 -0
- data/ext/dtprfs.c +119 -0
- data/ext/dtptri.c +78 -0
- data/ext/dtptrs.c +97 -0
- data/ext/dtpttf.c +75 -0
- data/ext/dtpttr.c +72 -0
- data/ext/dtrcon.c +78 -0
- data/ext/dtrevc.c +146 -0
- data/ext/dtrexc.c +112 -0
- data/ext/dtrrfs.c +119 -0
- data/ext/dtrsen.c +165 -0
- data/ext/dtrsna.c +133 -0
- data/ext/dtrsyl.c +112 -0
- data/ext/dtrti2.c +77 -0
- data/ext/dtrtri.c +77 -0
- data/ext/dtrtrs.c +95 -0
- data/ext/dtrttf.c +73 -0
- data/ext/dtrttp.c +69 -0
- data/ext/dtzrqf.c +79 -0
- data/ext/dtzrzf.c +97 -0
- data/ext/dzsum1.c +59 -0
- data/ext/extconf.rb +136 -0
- data/ext/f2c_minimal.h +36 -0
- data/ext/icmax1.c +59 -0
- data/ext/ieeeck.c +55 -0
- data/ext/ilaclc.c +61 -0
- data/ext/ilaclr.c +61 -0
- data/ext/iladiag.c +47 -0
- data/ext/iladlc.c +61 -0
- data/ext/iladlr.c +61 -0
- data/ext/ilaenv.c +71 -0
- data/ext/ilaprec.c +47 -0
- data/ext/ilaslc.c +61 -0
- data/ext/ilaslr.c +61 -0
- data/ext/ilatrans.c +47 -0
- data/ext/ilauplo.c +47 -0
- data/ext/ilaver.c +49 -0
- data/ext/ilazlc.c +61 -0
- data/ext/ilazlr.c +61 -0
- data/ext/iparmq.c +71 -0
- data/ext/izmax1.c +59 -0
- data/ext/lsamen.c +55 -0
- data/ext/rb_lapack.c +3279 -0
- data/ext/rb_lapack.h +18 -0
- data/ext/sbbcsd.c +283 -0
- data/ext/sbdsdc.c +153 -0
- data/ext/sbdsqr.c +178 -0
- data/ext/scsum1.c +59 -0
- data/ext/sdisna.c +71 -0
- data/ext/sgbbrd.c +150 -0
- data/ext/sgbcon.c +94 -0
- data/ext/sgbequ.c +94 -0
- data/ext/sgbequb.c +92 -0
- data/ext/sgbrfs.c +157 -0
- data/ext/sgbrfsx.c +245 -0
- data/ext/sgbsv.c +111 -0
- data/ext/sgbsvx.c +282 -0
- data/ext/sgbsvxx.c +285 -0
- data/ext/sgbtf2.c +89 -0
- data/ext/sgbtrf.c +89 -0
- data/ext/sgbtrs.c +107 -0
- data/ext/sgebak.c +97 -0
- data/ext/sgebal.c +87 -0
- data/ext/sgebd2.c +108 -0
- data/ext/sgebrd.c +123 -0
- data/ext/sgecon.c +74 -0
- data/ext/sgeequ.c +84 -0
- data/ext/sgeequb.c +84 -0
- data/ext/sgees.c +144 -0
- data/ext/sgeesx.c +166 -0
- data/ext/sgeev.c +133 -0
- data/ext/sgeevx.c +177 -0
- data/ext/sgegs.c +167 -0
- data/ext/sgegv.c +167 -0
- data/ext/sgehd2.c +88 -0
- data/ext/sgehrd.c +103 -0
- data/ext/sgejsv.c +155 -0
- data/ext/sgelq2.c +82 -0
- data/ext/sgelqf.c +99 -0
- data/ext/sgels.c +133 -0
- data/ext/sgelsd.c +145 -0
- data/ext/sgelss.c +144 -0
- data/ext/sgelsx.c +132 -0
- data/ext/sgelsy.c +159 -0
- data/ext/sgeql2.c +84 -0
- data/ext/sgeqlf.c +99 -0
- data/ext/sgeqp3.c +122 -0
- data/ext/sgeqpf.c +107 -0
- data/ext/sgeqr2.c +84 -0
- data/ext/sgeqr2p.c +84 -0
- data/ext/sgeqrf.c +99 -0
- data/ext/sgeqrfp.c +99 -0
- data/ext/sgerfs.c +149 -0
- data/ext/sgerfsx.c +215 -0
- data/ext/sgerq2.c +82 -0
- data/ext/sgerqf.c +99 -0
- data/ext/sgesc2.c +104 -0
- data/ext/sgesdd.c +128 -0
- data/ext/sgesv.c +103 -0
- data/ext/sgesvd.c +139 -0
- data/ext/sgesvj.c +152 -0
- data/ext/sgesvx.c +274 -0
- data/ext/sgesvxx.c +277 -0
- data/ext/sgetc2.c +85 -0
- data/ext/sgetf2.c +81 -0
- data/ext/sgetrf.c +81 -0
- data/ext/sgetri.c +99 -0
- data/ext/sgetrs.c +99 -0
- data/ext/sggbak.c +109 -0
- data/ext/sggbal.c +124 -0
- data/ext/sgges.c +194 -0
- data/ext/sggesx.c +227 -0
- data/ext/sggev.c +167 -0
- data/ext/sggevx.c +225 -0
- data/ext/sggglm.c +152 -0
- data/ext/sgghrd.c +163 -0
- data/ext/sgglse.c +167 -0
- data/ext/sggqrf.c +133 -0
- data/ext/sggrqf.c +137 -0
- data/ext/sggsvd.c +177 -0
- data/ext/sggsvp.c +167 -0
- data/ext/sgsvj0.c +178 -0
- data/ext/sgsvj1.c +182 -0
- data/ext/sgtcon.c +120 -0
- data/ext/sgtrfs.c +205 -0
- data/ext/sgtsv.c +138 -0
- data/ext/sgtsvx.c +252 -0
- data/ext/sgttrf.c +128 -0
- data/ext/sgttrs.c +133 -0
- data/ext/sgtts2.c +130 -0
- data/ext/shgeqz.c +209 -0
- data/ext/shsein.c +201 -0
- data/ext/shseqr.c +149 -0
- data/ext/sisnan.c +47 -0
- data/ext/sla_gbamv.c +125 -0
- data/ext/sla_gbrcond.c +138 -0
- data/ext/sla_gbrfsx_extended.c +287 -0
- data/ext/sla_gbrpvgrw.c +83 -0
- data/ext/sla_geamv.c +115 -0
- data/ext/sla_gercond.c +130 -0
- data/ext/sla_gerfsx_extended.c +279 -0
- data/ext/sla_lin_berr.c +80 -0
- data/ext/sla_porcond.c +118 -0
- data/ext/sla_porfsx_extended.c +267 -0
- data/ext/sla_porpvgrw.c +91 -0
- data/ext/sla_rpvgrw.c +75 -0
- data/ext/sla_syamv.c +112 -0
- data/ext/sla_syrcond.c +130 -0
- data/ext/sla_syrfsx_extended.c +279 -0
- data/ext/sla_syrpvgrw.c +103 -0
- data/ext/sla_wwaddw.c +98 -0
- data/ext/slabad.c +50 -0
- data/ext/slabrd.c +128 -0
- data/ext/slacn2.c +102 -0
- data/ext/slacon.c +79 -0
- data/ext/slacpy.c +73 -0
- data/ext/sladiv.c +62 -0
- data/ext/slae2.c +58 -0
- data/ext/slaebz.c +214 -0
- data/ext/slaed0.c +123 -0
- data/ext/slaed1.c +129 -0
- data/ext/slaed2.c +185 -0
- data/ext/slaed3.c +157 -0
- data/ext/slaed4.c +86 -0
- data/ext/slaed5.c +83 -0
- data/ext/slaed6.c +86 -0
- data/ext/slaed7.c +244 -0
- data/ext/slaed8.c +202 -0
- data/ext/slaed9.c +107 -0
- data/ext/slaeda.c +156 -0
- data/ext/slaein.c +139 -0
- data/ext/slaev2.c +64 -0
- data/ext/slaexc.c +114 -0
- data/ext/slag2.c +87 -0
- data/ext/slag2d.c +72 -0
- data/ext/slags2.c +86 -0
- data/ext/slagtf.c +136 -0
- data/ext/slagtm.c +128 -0
- data/ext/slagts.c +135 -0
- data/ext/slagv2.c +128 -0
- data/ext/slahqr.c +139 -0
- data/ext/slahr2.c +108 -0
- data/ext/slahrd.c +110 -0
- data/ext/slaic1.c +85 -0
- data/ext/slaln2.c +116 -0
- data/ext/slals0.c +197 -0
- data/ext/slalsa.c +266 -0
- data/ext/slalsd.c +138 -0
- data/ext/slamrg.c +76 -0
- data/ext/slaneg.c +79 -0
- data/ext/slangb.c +74 -0
- data/ext/slange.c +70 -0
- data/ext/slangt.c +83 -0
- data/ext/slanhs.c +66 -0
- data/ext/slansb.c +74 -0
- data/ext/slansf.c +74 -0
- data/ext/slansp.c +72 -0
- data/ext/slanst.c +71 -0
- data/ext/slansy.c +70 -0
- data/ext/slantb.c +78 -0
- data/ext/slantp.c +76 -0
- data/ext/slantr.c +78 -0
- data/ext/slanv2.c +78 -0
- data/ext/slapll.c +101 -0
- data/ext/slapmr.c +93 -0
- data/ext/slapmt.c +97 -0
- data/ext/slapy2.c +51 -0
- data/ext/slapy3.c +55 -0
- data/ext/slaqgb.c +113 -0
- data/ext/slaqge.c +105 -0
- data/ext/slaqp2.c +154 -0
- data/ext/slaqps.c +204 -0
- data/ext/slaqr0.c +149 -0
- data/ext/slaqr1.c +78 -0
- data/ext/slaqr2.c +178 -0
- data/ext/slaqr3.c +178 -0
- data/ext/slaqr4.c +149 -0
- data/ext/slaqr5.c +196 -0
- data/ext/slaqsb.c +97 -0
- data/ext/slaqsp.c +90 -0
- data/ext/slaqsy.c +93 -0
- data/ext/slaqtr.c +110 -0
- data/ext/slar1v.c +169 -0
- data/ext/slar2v.c +145 -0
- data/ext/slarf.c +97 -0
- data/ext/slarfb.c +119 -0
- data/ext/slarfg.c +79 -0
- data/ext/slarfgp.c +79 -0
- data/ext/slarft.c +101 -0
- data/ext/slarfx.c +89 -0
- data/ext/slargv.c +110 -0
- data/ext/slarnv.c +79 -0
- data/ext/slarra.c +120 -0
- data/ext/slarrb.c +174 -0
- data/ext/slarrc.c +92 -0
- data/ext/slarrd.c +186 -0
- data/ext/slarre.c +217 -0
- data/ext/slarrf.c +172 -0
- data/ext/slarrj.c +143 -0
- data/ext/slarrk.c +93 -0
- data/ext/slarrr.c +78 -0
- data/ext/slarrv.c +267 -0
- data/ext/slarscl2.c +78 -0
- data/ext/slartg.c +57 -0
- data/ext/slartgp.c +57 -0
- data/ext/slartgs.c +58 -0
- data/ext/slartv.c +126 -0
- data/ext/slaruv.c +75 -0
- data/ext/slarz.c +101 -0
- data/ext/slarzb.c +123 -0
- data/ext/slarzt.c +101 -0
- data/ext/slas2.c +58 -0
- data/ext/slascl.c +93 -0
- data/ext/slascl2.c +78 -0
- data/ext/slasd0.c +116 -0
- data/ext/slasd1.c +156 -0
- data/ext/slasd2.c +224 -0
- data/ext/slasd3.c +208 -0
- data/ext/slasd4.c +89 -0
- data/ext/slasd5.c +86 -0
- data/ext/slasd6.c +232 -0
- data/ext/slasd7.c +221 -0
- data/ext/slasd8.c +169 -0
- data/ext/slasda.c +217 -0
- data/ext/slasdq.c +182 -0
- data/ext/slasdt.c +78 -0
- data/ext/slaset.c +82 -0
- data/ext/slasq1.c +92 -0
- data/ext/slasq2.c +70 -0
- data/ext/slasq3.c +134 -0
- data/ext/slasq4.c +103 -0
- data/ext/slasq5.c +90 -0
- data/ext/slasq6.c +82 -0
- data/ext/slasr.c +106 -0
- data/ext/slasrt.c +70 -0
- data/ext/slassq.c +66 -0
- data/ext/slasv2.c +70 -0
- data/ext/slaswp.c +90 -0
- data/ext/slasy2.c +122 -0
- data/ext/slasyf.c +93 -0
- data/ext/slatbs.c +126 -0
- data/ext/slatdf.c +115 -0
- data/ext/slatps.c +120 -0
- data/ext/slatrd.c +101 -0
- data/ext/slatrs.c +122 -0
- data/ext/slatrz.c +83 -0
- data/ext/slatzm.c +127 -0
- data/ext/slauu2.c +73 -0
- data/ext/slauum.c +73 -0
- data/ext/sopgtr.c +87 -0
- data/ext/sopmtr.c +112 -0
- data/ext/sorbdb.c +228 -0
- data/ext/sorcsd.c +193 -0
- data/ext/sorg2l.c +88 -0
- data/ext/sorg2r.c +88 -0
- data/ext/sorgbr.c +111 -0
- data/ext/sorghr.c +107 -0
- data/ext/sorgl2.c +86 -0
- data/ext/sorglq.c +103 -0
- data/ext/sorgql.c +103 -0
- data/ext/sorgqr.c +103 -0
- data/ext/sorgr2.c +86 -0
- data/ext/sorgrq.c +103 -0
- data/ext/sorgtr.c +103 -0
- data/ext/sorm2l.c +110 -0
- data/ext/sorm2r.c +110 -0
- data/ext/sormbr.c +135 -0
- data/ext/sormhr.c +129 -0
- data/ext/sorml2.c +106 -0
- data/ext/sormlq.c +121 -0
- data/ext/sormql.c +125 -0
- data/ext/sormqr.c +125 -0
- data/ext/sormr2.c +106 -0
- data/ext/sormr3.c +110 -0
- data/ext/sormrq.c +121 -0
- data/ext/sormrz.c +125 -0
- data/ext/sormtr.c +125 -0
- data/ext/spbcon.c +78 -0
- data/ext/spbequ.c +79 -0
- data/ext/spbrfs.c +141 -0
- data/ext/spbstf.c +77 -0
- data/ext/spbsv.c +103 -0
- data/ext/spbsvx.c +197 -0
- data/ext/spbtf2.c +77 -0
- data/ext/spbtrf.c +77 -0
- data/ext/spbtrs.c +91 -0
- data/ext/spftrf.c +78 -0
- data/ext/spftri.c +78 -0
- data/ext/spftrs.c +93 -0
- data/ext/spocon.c +74 -0
- data/ext/spoequ.c +71 -0
- data/ext/spoequb.c +71 -0
- data/ext/sporfs.c +137 -0
- data/ext/sporfsx.c +202 -0
- data/ext/sposv.c +99 -0
- data/ext/sposvx.c +193 -0
- data/ext/sposvxx.c +231 -0
- data/ext/spotf2.c +73 -0
- data/ext/spotrf.c +73 -0
- data/ext/spotri.c +73 -0
- data/ext/spotrs.c +87 -0
- data/ext/sppcon.c +74 -0
- data/ext/sppequ.c +75 -0
- data/ext/spprfs.c +135 -0
- data/ext/sppsv.c +100 -0
- data/ext/sppsvx.c +187 -0
- data/ext/spptrf.c +74 -0
- data/ext/spptri.c +74 -0
- data/ext/spptrs.c +89 -0
- data/ext/spstf2.c +91 -0
- data/ext/spstrf.c +91 -0
- data/ext/sptcon.c +77 -0
- data/ext/spteqr.c +122 -0
- data/ext/sptrfs.c +150 -0
- data/ext/sptsv.c +115 -0
- data/ext/sptsvx.c +164 -0
- data/ext/spttrf.c +89 -0
- data/ext/spttrs.c +93 -0
- data/ext/sptts2.c +90 -0
- data/ext/srscl.c +75 -0
- data/ext/ssbev.c +103 -0
- data/ext/ssbevd.c +136 -0
- data/ext/ssbevx.c +153 -0
- data/ext/ssbgst.c +113 -0
- data/ext/ssbgv.c +133 -0
- data/ext/ssbgvd.c +166 -0
- data/ext/ssbgvx.c +193 -0
- data/ext/ssbtrd.c +126 -0
- data/ext/ssfrk.c +105 -0
- data/ext/sspcon.c +84 -0
- data/ext/sspev.c +98 -0
- data/ext/sspevd.c +131 -0
- data/ext/sspevx.c +137 -0
- data/ext/sspgst.c +90 -0
- data/ext/sspgv.c +125 -0
- data/ext/sspgvd.c +158 -0
- data/ext/sspgvx.c +164 -0
- data/ext/ssprfs.c +145 -0
- data/ext/sspsv.c +106 -0
- data/ext/sspsvx.c +159 -0
- data/ext/ssptrd.c +96 -0
- data/ext/ssptrf.c +80 -0
- data/ext/ssptri.c +85 -0
- data/ext/ssptrs.c +97 -0
- data/ext/sstebz.c +131 -0
- data/ext/sstedc.c +155 -0
- data/ext/sstegr.c +184 -0
- data/ext/sstein.c +130 -0
- data/ext/sstemr.c +189 -0
- data/ext/ssteqr.c +122 -0
- data/ext/ssterf.c +89 -0
- data/ext/sstev.c +107 -0
- data/ext/sstevd.c +140 -0
- data/ext/sstevr.c +184 -0
- data/ext/sstevx.c +154 -0
- data/ext/ssycon.c +86 -0
- data/ext/ssyconv.c +80 -0
- data/ext/ssyequb.c +78 -0
- data/ext/ssyev.c +103 -0
- data/ext/ssyevd.c +121 -0
- data/ext/ssyevr.c +168 -0
- data/ext/ssyevx.c +153 -0
- data/ext/ssygs2.c +91 -0
- data/ext/ssygst.c +91 -0
- data/ext/ssygv.c +133 -0
- data/ext/ssygvd.c +151 -0
- data/ext/ssygvx.c +187 -0
- data/ext/ssyrfs.c +149 -0
- data/ext/ssyrfsx.c +214 -0
- data/ext/ssysv.c +125 -0
- data/ext/ssysvx.c +179 -0
- data/ext/ssysvxx.c +254 -0
- data/ext/ssyswapr.c +78 -0
- data/ext/ssytd2.c +97 -0
- data/ext/ssytf2.c +81 -0
- data/ext/ssytrd.c +109 -0
- data/ext/ssytrf.c +93 -0
- data/ext/ssytri.c +88 -0
- data/ext/ssytri2.c +123 -0
- data/ext/ssytri2x.c +92 -0
- data/ext/ssytrs.c +99 -0
- data/ext/ssytrs2.c +102 -0
- data/ext/stbcon.c +82 -0
- data/ext/stbrfs.c +123 -0
- data/ext/stbtrs.c +99 -0
- data/ext/stfsm.c +108 -0
- data/ext/stftri.c +82 -0
- data/ext/stfttp.c +75 -0
- data/ext/stfttr.c +76 -0
- data/ext/stgevc.c +149 -0
- data/ext/stgex2.c +180 -0
- data/ext/stgexc.c +187 -0
- data/ext/stgsen.c +248 -0
- data/ext/stgsja.c +223 -0
- data/ext/stgsna.c +160 -0
- data/ext/stgsy2.c +178 -0
- data/ext/stgsyl.c +186 -0
- data/ext/stpcon.c +78 -0
- data/ext/stprfs.c +119 -0
- data/ext/stptri.c +78 -0
- data/ext/stptrs.c +97 -0
- data/ext/stpttf.c +75 -0
- data/ext/stpttr.c +72 -0
- data/ext/strcon.c +78 -0
- data/ext/strevc.c +146 -0
- data/ext/strexc.c +112 -0
- data/ext/strrfs.c +119 -0
- data/ext/strsen.c +165 -0
- data/ext/strsna.c +133 -0
- data/ext/strsyl.c +112 -0
- data/ext/strti2.c +77 -0
- data/ext/strtri.c +77 -0
- data/ext/strtrs.c +95 -0
- data/ext/strttf.c +73 -0
- data/ext/strttp.c +69 -0
- data/ext/stzrqf.c +79 -0
- data/ext/stzrzf.c +97 -0
- data/ext/xerbla.c +48 -0
- data/ext/xerbla_array.c +49 -0
- data/ext/zbbcsd.c +279 -0
- data/ext/zbdsqr.c +178 -0
- data/ext/zcgesv.c +114 -0
- data/ext/zcposv.c +110 -0
- data/ext/zdrscl.c +75 -0
- data/ext/zgbbrd.c +153 -0
- data/ext/zgbcon.c +94 -0
- data/ext/zgbequ.c +94 -0
- data/ext/zgbequb.c +92 -0
- data/ext/zgbrfs.c +157 -0
- data/ext/zgbrfsx.c +245 -0
- data/ext/zgbsv.c +111 -0
- data/ext/zgbsvx.c +282 -0
- data/ext/zgbsvxx.c +285 -0
- data/ext/zgbtf2.c +89 -0
- data/ext/zgbtrf.c +89 -0
- data/ext/zgbtrs.c +107 -0
- data/ext/zgebak.c +97 -0
- data/ext/zgebal.c +87 -0
- data/ext/zgebd2.c +108 -0
- data/ext/zgebrd.c +123 -0
- data/ext/zgecon.c +74 -0
- data/ext/zgeequ.c +84 -0
- data/ext/zgeequb.c +84 -0
- data/ext/zgees.c +138 -0
- data/ext/zgeesx.c +148 -0
- data/ext/zgeev.c +128 -0
- data/ext/zgeevx.c +169 -0
- data/ext/zgegs.c +162 -0
- data/ext/zgegv.c +167 -0
- data/ext/zgehd2.c +88 -0
- data/ext/zgehrd.c +103 -0
- data/ext/zgelq2.c +82 -0
- data/ext/zgelqf.c +99 -0
- data/ext/zgels.c +133 -0
- data/ext/zgelsd.c +150 -0
- data/ext/zgelss.c +147 -0
- data/ext/zgelsx.c +135 -0
- data/ext/zgelsy.c +162 -0
- data/ext/zgeql2.c +84 -0
- data/ext/zgeqlf.c +99 -0
- data/ext/zgeqp3.c +125 -0
- data/ext/zgeqpf.c +110 -0
- data/ext/zgeqr2.c +84 -0
- data/ext/zgeqr2p.c +84 -0
- data/ext/zgeqrf.c +99 -0
- data/ext/zgeqrfp.c +99 -0
- data/ext/zgerfs.c +149 -0
- data/ext/zgerfsx.c +215 -0
- data/ext/zgerq2.c +82 -0
- data/ext/zgerqf.c +99 -0
- data/ext/zgesc2.c +104 -0
- data/ext/zgesdd.c +131 -0
- data/ext/zgesv.c +103 -0
- data/ext/zgesvd.c +142 -0
- data/ext/zgesvx.c +274 -0
- data/ext/zgesvxx.c +277 -0
- data/ext/zgetc2.c +85 -0
- data/ext/zgetf2.c +81 -0
- data/ext/zgetrf.c +81 -0
- data/ext/zgetri.c +99 -0
- data/ext/zgetrs.c +99 -0
- data/ext/zggbak.c +109 -0
- data/ext/zggbal.c +124 -0
- data/ext/zgges.c +188 -0
- data/ext/zggesx.c +226 -0
- data/ext/zggev.c +167 -0
- data/ext/zggevx.c +222 -0
- data/ext/zggglm.c +152 -0
- data/ext/zgghrd.c +163 -0
- data/ext/zgglse.c +167 -0
- data/ext/zggqrf.c +133 -0
- data/ext/zggrqf.c +137 -0
- data/ext/zggsvd.c +180 -0
- data/ext/zggsvp.c +170 -0
- data/ext/zgtcon.c +117 -0
- data/ext/zgtrfs.c +205 -0
- data/ext/zgtsv.c +138 -0
- data/ext/zgtsvx.c +252 -0
- data/ext/zgttrf.c +128 -0
- data/ext/zgttrs.c +133 -0
- data/ext/zgtts2.c +130 -0
- data/ext/zhbev.c +106 -0
- data/ext/zhbevd.c +154 -0
- data/ext/zhbevx.c +156 -0
- data/ext/zhbgst.c +116 -0
- data/ext/zhbgv.c +136 -0
- data/ext/zhbgvd.c +184 -0
- data/ext/zhbgvx.c +185 -0
- data/ext/zhbtrd.c +126 -0
- data/ext/zhecon.c +83 -0
- data/ext/zheequb.c +78 -0
- data/ext/zheev.c +106 -0
- data/ext/zheevd.c +139 -0
- data/ext/zheevr.c +186 -0
- data/ext/zheevx.c +156 -0
- data/ext/zhegs2.c +91 -0
- data/ext/zhegst.c +91 -0
- data/ext/zhegv.c +136 -0
- data/ext/zhegvd.c +169 -0
- data/ext/zhegvx.c +186 -0
- data/ext/zherfs.c +149 -0
- data/ext/zherfsx.c +214 -0
- data/ext/zhesv.c +119 -0
- data/ext/zhesvx.c +179 -0
- data/ext/zhesvxx.c +254 -0
- data/ext/zhetd2.c +97 -0
- data/ext/zhetf2.c +81 -0
- data/ext/zhetrd.c +109 -0
- data/ext/zhetrf.c +93 -0
- data/ext/zhetri.c +88 -0
- data/ext/zhetrs.c +99 -0
- data/ext/zhetrs2.c +102 -0
- data/ext/zhfrk.c +105 -0
- data/ext/zhgeqz.c +204 -0
- data/ext/zhpcon.c +81 -0
- data/ext/zhpev.c +101 -0
- data/ext/zhpevd.c +149 -0
- data/ext/zhpevx.c +140 -0
- data/ext/zhpgst.c +90 -0
- data/ext/zhpgv.c +128 -0
- data/ext/zhpgvd.c +166 -0
- data/ext/zhpgvx.c +166 -0
- data/ext/zhprfs.c +145 -0
- data/ext/zhpsv.c +106 -0
- data/ext/zhpsvx.c +159 -0
- data/ext/zhptrd.c +96 -0
- data/ext/zhptrf.c +80 -0
- data/ext/zhptri.c +85 -0
- data/ext/zhptrs.c +97 -0
- data/ext/zhsein.c +181 -0
- data/ext/zhseqr.c +141 -0
- data/ext/zla_gbamv.c +123 -0
- data/ext/zla_gbrcond_c.c +138 -0
- data/ext/zla_gbrcond_x.c +134 -0
- data/ext/zla_gbrfsx_extended.c +291 -0
- data/ext/zla_gbrpvgrw.c +83 -0
- data/ext/zla_geamv.c +115 -0
- data/ext/zla_gercond_c.c +130 -0
- data/ext/zla_gercond_x.c +126 -0
- data/ext/zla_gerfsx_extended.c +277 -0
- data/ext/zla_heamv.c +112 -0
- data/ext/zla_hercond_c.c +130 -0
- data/ext/zla_hercond_x.c +126 -0
- data/ext/zla_herfsx_extended.c +279 -0
- data/ext/zla_herpvgrw.c +103 -0
- data/ext/zla_lin_berr.c +80 -0
- data/ext/zla_porcond_c.c +118 -0
- data/ext/zla_porcond_x.c +114 -0
- data/ext/zla_porfsx_extended.c +267 -0
- data/ext/zla_porpvgrw.c +91 -0
- data/ext/zla_rpvgrw.c +75 -0
- data/ext/zla_syamv.c +112 -0
- data/ext/zla_syrcond_c.c +130 -0
- data/ext/zla_syrcond_x.c +126 -0
- data/ext/zla_syrfsx_extended.c +279 -0
- data/ext/zla_syrpvgrw.c +103 -0
- data/ext/zla_wwaddw.c +98 -0
- data/ext/zlabrd.c +128 -0
- data/ext/zlacgv.c +71 -0
- data/ext/zlacn2.c +99 -0
- data/ext/zlacon.c +76 -0
- data/ext/zlacp2.c +73 -0
- data/ext/zlacpy.c +73 -0
- data/ext/zlacrm.c +86 -0
- data/ext/zlacrt.c +104 -0
- data/ext/zladiv.c +53 -0
- data/ext/zlaed0.c +130 -0
- data/ext/zlaed7.c +243 -0
- data/ext/zlaed8.c +194 -0
- data/ext/zlaein.c +109 -0
- data/ext/zlaesy.c +70 -0
- data/ext/zlaev2.c +67 -0
- data/ext/zlag2c.c +72 -0
- data/ext/zlags2.c +88 -0
- data/ext/zlagtm.c +128 -0
- data/ext/zlahef.c +93 -0
- data/ext/zlahqr.c +131 -0
- data/ext/zlahr2.c +108 -0
- data/ext/zlahrd.c +108 -0
- data/ext/zlaic1.c +86 -0
- data/ext/zlals0.c +197 -0
- data/ext/zlalsa.c +266 -0
- data/ext/zlalsd.c +141 -0
- data/ext/zlangb.c +72 -0
- data/ext/zlange.c +70 -0
- data/ext/zlangt.c +83 -0
- data/ext/zlanhb.c +74 -0
- data/ext/zlanhe.c +70 -0
- data/ext/zlanhf.c +76 -0
- data/ext/zlanhp.c +72 -0
- data/ext/zlanhs.c +66 -0
- data/ext/zlanht.c +71 -0
- data/ext/zlansb.c +74 -0
- data/ext/zlansp.c +72 -0
- data/ext/zlansy.c +70 -0
- data/ext/zlantb.c +78 -0
- data/ext/zlantp.c +76 -0
- data/ext/zlantr.c +78 -0
- data/ext/zlapll.c +101 -0
- data/ext/zlapmr.c +93 -0
- data/ext/zlapmt.c +97 -0
- data/ext/zlaqgb.c +113 -0
- data/ext/zlaqge.c +105 -0
- data/ext/zlaqhb.c +93 -0
- data/ext/zlaqhe.c +93 -0
- data/ext/zlaqhp.c +90 -0
- data/ext/zlaqp2.c +154 -0
- data/ext/zlaqps.c +204 -0
- data/ext/zlaqr0.c +149 -0
- data/ext/zlaqr1.c +72 -0
- data/ext/zlaqr2.c +170 -0
- data/ext/zlaqr3.c +170 -0
- data/ext/zlaqr4.c +143 -0
- data/ext/zlaqr5.c +175 -0
- data/ext/zlaqsb.c +97 -0
- data/ext/zlaqsp.c +90 -0
- data/ext/zlaqsy.c +93 -0
- data/ext/zlar1v.c +169 -0
- data/ext/zlar2v.c +145 -0
- data/ext/zlarcm.c +82 -0
- data/ext/zlarf.c +98 -0
- data/ext/zlarfb.c +119 -0
- data/ext/zlarfg.c +80 -0
- data/ext/zlarfgp.c +80 -0
- data/ext/zlarft.c +101 -0
- data/ext/zlarfx.c +90 -0
- data/ext/zlargv.c +110 -0
- data/ext/zlarnv.c +79 -0
- data/ext/zlarrv.c +267 -0
- data/ext/zlarscl2.c +78 -0
- data/ext/zlartg.c +59 -0
- data/ext/zlartv.c +126 -0
- data/ext/zlarz.c +102 -0
- data/ext/zlarzb.c +123 -0
- data/ext/zlarzt.c +101 -0
- data/ext/zlascl.c +93 -0
- data/ext/zlascl2.c +78 -0
- data/ext/zlaset.c +84 -0
- data/ext/zlasr.c +106 -0
- data/ext/zlassq.c +66 -0
- data/ext/zlaswp.c +90 -0
- data/ext/zlasyf.c +93 -0
- data/ext/zlat2c.c +72 -0
- data/ext/zlatbs.c +126 -0
- data/ext/zlatdf.c +115 -0
- data/ext/zlatps.c +120 -0
- data/ext/zlatrd.c +101 -0
- data/ext/zlatrs.c +122 -0
- data/ext/zlatrz.c +83 -0
- data/ext/zlatzm.c +128 -0
- data/ext/zlauu2.c +73 -0
- data/ext/zlauum.c +73 -0
- data/ext/zpbcon.c +78 -0
- data/ext/zpbequ.c +79 -0
- data/ext/zpbrfs.c +141 -0
- data/ext/zpbstf.c +77 -0
- data/ext/zpbsv.c +103 -0
- data/ext/zpbsvx.c +197 -0
- data/ext/zpbtf2.c +77 -0
- data/ext/zpbtrf.c +77 -0
- data/ext/zpbtrs.c +91 -0
- data/ext/zpftrf.c +78 -0
- data/ext/zpftri.c +78 -0
- data/ext/zpftrs.c +93 -0
- data/ext/zpocon.c +74 -0
- data/ext/zpoequ.c +71 -0
- data/ext/zpoequb.c +71 -0
- data/ext/zporfs.c +137 -0
- data/ext/zporfsx.c +202 -0
- data/ext/zposv.c +99 -0
- data/ext/zposvx.c +193 -0
- data/ext/zposvxx.c +231 -0
- data/ext/zpotf2.c +73 -0
- data/ext/zpotrf.c +73 -0
- data/ext/zpotri.c +73 -0
- data/ext/zpotrs.c +87 -0
- data/ext/zppcon.c +74 -0
- data/ext/zppequ.c +75 -0
- data/ext/zpprfs.c +135 -0
- data/ext/zppsv.c +100 -0
- data/ext/zppsvx.c +187 -0
- data/ext/zpptrf.c +74 -0
- data/ext/zpptri.c +74 -0
- data/ext/zpptrs.c +89 -0
- data/ext/zpstf2.c +91 -0
- data/ext/zpstrf.c +91 -0
- data/ext/zptcon.c +77 -0
- data/ext/zpteqr.c +122 -0
- data/ext/zptrfs.c +157 -0
- data/ext/zptsv.c +119 -0
- data/ext/zptsvx.c +167 -0
- data/ext/zpttrf.c +89 -0
- data/ext/zpttrs.c +97 -0
- data/ext/zptts2.c +94 -0
- data/ext/zrot.c +103 -0
- data/ext/zspcon.c +81 -0
- data/ext/zspmv.c +113 -0
- data/ext/zspr.c +92 -0
- data/ext/zsprfs.c +145 -0
- data/ext/zspsv.c +106 -0
- data/ext/zspsvx.c +159 -0
- data/ext/zsptrf.c +80 -0
- data/ext/zsptri.c +85 -0
- data/ext/zsptrs.c +97 -0
- data/ext/zstedc.c +173 -0
- data/ext/zstegr.c +184 -0
- data/ext/zstein.c +130 -0
- data/ext/zstemr.c +189 -0
- data/ext/zsteqr.c +122 -0
- data/ext/zsycon.c +83 -0
- data/ext/zsyconv.c +80 -0
- data/ext/zsyequb.c +78 -0
- data/ext/zsymv.c +111 -0
- data/ext/zsyr.c +91 -0
- data/ext/zsyrfs.c +149 -0
- data/ext/zsyrfsx.c +214 -0
- data/ext/zsysv.c +125 -0
- data/ext/zsysvx.c +179 -0
- data/ext/zsysvxx.c +254 -0
- data/ext/zsyswapr.c +78 -0
- data/ext/zsytf2.c +81 -0
- data/ext/zsytrf.c +93 -0
- data/ext/zsytri.c +88 -0
- data/ext/zsytri2.c +100 -0
- data/ext/zsytri2x.c +92 -0
- data/ext/zsytrs.c +99 -0
- data/ext/zsytrs2.c +102 -0
- data/ext/ztbcon.c +82 -0
- data/ext/ztbrfs.c +123 -0
- data/ext/ztbtrs.c +99 -0
- data/ext/ztfsm.c +107 -0
- data/ext/ztftri.c +82 -0
- data/ext/ztfttp.c +75 -0
- data/ext/ztfttr.c +76 -0
- data/ext/ztgevc.c +152 -0
- data/ext/ztgex2.c +167 -0
- data/ext/ztgexc.c +168 -0
- data/ext/ztgsen.c +240 -0
- data/ext/ztgsja.c +223 -0
- data/ext/ztgsna.c +160 -0
- data/ext/ztgsy2.c +172 -0
- data/ext/ztgsyl.c +186 -0
- data/ext/ztpcon.c +78 -0
- data/ext/ztprfs.c +119 -0
- data/ext/ztptri.c +78 -0
- data/ext/ztptrs.c +97 -0
- data/ext/ztpttf.c +75 -0
- data/ext/ztpttr.c +72 -0
- data/ext/ztrcon.c +78 -0
- data/ext/ztrevc.c +150 -0
- data/ext/ztrexc.c +107 -0
- data/ext/ztrrfs.c +119 -0
- data/ext/ztrsen.c +150 -0
- data/ext/ztrsna.c +133 -0
- data/ext/ztrsyl.c +112 -0
- data/ext/ztrti2.c +77 -0
- data/ext/ztrtri.c +77 -0
- data/ext/ztrtrs.c +95 -0
- data/ext/ztrttf.c +73 -0
- data/ext/ztrttp.c +69 -0
- data/ext/ztzrqf.c +79 -0
- data/ext/ztzrzf.c +97 -0
- data/ext/zunbdb.c +228 -0
- data/ext/zuncsd.c +200 -0
- data/ext/zung2l.c +88 -0
- data/ext/zung2r.c +88 -0
- data/ext/zungbr.c +111 -0
- data/ext/zunghr.c +107 -0
- data/ext/zungl2.c +86 -0
- data/ext/zunglq.c +103 -0
- data/ext/zungql.c +103 -0
- data/ext/zungqr.c +103 -0
- data/ext/zungr2.c +86 -0
- data/ext/zungrq.c +103 -0
- data/ext/zungtr.c +103 -0
- data/ext/zunm2l.c +110 -0
- data/ext/zunm2r.c +110 -0
- data/ext/zunmbr.c +135 -0
- data/ext/zunmhr.c +129 -0
- data/ext/zunml2.c +106 -0
- data/ext/zunmlq.c +121 -0
- data/ext/zunmql.c +125 -0
- data/ext/zunmqr.c +125 -0
- data/ext/zunmr2.c +106 -0
- data/ext/zunmr3.c +110 -0
- data/ext/zunmrq.c +121 -0
- data/ext/zunmrz.c +125 -0
- data/ext/zunmtr.c +125 -0
- data/ext/zupgtr.c +87 -0
- data/ext/zupmtr.c +112 -0
- data/lib/numru/lapack.rb +51 -0
- data/samples/dsyevr.rb +25 -0
- data/tests/eig/ge/test_gesdd.rb +90 -0
- data/tests/eig/ge/test_gesvd.rb +99 -0
- data/tests/eig/gg/test_ggev.rb +124 -0
- data/tests/eig/gg/test_ggsvd.rb +76 -0
- data/tests/eig/sb/test_sbev.rb +39 -0
- data/tests/lapack_test.rb +50 -0
- data/tests/lin/gb/test_gbsv.rb +46 -0
- data/tests/lin/gb/test_gbsvx.rb +56 -0
- data/tests/lin/ge/test_gels.rb +63 -0
- data/tests/lin/ge/test_gelsd.rb +54 -0
- data/tests/lin/ge/test_gelss.rb +73 -0
- data/tests/lin/ge/test_gelsy.rb +73 -0
- data/tests/lin/ge/test_gesv.rb +43 -0
- data/tests/lin/ge/test_gesvx.rb +52 -0
- data/tests/lin/gt/test_gtsv.rb +39 -0
- data/tests/test_all.rb +7 -0
- metadata +3513 -0
data/doc/str.html
ADDED
|
@@ -0,0 +1,1346 @@
|
|
|
1
|
+
<HTML>
|
|
2
|
+
<HEAD>
|
|
3
|
+
<TITLE>REAL routines for triangular (or in some cases quasi-triangular) matrix</TITLE>
|
|
4
|
+
</HEAD>
|
|
5
|
+
<BODY>
|
|
6
|
+
<A NAME="top"></A>
|
|
7
|
+
<H1>REAL routines for triangular (or in some cases quasi-triangular) matrix</H1>
|
|
8
|
+
<UL>
|
|
9
|
+
<LI><A HREF="#strcon">strcon</A></LI>
|
|
10
|
+
<LI><A HREF="#strevc">strevc</A></LI>
|
|
11
|
+
<LI><A HREF="#strexc">strexc</A></LI>
|
|
12
|
+
<LI><A HREF="#strrfs">strrfs</A></LI>
|
|
13
|
+
<LI><A HREF="#strsen">strsen</A></LI>
|
|
14
|
+
<LI><A HREF="#strsna">strsna</A></LI>
|
|
15
|
+
<LI><A HREF="#strsyl">strsyl</A></LI>
|
|
16
|
+
<LI><A HREF="#strti2">strti2</A></LI>
|
|
17
|
+
<LI><A HREF="#strtri">strtri</A></LI>
|
|
18
|
+
<LI><A HREF="#strtrs">strtrs</A></LI>
|
|
19
|
+
<LI><A HREF="#strttf">strttf</A></LI>
|
|
20
|
+
<LI><A HREF="#strttp">strttp</A></LI>
|
|
21
|
+
</UL>
|
|
22
|
+
|
|
23
|
+
<A NAME="strcon"></A>
|
|
24
|
+
<H2>strcon</H2>
|
|
25
|
+
<PRE>
|
|
26
|
+
USAGE:
|
|
27
|
+
rcond, info = NumRu::Lapack.strcon( norm, uplo, diag, a, [:usage => usage, :help => help])
|
|
28
|
+
|
|
29
|
+
|
|
30
|
+
FORTRAN MANUAL
|
|
31
|
+
SUBROUTINE STRCON( NORM, UPLO, DIAG, N, A, LDA, RCOND, WORK, IWORK, INFO )
|
|
32
|
+
|
|
33
|
+
* Purpose
|
|
34
|
+
* =======
|
|
35
|
+
*
|
|
36
|
+
* STRCON estimates the reciprocal of the condition number of a
|
|
37
|
+
* triangular matrix A, in either the 1-norm or the infinity-norm.
|
|
38
|
+
*
|
|
39
|
+
* The norm of A is computed and an estimate is obtained for
|
|
40
|
+
* norm(inv(A)), then the reciprocal of the condition number is
|
|
41
|
+
* computed as
|
|
42
|
+
* RCOND = 1 / ( norm(A) * norm(inv(A)) ).
|
|
43
|
+
*
|
|
44
|
+
|
|
45
|
+
* Arguments
|
|
46
|
+
* =========
|
|
47
|
+
*
|
|
48
|
+
* NORM (input) CHARACTER*1
|
|
49
|
+
* Specifies whether the 1-norm condition number or the
|
|
50
|
+
* infinity-norm condition number is required:
|
|
51
|
+
* = '1' or 'O': 1-norm;
|
|
52
|
+
* = 'I': Infinity-norm.
|
|
53
|
+
*
|
|
54
|
+
* UPLO (input) CHARACTER*1
|
|
55
|
+
* = 'U': A is upper triangular;
|
|
56
|
+
* = 'L': A is lower triangular.
|
|
57
|
+
*
|
|
58
|
+
* DIAG (input) CHARACTER*1
|
|
59
|
+
* = 'N': A is non-unit triangular;
|
|
60
|
+
* = 'U': A is unit triangular.
|
|
61
|
+
*
|
|
62
|
+
* N (input) INTEGER
|
|
63
|
+
* The order of the matrix A. N >= 0.
|
|
64
|
+
*
|
|
65
|
+
* A (input) REAL array, dimension (LDA,N)
|
|
66
|
+
* The triangular matrix A. If UPLO = 'U', the leading N-by-N
|
|
67
|
+
* upper triangular part of the array A contains the upper
|
|
68
|
+
* triangular matrix, and the strictly lower triangular part of
|
|
69
|
+
* A is not referenced. If UPLO = 'L', the leading N-by-N lower
|
|
70
|
+
* triangular part of the array A contains the lower triangular
|
|
71
|
+
* matrix, and the strictly upper triangular part of A is not
|
|
72
|
+
* referenced. If DIAG = 'U', the diagonal elements of A are
|
|
73
|
+
* also not referenced and are assumed to be 1.
|
|
74
|
+
*
|
|
75
|
+
* LDA (input) INTEGER
|
|
76
|
+
* The leading dimension of the array A. LDA >= max(1,N).
|
|
77
|
+
*
|
|
78
|
+
* RCOND (output) REAL
|
|
79
|
+
* The reciprocal of the condition number of the matrix A,
|
|
80
|
+
* computed as RCOND = 1/(norm(A) * norm(inv(A))).
|
|
81
|
+
*
|
|
82
|
+
* WORK (workspace) REAL array, dimension (3*N)
|
|
83
|
+
*
|
|
84
|
+
* IWORK (workspace) INTEGER array, dimension (N)
|
|
85
|
+
*
|
|
86
|
+
* INFO (output) INTEGER
|
|
87
|
+
* = 0: successful exit
|
|
88
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
89
|
+
*
|
|
90
|
+
|
|
91
|
+
* =====================================================================
|
|
92
|
+
*
|
|
93
|
+
|
|
94
|
+
|
|
95
|
+
</PRE>
|
|
96
|
+
<A HREF="#top">go to the page top</A>
|
|
97
|
+
|
|
98
|
+
<A NAME="strevc"></A>
|
|
99
|
+
<H2>strevc</H2>
|
|
100
|
+
<PRE>
|
|
101
|
+
USAGE:
|
|
102
|
+
m, info, select, vl, vr = NumRu::Lapack.strevc( side, howmny, select, t, vl, vr, [:usage => usage, :help => help])
|
|
103
|
+
|
|
104
|
+
|
|
105
|
+
FORTRAN MANUAL
|
|
106
|
+
SUBROUTINE STREVC( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR, LDVR, MM, M, WORK, INFO )
|
|
107
|
+
|
|
108
|
+
* Purpose
|
|
109
|
+
* =======
|
|
110
|
+
*
|
|
111
|
+
* STREVC computes some or all of the right and/or left eigenvectors of
|
|
112
|
+
* a real upper quasi-triangular matrix T.
|
|
113
|
+
* Matrices of this type are produced by the Schur factorization of
|
|
114
|
+
* a real general matrix: A = Q*T*Q**T, as computed by SHSEQR.
|
|
115
|
+
*
|
|
116
|
+
* The right eigenvector x and the left eigenvector y of T corresponding
|
|
117
|
+
* to an eigenvalue w are defined by:
|
|
118
|
+
*
|
|
119
|
+
* T*x = w*x, (y**H)*T = w*(y**H)
|
|
120
|
+
*
|
|
121
|
+
* where y**H denotes the conjugate transpose of y.
|
|
122
|
+
* The eigenvalues are not input to this routine, but are read directly
|
|
123
|
+
* from the diagonal blocks of T.
|
|
124
|
+
*
|
|
125
|
+
* This routine returns the matrices X and/or Y of right and left
|
|
126
|
+
* eigenvectors of T, or the products Q*X and/or Q*Y, where Q is an
|
|
127
|
+
* input matrix. If Q is the orthogonal factor that reduces a matrix
|
|
128
|
+
* A to Schur form T, then Q*X and Q*Y are the matrices of right and
|
|
129
|
+
* left eigenvectors of A.
|
|
130
|
+
*
|
|
131
|
+
|
|
132
|
+
* Arguments
|
|
133
|
+
* =========
|
|
134
|
+
*
|
|
135
|
+
* SIDE (input) CHARACTER*1
|
|
136
|
+
* = 'R': compute right eigenvectors only;
|
|
137
|
+
* = 'L': compute left eigenvectors only;
|
|
138
|
+
* = 'B': compute both right and left eigenvectors.
|
|
139
|
+
*
|
|
140
|
+
* HOWMNY (input) CHARACTER*1
|
|
141
|
+
* = 'A': compute all right and/or left eigenvectors;
|
|
142
|
+
* = 'B': compute all right and/or left eigenvectors,
|
|
143
|
+
* backtransformed by the matrices in VR and/or VL;
|
|
144
|
+
* = 'S': compute selected right and/or left eigenvectors,
|
|
145
|
+
* as indicated by the logical array SELECT.
|
|
146
|
+
*
|
|
147
|
+
* SELECT (input/output) LOGICAL array, dimension (N)
|
|
148
|
+
* If HOWMNY = 'S', SELECT specifies the eigenvectors to be
|
|
149
|
+
* computed.
|
|
150
|
+
* If w(j) is a real eigenvalue, the corresponding real
|
|
151
|
+
* eigenvector is computed if SELECT(j) is .TRUE..
|
|
152
|
+
* If w(j) and w(j+1) are the real and imaginary parts of a
|
|
153
|
+
* complex eigenvalue, the corresponding complex eigenvector is
|
|
154
|
+
* computed if either SELECT(j) or SELECT(j+1) is .TRUE., and
|
|
155
|
+
* on exit SELECT(j) is set to .TRUE. and SELECT(j+1) is set to
|
|
156
|
+
* .FALSE..
|
|
157
|
+
* Not referenced if HOWMNY = 'A' or 'B'.
|
|
158
|
+
*
|
|
159
|
+
* N (input) INTEGER
|
|
160
|
+
* The order of the matrix T. N >= 0.
|
|
161
|
+
*
|
|
162
|
+
* T (input) REAL array, dimension (LDT,N)
|
|
163
|
+
* The upper quasi-triangular matrix T in Schur canonical form.
|
|
164
|
+
*
|
|
165
|
+
* LDT (input) INTEGER
|
|
166
|
+
* The leading dimension of the array T. LDT >= max(1,N).
|
|
167
|
+
*
|
|
168
|
+
* VL (input/output) REAL array, dimension (LDVL,MM)
|
|
169
|
+
* On entry, if SIDE = 'L' or 'B' and HOWMNY = 'B', VL must
|
|
170
|
+
* contain an N-by-N matrix Q (usually the orthogonal matrix Q
|
|
171
|
+
* of Schur vectors returned by SHSEQR).
|
|
172
|
+
* On exit, if SIDE = 'L' or 'B', VL contains:
|
|
173
|
+
* if HOWMNY = 'A', the matrix Y of left eigenvectors of T;
|
|
174
|
+
* if HOWMNY = 'B', the matrix Q*Y;
|
|
175
|
+
* if HOWMNY = 'S', the left eigenvectors of T specified by
|
|
176
|
+
* SELECT, stored consecutively in the columns
|
|
177
|
+
* of VL, in the same order as their
|
|
178
|
+
* eigenvalues.
|
|
179
|
+
* A complex eigenvector corresponding to a complex eigenvalue
|
|
180
|
+
* is stored in two consecutive columns, the first holding the
|
|
181
|
+
* real part, and the second the imaginary part.
|
|
182
|
+
* Not referenced if SIDE = 'R'.
|
|
183
|
+
*
|
|
184
|
+
* LDVL (input) INTEGER
|
|
185
|
+
* The leading dimension of the array VL. LDVL >= 1, and if
|
|
186
|
+
* SIDE = 'L' or 'B', LDVL >= N.
|
|
187
|
+
*
|
|
188
|
+
* VR (input/output) REAL array, dimension (LDVR,MM)
|
|
189
|
+
* On entry, if SIDE = 'R' or 'B' and HOWMNY = 'B', VR must
|
|
190
|
+
* contain an N-by-N matrix Q (usually the orthogonal matrix Q
|
|
191
|
+
* of Schur vectors returned by SHSEQR).
|
|
192
|
+
* On exit, if SIDE = 'R' or 'B', VR contains:
|
|
193
|
+
* if HOWMNY = 'A', the matrix X of right eigenvectors of T;
|
|
194
|
+
* if HOWMNY = 'B', the matrix Q*X;
|
|
195
|
+
* if HOWMNY = 'S', the right eigenvectors of T specified by
|
|
196
|
+
* SELECT, stored consecutively in the columns
|
|
197
|
+
* of VR, in the same order as their
|
|
198
|
+
* eigenvalues.
|
|
199
|
+
* A complex eigenvector corresponding to a complex eigenvalue
|
|
200
|
+
* is stored in two consecutive columns, the first holding the
|
|
201
|
+
* real part and the second the imaginary part.
|
|
202
|
+
* Not referenced if SIDE = 'L'.
|
|
203
|
+
*
|
|
204
|
+
* LDVR (input) INTEGER
|
|
205
|
+
* The leading dimension of the array VR. LDVR >= 1, and if
|
|
206
|
+
* SIDE = 'R' or 'B', LDVR >= N.
|
|
207
|
+
*
|
|
208
|
+
* MM (input) INTEGER
|
|
209
|
+
* The number of columns in the arrays VL and/or VR. MM >= M.
|
|
210
|
+
*
|
|
211
|
+
* M (output) INTEGER
|
|
212
|
+
* The number of columns in the arrays VL and/or VR actually
|
|
213
|
+
* used to store the eigenvectors.
|
|
214
|
+
* If HOWMNY = 'A' or 'B', M is set to N.
|
|
215
|
+
* Each selected real eigenvector occupies one column and each
|
|
216
|
+
* selected complex eigenvector occupies two columns.
|
|
217
|
+
*
|
|
218
|
+
* WORK (workspace) REAL array, dimension (3*N)
|
|
219
|
+
*
|
|
220
|
+
* INFO (output) INTEGER
|
|
221
|
+
* = 0: successful exit
|
|
222
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
223
|
+
*
|
|
224
|
+
|
|
225
|
+
* Further Details
|
|
226
|
+
* ===============
|
|
227
|
+
*
|
|
228
|
+
* The algorithm used in this program is basically backward (forward)
|
|
229
|
+
* substitution, with scaling to make the the code robust against
|
|
230
|
+
* possible overflow.
|
|
231
|
+
*
|
|
232
|
+
* Each eigenvector is normalized so that the element of largest
|
|
233
|
+
* magnitude has magnitude 1; here the magnitude of a complex number
|
|
234
|
+
* (x,y) is taken to be |x| + |y|.
|
|
235
|
+
*
|
|
236
|
+
* =====================================================================
|
|
237
|
+
*
|
|
238
|
+
|
|
239
|
+
|
|
240
|
+
</PRE>
|
|
241
|
+
<A HREF="#top">go to the page top</A>
|
|
242
|
+
|
|
243
|
+
<A NAME="strexc"></A>
|
|
244
|
+
<H2>strexc</H2>
|
|
245
|
+
<PRE>
|
|
246
|
+
USAGE:
|
|
247
|
+
info, t, q, ifst, ilst = NumRu::Lapack.strexc( compq, t, q, ifst, ilst, [:usage => usage, :help => help])
|
|
248
|
+
|
|
249
|
+
|
|
250
|
+
FORTRAN MANUAL
|
|
251
|
+
SUBROUTINE STREXC( COMPQ, N, T, LDT, Q, LDQ, IFST, ILST, WORK, INFO )
|
|
252
|
+
|
|
253
|
+
* Purpose
|
|
254
|
+
* =======
|
|
255
|
+
*
|
|
256
|
+
* STREXC reorders the real Schur factorization of a real matrix
|
|
257
|
+
* A = Q*T*Q**T, so that the diagonal block of T with row index IFST is
|
|
258
|
+
* moved to row ILST.
|
|
259
|
+
*
|
|
260
|
+
* The real Schur form T is reordered by an orthogonal similarity
|
|
261
|
+
* transformation Z**T*T*Z, and optionally the matrix Q of Schur vectors
|
|
262
|
+
* is updated by postmultiplying it with Z.
|
|
263
|
+
*
|
|
264
|
+
* T must be in Schur canonical form (as returned by SHSEQR), that is,
|
|
265
|
+
* block upper triangular with 1-by-1 and 2-by-2 diagonal blocks; each
|
|
266
|
+
* 2-by-2 diagonal block has its diagonal elements equal and its
|
|
267
|
+
* off-diagonal elements of opposite sign.
|
|
268
|
+
*
|
|
269
|
+
|
|
270
|
+
* Arguments
|
|
271
|
+
* =========
|
|
272
|
+
*
|
|
273
|
+
* COMPQ (input) CHARACTER*1
|
|
274
|
+
* = 'V': update the matrix Q of Schur vectors;
|
|
275
|
+
* = 'N': do not update Q.
|
|
276
|
+
*
|
|
277
|
+
* N (input) INTEGER
|
|
278
|
+
* The order of the matrix T. N >= 0.
|
|
279
|
+
*
|
|
280
|
+
* T (input/output) REAL array, dimension (LDT,N)
|
|
281
|
+
* On entry, the upper quasi-triangular matrix T, in Schur
|
|
282
|
+
* Schur canonical form.
|
|
283
|
+
* On exit, the reordered upper quasi-triangular matrix, again
|
|
284
|
+
* in Schur canonical form.
|
|
285
|
+
*
|
|
286
|
+
* LDT (input) INTEGER
|
|
287
|
+
* The leading dimension of the array T. LDT >= max(1,N).
|
|
288
|
+
*
|
|
289
|
+
* Q (input/output) REAL array, dimension (LDQ,N)
|
|
290
|
+
* On entry, if COMPQ = 'V', the matrix Q of Schur vectors.
|
|
291
|
+
* On exit, if COMPQ = 'V', Q has been postmultiplied by the
|
|
292
|
+
* orthogonal transformation matrix Z which reorders T.
|
|
293
|
+
* If COMPQ = 'N', Q is not referenced.
|
|
294
|
+
*
|
|
295
|
+
* LDQ (input) INTEGER
|
|
296
|
+
* The leading dimension of the array Q. LDQ >= max(1,N).
|
|
297
|
+
*
|
|
298
|
+
* IFST (input/output) INTEGER
|
|
299
|
+
* ILST (input/output) INTEGER
|
|
300
|
+
* Specify the reordering of the diagonal blocks of T.
|
|
301
|
+
* The block with row index IFST is moved to row ILST, by a
|
|
302
|
+
* sequence of transpositions between adjacent blocks.
|
|
303
|
+
* On exit, if IFST pointed on entry to the second row of a
|
|
304
|
+
* 2-by-2 block, it is changed to point to the first row; ILST
|
|
305
|
+
* always points to the first row of the block in its final
|
|
306
|
+
* position (which may differ from its input value by +1 or -1).
|
|
307
|
+
* 1 <= IFST <= N; 1 <= ILST <= N.
|
|
308
|
+
*
|
|
309
|
+
* WORK (workspace) REAL array, dimension (N)
|
|
310
|
+
*
|
|
311
|
+
* INFO (output) INTEGER
|
|
312
|
+
* = 0: successful exit
|
|
313
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
314
|
+
* = 1: two adjacent blocks were too close to swap (the problem
|
|
315
|
+
* is very ill-conditioned); T may have been partially
|
|
316
|
+
* reordered, and ILST points to the first row of the
|
|
317
|
+
* current position of the block being moved.
|
|
318
|
+
*
|
|
319
|
+
|
|
320
|
+
* =====================================================================
|
|
321
|
+
*
|
|
322
|
+
|
|
323
|
+
|
|
324
|
+
</PRE>
|
|
325
|
+
<A HREF="#top">go to the page top</A>
|
|
326
|
+
|
|
327
|
+
<A NAME="strrfs"></A>
|
|
328
|
+
<H2>strrfs</H2>
|
|
329
|
+
<PRE>
|
|
330
|
+
USAGE:
|
|
331
|
+
ferr, berr, info = NumRu::Lapack.strrfs( uplo, trans, diag, a, b, x, [:usage => usage, :help => help])
|
|
332
|
+
|
|
333
|
+
|
|
334
|
+
FORTRAN MANUAL
|
|
335
|
+
SUBROUTINE STRRFS( UPLO, TRANS, DIAG, N, NRHS, A, LDA, B, LDB, X, LDX, FERR, BERR, WORK, IWORK, INFO )
|
|
336
|
+
|
|
337
|
+
* Purpose
|
|
338
|
+
* =======
|
|
339
|
+
*
|
|
340
|
+
* STRRFS provides error bounds and backward error estimates for the
|
|
341
|
+
* solution to a system of linear equations with a triangular
|
|
342
|
+
* coefficient matrix.
|
|
343
|
+
*
|
|
344
|
+
* The solution matrix X must be computed by STRTRS or some other
|
|
345
|
+
* means before entering this routine. STRRFS does not do iterative
|
|
346
|
+
* refinement because doing so cannot improve the backward error.
|
|
347
|
+
*
|
|
348
|
+
|
|
349
|
+
* Arguments
|
|
350
|
+
* =========
|
|
351
|
+
*
|
|
352
|
+
* UPLO (input) CHARACTER*1
|
|
353
|
+
* = 'U': A is upper triangular;
|
|
354
|
+
* = 'L': A is lower triangular.
|
|
355
|
+
*
|
|
356
|
+
* TRANS (input) CHARACTER*1
|
|
357
|
+
* Specifies the form of the system of equations:
|
|
358
|
+
* = 'N': A * X = B (No transpose)
|
|
359
|
+
* = 'T': A**T * X = B (Transpose)
|
|
360
|
+
* = 'C': A**H * X = B (Conjugate transpose = Transpose)
|
|
361
|
+
*
|
|
362
|
+
* DIAG (input) CHARACTER*1
|
|
363
|
+
* = 'N': A is non-unit triangular;
|
|
364
|
+
* = 'U': A is unit triangular.
|
|
365
|
+
*
|
|
366
|
+
* N (input) INTEGER
|
|
367
|
+
* The order of the matrix A. N >= 0.
|
|
368
|
+
*
|
|
369
|
+
* NRHS (input) INTEGER
|
|
370
|
+
* The number of right hand sides, i.e., the number of columns
|
|
371
|
+
* of the matrices B and X. NRHS >= 0.
|
|
372
|
+
*
|
|
373
|
+
* A (input) REAL array, dimension (LDA,N)
|
|
374
|
+
* The triangular matrix A. If UPLO = 'U', the leading N-by-N
|
|
375
|
+
* upper triangular part of the array A contains the upper
|
|
376
|
+
* triangular matrix, and the strictly lower triangular part of
|
|
377
|
+
* A is not referenced. If UPLO = 'L', the leading N-by-N lower
|
|
378
|
+
* triangular part of the array A contains the lower triangular
|
|
379
|
+
* matrix, and the strictly upper triangular part of A is not
|
|
380
|
+
* referenced. If DIAG = 'U', the diagonal elements of A are
|
|
381
|
+
* also not referenced and are assumed to be 1.
|
|
382
|
+
*
|
|
383
|
+
* LDA (input) INTEGER
|
|
384
|
+
* The leading dimension of the array A. LDA >= max(1,N).
|
|
385
|
+
*
|
|
386
|
+
* B (input) REAL array, dimension (LDB,NRHS)
|
|
387
|
+
* The right hand side matrix B.
|
|
388
|
+
*
|
|
389
|
+
* LDB (input) INTEGER
|
|
390
|
+
* The leading dimension of the array B. LDB >= max(1,N).
|
|
391
|
+
*
|
|
392
|
+
* X (input) REAL array, dimension (LDX,NRHS)
|
|
393
|
+
* The solution matrix X.
|
|
394
|
+
*
|
|
395
|
+
* LDX (input) INTEGER
|
|
396
|
+
* The leading dimension of the array X. LDX >= max(1,N).
|
|
397
|
+
*
|
|
398
|
+
* FERR (output) REAL array, dimension (NRHS)
|
|
399
|
+
* The estimated forward error bound for each solution vector
|
|
400
|
+
* X(j) (the j-th column of the solution matrix X).
|
|
401
|
+
* If XTRUE is the true solution corresponding to X(j), FERR(j)
|
|
402
|
+
* is an estimated upper bound for the magnitude of the largest
|
|
403
|
+
* element in (X(j) - XTRUE) divided by the magnitude of the
|
|
404
|
+
* largest element in X(j). The estimate is as reliable as
|
|
405
|
+
* the estimate for RCOND, and is almost always a slight
|
|
406
|
+
* overestimate of the true error.
|
|
407
|
+
*
|
|
408
|
+
* BERR (output) REAL array, dimension (NRHS)
|
|
409
|
+
* The componentwise relative backward error of each solution
|
|
410
|
+
* vector X(j) (i.e., the smallest relative change in
|
|
411
|
+
* any element of A or B that makes X(j) an exact solution).
|
|
412
|
+
*
|
|
413
|
+
* WORK (workspace) REAL array, dimension (3*N)
|
|
414
|
+
*
|
|
415
|
+
* IWORK (workspace) INTEGER array, dimension (N)
|
|
416
|
+
*
|
|
417
|
+
* INFO (output) INTEGER
|
|
418
|
+
* = 0: successful exit
|
|
419
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
420
|
+
*
|
|
421
|
+
|
|
422
|
+
* =====================================================================
|
|
423
|
+
*
|
|
424
|
+
|
|
425
|
+
|
|
426
|
+
</PRE>
|
|
427
|
+
<A HREF="#top">go to the page top</A>
|
|
428
|
+
|
|
429
|
+
<A NAME="strsen"></A>
|
|
430
|
+
<H2>strsen</H2>
|
|
431
|
+
<PRE>
|
|
432
|
+
USAGE:
|
|
433
|
+
wr, wi, m, s, sep, work, info, t, q = NumRu::Lapack.strsen( job, compq, select, t, q, liwork, [:lwork => lwork, :usage => usage, :help => help])
|
|
434
|
+
|
|
435
|
+
|
|
436
|
+
FORTRAN MANUAL
|
|
437
|
+
SUBROUTINE STRSEN( JOB, COMPQ, SELECT, N, T, LDT, Q, LDQ, WR, WI, M, S, SEP, WORK, LWORK, IWORK, LIWORK, INFO )
|
|
438
|
+
|
|
439
|
+
* Purpose
|
|
440
|
+
* =======
|
|
441
|
+
*
|
|
442
|
+
* STRSEN reorders the real Schur factorization of a real matrix
|
|
443
|
+
* A = Q*T*Q**T, so that a selected cluster of eigenvalues appears in
|
|
444
|
+
* the leading diagonal blocks of the upper quasi-triangular matrix T,
|
|
445
|
+
* and the leading columns of Q form an orthonormal basis of the
|
|
446
|
+
* corresponding right invariant subspace.
|
|
447
|
+
*
|
|
448
|
+
* Optionally the routine computes the reciprocal condition numbers of
|
|
449
|
+
* the cluster of eigenvalues and/or the invariant subspace.
|
|
450
|
+
*
|
|
451
|
+
* T must be in Schur canonical form (as returned by SHSEQR), that is,
|
|
452
|
+
* block upper triangular with 1-by-1 and 2-by-2 diagonal blocks; each
|
|
453
|
+
* 2-by-2 diagonal block has its diagonal elemnts equal and its
|
|
454
|
+
* off-diagonal elements of opposite sign.
|
|
455
|
+
*
|
|
456
|
+
|
|
457
|
+
* Arguments
|
|
458
|
+
* =========
|
|
459
|
+
*
|
|
460
|
+
* JOB (input) CHARACTER*1
|
|
461
|
+
* Specifies whether condition numbers are required for the
|
|
462
|
+
* cluster of eigenvalues (S) or the invariant subspace (SEP):
|
|
463
|
+
* = 'N': none;
|
|
464
|
+
* = 'E': for eigenvalues only (S);
|
|
465
|
+
* = 'V': for invariant subspace only (SEP);
|
|
466
|
+
* = 'B': for both eigenvalues and invariant subspace (S and
|
|
467
|
+
* SEP).
|
|
468
|
+
*
|
|
469
|
+
* COMPQ (input) CHARACTER*1
|
|
470
|
+
* = 'V': update the matrix Q of Schur vectors;
|
|
471
|
+
* = 'N': do not update Q.
|
|
472
|
+
*
|
|
473
|
+
* SELECT (input) LOGICAL array, dimension (N)
|
|
474
|
+
* SELECT specifies the eigenvalues in the selected cluster. To
|
|
475
|
+
* select a real eigenvalue w(j), SELECT(j) must be set to
|
|
476
|
+
* .TRUE.. To select a complex conjugate pair of eigenvalues
|
|
477
|
+
* w(j) and w(j+1), corresponding to a 2-by-2 diagonal block,
|
|
478
|
+
* either SELECT(j) or SELECT(j+1) or both must be set to
|
|
479
|
+
* .TRUE.; a complex conjugate pair of eigenvalues must be
|
|
480
|
+
* either both included in the cluster or both excluded.
|
|
481
|
+
*
|
|
482
|
+
* N (input) INTEGER
|
|
483
|
+
* The order of the matrix T. N >= 0.
|
|
484
|
+
*
|
|
485
|
+
* T (input/output) REAL array, dimension (LDT,N)
|
|
486
|
+
* On entry, the upper quasi-triangular matrix T, in Schur
|
|
487
|
+
* canonical form.
|
|
488
|
+
* On exit, T is overwritten by the reordered matrix T, again in
|
|
489
|
+
* Schur canonical form, with the selected eigenvalues in the
|
|
490
|
+
* leading diagonal blocks.
|
|
491
|
+
*
|
|
492
|
+
* LDT (input) INTEGER
|
|
493
|
+
* The leading dimension of the array T. LDT >= max(1,N).
|
|
494
|
+
*
|
|
495
|
+
* Q (input/output) REAL array, dimension (LDQ,N)
|
|
496
|
+
* On entry, if COMPQ = 'V', the matrix Q of Schur vectors.
|
|
497
|
+
* On exit, if COMPQ = 'V', Q has been postmultiplied by the
|
|
498
|
+
* orthogonal transformation matrix which reorders T; the
|
|
499
|
+
* leading M columns of Q form an orthonormal basis for the
|
|
500
|
+
* specified invariant subspace.
|
|
501
|
+
* If COMPQ = 'N', Q is not referenced.
|
|
502
|
+
*
|
|
503
|
+
* LDQ (input) INTEGER
|
|
504
|
+
* The leading dimension of the array Q.
|
|
505
|
+
* LDQ >= 1; and if COMPQ = 'V', LDQ >= N.
|
|
506
|
+
*
|
|
507
|
+
* WR (output) REAL array, dimension (N)
|
|
508
|
+
* WI (output) REAL array, dimension (N)
|
|
509
|
+
* The real and imaginary parts, respectively, of the reordered
|
|
510
|
+
* eigenvalues of T. The eigenvalues are stored in the same
|
|
511
|
+
* order as on the diagonal of T, with WR(i) = T(i,i) and, if
|
|
512
|
+
* T(i:i+1,i:i+1) is a 2-by-2 diagonal block, WI(i) > 0 and
|
|
513
|
+
* WI(i+1) = -WI(i). Note that if a complex eigenvalue is
|
|
514
|
+
* sufficiently ill-conditioned, then its value may differ
|
|
515
|
+
* significantly from its value before reordering.
|
|
516
|
+
*
|
|
517
|
+
* M (output) INTEGER
|
|
518
|
+
* The dimension of the specified invariant subspace.
|
|
519
|
+
* 0 < = M <= N.
|
|
520
|
+
*
|
|
521
|
+
* S (output) REAL
|
|
522
|
+
* If JOB = 'E' or 'B', S is a lower bound on the reciprocal
|
|
523
|
+
* condition number for the selected cluster of eigenvalues.
|
|
524
|
+
* S cannot underestimate the true reciprocal condition number
|
|
525
|
+
* by more than a factor of sqrt(N). If M = 0 or N, S = 1.
|
|
526
|
+
* If JOB = 'N' or 'V', S is not referenced.
|
|
527
|
+
*
|
|
528
|
+
* SEP (output) REAL
|
|
529
|
+
* If JOB = 'V' or 'B', SEP is the estimated reciprocal
|
|
530
|
+
* condition number of the specified invariant subspace. If
|
|
531
|
+
* M = 0 or N, SEP = norm(T).
|
|
532
|
+
* If JOB = 'N' or 'E', SEP is not referenced.
|
|
533
|
+
*
|
|
534
|
+
* WORK (workspace/output) REAL array, dimension (MAX(1,LWORK))
|
|
535
|
+
* On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
|
|
536
|
+
*
|
|
537
|
+
* LWORK (input) INTEGER
|
|
538
|
+
* The dimension of the array WORK.
|
|
539
|
+
* If JOB = 'N', LWORK >= max(1,N);
|
|
540
|
+
* if JOB = 'E', LWORK >= max(1,M*(N-M));
|
|
541
|
+
* if JOB = 'V' or 'B', LWORK >= max(1,2*M*(N-M)).
|
|
542
|
+
*
|
|
543
|
+
* If LWORK = -1, then a workspace query is assumed; the routine
|
|
544
|
+
* only calculates the optimal size of the WORK array, returns
|
|
545
|
+
* this value as the first entry of the WORK array, and no error
|
|
546
|
+
* message related to LWORK is issued by XERBLA.
|
|
547
|
+
*
|
|
548
|
+
* IWORK (workspace) INTEGER array, dimension (MAX(1,LIWORK))
|
|
549
|
+
* On exit, if INFO = 0, IWORK(1) returns the optimal LIWORK.
|
|
550
|
+
*
|
|
551
|
+
* LIWORK (input) INTEGER
|
|
552
|
+
* The dimension of the array IWORK.
|
|
553
|
+
* If JOB = 'N' or 'E', LIWORK >= 1;
|
|
554
|
+
* if JOB = 'V' or 'B', LIWORK >= max(1,M*(N-M)).
|
|
555
|
+
*
|
|
556
|
+
* If LIWORK = -1, then a workspace query is assumed; the
|
|
557
|
+
* routine only calculates the optimal size of the IWORK array,
|
|
558
|
+
* returns this value as the first entry of the IWORK array, and
|
|
559
|
+
* no error message related to LIWORK is issued by XERBLA.
|
|
560
|
+
*
|
|
561
|
+
* INFO (output) INTEGER
|
|
562
|
+
* = 0: successful exit
|
|
563
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
564
|
+
* = 1: reordering of T failed because some eigenvalues are too
|
|
565
|
+
* close to separate (the problem is very ill-conditioned);
|
|
566
|
+
* T may have been partially reordered, and WR and WI
|
|
567
|
+
* contain the eigenvalues in the same order as in T; S and
|
|
568
|
+
* SEP (if requested) are set to zero.
|
|
569
|
+
*
|
|
570
|
+
|
|
571
|
+
* Further Details
|
|
572
|
+
* ===============
|
|
573
|
+
*
|
|
574
|
+
* STRSEN first collects the selected eigenvalues by computing an
|
|
575
|
+
* orthogonal transformation Z to move them to the top left corner of T.
|
|
576
|
+
* In other words, the selected eigenvalues are the eigenvalues of T11
|
|
577
|
+
* in:
|
|
578
|
+
*
|
|
579
|
+
* Z'*T*Z = ( T11 T12 ) n1
|
|
580
|
+
* ( 0 T22 ) n2
|
|
581
|
+
* n1 n2
|
|
582
|
+
*
|
|
583
|
+
* where N = n1+n2 and Z' means the transpose of Z. The first n1 columns
|
|
584
|
+
* of Z span the specified invariant subspace of T.
|
|
585
|
+
*
|
|
586
|
+
* If T has been obtained from the real Schur factorization of a matrix
|
|
587
|
+
* A = Q*T*Q', then the reordered real Schur factorization of A is given
|
|
588
|
+
* by A = (Q*Z)*(Z'*T*Z)*(Q*Z)', and the first n1 columns of Q*Z span
|
|
589
|
+
* the corresponding invariant subspace of A.
|
|
590
|
+
*
|
|
591
|
+
* The reciprocal condition number of the average of the eigenvalues of
|
|
592
|
+
* T11 may be returned in S. S lies between 0 (very badly conditioned)
|
|
593
|
+
* and 1 (very well conditioned). It is computed as follows. First we
|
|
594
|
+
* compute R so that
|
|
595
|
+
*
|
|
596
|
+
* P = ( I R ) n1
|
|
597
|
+
* ( 0 0 ) n2
|
|
598
|
+
* n1 n2
|
|
599
|
+
*
|
|
600
|
+
* is the projector on the invariant subspace associated with T11.
|
|
601
|
+
* R is the solution of the Sylvester equation:
|
|
602
|
+
*
|
|
603
|
+
* T11*R - R*T22 = T12.
|
|
604
|
+
*
|
|
605
|
+
* Let F-norm(M) denote the Frobenius-norm of M and 2-norm(M) denote
|
|
606
|
+
* the two-norm of M. Then S is computed as the lower bound
|
|
607
|
+
*
|
|
608
|
+
* (1 + F-norm(R)**2)**(-1/2)
|
|
609
|
+
*
|
|
610
|
+
* on the reciprocal of 2-norm(P), the true reciprocal condition number.
|
|
611
|
+
* S cannot underestimate 1 / 2-norm(P) by more than a factor of
|
|
612
|
+
* sqrt(N).
|
|
613
|
+
*
|
|
614
|
+
* An approximate error bound for the computed average of the
|
|
615
|
+
* eigenvalues of T11 is
|
|
616
|
+
*
|
|
617
|
+
* EPS * norm(T) / S
|
|
618
|
+
*
|
|
619
|
+
* where EPS is the machine precision.
|
|
620
|
+
*
|
|
621
|
+
* The reciprocal condition number of the right invariant subspace
|
|
622
|
+
* spanned by the first n1 columns of Z (or of Q*Z) is returned in SEP.
|
|
623
|
+
* SEP is defined as the separation of T11 and T22:
|
|
624
|
+
*
|
|
625
|
+
* sep( T11, T22 ) = sigma-min( C )
|
|
626
|
+
*
|
|
627
|
+
* where sigma-min(C) is the smallest singular value of the
|
|
628
|
+
* n1*n2-by-n1*n2 matrix
|
|
629
|
+
*
|
|
630
|
+
* C = kprod( I(n2), T11 ) - kprod( transpose(T22), I(n1) )
|
|
631
|
+
*
|
|
632
|
+
* I(m) is an m by m identity matrix, and kprod denotes the Kronecker
|
|
633
|
+
* product. We estimate sigma-min(C) by the reciprocal of an estimate of
|
|
634
|
+
* the 1-norm of inverse(C). The true reciprocal 1-norm of inverse(C)
|
|
635
|
+
* cannot differ from sigma-min(C) by more than a factor of sqrt(n1*n2).
|
|
636
|
+
*
|
|
637
|
+
* When SEP is small, small changes in T can cause large changes in
|
|
638
|
+
* the invariant subspace. An approximate bound on the maximum angular
|
|
639
|
+
* error in the computed right invariant subspace is
|
|
640
|
+
*
|
|
641
|
+
* EPS * norm(T) / SEP
|
|
642
|
+
*
|
|
643
|
+
* =====================================================================
|
|
644
|
+
*
|
|
645
|
+
|
|
646
|
+
|
|
647
|
+
</PRE>
|
|
648
|
+
<A HREF="#top">go to the page top</A>
|
|
649
|
+
|
|
650
|
+
<A NAME="strsna"></A>
|
|
651
|
+
<H2>strsna</H2>
|
|
652
|
+
<PRE>
|
|
653
|
+
USAGE:
|
|
654
|
+
s, sep, m, info = NumRu::Lapack.strsna( job, howmny, select, t, vl, vr, [:usage => usage, :help => help])
|
|
655
|
+
|
|
656
|
+
|
|
657
|
+
FORTRAN MANUAL
|
|
658
|
+
SUBROUTINE STRSNA( JOB, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR, LDVR, S, SEP, MM, M, WORK, LDWORK, IWORK, INFO )
|
|
659
|
+
|
|
660
|
+
* Purpose
|
|
661
|
+
* =======
|
|
662
|
+
*
|
|
663
|
+
* STRSNA estimates reciprocal condition numbers for specified
|
|
664
|
+
* eigenvalues and/or right eigenvectors of a real upper
|
|
665
|
+
* quasi-triangular matrix T (or of any matrix Q*T*Q**T with Q
|
|
666
|
+
* orthogonal).
|
|
667
|
+
*
|
|
668
|
+
* T must be in Schur canonical form (as returned by SHSEQR), that is,
|
|
669
|
+
* block upper triangular with 1-by-1 and 2-by-2 diagonal blocks; each
|
|
670
|
+
* 2-by-2 diagonal block has its diagonal elements equal and its
|
|
671
|
+
* off-diagonal elements of opposite sign.
|
|
672
|
+
*
|
|
673
|
+
|
|
674
|
+
* Arguments
|
|
675
|
+
* =========
|
|
676
|
+
*
|
|
677
|
+
* JOB (input) CHARACTER*1
|
|
678
|
+
* Specifies whether condition numbers are required for
|
|
679
|
+
* eigenvalues (S) or eigenvectors (SEP):
|
|
680
|
+
* = 'E': for eigenvalues only (S);
|
|
681
|
+
* = 'V': for eigenvectors only (SEP);
|
|
682
|
+
* = 'B': for both eigenvalues and eigenvectors (S and SEP).
|
|
683
|
+
*
|
|
684
|
+
* HOWMNY (input) CHARACTER*1
|
|
685
|
+
* = 'A': compute condition numbers for all eigenpairs;
|
|
686
|
+
* = 'S': compute condition numbers for selected eigenpairs
|
|
687
|
+
* specified by the array SELECT.
|
|
688
|
+
*
|
|
689
|
+
* SELECT (input) LOGICAL array, dimension (N)
|
|
690
|
+
* If HOWMNY = 'S', SELECT specifies the eigenpairs for which
|
|
691
|
+
* condition numbers are required. To select condition numbers
|
|
692
|
+
* for the eigenpair corresponding to a real eigenvalue w(j),
|
|
693
|
+
* SELECT(j) must be set to .TRUE.. To select condition numbers
|
|
694
|
+
* corresponding to a complex conjugate pair of eigenvalues w(j)
|
|
695
|
+
* and w(j+1), either SELECT(j) or SELECT(j+1) or both, must be
|
|
696
|
+
* set to .TRUE..
|
|
697
|
+
* If HOWMNY = 'A', SELECT is not referenced.
|
|
698
|
+
*
|
|
699
|
+
* N (input) INTEGER
|
|
700
|
+
* The order of the matrix T. N >= 0.
|
|
701
|
+
*
|
|
702
|
+
* T (input) REAL array, dimension (LDT,N)
|
|
703
|
+
* The upper quasi-triangular matrix T, in Schur canonical form.
|
|
704
|
+
*
|
|
705
|
+
* LDT (input) INTEGER
|
|
706
|
+
* The leading dimension of the array T. LDT >= max(1,N).
|
|
707
|
+
*
|
|
708
|
+
* VL (input) REAL array, dimension (LDVL,M)
|
|
709
|
+
* If JOB = 'E' or 'B', VL must contain left eigenvectors of T
|
|
710
|
+
* (or of any Q*T*Q**T with Q orthogonal), corresponding to the
|
|
711
|
+
* eigenpairs specified by HOWMNY and SELECT. The eigenvectors
|
|
712
|
+
* must be stored in consecutive columns of VL, as returned by
|
|
713
|
+
* SHSEIN or STREVC.
|
|
714
|
+
* If JOB = 'V', VL is not referenced.
|
|
715
|
+
*
|
|
716
|
+
* LDVL (input) INTEGER
|
|
717
|
+
* The leading dimension of the array VL.
|
|
718
|
+
* LDVL >= 1; and if JOB = 'E' or 'B', LDVL >= N.
|
|
719
|
+
*
|
|
720
|
+
* VR (input) REAL array, dimension (LDVR,M)
|
|
721
|
+
* If JOB = 'E' or 'B', VR must contain right eigenvectors of T
|
|
722
|
+
* (or of any Q*T*Q**T with Q orthogonal), corresponding to the
|
|
723
|
+
* eigenpairs specified by HOWMNY and SELECT. The eigenvectors
|
|
724
|
+
* must be stored in consecutive columns of VR, as returned by
|
|
725
|
+
* SHSEIN or STREVC.
|
|
726
|
+
* If JOB = 'V', VR is not referenced.
|
|
727
|
+
*
|
|
728
|
+
* LDVR (input) INTEGER
|
|
729
|
+
* The leading dimension of the array VR.
|
|
730
|
+
* LDVR >= 1; and if JOB = 'E' or 'B', LDVR >= N.
|
|
731
|
+
*
|
|
732
|
+
* S (output) REAL array, dimension (MM)
|
|
733
|
+
* If JOB = 'E' or 'B', the reciprocal condition numbers of the
|
|
734
|
+
* selected eigenvalues, stored in consecutive elements of the
|
|
735
|
+
* array. For a complex conjugate pair of eigenvalues two
|
|
736
|
+
* consecutive elements of S are set to the same value. Thus
|
|
737
|
+
* S(j), SEP(j), and the j-th columns of VL and VR all
|
|
738
|
+
* correspond to the same eigenpair (but not in general the
|
|
739
|
+
* j-th eigenpair, unless all eigenpairs are selected).
|
|
740
|
+
* If JOB = 'V', S is not referenced.
|
|
741
|
+
*
|
|
742
|
+
* SEP (output) REAL array, dimension (MM)
|
|
743
|
+
* If JOB = 'V' or 'B', the estimated reciprocal condition
|
|
744
|
+
* numbers of the selected eigenvectors, stored in consecutive
|
|
745
|
+
* elements of the array. For a complex eigenvector two
|
|
746
|
+
* consecutive elements of SEP are set to the same value. If
|
|
747
|
+
* the eigenvalues cannot be reordered to compute SEP(j), SEP(j)
|
|
748
|
+
* is set to 0; this can only occur when the true value would be
|
|
749
|
+
* very small anyway.
|
|
750
|
+
* If JOB = 'E', SEP is not referenced.
|
|
751
|
+
*
|
|
752
|
+
* MM (input) INTEGER
|
|
753
|
+
* The number of elements in the arrays S (if JOB = 'E' or 'B')
|
|
754
|
+
* and/or SEP (if JOB = 'V' or 'B'). MM >= M.
|
|
755
|
+
*
|
|
756
|
+
* M (output) INTEGER
|
|
757
|
+
* The number of elements of the arrays S and/or SEP actually
|
|
758
|
+
* used to store the estimated condition numbers.
|
|
759
|
+
* If HOWMNY = 'A', M is set to N.
|
|
760
|
+
*
|
|
761
|
+
* WORK (workspace) REAL array, dimension (LDWORK,N+6)
|
|
762
|
+
* If JOB = 'E', WORK is not referenced.
|
|
763
|
+
*
|
|
764
|
+
* LDWORK (input) INTEGER
|
|
765
|
+
* The leading dimension of the array WORK.
|
|
766
|
+
* LDWORK >= 1; and if JOB = 'V' or 'B', LDWORK >= N.
|
|
767
|
+
*
|
|
768
|
+
* IWORK (workspace) INTEGER array, dimension (2*(N-1))
|
|
769
|
+
* If JOB = 'E', IWORK is not referenced.
|
|
770
|
+
*
|
|
771
|
+
* INFO (output) INTEGER
|
|
772
|
+
* = 0: successful exit
|
|
773
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
774
|
+
*
|
|
775
|
+
|
|
776
|
+
* Further Details
|
|
777
|
+
* ===============
|
|
778
|
+
*
|
|
779
|
+
* The reciprocal of the condition number of an eigenvalue lambda is
|
|
780
|
+
* defined as
|
|
781
|
+
*
|
|
782
|
+
* S(lambda) = |v'*u| / (norm(u)*norm(v))
|
|
783
|
+
*
|
|
784
|
+
* where u and v are the right and left eigenvectors of T corresponding
|
|
785
|
+
* to lambda; v' denotes the conjugate-transpose of v, and norm(u)
|
|
786
|
+
* denotes the Euclidean norm. These reciprocal condition numbers always
|
|
787
|
+
* lie between zero (very badly conditioned) and one (very well
|
|
788
|
+
* conditioned). If n = 1, S(lambda) is defined to be 1.
|
|
789
|
+
*
|
|
790
|
+
* An approximate error bound for a computed eigenvalue W(i) is given by
|
|
791
|
+
*
|
|
792
|
+
* EPS * norm(T) / S(i)
|
|
793
|
+
*
|
|
794
|
+
* where EPS is the machine precision.
|
|
795
|
+
*
|
|
796
|
+
* The reciprocal of the condition number of the right eigenvector u
|
|
797
|
+
* corresponding to lambda is defined as follows. Suppose
|
|
798
|
+
*
|
|
799
|
+
* T = ( lambda c )
|
|
800
|
+
* ( 0 T22 )
|
|
801
|
+
*
|
|
802
|
+
* Then the reciprocal condition number is
|
|
803
|
+
*
|
|
804
|
+
* SEP( lambda, T22 ) = sigma-min( T22 - lambda*I )
|
|
805
|
+
*
|
|
806
|
+
* where sigma-min denotes the smallest singular value. We approximate
|
|
807
|
+
* the smallest singular value by the reciprocal of an estimate of the
|
|
808
|
+
* one-norm of the inverse of T22 - lambda*I. If n = 1, SEP(1) is
|
|
809
|
+
* defined to be abs(T(1,1)).
|
|
810
|
+
*
|
|
811
|
+
* An approximate error bound for a computed right eigenvector VR(i)
|
|
812
|
+
* is given by
|
|
813
|
+
*
|
|
814
|
+
* EPS * norm(T) / SEP(i)
|
|
815
|
+
*
|
|
816
|
+
* =====================================================================
|
|
817
|
+
*
|
|
818
|
+
|
|
819
|
+
|
|
820
|
+
</PRE>
|
|
821
|
+
<A HREF="#top">go to the page top</A>
|
|
822
|
+
|
|
823
|
+
<A NAME="strsyl"></A>
|
|
824
|
+
<H2>strsyl</H2>
|
|
825
|
+
<PRE>
|
|
826
|
+
USAGE:
|
|
827
|
+
scale, info, c = NumRu::Lapack.strsyl( trana, tranb, isgn, a, b, c, [:usage => usage, :help => help])
|
|
828
|
+
|
|
829
|
+
|
|
830
|
+
FORTRAN MANUAL
|
|
831
|
+
SUBROUTINE STRSYL( TRANA, TRANB, ISGN, M, N, A, LDA, B, LDB, C, LDC, SCALE, INFO )
|
|
832
|
+
|
|
833
|
+
* Purpose
|
|
834
|
+
* =======
|
|
835
|
+
*
|
|
836
|
+
* STRSYL solves the real Sylvester matrix equation:
|
|
837
|
+
*
|
|
838
|
+
* op(A)*X + X*op(B) = scale*C or
|
|
839
|
+
* op(A)*X - X*op(B) = scale*C,
|
|
840
|
+
*
|
|
841
|
+
* where op(A) = A or A**T, and A and B are both upper quasi-
|
|
842
|
+
* triangular. A is M-by-M and B is N-by-N; the right hand side C and
|
|
843
|
+
* the solution X are M-by-N; and scale is an output scale factor, set
|
|
844
|
+
* <= 1 to avoid overflow in X.
|
|
845
|
+
*
|
|
846
|
+
* A and B must be in Schur canonical form (as returned by SHSEQR), that
|
|
847
|
+
* is, block upper triangular with 1-by-1 and 2-by-2 diagonal blocks;
|
|
848
|
+
* each 2-by-2 diagonal block has its diagonal elements equal and its
|
|
849
|
+
* off-diagonal elements of opposite sign.
|
|
850
|
+
*
|
|
851
|
+
|
|
852
|
+
* Arguments
|
|
853
|
+
* =========
|
|
854
|
+
*
|
|
855
|
+
* TRANA (input) CHARACTER*1
|
|
856
|
+
* Specifies the option op(A):
|
|
857
|
+
* = 'N': op(A) = A (No transpose)
|
|
858
|
+
* = 'T': op(A) = A**T (Transpose)
|
|
859
|
+
* = 'C': op(A) = A**H (Conjugate transpose = Transpose)
|
|
860
|
+
*
|
|
861
|
+
* TRANB (input) CHARACTER*1
|
|
862
|
+
* Specifies the option op(B):
|
|
863
|
+
* = 'N': op(B) = B (No transpose)
|
|
864
|
+
* = 'T': op(B) = B**T (Transpose)
|
|
865
|
+
* = 'C': op(B) = B**H (Conjugate transpose = Transpose)
|
|
866
|
+
*
|
|
867
|
+
* ISGN (input) INTEGER
|
|
868
|
+
* Specifies the sign in the equation:
|
|
869
|
+
* = +1: solve op(A)*X + X*op(B) = scale*C
|
|
870
|
+
* = -1: solve op(A)*X - X*op(B) = scale*C
|
|
871
|
+
*
|
|
872
|
+
* M (input) INTEGER
|
|
873
|
+
* The order of the matrix A, and the number of rows in the
|
|
874
|
+
* matrices X and C. M >= 0.
|
|
875
|
+
*
|
|
876
|
+
* N (input) INTEGER
|
|
877
|
+
* The order of the matrix B, and the number of columns in the
|
|
878
|
+
* matrices X and C. N >= 0.
|
|
879
|
+
*
|
|
880
|
+
* A (input) REAL array, dimension (LDA,M)
|
|
881
|
+
* The upper quasi-triangular matrix A, in Schur canonical form.
|
|
882
|
+
*
|
|
883
|
+
* LDA (input) INTEGER
|
|
884
|
+
* The leading dimension of the array A. LDA >= max(1,M).
|
|
885
|
+
*
|
|
886
|
+
* B (input) REAL array, dimension (LDB,N)
|
|
887
|
+
* The upper quasi-triangular matrix B, in Schur canonical form.
|
|
888
|
+
*
|
|
889
|
+
* LDB (input) INTEGER
|
|
890
|
+
* The leading dimension of the array B. LDB >= max(1,N).
|
|
891
|
+
*
|
|
892
|
+
* C (input/output) REAL array, dimension (LDC,N)
|
|
893
|
+
* On entry, the M-by-N right hand side matrix C.
|
|
894
|
+
* On exit, C is overwritten by the solution matrix X.
|
|
895
|
+
*
|
|
896
|
+
* LDC (input) INTEGER
|
|
897
|
+
* The leading dimension of the array C. LDC >= max(1,M)
|
|
898
|
+
*
|
|
899
|
+
* SCALE (output) REAL
|
|
900
|
+
* The scale factor, scale, set <= 1 to avoid overflow in X.
|
|
901
|
+
*
|
|
902
|
+
* INFO (output) INTEGER
|
|
903
|
+
* = 0: successful exit
|
|
904
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
905
|
+
* = 1: A and B have common or very close eigenvalues; perturbed
|
|
906
|
+
* values were used to solve the equation (but the matrices
|
|
907
|
+
* A and B are unchanged).
|
|
908
|
+
*
|
|
909
|
+
|
|
910
|
+
* =====================================================================
|
|
911
|
+
*
|
|
912
|
+
|
|
913
|
+
|
|
914
|
+
</PRE>
|
|
915
|
+
<A HREF="#top">go to the page top</A>
|
|
916
|
+
|
|
917
|
+
<A NAME="strti2"></A>
|
|
918
|
+
<H2>strti2</H2>
|
|
919
|
+
<PRE>
|
|
920
|
+
USAGE:
|
|
921
|
+
info, a = NumRu::Lapack.strti2( uplo, diag, a, [:usage => usage, :help => help])
|
|
922
|
+
|
|
923
|
+
|
|
924
|
+
FORTRAN MANUAL
|
|
925
|
+
SUBROUTINE STRTI2( UPLO, DIAG, N, A, LDA, INFO )
|
|
926
|
+
|
|
927
|
+
* Purpose
|
|
928
|
+
* =======
|
|
929
|
+
*
|
|
930
|
+
* STRTI2 computes the inverse of a real upper or lower triangular
|
|
931
|
+
* matrix.
|
|
932
|
+
*
|
|
933
|
+
* This is the Level 2 BLAS version of the algorithm.
|
|
934
|
+
*
|
|
935
|
+
|
|
936
|
+
* Arguments
|
|
937
|
+
* =========
|
|
938
|
+
*
|
|
939
|
+
* UPLO (input) CHARACTER*1
|
|
940
|
+
* Specifies whether the matrix A is upper or lower triangular.
|
|
941
|
+
* = 'U': Upper triangular
|
|
942
|
+
* = 'L': Lower triangular
|
|
943
|
+
*
|
|
944
|
+
* DIAG (input) CHARACTER*1
|
|
945
|
+
* Specifies whether or not the matrix A is unit triangular.
|
|
946
|
+
* = 'N': Non-unit triangular
|
|
947
|
+
* = 'U': Unit triangular
|
|
948
|
+
*
|
|
949
|
+
* N (input) INTEGER
|
|
950
|
+
* The order of the matrix A. N >= 0.
|
|
951
|
+
*
|
|
952
|
+
* A (input/output) REAL array, dimension (LDA,N)
|
|
953
|
+
* On entry, the triangular matrix A. If UPLO = 'U', the
|
|
954
|
+
* leading n by n upper triangular part of the array A contains
|
|
955
|
+
* the upper triangular matrix, and the strictly lower
|
|
956
|
+
* triangular part of A is not referenced. If UPLO = 'L', the
|
|
957
|
+
* leading n by n lower triangular part of the array A contains
|
|
958
|
+
* the lower triangular matrix, and the strictly upper
|
|
959
|
+
* triangular part of A is not referenced. If DIAG = 'U', the
|
|
960
|
+
* diagonal elements of A are also not referenced and are
|
|
961
|
+
* assumed to be 1.
|
|
962
|
+
*
|
|
963
|
+
* On exit, the (triangular) inverse of the original matrix, in
|
|
964
|
+
* the same storage format.
|
|
965
|
+
*
|
|
966
|
+
* LDA (input) INTEGER
|
|
967
|
+
* The leading dimension of the array A. LDA >= max(1,N).
|
|
968
|
+
*
|
|
969
|
+
* INFO (output) INTEGER
|
|
970
|
+
* = 0: successful exit
|
|
971
|
+
* < 0: if INFO = -k, the k-th argument had an illegal value
|
|
972
|
+
*
|
|
973
|
+
|
|
974
|
+
* =====================================================================
|
|
975
|
+
*
|
|
976
|
+
|
|
977
|
+
|
|
978
|
+
</PRE>
|
|
979
|
+
<A HREF="#top">go to the page top</A>
|
|
980
|
+
|
|
981
|
+
<A NAME="strtri"></A>
|
|
982
|
+
<H2>strtri</H2>
|
|
983
|
+
<PRE>
|
|
984
|
+
USAGE:
|
|
985
|
+
info, a = NumRu::Lapack.strtri( uplo, diag, a, [:usage => usage, :help => help])
|
|
986
|
+
|
|
987
|
+
|
|
988
|
+
FORTRAN MANUAL
|
|
989
|
+
SUBROUTINE STRTRI( UPLO, DIAG, N, A, LDA, INFO )
|
|
990
|
+
|
|
991
|
+
* Purpose
|
|
992
|
+
* =======
|
|
993
|
+
*
|
|
994
|
+
* STRTRI computes the inverse of a real upper or lower triangular
|
|
995
|
+
* matrix A.
|
|
996
|
+
*
|
|
997
|
+
* This is the Level 3 BLAS version of the algorithm.
|
|
998
|
+
*
|
|
999
|
+
|
|
1000
|
+
* Arguments
|
|
1001
|
+
* =========
|
|
1002
|
+
*
|
|
1003
|
+
* UPLO (input) CHARACTER*1
|
|
1004
|
+
* = 'U': A is upper triangular;
|
|
1005
|
+
* = 'L': A is lower triangular.
|
|
1006
|
+
*
|
|
1007
|
+
* DIAG (input) CHARACTER*1
|
|
1008
|
+
* = 'N': A is non-unit triangular;
|
|
1009
|
+
* = 'U': A is unit triangular.
|
|
1010
|
+
*
|
|
1011
|
+
* N (input) INTEGER
|
|
1012
|
+
* The order of the matrix A. N >= 0.
|
|
1013
|
+
*
|
|
1014
|
+
* A (input/output) REAL array, dimension (LDA,N)
|
|
1015
|
+
* On entry, the triangular matrix A. If UPLO = 'U', the
|
|
1016
|
+
* leading N-by-N upper triangular part of the array A contains
|
|
1017
|
+
* the upper triangular matrix, and the strictly lower
|
|
1018
|
+
* triangular part of A is not referenced. If UPLO = 'L', the
|
|
1019
|
+
* leading N-by-N lower triangular part of the array A contains
|
|
1020
|
+
* the lower triangular matrix, and the strictly upper
|
|
1021
|
+
* triangular part of A is not referenced. If DIAG = 'U', the
|
|
1022
|
+
* diagonal elements of A are also not referenced and are
|
|
1023
|
+
* assumed to be 1.
|
|
1024
|
+
* On exit, the (triangular) inverse of the original matrix, in
|
|
1025
|
+
* the same storage format.
|
|
1026
|
+
*
|
|
1027
|
+
* LDA (input) INTEGER
|
|
1028
|
+
* The leading dimension of the array A. LDA >= max(1,N).
|
|
1029
|
+
*
|
|
1030
|
+
* INFO (output) INTEGER
|
|
1031
|
+
* = 0: successful exit
|
|
1032
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1033
|
+
* > 0: if INFO = i, A(i,i) is exactly zero. The triangular
|
|
1034
|
+
* matrix is singular and its inverse can not be computed.
|
|
1035
|
+
*
|
|
1036
|
+
|
|
1037
|
+
* =====================================================================
|
|
1038
|
+
*
|
|
1039
|
+
|
|
1040
|
+
|
|
1041
|
+
</PRE>
|
|
1042
|
+
<A HREF="#top">go to the page top</A>
|
|
1043
|
+
|
|
1044
|
+
<A NAME="strtrs"></A>
|
|
1045
|
+
<H2>strtrs</H2>
|
|
1046
|
+
<PRE>
|
|
1047
|
+
USAGE:
|
|
1048
|
+
info, b = NumRu::Lapack.strtrs( uplo, trans, diag, a, b, [:usage => usage, :help => help])
|
|
1049
|
+
|
|
1050
|
+
|
|
1051
|
+
FORTRAN MANUAL
|
|
1052
|
+
SUBROUTINE STRTRS( UPLO, TRANS, DIAG, N, NRHS, A, LDA, B, LDB, INFO )
|
|
1053
|
+
|
|
1054
|
+
* Purpose
|
|
1055
|
+
* =======
|
|
1056
|
+
*
|
|
1057
|
+
* STRTRS solves a triangular system of the form
|
|
1058
|
+
*
|
|
1059
|
+
* A * X = B or A**T * X = B,
|
|
1060
|
+
*
|
|
1061
|
+
* where A is a triangular matrix of order N, and B is an N-by-NRHS
|
|
1062
|
+
* matrix. A check is made to verify that A is nonsingular.
|
|
1063
|
+
*
|
|
1064
|
+
|
|
1065
|
+
* Arguments
|
|
1066
|
+
* =========
|
|
1067
|
+
*
|
|
1068
|
+
* UPLO (input) CHARACTER*1
|
|
1069
|
+
* = 'U': A is upper triangular;
|
|
1070
|
+
* = 'L': A is lower triangular.
|
|
1071
|
+
*
|
|
1072
|
+
* TRANS (input) CHARACTER*1
|
|
1073
|
+
* Specifies the form of the system of equations:
|
|
1074
|
+
* = 'N': A * X = B (No transpose)
|
|
1075
|
+
* = 'T': A**T * X = B (Transpose)
|
|
1076
|
+
* = 'C': A**H * X = B (Conjugate transpose = Transpose)
|
|
1077
|
+
*
|
|
1078
|
+
* DIAG (input) CHARACTER*1
|
|
1079
|
+
* = 'N': A is non-unit triangular;
|
|
1080
|
+
* = 'U': A is unit triangular.
|
|
1081
|
+
*
|
|
1082
|
+
* N (input) INTEGER
|
|
1083
|
+
* The order of the matrix A. N >= 0.
|
|
1084
|
+
*
|
|
1085
|
+
* NRHS (input) INTEGER
|
|
1086
|
+
* The number of right hand sides, i.e., the number of columns
|
|
1087
|
+
* of the matrix B. NRHS >= 0.
|
|
1088
|
+
*
|
|
1089
|
+
* A (input) REAL array, dimension (LDA,N)
|
|
1090
|
+
* The triangular matrix A. If UPLO = 'U', the leading N-by-N
|
|
1091
|
+
* upper triangular part of the array A contains the upper
|
|
1092
|
+
* triangular matrix, and the strictly lower triangular part of
|
|
1093
|
+
* A is not referenced. If UPLO = 'L', the leading N-by-N lower
|
|
1094
|
+
* triangular part of the array A contains the lower triangular
|
|
1095
|
+
* matrix, and the strictly upper triangular part of A is not
|
|
1096
|
+
* referenced. If DIAG = 'U', the diagonal elements of A are
|
|
1097
|
+
* also not referenced and are assumed to be 1.
|
|
1098
|
+
*
|
|
1099
|
+
* LDA (input) INTEGER
|
|
1100
|
+
* The leading dimension of the array A. LDA >= max(1,N).
|
|
1101
|
+
*
|
|
1102
|
+
* B (input/output) REAL array, dimension (LDB,NRHS)
|
|
1103
|
+
* On entry, the right hand side matrix B.
|
|
1104
|
+
* On exit, if INFO = 0, the solution matrix X.
|
|
1105
|
+
*
|
|
1106
|
+
* LDB (input) INTEGER
|
|
1107
|
+
* The leading dimension of the array B. LDB >= max(1,N).
|
|
1108
|
+
*
|
|
1109
|
+
* INFO (output) INTEGER
|
|
1110
|
+
* = 0: successful exit
|
|
1111
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1112
|
+
* > 0: if INFO = i, the i-th diagonal element of A is zero,
|
|
1113
|
+
* indicating that the matrix is singular and the solutions
|
|
1114
|
+
* X have not been computed.
|
|
1115
|
+
*
|
|
1116
|
+
|
|
1117
|
+
* =====================================================================
|
|
1118
|
+
*
|
|
1119
|
+
|
|
1120
|
+
|
|
1121
|
+
</PRE>
|
|
1122
|
+
<A HREF="#top">go to the page top</A>
|
|
1123
|
+
|
|
1124
|
+
<A NAME="strttf"></A>
|
|
1125
|
+
<H2>strttf</H2>
|
|
1126
|
+
<PRE>
|
|
1127
|
+
USAGE:
|
|
1128
|
+
arf, info = NumRu::Lapack.strttf( transr, uplo, a, [:usage => usage, :help => help])
|
|
1129
|
+
|
|
1130
|
+
|
|
1131
|
+
FORTRAN MANUAL
|
|
1132
|
+
SUBROUTINE STRTTF( TRANSR, UPLO, N, A, LDA, ARF, INFO )
|
|
1133
|
+
|
|
1134
|
+
* Purpose
|
|
1135
|
+
* =======
|
|
1136
|
+
*
|
|
1137
|
+
* STRTTF copies a triangular matrix A from standard full format (TR)
|
|
1138
|
+
* to rectangular full packed format (TF) .
|
|
1139
|
+
*
|
|
1140
|
+
|
|
1141
|
+
* Arguments
|
|
1142
|
+
* =========
|
|
1143
|
+
*
|
|
1144
|
+
* TRANSR (input) CHARACTER*1
|
|
1145
|
+
* = 'N': ARF in Normal form is wanted;
|
|
1146
|
+
* = 'T': ARF in Transpose form is wanted.
|
|
1147
|
+
*
|
|
1148
|
+
* UPLO (input) CHARACTER*1
|
|
1149
|
+
* = 'U': Upper triangle of A is stored;
|
|
1150
|
+
* = 'L': Lower triangle of A is stored.
|
|
1151
|
+
*
|
|
1152
|
+
* N (input) INTEGER
|
|
1153
|
+
* The order of the matrix A. N >= 0.
|
|
1154
|
+
*
|
|
1155
|
+
* A (input) REAL array, dimension (LDA,N).
|
|
1156
|
+
* On entry, the triangular matrix A. If UPLO = 'U', the
|
|
1157
|
+
* leading N-by-N upper triangular part of the array A contains
|
|
1158
|
+
* the upper triangular matrix, and the strictly lower
|
|
1159
|
+
* triangular part of A is not referenced. If UPLO = 'L', the
|
|
1160
|
+
* leading N-by-N lower triangular part of the array A contains
|
|
1161
|
+
* the lower triangular matrix, and the strictly upper
|
|
1162
|
+
* triangular part of A is not referenced.
|
|
1163
|
+
*
|
|
1164
|
+
* LDA (input) INTEGER
|
|
1165
|
+
* The leading dimension of the matrix A. LDA >= max(1,N).
|
|
1166
|
+
*
|
|
1167
|
+
* ARF (output) REAL array, dimension (NT).
|
|
1168
|
+
* NT=N*(N+1)/2. On exit, the triangular matrix A in RFP format.
|
|
1169
|
+
*
|
|
1170
|
+
* INFO (output) INTEGER
|
|
1171
|
+
* = 0: successful exit
|
|
1172
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1173
|
+
*
|
|
1174
|
+
|
|
1175
|
+
* Further Details
|
|
1176
|
+
* ===============
|
|
1177
|
+
*
|
|
1178
|
+
* We first consider Rectangular Full Packed (RFP) Format when N is
|
|
1179
|
+
* even. We give an example where N = 6.
|
|
1180
|
+
*
|
|
1181
|
+
* AP is Upper AP is Lower
|
|
1182
|
+
*
|
|
1183
|
+
* 00 01 02 03 04 05 00
|
|
1184
|
+
* 11 12 13 14 15 10 11
|
|
1185
|
+
* 22 23 24 25 20 21 22
|
|
1186
|
+
* 33 34 35 30 31 32 33
|
|
1187
|
+
* 44 45 40 41 42 43 44
|
|
1188
|
+
* 55 50 51 52 53 54 55
|
|
1189
|
+
*
|
|
1190
|
+
*
|
|
1191
|
+
* Let TRANSR = 'N'. RFP holds AP as follows:
|
|
1192
|
+
* For UPLO = 'U' the upper trapezoid A(0:5,0:2) consists of the last
|
|
1193
|
+
* three columns of AP upper. The lower triangle A(4:6,0:2) consists of
|
|
1194
|
+
* the transpose of the first three columns of AP upper.
|
|
1195
|
+
* For UPLO = 'L' the lower trapezoid A(1:6,0:2) consists of the first
|
|
1196
|
+
* three columns of AP lower. The upper triangle A(0:2,0:2) consists of
|
|
1197
|
+
* the transpose of the last three columns of AP lower.
|
|
1198
|
+
* This covers the case N even and TRANSR = 'N'.
|
|
1199
|
+
*
|
|
1200
|
+
* RFP A RFP A
|
|
1201
|
+
*
|
|
1202
|
+
* 03 04 05 33 43 53
|
|
1203
|
+
* 13 14 15 00 44 54
|
|
1204
|
+
* 23 24 25 10 11 55
|
|
1205
|
+
* 33 34 35 20 21 22
|
|
1206
|
+
* 00 44 45 30 31 32
|
|
1207
|
+
* 01 11 55 40 41 42
|
|
1208
|
+
* 02 12 22 50 51 52
|
|
1209
|
+
*
|
|
1210
|
+
* Now let TRANSR = 'T'. RFP A in both UPLO cases is just the
|
|
1211
|
+
* transpose of RFP A above. One therefore gets:
|
|
1212
|
+
*
|
|
1213
|
+
*
|
|
1214
|
+
* RFP A RFP A
|
|
1215
|
+
*
|
|
1216
|
+
* 03 13 23 33 00 01 02 33 00 10 20 30 40 50
|
|
1217
|
+
* 04 14 24 34 44 11 12 43 44 11 21 31 41 51
|
|
1218
|
+
* 05 15 25 35 45 55 22 53 54 55 22 32 42 52
|
|
1219
|
+
*
|
|
1220
|
+
*
|
|
1221
|
+
* We then consider Rectangular Full Packed (RFP) Format when N is
|
|
1222
|
+
* odd. We give an example where N = 5.
|
|
1223
|
+
*
|
|
1224
|
+
* AP is Upper AP is Lower
|
|
1225
|
+
*
|
|
1226
|
+
* 00 01 02 03 04 00
|
|
1227
|
+
* 11 12 13 14 10 11
|
|
1228
|
+
* 22 23 24 20 21 22
|
|
1229
|
+
* 33 34 30 31 32 33
|
|
1230
|
+
* 44 40 41 42 43 44
|
|
1231
|
+
*
|
|
1232
|
+
*
|
|
1233
|
+
* Let TRANSR = 'N'. RFP holds AP as follows:
|
|
1234
|
+
* For UPLO = 'U' the upper trapezoid A(0:4,0:2) consists of the last
|
|
1235
|
+
* three columns of AP upper. The lower triangle A(3:4,0:1) consists of
|
|
1236
|
+
* the transpose of the first two columns of AP upper.
|
|
1237
|
+
* For UPLO = 'L' the lower trapezoid A(0:4,0:2) consists of the first
|
|
1238
|
+
* three columns of AP lower. The upper triangle A(0:1,1:2) consists of
|
|
1239
|
+
* the transpose of the last two columns of AP lower.
|
|
1240
|
+
* This covers the case N odd and TRANSR = 'N'.
|
|
1241
|
+
*
|
|
1242
|
+
* RFP A RFP A
|
|
1243
|
+
*
|
|
1244
|
+
* 02 03 04 00 33 43
|
|
1245
|
+
* 12 13 14 10 11 44
|
|
1246
|
+
* 22 23 24 20 21 22
|
|
1247
|
+
* 00 33 34 30 31 32
|
|
1248
|
+
* 01 11 44 40 41 42
|
|
1249
|
+
*
|
|
1250
|
+
* Now let TRANSR = 'T'. RFP A in both UPLO cases is just the
|
|
1251
|
+
* transpose of RFP A above. One therefore gets:
|
|
1252
|
+
*
|
|
1253
|
+
* RFP A RFP A
|
|
1254
|
+
*
|
|
1255
|
+
* 02 12 22 00 01 00 10 20 30 40 50
|
|
1256
|
+
* 03 13 23 33 11 33 11 21 31 41 51
|
|
1257
|
+
* 04 14 24 34 44 43 44 22 32 42 52
|
|
1258
|
+
*
|
|
1259
|
+
* Reference
|
|
1260
|
+
* =========
|
|
1261
|
+
*
|
|
1262
|
+
* =====================================================================
|
|
1263
|
+
*
|
|
1264
|
+
* ..
|
|
1265
|
+
* .. Local Scalars ..
|
|
1266
|
+
LOGICAL LOWER, NISODD, NORMALTRANSR
|
|
1267
|
+
INTEGER I, IJ, J, K, L, N1, N2, NT, NX2, NP1X2
|
|
1268
|
+
* ..
|
|
1269
|
+
* .. External Functions ..
|
|
1270
|
+
LOGICAL LSAME
|
|
1271
|
+
EXTERNAL LSAME
|
|
1272
|
+
* ..
|
|
1273
|
+
* .. External Subroutines ..
|
|
1274
|
+
EXTERNAL XERBLA
|
|
1275
|
+
* ..
|
|
1276
|
+
* .. Intrinsic Functions ..
|
|
1277
|
+
INTRINSIC MAX, MOD
|
|
1278
|
+
* ..
|
|
1279
|
+
|
|
1280
|
+
|
|
1281
|
+
</PRE>
|
|
1282
|
+
<A HREF="#top">go to the page top</A>
|
|
1283
|
+
|
|
1284
|
+
<A NAME="strttp"></A>
|
|
1285
|
+
<H2>strttp</H2>
|
|
1286
|
+
<PRE>
|
|
1287
|
+
USAGE:
|
|
1288
|
+
ap, info = NumRu::Lapack.strttp( uplo, a, [:usage => usage, :help => help])
|
|
1289
|
+
|
|
1290
|
+
|
|
1291
|
+
FORTRAN MANUAL
|
|
1292
|
+
SUBROUTINE STRTTP( UPLO, N, A, LDA, AP, INFO )
|
|
1293
|
+
|
|
1294
|
+
* Purpose
|
|
1295
|
+
* =======
|
|
1296
|
+
*
|
|
1297
|
+
* STRTTP copies a triangular matrix A from full format (TR) to standard
|
|
1298
|
+
* packed format (TP).
|
|
1299
|
+
*
|
|
1300
|
+
|
|
1301
|
+
* Arguments
|
|
1302
|
+
* =========
|
|
1303
|
+
*
|
|
1304
|
+
* UPLO (input) CHARACTER*1
|
|
1305
|
+
* = 'U': A is upper triangular.
|
|
1306
|
+
* = 'L': A is lower triangular.
|
|
1307
|
+
*
|
|
1308
|
+
* N (input) INTEGER
|
|
1309
|
+
* The order of the matrices AP and A. N >= 0.
|
|
1310
|
+
*
|
|
1311
|
+
* A (input) REAL array, dimension (LDA,N)
|
|
1312
|
+
* On exit, the triangular matrix A. If UPLO = 'U', the leading
|
|
1313
|
+
* N-by-N upper triangular part of A contains the upper
|
|
1314
|
+
* triangular part of the matrix A, and the strictly lower
|
|
1315
|
+
* triangular part of A is not referenced. If UPLO = 'L', the
|
|
1316
|
+
* leading N-by-N lower triangular part of A contains the lower
|
|
1317
|
+
* triangular part of the matrix A, and the strictly upper
|
|
1318
|
+
* triangular part of A is not referenced.
|
|
1319
|
+
*
|
|
1320
|
+
* LDA (input) INTEGER
|
|
1321
|
+
* The leading dimension of the array A. LDA >= max(1,N).
|
|
1322
|
+
*
|
|
1323
|
+
* AP (output) REAL array, dimension (N*(N+1)/2
|
|
1324
|
+
* On exit, the upper or lower triangular matrix A, packed
|
|
1325
|
+
* columnwise in a linear array. The j-th column of A is stored
|
|
1326
|
+
* in the array AP as follows:
|
|
1327
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
1328
|
+
* if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.
|
|
1329
|
+
*
|
|
1330
|
+
* INFO (output) INTEGER
|
|
1331
|
+
* = 0: successful exit
|
|
1332
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1333
|
+
*
|
|
1334
|
+
|
|
1335
|
+
* =====================================================================
|
|
1336
|
+
*
|
|
1337
|
+
|
|
1338
|
+
|
|
1339
|
+
</PRE>
|
|
1340
|
+
<A HREF="#top">go to the page top</A>
|
|
1341
|
+
|
|
1342
|
+
<HR />
|
|
1343
|
+
<A HREF="s.html">back to matrix types</A><BR>
|
|
1344
|
+
<A HREF="s.html">back to data types</A>
|
|
1345
|
+
</BODY>
|
|
1346
|
+
</HTML>
|