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/ssp.html
ADDED
|
@@ -0,0 +1,1647 @@
|
|
|
1
|
+
<HTML>
|
|
2
|
+
<HEAD>
|
|
3
|
+
<TITLE>REAL routines for symmetric, packed storage matrix</TITLE>
|
|
4
|
+
</HEAD>
|
|
5
|
+
<BODY>
|
|
6
|
+
<A NAME="top"></A>
|
|
7
|
+
<H1>REAL routines for symmetric, packed storage matrix</H1>
|
|
8
|
+
<UL>
|
|
9
|
+
<LI><A HREF="#sspcon">sspcon</A></LI>
|
|
10
|
+
<LI><A HREF="#sspev">sspev</A></LI>
|
|
11
|
+
<LI><A HREF="#sspevd">sspevd</A></LI>
|
|
12
|
+
<LI><A HREF="#sspevx">sspevx</A></LI>
|
|
13
|
+
<LI><A HREF="#sspgst">sspgst</A></LI>
|
|
14
|
+
<LI><A HREF="#sspgv">sspgv</A></LI>
|
|
15
|
+
<LI><A HREF="#sspgvd">sspgvd</A></LI>
|
|
16
|
+
<LI><A HREF="#sspgvx">sspgvx</A></LI>
|
|
17
|
+
<LI><A HREF="#ssprfs">ssprfs</A></LI>
|
|
18
|
+
<LI><A HREF="#sspsv">sspsv</A></LI>
|
|
19
|
+
<LI><A HREF="#sspsvx">sspsvx</A></LI>
|
|
20
|
+
<LI><A HREF="#ssptrd">ssptrd</A></LI>
|
|
21
|
+
<LI><A HREF="#ssptrf">ssptrf</A></LI>
|
|
22
|
+
<LI><A HREF="#ssptri">ssptri</A></LI>
|
|
23
|
+
<LI><A HREF="#ssptrs">ssptrs</A></LI>
|
|
24
|
+
</UL>
|
|
25
|
+
|
|
26
|
+
<A NAME="sspcon"></A>
|
|
27
|
+
<H2>sspcon</H2>
|
|
28
|
+
<PRE>
|
|
29
|
+
USAGE:
|
|
30
|
+
rcond, info = NumRu::Lapack.sspcon( uplo, ap, ipiv, anorm, [:usage => usage, :help => help])
|
|
31
|
+
|
|
32
|
+
|
|
33
|
+
FORTRAN MANUAL
|
|
34
|
+
SUBROUTINE SSPCON( UPLO, N, AP, IPIV, ANORM, RCOND, WORK, IWORK, INFO )
|
|
35
|
+
|
|
36
|
+
* Purpose
|
|
37
|
+
* =======
|
|
38
|
+
*
|
|
39
|
+
* SSPCON estimates the reciprocal of the condition number (in the
|
|
40
|
+
* 1-norm) of a real symmetric packed matrix A using the factorization
|
|
41
|
+
* A = U*D*U**T or A = L*D*L**T computed by SSPTRF.
|
|
42
|
+
*
|
|
43
|
+
* An estimate is obtained for norm(inv(A)), and the reciprocal of the
|
|
44
|
+
* condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
|
|
45
|
+
*
|
|
46
|
+
|
|
47
|
+
* Arguments
|
|
48
|
+
* =========
|
|
49
|
+
*
|
|
50
|
+
* UPLO (input) CHARACTER*1
|
|
51
|
+
* Specifies whether the details of the factorization are stored
|
|
52
|
+
* as an upper or lower triangular matrix.
|
|
53
|
+
* = 'U': Upper triangular, form is A = U*D*U**T;
|
|
54
|
+
* = 'L': Lower triangular, form is A = L*D*L**T.
|
|
55
|
+
*
|
|
56
|
+
* N (input) INTEGER
|
|
57
|
+
* The order of the matrix A. N >= 0.
|
|
58
|
+
*
|
|
59
|
+
* AP (input) REAL array, dimension (N*(N+1)/2)
|
|
60
|
+
* The block diagonal matrix D and the multipliers used to
|
|
61
|
+
* obtain the factor U or L as computed by SSPTRF, stored as a
|
|
62
|
+
* packed triangular matrix.
|
|
63
|
+
*
|
|
64
|
+
* IPIV (input) INTEGER array, dimension (N)
|
|
65
|
+
* Details of the interchanges and the block structure of D
|
|
66
|
+
* as determined by SSPTRF.
|
|
67
|
+
*
|
|
68
|
+
* ANORM (input) REAL
|
|
69
|
+
* The 1-norm of the original matrix A.
|
|
70
|
+
*
|
|
71
|
+
* RCOND (output) REAL
|
|
72
|
+
* The reciprocal of the condition number of the matrix A,
|
|
73
|
+
* computed as RCOND = 1/(ANORM * AINVNM), where AINVNM is an
|
|
74
|
+
* estimate of the 1-norm of inv(A) computed in this routine.
|
|
75
|
+
*
|
|
76
|
+
* WORK (workspace) REAL array, dimension (2*N)
|
|
77
|
+
*
|
|
78
|
+
* IWORK (workspace) INTEGER array, dimension (N)
|
|
79
|
+
*
|
|
80
|
+
* INFO (output) INTEGER
|
|
81
|
+
* = 0: successful exit
|
|
82
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
83
|
+
*
|
|
84
|
+
|
|
85
|
+
* =====================================================================
|
|
86
|
+
*
|
|
87
|
+
|
|
88
|
+
|
|
89
|
+
</PRE>
|
|
90
|
+
<A HREF="#top">go to the page top</A>
|
|
91
|
+
|
|
92
|
+
<A NAME="sspev"></A>
|
|
93
|
+
<H2>sspev</H2>
|
|
94
|
+
<PRE>
|
|
95
|
+
USAGE:
|
|
96
|
+
w, z, info, ap = NumRu::Lapack.sspev( jobz, uplo, ap, [:usage => usage, :help => help])
|
|
97
|
+
|
|
98
|
+
|
|
99
|
+
FORTRAN MANUAL
|
|
100
|
+
SUBROUTINE SSPEV( JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, INFO )
|
|
101
|
+
|
|
102
|
+
* Purpose
|
|
103
|
+
* =======
|
|
104
|
+
*
|
|
105
|
+
* SSPEV computes all the eigenvalues and, optionally, eigenvectors of a
|
|
106
|
+
* real symmetric matrix A in packed storage.
|
|
107
|
+
*
|
|
108
|
+
|
|
109
|
+
* Arguments
|
|
110
|
+
* =========
|
|
111
|
+
*
|
|
112
|
+
* JOBZ (input) CHARACTER*1
|
|
113
|
+
* = 'N': Compute eigenvalues only;
|
|
114
|
+
* = 'V': Compute eigenvalues and eigenvectors.
|
|
115
|
+
*
|
|
116
|
+
* UPLO (input) CHARACTER*1
|
|
117
|
+
* = 'U': Upper triangle of A is stored;
|
|
118
|
+
* = 'L': Lower triangle of A is stored.
|
|
119
|
+
*
|
|
120
|
+
* N (input) INTEGER
|
|
121
|
+
* The order of the matrix A. N >= 0.
|
|
122
|
+
*
|
|
123
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
124
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
125
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
126
|
+
* is stored in the array AP as follows:
|
|
127
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
128
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
129
|
+
*
|
|
130
|
+
* On exit, AP is overwritten by values generated during the
|
|
131
|
+
* reduction to tridiagonal form. If UPLO = 'U', the diagonal
|
|
132
|
+
* and first superdiagonal of the tridiagonal matrix T overwrite
|
|
133
|
+
* the corresponding elements of A, and if UPLO = 'L', the
|
|
134
|
+
* diagonal and first subdiagonal of T overwrite the
|
|
135
|
+
* corresponding elements of A.
|
|
136
|
+
*
|
|
137
|
+
* W (output) REAL array, dimension (N)
|
|
138
|
+
* If INFO = 0, the eigenvalues in ascending order.
|
|
139
|
+
*
|
|
140
|
+
* Z (output) REAL array, dimension (LDZ, N)
|
|
141
|
+
* If JOBZ = 'V', then if INFO = 0, Z contains the orthonormal
|
|
142
|
+
* eigenvectors of the matrix A, with the i-th column of Z
|
|
143
|
+
* holding the eigenvector associated with W(i).
|
|
144
|
+
* If JOBZ = 'N', then Z is not referenced.
|
|
145
|
+
*
|
|
146
|
+
* LDZ (input) INTEGER
|
|
147
|
+
* The leading dimension of the array Z. LDZ >= 1, and if
|
|
148
|
+
* JOBZ = 'V', LDZ >= max(1,N).
|
|
149
|
+
*
|
|
150
|
+
* WORK (workspace) REAL array, dimension (3*N)
|
|
151
|
+
*
|
|
152
|
+
* INFO (output) INTEGER
|
|
153
|
+
* = 0: successful exit.
|
|
154
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value.
|
|
155
|
+
* > 0: if INFO = i, the algorithm failed to converge; i
|
|
156
|
+
* off-diagonal elements of an intermediate tridiagonal
|
|
157
|
+
* form did not converge to zero.
|
|
158
|
+
*
|
|
159
|
+
|
|
160
|
+
* =====================================================================
|
|
161
|
+
*
|
|
162
|
+
|
|
163
|
+
|
|
164
|
+
</PRE>
|
|
165
|
+
<A HREF="#top">go to the page top</A>
|
|
166
|
+
|
|
167
|
+
<A NAME="sspevd"></A>
|
|
168
|
+
<H2>sspevd</H2>
|
|
169
|
+
<PRE>
|
|
170
|
+
USAGE:
|
|
171
|
+
w, z, work, iwork, info, ap = NumRu::Lapack.sspevd( jobz, uplo, ap, [:lwork => lwork, :liwork => liwork, :usage => usage, :help => help])
|
|
172
|
+
|
|
173
|
+
|
|
174
|
+
FORTRAN MANUAL
|
|
175
|
+
SUBROUTINE SSPEVD( JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, LWORK, IWORK, LIWORK, INFO )
|
|
176
|
+
|
|
177
|
+
* Purpose
|
|
178
|
+
* =======
|
|
179
|
+
*
|
|
180
|
+
* SSPEVD computes all the eigenvalues and, optionally, eigenvectors
|
|
181
|
+
* of a real symmetric matrix A in packed storage. If eigenvectors are
|
|
182
|
+
* desired, it uses a divide and conquer algorithm.
|
|
183
|
+
*
|
|
184
|
+
* The divide and conquer algorithm makes very mild assumptions about
|
|
185
|
+
* floating point arithmetic. It will work on machines with a guard
|
|
186
|
+
* digit in add/subtract, or on those binary machines without guard
|
|
187
|
+
* digits which subtract like the Cray X-MP, Cray Y-MP, Cray C-90, or
|
|
188
|
+
* Cray-2. It could conceivably fail on hexadecimal or decimal machines
|
|
189
|
+
* without guard digits, but we know of none.
|
|
190
|
+
*
|
|
191
|
+
|
|
192
|
+
* Arguments
|
|
193
|
+
* =========
|
|
194
|
+
*
|
|
195
|
+
* JOBZ (input) CHARACTER*1
|
|
196
|
+
* = 'N': Compute eigenvalues only;
|
|
197
|
+
* = 'V': Compute eigenvalues and eigenvectors.
|
|
198
|
+
*
|
|
199
|
+
* UPLO (input) CHARACTER*1
|
|
200
|
+
* = 'U': Upper triangle of A is stored;
|
|
201
|
+
* = 'L': Lower triangle of A is stored.
|
|
202
|
+
*
|
|
203
|
+
* N (input) INTEGER
|
|
204
|
+
* The order of the matrix A. N >= 0.
|
|
205
|
+
*
|
|
206
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
207
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
208
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
209
|
+
* is stored in the array AP as follows:
|
|
210
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
211
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
212
|
+
*
|
|
213
|
+
* On exit, AP is overwritten by values generated during the
|
|
214
|
+
* reduction to tridiagonal form. If UPLO = 'U', the diagonal
|
|
215
|
+
* and first superdiagonal of the tridiagonal matrix T overwrite
|
|
216
|
+
* the corresponding elements of A, and if UPLO = 'L', the
|
|
217
|
+
* diagonal and first subdiagonal of T overwrite the
|
|
218
|
+
* corresponding elements of A.
|
|
219
|
+
*
|
|
220
|
+
* W (output) REAL array, dimension (N)
|
|
221
|
+
* If INFO = 0, the eigenvalues in ascending order.
|
|
222
|
+
*
|
|
223
|
+
* Z (output) REAL array, dimension (LDZ, N)
|
|
224
|
+
* If JOBZ = 'V', then if INFO = 0, Z contains the orthonormal
|
|
225
|
+
* eigenvectors of the matrix A, with the i-th column of Z
|
|
226
|
+
* holding the eigenvector associated with W(i).
|
|
227
|
+
* If JOBZ = 'N', then Z is not referenced.
|
|
228
|
+
*
|
|
229
|
+
* LDZ (input) INTEGER
|
|
230
|
+
* The leading dimension of the array Z. LDZ >= 1, and if
|
|
231
|
+
* JOBZ = 'V', LDZ >= max(1,N).
|
|
232
|
+
*
|
|
233
|
+
* WORK (workspace/output) REAL array, dimension (MAX(1,LWORK))
|
|
234
|
+
* On exit, if INFO = 0, WORK(1) returns the required LWORK.
|
|
235
|
+
*
|
|
236
|
+
* LWORK (input) INTEGER
|
|
237
|
+
* The dimension of the array WORK.
|
|
238
|
+
* If N <= 1, LWORK must be at least 1.
|
|
239
|
+
* If JOBZ = 'N' and N > 1, LWORK must be at least 2*N.
|
|
240
|
+
* If JOBZ = 'V' and N > 1, LWORK must be at least
|
|
241
|
+
* 1 + 6*N + N**2.
|
|
242
|
+
*
|
|
243
|
+
* If LWORK = -1, then a workspace query is assumed; the routine
|
|
244
|
+
* only calculates the required sizes of the WORK and IWORK
|
|
245
|
+
* arrays, returns these values as the first entries of the WORK
|
|
246
|
+
* and IWORK arrays, and no error message related to LWORK or
|
|
247
|
+
* LIWORK is issued by XERBLA.
|
|
248
|
+
*
|
|
249
|
+
* IWORK (workspace/output) INTEGER array, dimension (MAX(1,LIWORK))
|
|
250
|
+
* On exit, if INFO = 0, IWORK(1) returns the required LIWORK.
|
|
251
|
+
*
|
|
252
|
+
* LIWORK (input) INTEGER
|
|
253
|
+
* The dimension of the array IWORK.
|
|
254
|
+
* If JOBZ = 'N' or N <= 1, LIWORK must be at least 1.
|
|
255
|
+
* If JOBZ = 'V' and N > 1, LIWORK must be at least 3 + 5*N.
|
|
256
|
+
*
|
|
257
|
+
* If LIWORK = -1, then a workspace query is assumed; the
|
|
258
|
+
* routine only calculates the required sizes of the WORK and
|
|
259
|
+
* IWORK arrays, returns these values as the first entries of
|
|
260
|
+
* the WORK and IWORK arrays, and no error message related to
|
|
261
|
+
* LWORK or LIWORK is issued by XERBLA.
|
|
262
|
+
*
|
|
263
|
+
* INFO (output) INTEGER
|
|
264
|
+
* = 0: successful exit
|
|
265
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value.
|
|
266
|
+
* > 0: if INFO = i, the algorithm failed to converge; i
|
|
267
|
+
* off-diagonal elements of an intermediate tridiagonal
|
|
268
|
+
* form did not converge to zero.
|
|
269
|
+
*
|
|
270
|
+
|
|
271
|
+
* =====================================================================
|
|
272
|
+
*
|
|
273
|
+
|
|
274
|
+
|
|
275
|
+
</PRE>
|
|
276
|
+
<A HREF="#top">go to the page top</A>
|
|
277
|
+
|
|
278
|
+
<A NAME="sspevx"></A>
|
|
279
|
+
<H2>sspevx</H2>
|
|
280
|
+
<PRE>
|
|
281
|
+
USAGE:
|
|
282
|
+
m, w, z, ifail, info, ap = NumRu::Lapack.sspevx( jobz, range, uplo, ap, vl, vu, il, iu, abstol, [:usage => usage, :help => help])
|
|
283
|
+
|
|
284
|
+
|
|
285
|
+
FORTRAN MANUAL
|
|
286
|
+
SUBROUTINE SSPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, IWORK, IFAIL, INFO )
|
|
287
|
+
|
|
288
|
+
* Purpose
|
|
289
|
+
* =======
|
|
290
|
+
*
|
|
291
|
+
* SSPEVX computes selected eigenvalues and, optionally, eigenvectors
|
|
292
|
+
* of a real symmetric matrix A in packed storage. Eigenvalues/vectors
|
|
293
|
+
* can be selected by specifying either a range of values or a range of
|
|
294
|
+
* indices for the desired eigenvalues.
|
|
295
|
+
*
|
|
296
|
+
|
|
297
|
+
* Arguments
|
|
298
|
+
* =========
|
|
299
|
+
*
|
|
300
|
+
* JOBZ (input) CHARACTER*1
|
|
301
|
+
* = 'N': Compute eigenvalues only;
|
|
302
|
+
* = 'V': Compute eigenvalues and eigenvectors.
|
|
303
|
+
*
|
|
304
|
+
* RANGE (input) CHARACTER*1
|
|
305
|
+
* = 'A': all eigenvalues will be found;
|
|
306
|
+
* = 'V': all eigenvalues in the half-open interval (VL,VU]
|
|
307
|
+
* will be found;
|
|
308
|
+
* = 'I': the IL-th through IU-th eigenvalues will be found.
|
|
309
|
+
*
|
|
310
|
+
* UPLO (input) CHARACTER*1
|
|
311
|
+
* = 'U': Upper triangle of A is stored;
|
|
312
|
+
* = 'L': Lower triangle of A is stored.
|
|
313
|
+
*
|
|
314
|
+
* N (input) INTEGER
|
|
315
|
+
* The order of the matrix A. N >= 0.
|
|
316
|
+
*
|
|
317
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
318
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
319
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
320
|
+
* is stored in the array AP as follows:
|
|
321
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
322
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
323
|
+
*
|
|
324
|
+
* On exit, AP is overwritten by values generated during the
|
|
325
|
+
* reduction to tridiagonal form. If UPLO = 'U', the diagonal
|
|
326
|
+
* and first superdiagonal of the tridiagonal matrix T overwrite
|
|
327
|
+
* the corresponding elements of A, and if UPLO = 'L', the
|
|
328
|
+
* diagonal and first subdiagonal of T overwrite the
|
|
329
|
+
* corresponding elements of A.
|
|
330
|
+
*
|
|
331
|
+
* VL (input) REAL
|
|
332
|
+
* VU (input) REAL
|
|
333
|
+
* If RANGE='V', the lower and upper bounds of the interval to
|
|
334
|
+
* be searched for eigenvalues. VL < VU.
|
|
335
|
+
* Not referenced if RANGE = 'A' or 'I'.
|
|
336
|
+
*
|
|
337
|
+
* IL (input) INTEGER
|
|
338
|
+
* IU (input) INTEGER
|
|
339
|
+
* If RANGE='I', the indices (in ascending order) of the
|
|
340
|
+
* smallest and largest eigenvalues to be returned.
|
|
341
|
+
* 1 <= IL <= IU <= N, if N > 0; IL = 1 and IU = 0 if N = 0.
|
|
342
|
+
* Not referenced if RANGE = 'A' or 'V'.
|
|
343
|
+
*
|
|
344
|
+
* ABSTOL (input) REAL
|
|
345
|
+
* The absolute error tolerance for the eigenvalues.
|
|
346
|
+
* An approximate eigenvalue is accepted as converged
|
|
347
|
+
* when it is determined to lie in an interval [a,b]
|
|
348
|
+
* of width less than or equal to
|
|
349
|
+
*
|
|
350
|
+
* ABSTOL + EPS * max( |a|,|b| ) ,
|
|
351
|
+
*
|
|
352
|
+
* where EPS is the machine precision. If ABSTOL is less than
|
|
353
|
+
* or equal to zero, then EPS*|T| will be used in its place,
|
|
354
|
+
* where |T| is the 1-norm of the tridiagonal matrix obtained
|
|
355
|
+
* by reducing AP to tridiagonal form.
|
|
356
|
+
*
|
|
357
|
+
* Eigenvalues will be computed most accurately when ABSTOL is
|
|
358
|
+
* set to twice the underflow threshold 2*SLAMCH('S'), not zero.
|
|
359
|
+
* If this routine returns with INFO>0, indicating that some
|
|
360
|
+
* eigenvectors did not converge, try setting ABSTOL to
|
|
361
|
+
* 2*SLAMCH('S').
|
|
362
|
+
*
|
|
363
|
+
* See "Computing Small Singular Values of Bidiagonal Matrices
|
|
364
|
+
* with Guaranteed High Relative Accuracy," by Demmel and
|
|
365
|
+
* Kahan, LAPACK Working Note #3.
|
|
366
|
+
*
|
|
367
|
+
* M (output) INTEGER
|
|
368
|
+
* The total number of eigenvalues found. 0 <= M <= N.
|
|
369
|
+
* If RANGE = 'A', M = N, and if RANGE = 'I', M = IU-IL+1.
|
|
370
|
+
*
|
|
371
|
+
* W (output) REAL array, dimension (N)
|
|
372
|
+
* If INFO = 0, the selected eigenvalues in ascending order.
|
|
373
|
+
*
|
|
374
|
+
* Z (output) REAL array, dimension (LDZ, max(1,M))
|
|
375
|
+
* If JOBZ = 'V', then if INFO = 0, the first M columns of Z
|
|
376
|
+
* contain the orthonormal eigenvectors of the matrix A
|
|
377
|
+
* corresponding to the selected eigenvalues, with the i-th
|
|
378
|
+
* column of Z holding the eigenvector associated with W(i).
|
|
379
|
+
* If an eigenvector fails to converge, then that column of Z
|
|
380
|
+
* contains the latest approximation to the eigenvector, and the
|
|
381
|
+
* index of the eigenvector is returned in IFAIL.
|
|
382
|
+
* If JOBZ = 'N', then Z is not referenced.
|
|
383
|
+
* Note: the user must ensure that at least max(1,M) columns are
|
|
384
|
+
* supplied in the array Z; if RANGE = 'V', the exact value of M
|
|
385
|
+
* is not known in advance and an upper bound must be used.
|
|
386
|
+
*
|
|
387
|
+
* LDZ (input) INTEGER
|
|
388
|
+
* The leading dimension of the array Z. LDZ >= 1, and if
|
|
389
|
+
* JOBZ = 'V', LDZ >= max(1,N).
|
|
390
|
+
*
|
|
391
|
+
* WORK (workspace) REAL array, dimension (8*N)
|
|
392
|
+
*
|
|
393
|
+
* IWORK (workspace) INTEGER array, dimension (5*N)
|
|
394
|
+
*
|
|
395
|
+
* IFAIL (output) INTEGER array, dimension (N)
|
|
396
|
+
* If JOBZ = 'V', then if INFO = 0, the first M elements of
|
|
397
|
+
* IFAIL are zero. If INFO > 0, then IFAIL contains the
|
|
398
|
+
* indices of the eigenvectors that failed to converge.
|
|
399
|
+
* If JOBZ = 'N', then IFAIL is not referenced.
|
|
400
|
+
*
|
|
401
|
+
* INFO (output) INTEGER
|
|
402
|
+
* = 0: successful exit
|
|
403
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
404
|
+
* > 0: if INFO = i, then i eigenvectors failed to converge.
|
|
405
|
+
* Their indices are stored in array IFAIL.
|
|
406
|
+
*
|
|
407
|
+
|
|
408
|
+
* =====================================================================
|
|
409
|
+
*
|
|
410
|
+
|
|
411
|
+
|
|
412
|
+
</PRE>
|
|
413
|
+
<A HREF="#top">go to the page top</A>
|
|
414
|
+
|
|
415
|
+
<A NAME="sspgst"></A>
|
|
416
|
+
<H2>sspgst</H2>
|
|
417
|
+
<PRE>
|
|
418
|
+
USAGE:
|
|
419
|
+
info, ap = NumRu::Lapack.sspgst( itype, uplo, n, ap, bp, [:usage => usage, :help => help])
|
|
420
|
+
|
|
421
|
+
|
|
422
|
+
FORTRAN MANUAL
|
|
423
|
+
SUBROUTINE SSPGST( ITYPE, UPLO, N, AP, BP, INFO )
|
|
424
|
+
|
|
425
|
+
* Purpose
|
|
426
|
+
* =======
|
|
427
|
+
*
|
|
428
|
+
* SSPGST reduces a real symmetric-definite generalized eigenproblem
|
|
429
|
+
* to standard form, using packed storage.
|
|
430
|
+
*
|
|
431
|
+
* If ITYPE = 1, the problem is A*x = lambda*B*x,
|
|
432
|
+
* and A is overwritten by inv(U**T)*A*inv(U) or inv(L)*A*inv(L**T)
|
|
433
|
+
*
|
|
434
|
+
* If ITYPE = 2 or 3, the problem is A*B*x = lambda*x or
|
|
435
|
+
* B*A*x = lambda*x, and A is overwritten by U*A*U**T or L**T*A*L.
|
|
436
|
+
*
|
|
437
|
+
* B must have been previously factorized as U**T*U or L*L**T by SPPTRF.
|
|
438
|
+
*
|
|
439
|
+
|
|
440
|
+
* Arguments
|
|
441
|
+
* =========
|
|
442
|
+
*
|
|
443
|
+
* ITYPE (input) INTEGER
|
|
444
|
+
* = 1: compute inv(U**T)*A*inv(U) or inv(L)*A*inv(L**T);
|
|
445
|
+
* = 2 or 3: compute U*A*U**T or L**T*A*L.
|
|
446
|
+
*
|
|
447
|
+
* UPLO (input) CHARACTER*1
|
|
448
|
+
* = 'U': Upper triangle of A is stored and B is factored as
|
|
449
|
+
* U**T*U;
|
|
450
|
+
* = 'L': Lower triangle of A is stored and B is factored as
|
|
451
|
+
* L*L**T.
|
|
452
|
+
*
|
|
453
|
+
* N (input) INTEGER
|
|
454
|
+
* The order of the matrices A and B. N >= 0.
|
|
455
|
+
*
|
|
456
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
457
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
458
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
459
|
+
* is stored in the array AP as follows:
|
|
460
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
461
|
+
* if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.
|
|
462
|
+
*
|
|
463
|
+
* On exit, if INFO = 0, the transformed matrix, stored in the
|
|
464
|
+
* same format as A.
|
|
465
|
+
*
|
|
466
|
+
* BP (input) REAL array, dimension (N*(N+1)/2)
|
|
467
|
+
* The triangular factor from the Cholesky factorization of B,
|
|
468
|
+
* stored in the same format as A, as returned by SPPTRF.
|
|
469
|
+
*
|
|
470
|
+
* INFO (output) INTEGER
|
|
471
|
+
* = 0: successful exit
|
|
472
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
473
|
+
*
|
|
474
|
+
|
|
475
|
+
* =====================================================================
|
|
476
|
+
*
|
|
477
|
+
|
|
478
|
+
|
|
479
|
+
</PRE>
|
|
480
|
+
<A HREF="#top">go to the page top</A>
|
|
481
|
+
|
|
482
|
+
<A NAME="sspgv"></A>
|
|
483
|
+
<H2>sspgv</H2>
|
|
484
|
+
<PRE>
|
|
485
|
+
USAGE:
|
|
486
|
+
w, z, info, ap, bp = NumRu::Lapack.sspgv( itype, jobz, uplo, ap, bp, [:usage => usage, :help => help])
|
|
487
|
+
|
|
488
|
+
|
|
489
|
+
FORTRAN MANUAL
|
|
490
|
+
SUBROUTINE SSPGV( ITYPE, JOBZ, UPLO, N, AP, BP, W, Z, LDZ, WORK, INFO )
|
|
491
|
+
|
|
492
|
+
* Purpose
|
|
493
|
+
* =======
|
|
494
|
+
*
|
|
495
|
+
* SSPGV computes all the eigenvalues and, optionally, the eigenvectors
|
|
496
|
+
* of a real generalized symmetric-definite eigenproblem, of the form
|
|
497
|
+
* A*x=(lambda)*B*x, A*Bx=(lambda)*x, or B*A*x=(lambda)*x.
|
|
498
|
+
* Here A and B are assumed to be symmetric, stored in packed format,
|
|
499
|
+
* and B is also positive definite.
|
|
500
|
+
*
|
|
501
|
+
|
|
502
|
+
* Arguments
|
|
503
|
+
* =========
|
|
504
|
+
*
|
|
505
|
+
* ITYPE (input) INTEGER
|
|
506
|
+
* Specifies the problem type to be solved:
|
|
507
|
+
* = 1: A*x = (lambda)*B*x
|
|
508
|
+
* = 2: A*B*x = (lambda)*x
|
|
509
|
+
* = 3: B*A*x = (lambda)*x
|
|
510
|
+
*
|
|
511
|
+
* JOBZ (input) CHARACTER*1
|
|
512
|
+
* = 'N': Compute eigenvalues only;
|
|
513
|
+
* = 'V': Compute eigenvalues and eigenvectors.
|
|
514
|
+
*
|
|
515
|
+
* UPLO (input) CHARACTER*1
|
|
516
|
+
* = 'U': Upper triangles of A and B are stored;
|
|
517
|
+
* = 'L': Lower triangles of A and B are stored.
|
|
518
|
+
*
|
|
519
|
+
* N (input) INTEGER
|
|
520
|
+
* The order of the matrices A and B. N >= 0.
|
|
521
|
+
*
|
|
522
|
+
* AP (input/output) REAL array, dimension
|
|
523
|
+
* (N*(N+1)/2)
|
|
524
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
525
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
526
|
+
* is stored in the array AP as follows:
|
|
527
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
528
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
529
|
+
*
|
|
530
|
+
* On exit, the contents of AP are destroyed.
|
|
531
|
+
*
|
|
532
|
+
* BP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
533
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
534
|
+
* B, packed columnwise in a linear array. The j-th column of B
|
|
535
|
+
* is stored in the array BP as follows:
|
|
536
|
+
* if UPLO = 'U', BP(i + (j-1)*j/2) = B(i,j) for 1<=i<=j;
|
|
537
|
+
* if UPLO = 'L', BP(i + (j-1)*(2*n-j)/2) = B(i,j) for j<=i<=n.
|
|
538
|
+
*
|
|
539
|
+
* On exit, the triangular factor U or L from the Cholesky
|
|
540
|
+
* factorization B = U**T*U or B = L*L**T, in the same storage
|
|
541
|
+
* format as B.
|
|
542
|
+
*
|
|
543
|
+
* W (output) REAL array, dimension (N)
|
|
544
|
+
* If INFO = 0, the eigenvalues in ascending order.
|
|
545
|
+
*
|
|
546
|
+
* Z (output) REAL array, dimension (LDZ, N)
|
|
547
|
+
* If JOBZ = 'V', then if INFO = 0, Z contains the matrix Z of
|
|
548
|
+
* eigenvectors. The eigenvectors are normalized as follows:
|
|
549
|
+
* if ITYPE = 1 or 2, Z**T*B*Z = I;
|
|
550
|
+
* if ITYPE = 3, Z**T*inv(B)*Z = I.
|
|
551
|
+
* If JOBZ = 'N', then Z is not referenced.
|
|
552
|
+
*
|
|
553
|
+
* LDZ (input) INTEGER
|
|
554
|
+
* The leading dimension of the array Z. LDZ >= 1, and if
|
|
555
|
+
* JOBZ = 'V', LDZ >= max(1,N).
|
|
556
|
+
*
|
|
557
|
+
* WORK (workspace) REAL array, dimension (3*N)
|
|
558
|
+
*
|
|
559
|
+
* INFO (output) INTEGER
|
|
560
|
+
* = 0: successful exit
|
|
561
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
562
|
+
* > 0: SPPTRF or SSPEV returned an error code:
|
|
563
|
+
* <= N: if INFO = i, SSPEV failed to converge;
|
|
564
|
+
* i off-diagonal elements of an intermediate
|
|
565
|
+
* tridiagonal form did not converge to zero.
|
|
566
|
+
* > N: if INFO = n + i, for 1 <= i <= n, then the leading
|
|
567
|
+
* minor of order i of B is not positive definite.
|
|
568
|
+
* The factorization of B could not be completed and
|
|
569
|
+
* no eigenvalues or eigenvectors were computed.
|
|
570
|
+
*
|
|
571
|
+
|
|
572
|
+
* =====================================================================
|
|
573
|
+
*
|
|
574
|
+
* .. Local Scalars ..
|
|
575
|
+
LOGICAL UPPER, WANTZ
|
|
576
|
+
CHARACTER TRANS
|
|
577
|
+
INTEGER J, NEIG
|
|
578
|
+
* ..
|
|
579
|
+
* .. External Functions ..
|
|
580
|
+
LOGICAL LSAME
|
|
581
|
+
EXTERNAL LSAME
|
|
582
|
+
* ..
|
|
583
|
+
* .. External Subroutines ..
|
|
584
|
+
EXTERNAL SPPTRF, SSPEV, SSPGST, STPMV, STPSV, XERBLA
|
|
585
|
+
* ..
|
|
586
|
+
|
|
587
|
+
|
|
588
|
+
</PRE>
|
|
589
|
+
<A HREF="#top">go to the page top</A>
|
|
590
|
+
|
|
591
|
+
<A NAME="sspgvd"></A>
|
|
592
|
+
<H2>sspgvd</H2>
|
|
593
|
+
<PRE>
|
|
594
|
+
USAGE:
|
|
595
|
+
w, z, work, iwork, info, ap, bp = NumRu::Lapack.sspgvd( itype, jobz, uplo, ap, bp, [:lwork => lwork, :liwork => liwork, :usage => usage, :help => help])
|
|
596
|
+
|
|
597
|
+
|
|
598
|
+
FORTRAN MANUAL
|
|
599
|
+
SUBROUTINE SSPGVD( ITYPE, JOBZ, UPLO, N, AP, BP, W, Z, LDZ, WORK, LWORK, IWORK, LIWORK, INFO )
|
|
600
|
+
|
|
601
|
+
* Purpose
|
|
602
|
+
* =======
|
|
603
|
+
*
|
|
604
|
+
* SSPGVD computes all the eigenvalues, and optionally, the eigenvectors
|
|
605
|
+
* of a real generalized symmetric-definite eigenproblem, of the form
|
|
606
|
+
* A*x=(lambda)*B*x, A*Bx=(lambda)*x, or B*A*x=(lambda)*x. Here A and
|
|
607
|
+
* B are assumed to be symmetric, stored in packed format, and B is also
|
|
608
|
+
* positive definite.
|
|
609
|
+
* If eigenvectors are desired, it uses a divide and conquer algorithm.
|
|
610
|
+
*
|
|
611
|
+
* The divide and conquer algorithm makes very mild assumptions about
|
|
612
|
+
* floating point arithmetic. It will work on machines with a guard
|
|
613
|
+
* digit in add/subtract, or on those binary machines without guard
|
|
614
|
+
* digits which subtract like the Cray X-MP, Cray Y-MP, Cray C-90, or
|
|
615
|
+
* Cray-2. It could conceivably fail on hexadecimal or decimal machines
|
|
616
|
+
* without guard digits, but we know of none.
|
|
617
|
+
*
|
|
618
|
+
|
|
619
|
+
* Arguments
|
|
620
|
+
* =========
|
|
621
|
+
*
|
|
622
|
+
* ITYPE (input) INTEGER
|
|
623
|
+
* Specifies the problem type to be solved:
|
|
624
|
+
* = 1: A*x = (lambda)*B*x
|
|
625
|
+
* = 2: A*B*x = (lambda)*x
|
|
626
|
+
* = 3: B*A*x = (lambda)*x
|
|
627
|
+
*
|
|
628
|
+
* JOBZ (input) CHARACTER*1
|
|
629
|
+
* = 'N': Compute eigenvalues only;
|
|
630
|
+
* = 'V': Compute eigenvalues and eigenvectors.
|
|
631
|
+
*
|
|
632
|
+
* UPLO (input) CHARACTER*1
|
|
633
|
+
* = 'U': Upper triangles of A and B are stored;
|
|
634
|
+
* = 'L': Lower triangles of A and B are stored.
|
|
635
|
+
*
|
|
636
|
+
* N (input) INTEGER
|
|
637
|
+
* The order of the matrices A and B. N >= 0.
|
|
638
|
+
*
|
|
639
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
640
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
641
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
642
|
+
* is stored in the array AP as follows:
|
|
643
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
644
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
645
|
+
*
|
|
646
|
+
* On exit, the contents of AP are destroyed.
|
|
647
|
+
*
|
|
648
|
+
* BP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
649
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
650
|
+
* B, packed columnwise in a linear array. The j-th column of B
|
|
651
|
+
* is stored in the array BP as follows:
|
|
652
|
+
* if UPLO = 'U', BP(i + (j-1)*j/2) = B(i,j) for 1<=i<=j;
|
|
653
|
+
* if UPLO = 'L', BP(i + (j-1)*(2*n-j)/2) = B(i,j) for j<=i<=n.
|
|
654
|
+
*
|
|
655
|
+
* On exit, the triangular factor U or L from the Cholesky
|
|
656
|
+
* factorization B = U**T*U or B = L*L**T, in the same storage
|
|
657
|
+
* format as B.
|
|
658
|
+
*
|
|
659
|
+
* W (output) REAL array, dimension (N)
|
|
660
|
+
* If INFO = 0, the eigenvalues in ascending order.
|
|
661
|
+
*
|
|
662
|
+
* Z (output) REAL array, dimension (LDZ, N)
|
|
663
|
+
* If JOBZ = 'V', then if INFO = 0, Z contains the matrix Z of
|
|
664
|
+
* eigenvectors. The eigenvectors are normalized as follows:
|
|
665
|
+
* if ITYPE = 1 or 2, Z**T*B*Z = I;
|
|
666
|
+
* if ITYPE = 3, Z**T*inv(B)*Z = I.
|
|
667
|
+
* If JOBZ = 'N', then Z is not referenced.
|
|
668
|
+
*
|
|
669
|
+
* LDZ (input) INTEGER
|
|
670
|
+
* The leading dimension of the array Z. LDZ >= 1, and if
|
|
671
|
+
* JOBZ = 'V', LDZ >= max(1,N).
|
|
672
|
+
*
|
|
673
|
+
* WORK (workspace/output) REAL array, dimension (MAX(1,LWORK))
|
|
674
|
+
* On exit, if INFO = 0, WORK(1) returns the required LWORK.
|
|
675
|
+
*
|
|
676
|
+
* LWORK (input) INTEGER
|
|
677
|
+
* The dimension of the array WORK.
|
|
678
|
+
* If N <= 1, LWORK >= 1.
|
|
679
|
+
* If JOBZ = 'N' and N > 1, LWORK >= 2*N.
|
|
680
|
+
* If JOBZ = 'V' and N > 1, LWORK >= 1 + 6*N + 2*N**2.
|
|
681
|
+
*
|
|
682
|
+
* If LWORK = -1, then a workspace query is assumed; the routine
|
|
683
|
+
* only calculates the required sizes of the WORK and IWORK
|
|
684
|
+
* arrays, returns these values as the first entries of the WORK
|
|
685
|
+
* and IWORK arrays, and no error message related to LWORK or
|
|
686
|
+
* LIWORK is issued by XERBLA.
|
|
687
|
+
*
|
|
688
|
+
* IWORK (workspace/output) INTEGER array, dimension (MAX(1,LIWORK))
|
|
689
|
+
* On exit, if INFO = 0, IWORK(1) returns the required LIWORK.
|
|
690
|
+
*
|
|
691
|
+
* LIWORK (input) INTEGER
|
|
692
|
+
* The dimension of the array IWORK.
|
|
693
|
+
* If JOBZ = 'N' or N <= 1, LIWORK >= 1.
|
|
694
|
+
* If JOBZ = 'V' and N > 1, LIWORK >= 3 + 5*N.
|
|
695
|
+
*
|
|
696
|
+
* If LIWORK = -1, then a workspace query is assumed; the
|
|
697
|
+
* routine only calculates the required sizes of the WORK and
|
|
698
|
+
* IWORK arrays, returns these values as the first entries of
|
|
699
|
+
* the WORK and IWORK arrays, and no error message related to
|
|
700
|
+
* LWORK or LIWORK is issued by XERBLA.
|
|
701
|
+
*
|
|
702
|
+
* INFO (output) INTEGER
|
|
703
|
+
* = 0: successful exit
|
|
704
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
705
|
+
* > 0: SPPTRF or SSPEVD returned an error code:
|
|
706
|
+
* <= N: if INFO = i, SSPEVD failed to converge;
|
|
707
|
+
* i off-diagonal elements of an intermediate
|
|
708
|
+
* tridiagonal form did not converge to zero;
|
|
709
|
+
* > N: if INFO = N + i, for 1 <= i <= N, then the leading
|
|
710
|
+
* minor of order i of B is not positive definite.
|
|
711
|
+
* The factorization of B could not be completed and
|
|
712
|
+
* no eigenvalues or eigenvectors were computed.
|
|
713
|
+
*
|
|
714
|
+
|
|
715
|
+
* Further Details
|
|
716
|
+
* ===============
|
|
717
|
+
*
|
|
718
|
+
* Based on contributions by
|
|
719
|
+
* Mark Fahey, Department of Mathematics, Univ. of Kentucky, USA
|
|
720
|
+
*
|
|
721
|
+
* =====================================================================
|
|
722
|
+
*
|
|
723
|
+
|
|
724
|
+
|
|
725
|
+
</PRE>
|
|
726
|
+
<A HREF="#top">go to the page top</A>
|
|
727
|
+
|
|
728
|
+
<A NAME="sspgvx"></A>
|
|
729
|
+
<H2>sspgvx</H2>
|
|
730
|
+
<PRE>
|
|
731
|
+
USAGE:
|
|
732
|
+
m, w, z, ifail, info, ap, bp = NumRu::Lapack.sspgvx( itype, jobz, range, uplo, ap, bp, vl, vu, il, iu, abstol, [:usage => usage, :help => help])
|
|
733
|
+
|
|
734
|
+
|
|
735
|
+
FORTRAN MANUAL
|
|
736
|
+
SUBROUTINE SSPGVX( ITYPE, JOBZ, RANGE, UPLO, N, AP, BP, VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, IWORK, IFAIL, INFO )
|
|
737
|
+
|
|
738
|
+
* Purpose
|
|
739
|
+
* =======
|
|
740
|
+
*
|
|
741
|
+
* SSPGVX computes selected eigenvalues, and optionally, eigenvectors
|
|
742
|
+
* of a real generalized symmetric-definite eigenproblem, of the form
|
|
743
|
+
* A*x=(lambda)*B*x, A*Bx=(lambda)*x, or B*A*x=(lambda)*x. Here A
|
|
744
|
+
* and B are assumed to be symmetric, stored in packed storage, and B
|
|
745
|
+
* is also positive definite. Eigenvalues and eigenvectors can be
|
|
746
|
+
* selected by specifying either a range of values or a range of indices
|
|
747
|
+
* for the desired eigenvalues.
|
|
748
|
+
*
|
|
749
|
+
|
|
750
|
+
* Arguments
|
|
751
|
+
* =========
|
|
752
|
+
*
|
|
753
|
+
* ITYPE (input) INTEGER
|
|
754
|
+
* Specifies the problem type to be solved:
|
|
755
|
+
* = 1: A*x = (lambda)*B*x
|
|
756
|
+
* = 2: A*B*x = (lambda)*x
|
|
757
|
+
* = 3: B*A*x = (lambda)*x
|
|
758
|
+
*
|
|
759
|
+
* JOBZ (input) CHARACTER*1
|
|
760
|
+
* = 'N': Compute eigenvalues only;
|
|
761
|
+
* = 'V': Compute eigenvalues and eigenvectors.
|
|
762
|
+
*
|
|
763
|
+
* RANGE (input) CHARACTER*1
|
|
764
|
+
* = 'A': all eigenvalues will be found.
|
|
765
|
+
* = 'V': all eigenvalues in the half-open interval (VL,VU]
|
|
766
|
+
* will be found.
|
|
767
|
+
* = 'I': the IL-th through IU-th eigenvalues will be found.
|
|
768
|
+
*
|
|
769
|
+
* UPLO (input) CHARACTER*1
|
|
770
|
+
* = 'U': Upper triangle of A and B are stored;
|
|
771
|
+
* = 'L': Lower triangle of A and B are stored.
|
|
772
|
+
*
|
|
773
|
+
* N (input) INTEGER
|
|
774
|
+
* The order of the matrix pencil (A,B). N >= 0.
|
|
775
|
+
*
|
|
776
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
777
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
778
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
779
|
+
* is stored in the array AP as follows:
|
|
780
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
781
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
782
|
+
*
|
|
783
|
+
* On exit, the contents of AP are destroyed.
|
|
784
|
+
*
|
|
785
|
+
* BP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
786
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
787
|
+
* B, packed columnwise in a linear array. The j-th column of B
|
|
788
|
+
* is stored in the array BP as follows:
|
|
789
|
+
* if UPLO = 'U', BP(i + (j-1)*j/2) = B(i,j) for 1<=i<=j;
|
|
790
|
+
* if UPLO = 'L', BP(i + (j-1)*(2*n-j)/2) = B(i,j) for j<=i<=n.
|
|
791
|
+
*
|
|
792
|
+
* On exit, the triangular factor U or L from the Cholesky
|
|
793
|
+
* factorization B = U**T*U or B = L*L**T, in the same storage
|
|
794
|
+
* format as B.
|
|
795
|
+
*
|
|
796
|
+
* VL (input) REAL
|
|
797
|
+
* VU (input) REAL
|
|
798
|
+
* If RANGE='V', the lower and upper bounds of the interval to
|
|
799
|
+
* be searched for eigenvalues. VL < VU.
|
|
800
|
+
* Not referenced if RANGE = 'A' or 'I'.
|
|
801
|
+
*
|
|
802
|
+
* IL (input) INTEGER
|
|
803
|
+
* IU (input) INTEGER
|
|
804
|
+
* If RANGE='I', the indices (in ascending order) of the
|
|
805
|
+
* smallest and largest eigenvalues to be returned.
|
|
806
|
+
* 1 <= IL <= IU <= N, if N > 0; IL = 1 and IU = 0 if N = 0.
|
|
807
|
+
* Not referenced if RANGE = 'A' or 'V'.
|
|
808
|
+
*
|
|
809
|
+
* ABSTOL (input) REAL
|
|
810
|
+
* The absolute error tolerance for the eigenvalues.
|
|
811
|
+
* An approximate eigenvalue is accepted as converged
|
|
812
|
+
* when it is determined to lie in an interval [a,b]
|
|
813
|
+
* of width less than or equal to
|
|
814
|
+
*
|
|
815
|
+
* ABSTOL + EPS * max( |a|,|b| ) ,
|
|
816
|
+
*
|
|
817
|
+
* where EPS is the machine precision. If ABSTOL is less than
|
|
818
|
+
* or equal to zero, then EPS*|T| will be used in its place,
|
|
819
|
+
* where |T| is the 1-norm of the tridiagonal matrix obtained
|
|
820
|
+
* by reducing A to tridiagonal form.
|
|
821
|
+
*
|
|
822
|
+
* Eigenvalues will be computed most accurately when ABSTOL is
|
|
823
|
+
* set to twice the underflow threshold 2*SLAMCH('S'), not zero.
|
|
824
|
+
* If this routine returns with INFO>0, indicating that some
|
|
825
|
+
* eigenvectors did not converge, try setting ABSTOL to
|
|
826
|
+
* 2*SLAMCH('S').
|
|
827
|
+
*
|
|
828
|
+
* M (output) INTEGER
|
|
829
|
+
* The total number of eigenvalues found. 0 <= M <= N.
|
|
830
|
+
* If RANGE = 'A', M = N, and if RANGE = 'I', M = IU-IL+1.
|
|
831
|
+
*
|
|
832
|
+
* W (output) REAL array, dimension (N)
|
|
833
|
+
* On normal exit, the first M elements contain the selected
|
|
834
|
+
* eigenvalues in ascending order.
|
|
835
|
+
*
|
|
836
|
+
* Z (output) REAL array, dimension (LDZ, max(1,M))
|
|
837
|
+
* If JOBZ = 'N', then Z is not referenced.
|
|
838
|
+
* If JOBZ = 'V', then if INFO = 0, the first M columns of Z
|
|
839
|
+
* contain the orthonormal eigenvectors of the matrix A
|
|
840
|
+
* corresponding to the selected eigenvalues, with the i-th
|
|
841
|
+
* column of Z holding the eigenvector associated with W(i).
|
|
842
|
+
* The eigenvectors are normalized as follows:
|
|
843
|
+
* if ITYPE = 1 or 2, Z**T*B*Z = I;
|
|
844
|
+
* if ITYPE = 3, Z**T*inv(B)*Z = I.
|
|
845
|
+
*
|
|
846
|
+
* If an eigenvector fails to converge, then that column of Z
|
|
847
|
+
* contains the latest approximation to the eigenvector, and the
|
|
848
|
+
* index of the eigenvector is returned in IFAIL.
|
|
849
|
+
* Note: the user must ensure that at least max(1,M) columns are
|
|
850
|
+
* supplied in the array Z; if RANGE = 'V', the exact value of M
|
|
851
|
+
* is not known in advance and an upper bound must be used.
|
|
852
|
+
*
|
|
853
|
+
* LDZ (input) INTEGER
|
|
854
|
+
* The leading dimension of the array Z. LDZ >= 1, and if
|
|
855
|
+
* JOBZ = 'V', LDZ >= max(1,N).
|
|
856
|
+
*
|
|
857
|
+
* WORK (workspace) REAL array, dimension (8*N)
|
|
858
|
+
*
|
|
859
|
+
* IWORK (workspace) INTEGER array, dimension (5*N)
|
|
860
|
+
*
|
|
861
|
+
* IFAIL (output) INTEGER array, dimension (N)
|
|
862
|
+
* If JOBZ = 'V', then if INFO = 0, the first M elements of
|
|
863
|
+
* IFAIL are zero. If INFO > 0, then IFAIL contains the
|
|
864
|
+
* indices of the eigenvectors that failed to converge.
|
|
865
|
+
* If JOBZ = 'N', then IFAIL is not referenced.
|
|
866
|
+
*
|
|
867
|
+
* INFO (output) INTEGER
|
|
868
|
+
* = 0: successful exit
|
|
869
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
870
|
+
* > 0: SPPTRF or SSPEVX returned an error code:
|
|
871
|
+
* <= N: if INFO = i, SSPEVX failed to converge;
|
|
872
|
+
* i eigenvectors failed to converge. Their indices
|
|
873
|
+
* are stored in array IFAIL.
|
|
874
|
+
* > N: if INFO = N + i, for 1 <= i <= N, then the leading
|
|
875
|
+
* minor of order i of B is not positive definite.
|
|
876
|
+
* The factorization of B could not be completed and
|
|
877
|
+
* no eigenvalues or eigenvectors were computed.
|
|
878
|
+
*
|
|
879
|
+
|
|
880
|
+
* Further Details
|
|
881
|
+
* ===============
|
|
882
|
+
*
|
|
883
|
+
* Based on contributions by
|
|
884
|
+
* Mark Fahey, Department of Mathematics, Univ. of Kentucky, USA
|
|
885
|
+
*
|
|
886
|
+
* =====================================================================
|
|
887
|
+
*
|
|
888
|
+
* .. Local Scalars ..
|
|
889
|
+
LOGICAL ALLEIG, INDEIG, UPPER, VALEIG, WANTZ
|
|
890
|
+
CHARACTER TRANS
|
|
891
|
+
INTEGER J
|
|
892
|
+
* ..
|
|
893
|
+
* .. External Functions ..
|
|
894
|
+
LOGICAL LSAME
|
|
895
|
+
EXTERNAL LSAME
|
|
896
|
+
* ..
|
|
897
|
+
* .. External Subroutines ..
|
|
898
|
+
EXTERNAL SPPTRF, SSPEVX, SSPGST, STPMV, STPSV, XERBLA
|
|
899
|
+
* ..
|
|
900
|
+
* .. Intrinsic Functions ..
|
|
901
|
+
INTRINSIC MIN
|
|
902
|
+
* ..
|
|
903
|
+
|
|
904
|
+
|
|
905
|
+
</PRE>
|
|
906
|
+
<A HREF="#top">go to the page top</A>
|
|
907
|
+
|
|
908
|
+
<A NAME="ssprfs"></A>
|
|
909
|
+
<H2>ssprfs</H2>
|
|
910
|
+
<PRE>
|
|
911
|
+
USAGE:
|
|
912
|
+
ferr, berr, info, x = NumRu::Lapack.ssprfs( uplo, ap, afp, ipiv, b, x, [:usage => usage, :help => help])
|
|
913
|
+
|
|
914
|
+
|
|
915
|
+
FORTRAN MANUAL
|
|
916
|
+
SUBROUTINE SSPRFS( UPLO, N, NRHS, AP, AFP, IPIV, B, LDB, X, LDX, FERR, BERR, WORK, IWORK, INFO )
|
|
917
|
+
|
|
918
|
+
* Purpose
|
|
919
|
+
* =======
|
|
920
|
+
*
|
|
921
|
+
* SSPRFS improves the computed solution to a system of linear
|
|
922
|
+
* equations when the coefficient matrix is symmetric indefinite
|
|
923
|
+
* and packed, and provides error bounds and backward error estimates
|
|
924
|
+
* for the solution.
|
|
925
|
+
*
|
|
926
|
+
|
|
927
|
+
* Arguments
|
|
928
|
+
* =========
|
|
929
|
+
*
|
|
930
|
+
* UPLO (input) CHARACTER*1
|
|
931
|
+
* = 'U': Upper triangle of A is stored;
|
|
932
|
+
* = 'L': Lower triangle of A is stored.
|
|
933
|
+
*
|
|
934
|
+
* N (input) INTEGER
|
|
935
|
+
* The order of the matrix A. N >= 0.
|
|
936
|
+
*
|
|
937
|
+
* NRHS (input) INTEGER
|
|
938
|
+
* The number of right hand sides, i.e., the number of columns
|
|
939
|
+
* of the matrices B and X. NRHS >= 0.
|
|
940
|
+
*
|
|
941
|
+
* AP (input) REAL array, dimension (N*(N+1)/2)
|
|
942
|
+
* The upper or lower triangle of the symmetric matrix A, packed
|
|
943
|
+
* columnwise in a linear array. The j-th column of A is stored
|
|
944
|
+
* in the array AP as follows:
|
|
945
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
946
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
947
|
+
*
|
|
948
|
+
* AFP (input) REAL array, dimension (N*(N+1)/2)
|
|
949
|
+
* The factored form of the matrix A. AFP contains the block
|
|
950
|
+
* diagonal matrix D and the multipliers used to obtain the
|
|
951
|
+
* factor U or L from the factorization A = U*D*U**T or
|
|
952
|
+
* A = L*D*L**T as computed by SSPTRF, stored as a packed
|
|
953
|
+
* triangular matrix.
|
|
954
|
+
*
|
|
955
|
+
* IPIV (input) INTEGER array, dimension (N)
|
|
956
|
+
* Details of the interchanges and the block structure of D
|
|
957
|
+
* as determined by SSPTRF.
|
|
958
|
+
*
|
|
959
|
+
* B (input) REAL array, dimension (LDB,NRHS)
|
|
960
|
+
* The right hand side matrix B.
|
|
961
|
+
*
|
|
962
|
+
* LDB (input) INTEGER
|
|
963
|
+
* The leading dimension of the array B. LDB >= max(1,N).
|
|
964
|
+
*
|
|
965
|
+
* X (input/output) REAL array, dimension (LDX,NRHS)
|
|
966
|
+
* On entry, the solution matrix X, as computed by SSPTRS.
|
|
967
|
+
* On exit, the improved solution matrix X.
|
|
968
|
+
*
|
|
969
|
+
* LDX (input) INTEGER
|
|
970
|
+
* The leading dimension of the array X. LDX >= max(1,N).
|
|
971
|
+
*
|
|
972
|
+
* FERR (output) REAL array, dimension (NRHS)
|
|
973
|
+
* The estimated forward error bound for each solution vector
|
|
974
|
+
* X(j) (the j-th column of the solution matrix X).
|
|
975
|
+
* If XTRUE is the true solution corresponding to X(j), FERR(j)
|
|
976
|
+
* is an estimated upper bound for the magnitude of the largest
|
|
977
|
+
* element in (X(j) - XTRUE) divided by the magnitude of the
|
|
978
|
+
* largest element in X(j). The estimate is as reliable as
|
|
979
|
+
* the estimate for RCOND, and is almost always a slight
|
|
980
|
+
* overestimate of the true error.
|
|
981
|
+
*
|
|
982
|
+
* BERR (output) REAL array, dimension (NRHS)
|
|
983
|
+
* The componentwise relative backward error of each solution
|
|
984
|
+
* vector X(j) (i.e., the smallest relative change in
|
|
985
|
+
* any element of A or B that makes X(j) an exact solution).
|
|
986
|
+
*
|
|
987
|
+
* WORK (workspace) REAL array, dimension (3*N)
|
|
988
|
+
*
|
|
989
|
+
* IWORK (workspace) INTEGER array, dimension (N)
|
|
990
|
+
*
|
|
991
|
+
* INFO (output) INTEGER
|
|
992
|
+
* = 0: successful exit
|
|
993
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
994
|
+
*
|
|
995
|
+
* Internal Parameters
|
|
996
|
+
* ===================
|
|
997
|
+
*
|
|
998
|
+
* ITMAX is the maximum number of steps of iterative refinement.
|
|
999
|
+
*
|
|
1000
|
+
|
|
1001
|
+
* =====================================================================
|
|
1002
|
+
*
|
|
1003
|
+
|
|
1004
|
+
|
|
1005
|
+
</PRE>
|
|
1006
|
+
<A HREF="#top">go to the page top</A>
|
|
1007
|
+
|
|
1008
|
+
<A NAME="sspsv"></A>
|
|
1009
|
+
<H2>sspsv</H2>
|
|
1010
|
+
<PRE>
|
|
1011
|
+
USAGE:
|
|
1012
|
+
ipiv, info, ap, b = NumRu::Lapack.sspsv( uplo, ap, b, [:usage => usage, :help => help])
|
|
1013
|
+
|
|
1014
|
+
|
|
1015
|
+
FORTRAN MANUAL
|
|
1016
|
+
SUBROUTINE SSPSV( UPLO, N, NRHS, AP, IPIV, B, LDB, INFO )
|
|
1017
|
+
|
|
1018
|
+
* Purpose
|
|
1019
|
+
* =======
|
|
1020
|
+
*
|
|
1021
|
+
* SSPSV computes the solution to a real system of linear equations
|
|
1022
|
+
* A * X = B,
|
|
1023
|
+
* where A is an N-by-N symmetric matrix stored in packed format and X
|
|
1024
|
+
* and B are N-by-NRHS matrices.
|
|
1025
|
+
*
|
|
1026
|
+
* The diagonal pivoting method is used to factor A as
|
|
1027
|
+
* A = U * D * U**T, if UPLO = 'U', or
|
|
1028
|
+
* A = L * D * L**T, if UPLO = 'L',
|
|
1029
|
+
* where U (or L) is a product of permutation and unit upper (lower)
|
|
1030
|
+
* triangular matrices, D is symmetric and block diagonal with 1-by-1
|
|
1031
|
+
* and 2-by-2 diagonal blocks. The factored form of A is then used to
|
|
1032
|
+
* solve the system of equations A * X = B.
|
|
1033
|
+
*
|
|
1034
|
+
|
|
1035
|
+
* Arguments
|
|
1036
|
+
* =========
|
|
1037
|
+
*
|
|
1038
|
+
* UPLO (input) CHARACTER*1
|
|
1039
|
+
* = 'U': Upper triangle of A is stored;
|
|
1040
|
+
* = 'L': Lower triangle of A is stored.
|
|
1041
|
+
*
|
|
1042
|
+
* N (input) INTEGER
|
|
1043
|
+
* The number of linear equations, i.e., the order of the
|
|
1044
|
+
* matrix A. N >= 0.
|
|
1045
|
+
*
|
|
1046
|
+
* NRHS (input) INTEGER
|
|
1047
|
+
* The number of right hand sides, i.e., the number of columns
|
|
1048
|
+
* of the matrix B. NRHS >= 0.
|
|
1049
|
+
*
|
|
1050
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
1051
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
1052
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
1053
|
+
* is stored in the array AP as follows:
|
|
1054
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
1055
|
+
* if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.
|
|
1056
|
+
* See below for further details.
|
|
1057
|
+
*
|
|
1058
|
+
* On exit, the block diagonal matrix D and the multipliers used
|
|
1059
|
+
* to obtain the factor U or L from the factorization
|
|
1060
|
+
* A = U*D*U**T or A = L*D*L**T as computed by SSPTRF, stored as
|
|
1061
|
+
* a packed triangular matrix in the same storage format as A.
|
|
1062
|
+
*
|
|
1063
|
+
* IPIV (output) INTEGER array, dimension (N)
|
|
1064
|
+
* Details of the interchanges and the block structure of D, as
|
|
1065
|
+
* determined by SSPTRF. If IPIV(k) > 0, then rows and columns
|
|
1066
|
+
* k and IPIV(k) were interchanged, and D(k,k) is a 1-by-1
|
|
1067
|
+
* diagonal block. If UPLO = 'U' and IPIV(k) = IPIV(k-1) < 0,
|
|
1068
|
+
* then rows and columns k-1 and -IPIV(k) were interchanged and
|
|
1069
|
+
* D(k-1:k,k-1:k) is a 2-by-2 diagonal block. If UPLO = 'L' and
|
|
1070
|
+
* IPIV(k) = IPIV(k+1) < 0, then rows and columns k+1 and
|
|
1071
|
+
* -IPIV(k) were interchanged and D(k:k+1,k:k+1) is a 2-by-2
|
|
1072
|
+
* diagonal block.
|
|
1073
|
+
*
|
|
1074
|
+
* B (input/output) REAL array, dimension (LDB,NRHS)
|
|
1075
|
+
* On entry, the N-by-NRHS right hand side matrix B.
|
|
1076
|
+
* On exit, if INFO = 0, the N-by-NRHS solution matrix X.
|
|
1077
|
+
*
|
|
1078
|
+
* LDB (input) INTEGER
|
|
1079
|
+
* The leading dimension of the array B. LDB >= max(1,N).
|
|
1080
|
+
*
|
|
1081
|
+
* INFO (output) INTEGER
|
|
1082
|
+
* = 0: successful exit
|
|
1083
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1084
|
+
* > 0: if INFO = i, D(i,i) is exactly zero. The factorization
|
|
1085
|
+
* has been completed, but the block diagonal matrix D is
|
|
1086
|
+
* exactly singular, so the solution could not be
|
|
1087
|
+
* computed.
|
|
1088
|
+
*
|
|
1089
|
+
|
|
1090
|
+
* Further Details
|
|
1091
|
+
* ===============
|
|
1092
|
+
*
|
|
1093
|
+
* The packed storage scheme is illustrated by the following example
|
|
1094
|
+
* when N = 4, UPLO = 'U':
|
|
1095
|
+
*
|
|
1096
|
+
* Two-dimensional storage of the symmetric matrix A:
|
|
1097
|
+
*
|
|
1098
|
+
* a11 a12 a13 a14
|
|
1099
|
+
* a22 a23 a24
|
|
1100
|
+
* a33 a34 (aij = aji)
|
|
1101
|
+
* a44
|
|
1102
|
+
*
|
|
1103
|
+
* Packed storage of the upper triangle of A:
|
|
1104
|
+
*
|
|
1105
|
+
* AP = [ a11, a12, a22, a13, a23, a33, a14, a24, a34, a44 ]
|
|
1106
|
+
*
|
|
1107
|
+
* =====================================================================
|
|
1108
|
+
*
|
|
1109
|
+
* .. External Functions ..
|
|
1110
|
+
LOGICAL LSAME
|
|
1111
|
+
EXTERNAL LSAME
|
|
1112
|
+
* ..
|
|
1113
|
+
* .. External Subroutines ..
|
|
1114
|
+
EXTERNAL SSPTRF, SSPTRS, XERBLA
|
|
1115
|
+
* ..
|
|
1116
|
+
* .. Intrinsic Functions ..
|
|
1117
|
+
INTRINSIC MAX
|
|
1118
|
+
* ..
|
|
1119
|
+
|
|
1120
|
+
|
|
1121
|
+
</PRE>
|
|
1122
|
+
<A HREF="#top">go to the page top</A>
|
|
1123
|
+
|
|
1124
|
+
<A NAME="sspsvx"></A>
|
|
1125
|
+
<H2>sspsvx</H2>
|
|
1126
|
+
<PRE>
|
|
1127
|
+
USAGE:
|
|
1128
|
+
x, rcond, ferr, berr, info, afp, ipiv = NumRu::Lapack.sspsvx( fact, uplo, ap, afp, ipiv, b, [:usage => usage, :help => help])
|
|
1129
|
+
|
|
1130
|
+
|
|
1131
|
+
FORTRAN MANUAL
|
|
1132
|
+
SUBROUTINE SSPSVX( FACT, UPLO, N, NRHS, AP, AFP, IPIV, B, LDB, X, LDX, RCOND, FERR, BERR, WORK, IWORK, INFO )
|
|
1133
|
+
|
|
1134
|
+
* Purpose
|
|
1135
|
+
* =======
|
|
1136
|
+
*
|
|
1137
|
+
* SSPSVX uses the diagonal pivoting factorization A = U*D*U**T or
|
|
1138
|
+
* A = L*D*L**T to compute the solution to a real system of linear
|
|
1139
|
+
* equations A * X = B, where A is an N-by-N symmetric matrix stored
|
|
1140
|
+
* in packed format and X and B are N-by-NRHS matrices.
|
|
1141
|
+
*
|
|
1142
|
+
* Error bounds on the solution and a condition estimate are also
|
|
1143
|
+
* provided.
|
|
1144
|
+
*
|
|
1145
|
+
* Description
|
|
1146
|
+
* ===========
|
|
1147
|
+
*
|
|
1148
|
+
* The following steps are performed:
|
|
1149
|
+
*
|
|
1150
|
+
* 1. If FACT = 'N', the diagonal pivoting method is used to factor A as
|
|
1151
|
+
* A = U * D * U**T, if UPLO = 'U', or
|
|
1152
|
+
* A = L * D * L**T, if UPLO = 'L',
|
|
1153
|
+
* where U (or L) is a product of permutation and unit upper (lower)
|
|
1154
|
+
* triangular matrices and D is symmetric and block diagonal with
|
|
1155
|
+
* 1-by-1 and 2-by-2 diagonal blocks.
|
|
1156
|
+
*
|
|
1157
|
+
* 2. If some D(i,i)=0, so that D is exactly singular, then the routine
|
|
1158
|
+
* returns with INFO = i. Otherwise, the factored form of A is used
|
|
1159
|
+
* to estimate the condition number of the matrix A. If the
|
|
1160
|
+
* reciprocal of the condition number is less than machine precision,
|
|
1161
|
+
* INFO = N+1 is returned as a warning, but the routine still goes on
|
|
1162
|
+
* to solve for X and compute error bounds as described below.
|
|
1163
|
+
*
|
|
1164
|
+
* 3. The system of equations is solved for X using the factored form
|
|
1165
|
+
* of A.
|
|
1166
|
+
*
|
|
1167
|
+
* 4. Iterative refinement is applied to improve the computed solution
|
|
1168
|
+
* matrix and calculate error bounds and backward error estimates
|
|
1169
|
+
* for it.
|
|
1170
|
+
*
|
|
1171
|
+
|
|
1172
|
+
* Arguments
|
|
1173
|
+
* =========
|
|
1174
|
+
*
|
|
1175
|
+
* FACT (input) CHARACTER*1
|
|
1176
|
+
* Specifies whether or not the factored form of A has been
|
|
1177
|
+
* supplied on entry.
|
|
1178
|
+
* = 'F': On entry, AFP and IPIV contain the factored form of
|
|
1179
|
+
* A. AP, AFP and IPIV will not be modified.
|
|
1180
|
+
* = 'N': The matrix A will be copied to AFP and factored.
|
|
1181
|
+
*
|
|
1182
|
+
* UPLO (input) CHARACTER*1
|
|
1183
|
+
* = 'U': Upper triangle of A is stored;
|
|
1184
|
+
* = 'L': Lower triangle of A is stored.
|
|
1185
|
+
*
|
|
1186
|
+
* N (input) INTEGER
|
|
1187
|
+
* The number of linear equations, i.e., the order of the
|
|
1188
|
+
* matrix A. N >= 0.
|
|
1189
|
+
*
|
|
1190
|
+
* NRHS (input) INTEGER
|
|
1191
|
+
* The number of right hand sides, i.e., the number of columns
|
|
1192
|
+
* of the matrices B and X. NRHS >= 0.
|
|
1193
|
+
*
|
|
1194
|
+
* AP (input) REAL array, dimension (N*(N+1)/2)
|
|
1195
|
+
* The upper or lower triangle of the symmetric matrix A, packed
|
|
1196
|
+
* columnwise in a linear array. The j-th column of A is stored
|
|
1197
|
+
* in the array AP as follows:
|
|
1198
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
1199
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
1200
|
+
* See below for further details.
|
|
1201
|
+
*
|
|
1202
|
+
* AFP (input or output) REAL array, dimension
|
|
1203
|
+
* (N*(N+1)/2)
|
|
1204
|
+
* If FACT = 'F', then AFP is an input argument and on entry
|
|
1205
|
+
* contains the block diagonal matrix D and the multipliers used
|
|
1206
|
+
* to obtain the factor U or L from the factorization
|
|
1207
|
+
* A = U*D*U**T or A = L*D*L**T as computed by SSPTRF, stored as
|
|
1208
|
+
* a packed triangular matrix in the same storage format as A.
|
|
1209
|
+
*
|
|
1210
|
+
* If FACT = 'N', then AFP is an output argument and on exit
|
|
1211
|
+
* contains the block diagonal matrix D and the multipliers used
|
|
1212
|
+
* to obtain the factor U or L from the factorization
|
|
1213
|
+
* A = U*D*U**T or A = L*D*L**T as computed by SSPTRF, stored as
|
|
1214
|
+
* a packed triangular matrix in the same storage format as A.
|
|
1215
|
+
*
|
|
1216
|
+
* IPIV (input or output) INTEGER array, dimension (N)
|
|
1217
|
+
* If FACT = 'F', then IPIV is an input argument and on entry
|
|
1218
|
+
* contains details of the interchanges and the block structure
|
|
1219
|
+
* of D, as determined by SSPTRF.
|
|
1220
|
+
* If IPIV(k) > 0, then rows and columns k and IPIV(k) were
|
|
1221
|
+
* interchanged and D(k,k) is a 1-by-1 diagonal block.
|
|
1222
|
+
* If UPLO = 'U' and IPIV(k) = IPIV(k-1) < 0, then rows and
|
|
1223
|
+
* columns k-1 and -IPIV(k) were interchanged and D(k-1:k,k-1:k)
|
|
1224
|
+
* is a 2-by-2 diagonal block. If UPLO = 'L' and IPIV(k) =
|
|
1225
|
+
* IPIV(k+1) < 0, then rows and columns k+1 and -IPIV(k) were
|
|
1226
|
+
* interchanged and D(k:k+1,k:k+1) is a 2-by-2 diagonal block.
|
|
1227
|
+
*
|
|
1228
|
+
* If FACT = 'N', then IPIV is an output argument and on exit
|
|
1229
|
+
* contains details of the interchanges and the block structure
|
|
1230
|
+
* of D, as determined by SSPTRF.
|
|
1231
|
+
*
|
|
1232
|
+
* B (input) REAL array, dimension (LDB,NRHS)
|
|
1233
|
+
* The N-by-NRHS right hand side matrix B.
|
|
1234
|
+
*
|
|
1235
|
+
* LDB (input) INTEGER
|
|
1236
|
+
* The leading dimension of the array B. LDB >= max(1,N).
|
|
1237
|
+
*
|
|
1238
|
+
* X (output) REAL array, dimension (LDX,NRHS)
|
|
1239
|
+
* If INFO = 0 or INFO = N+1, the N-by-NRHS solution matrix X.
|
|
1240
|
+
*
|
|
1241
|
+
* LDX (input) INTEGER
|
|
1242
|
+
* The leading dimension of the array X. LDX >= max(1,N).
|
|
1243
|
+
*
|
|
1244
|
+
* RCOND (output) REAL
|
|
1245
|
+
* The estimate of the reciprocal condition number of the matrix
|
|
1246
|
+
* A. If RCOND is less than the machine precision (in
|
|
1247
|
+
* particular, if RCOND = 0), the matrix is singular to working
|
|
1248
|
+
* precision. This condition is indicated by a return code of
|
|
1249
|
+
* INFO > 0.
|
|
1250
|
+
*
|
|
1251
|
+
* FERR (output) REAL array, dimension (NRHS)
|
|
1252
|
+
* The estimated forward error bound for each solution vector
|
|
1253
|
+
* X(j) (the j-th column of the solution matrix X).
|
|
1254
|
+
* If XTRUE is the true solution corresponding to X(j), FERR(j)
|
|
1255
|
+
* is an estimated upper bound for the magnitude of the largest
|
|
1256
|
+
* element in (X(j) - XTRUE) divided by the magnitude of the
|
|
1257
|
+
* largest element in X(j). The estimate is as reliable as
|
|
1258
|
+
* the estimate for RCOND, and is almost always a slight
|
|
1259
|
+
* overestimate of the true error.
|
|
1260
|
+
*
|
|
1261
|
+
* BERR (output) REAL array, dimension (NRHS)
|
|
1262
|
+
* The componentwise relative backward error of each solution
|
|
1263
|
+
* vector X(j) (i.e., the smallest relative change in
|
|
1264
|
+
* any element of A or B that makes X(j) an exact solution).
|
|
1265
|
+
*
|
|
1266
|
+
* WORK (workspace) REAL array, dimension (3*N)
|
|
1267
|
+
*
|
|
1268
|
+
* IWORK (workspace) INTEGER array, dimension (N)
|
|
1269
|
+
*
|
|
1270
|
+
* INFO (output) INTEGER
|
|
1271
|
+
* = 0: successful exit
|
|
1272
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1273
|
+
* > 0: if INFO = i, and i is
|
|
1274
|
+
* <= N: D(i,i) is exactly zero. The factorization
|
|
1275
|
+
* has been completed but the factor D is exactly
|
|
1276
|
+
* singular, so the solution and error bounds could
|
|
1277
|
+
* not be computed. RCOND = 0 is returned.
|
|
1278
|
+
* = N+1: D is nonsingular, but RCOND is less than machine
|
|
1279
|
+
* precision, meaning that the matrix is singular
|
|
1280
|
+
* to working precision. Nevertheless, the
|
|
1281
|
+
* solution and error bounds are computed because
|
|
1282
|
+
* there are a number of situations where the
|
|
1283
|
+
* computed solution can be more accurate than the
|
|
1284
|
+
* value of RCOND would suggest.
|
|
1285
|
+
*
|
|
1286
|
+
|
|
1287
|
+
* Further Details
|
|
1288
|
+
* ===============
|
|
1289
|
+
*
|
|
1290
|
+
* The packed storage scheme is illustrated by the following example
|
|
1291
|
+
* when N = 4, UPLO = 'U':
|
|
1292
|
+
*
|
|
1293
|
+
* Two-dimensional storage of the symmetric matrix A:
|
|
1294
|
+
*
|
|
1295
|
+
* a11 a12 a13 a14
|
|
1296
|
+
* a22 a23 a24
|
|
1297
|
+
* a33 a34 (aij = aji)
|
|
1298
|
+
* a44
|
|
1299
|
+
*
|
|
1300
|
+
* Packed storage of the upper triangle of A:
|
|
1301
|
+
*
|
|
1302
|
+
* AP = [ a11, a12, a22, a13, a23, a33, a14, a24, a34, a44 ]
|
|
1303
|
+
*
|
|
1304
|
+
* =====================================================================
|
|
1305
|
+
*
|
|
1306
|
+
|
|
1307
|
+
|
|
1308
|
+
</PRE>
|
|
1309
|
+
<A HREF="#top">go to the page top</A>
|
|
1310
|
+
|
|
1311
|
+
<A NAME="ssptrd"></A>
|
|
1312
|
+
<H2>ssptrd</H2>
|
|
1313
|
+
<PRE>
|
|
1314
|
+
USAGE:
|
|
1315
|
+
d, e, tau, info, ap = NumRu::Lapack.ssptrd( uplo, ap, [:usage => usage, :help => help])
|
|
1316
|
+
|
|
1317
|
+
|
|
1318
|
+
FORTRAN MANUAL
|
|
1319
|
+
SUBROUTINE SSPTRD( UPLO, N, AP, D, E, TAU, INFO )
|
|
1320
|
+
|
|
1321
|
+
* Purpose
|
|
1322
|
+
* =======
|
|
1323
|
+
*
|
|
1324
|
+
* SSPTRD reduces a real symmetric matrix A stored in packed form to
|
|
1325
|
+
* symmetric tridiagonal form T by an orthogonal similarity
|
|
1326
|
+
* transformation: Q**T * A * Q = T.
|
|
1327
|
+
*
|
|
1328
|
+
|
|
1329
|
+
* Arguments
|
|
1330
|
+
* =========
|
|
1331
|
+
*
|
|
1332
|
+
* UPLO (input) CHARACTER*1
|
|
1333
|
+
* = 'U': Upper triangle of A is stored;
|
|
1334
|
+
* = 'L': Lower triangle of A is stored.
|
|
1335
|
+
*
|
|
1336
|
+
* N (input) INTEGER
|
|
1337
|
+
* The order of the matrix A. N >= 0.
|
|
1338
|
+
*
|
|
1339
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
1340
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
1341
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
1342
|
+
* is stored in the array AP as follows:
|
|
1343
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
1344
|
+
* if UPLO = 'L', AP(i + (j-1)*(2*n-j)/2) = A(i,j) for j<=i<=n.
|
|
1345
|
+
* On exit, if UPLO = 'U', the diagonal and first superdiagonal
|
|
1346
|
+
* of A are overwritten by the corresponding elements of the
|
|
1347
|
+
* tridiagonal matrix T, and the elements above the first
|
|
1348
|
+
* superdiagonal, with the array TAU, represent the orthogonal
|
|
1349
|
+
* matrix Q as a product of elementary reflectors; if UPLO
|
|
1350
|
+
* = 'L', the diagonal and first subdiagonal of A are over-
|
|
1351
|
+
* written by the corresponding elements of the tridiagonal
|
|
1352
|
+
* matrix T, and the elements below the first subdiagonal, with
|
|
1353
|
+
* the array TAU, represent the orthogonal matrix Q as a product
|
|
1354
|
+
* of elementary reflectors. See Further Details.
|
|
1355
|
+
*
|
|
1356
|
+
* D (output) REAL array, dimension (N)
|
|
1357
|
+
* The diagonal elements of the tridiagonal matrix T:
|
|
1358
|
+
* D(i) = A(i,i).
|
|
1359
|
+
*
|
|
1360
|
+
* E (output) REAL array, dimension (N-1)
|
|
1361
|
+
* The off-diagonal elements of the tridiagonal matrix T:
|
|
1362
|
+
* E(i) = A(i,i+1) if UPLO = 'U', E(i) = A(i+1,i) if UPLO = 'L'.
|
|
1363
|
+
*
|
|
1364
|
+
* TAU (output) REAL array, dimension (N-1)
|
|
1365
|
+
* The scalar factors of the elementary reflectors (see Further
|
|
1366
|
+
* Details).
|
|
1367
|
+
*
|
|
1368
|
+
* INFO (output) INTEGER
|
|
1369
|
+
* = 0: successful exit
|
|
1370
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1371
|
+
*
|
|
1372
|
+
|
|
1373
|
+
* Further Details
|
|
1374
|
+
* ===============
|
|
1375
|
+
*
|
|
1376
|
+
* If UPLO = 'U', the matrix Q is represented as a product of elementary
|
|
1377
|
+
* reflectors
|
|
1378
|
+
*
|
|
1379
|
+
* Q = H(n-1) . . . H(2) H(1).
|
|
1380
|
+
*
|
|
1381
|
+
* Each H(i) has the form
|
|
1382
|
+
*
|
|
1383
|
+
* H(i) = I - tau * v * v'
|
|
1384
|
+
*
|
|
1385
|
+
* where tau is a real scalar, and v is a real vector with
|
|
1386
|
+
* v(i+1:n) = 0 and v(i) = 1; v(1:i-1) is stored on exit in AP,
|
|
1387
|
+
* overwriting A(1:i-1,i+1), and tau is stored in TAU(i).
|
|
1388
|
+
*
|
|
1389
|
+
* If UPLO = 'L', the matrix Q is represented as a product of elementary
|
|
1390
|
+
* reflectors
|
|
1391
|
+
*
|
|
1392
|
+
* Q = H(1) H(2) . . . H(n-1).
|
|
1393
|
+
*
|
|
1394
|
+
* Each H(i) has the form
|
|
1395
|
+
*
|
|
1396
|
+
* H(i) = I - tau * v * v'
|
|
1397
|
+
*
|
|
1398
|
+
* where tau is a real scalar, and v is a real vector with
|
|
1399
|
+
* v(1:i) = 0 and v(i+1) = 1; v(i+2:n) is stored on exit in AP,
|
|
1400
|
+
* overwriting A(i+2:n,i), and tau is stored in TAU(i).
|
|
1401
|
+
*
|
|
1402
|
+
* =====================================================================
|
|
1403
|
+
*
|
|
1404
|
+
|
|
1405
|
+
|
|
1406
|
+
</PRE>
|
|
1407
|
+
<A HREF="#top">go to the page top</A>
|
|
1408
|
+
|
|
1409
|
+
<A NAME="ssptrf"></A>
|
|
1410
|
+
<H2>ssptrf</H2>
|
|
1411
|
+
<PRE>
|
|
1412
|
+
USAGE:
|
|
1413
|
+
ipiv, info, ap = NumRu::Lapack.ssptrf( uplo, ap, [:usage => usage, :help => help])
|
|
1414
|
+
|
|
1415
|
+
|
|
1416
|
+
FORTRAN MANUAL
|
|
1417
|
+
SUBROUTINE SSPTRF( UPLO, N, AP, IPIV, INFO )
|
|
1418
|
+
|
|
1419
|
+
* Purpose
|
|
1420
|
+
* =======
|
|
1421
|
+
*
|
|
1422
|
+
* SSPTRF computes the factorization of a real symmetric matrix A stored
|
|
1423
|
+
* in packed format using the Bunch-Kaufman diagonal pivoting method:
|
|
1424
|
+
*
|
|
1425
|
+
* A = U*D*U**T or A = L*D*L**T
|
|
1426
|
+
*
|
|
1427
|
+
* where U (or L) is a product of permutation and unit upper (lower)
|
|
1428
|
+
* triangular matrices, and D is symmetric and block diagonal with
|
|
1429
|
+
* 1-by-1 and 2-by-2 diagonal blocks.
|
|
1430
|
+
*
|
|
1431
|
+
|
|
1432
|
+
* Arguments
|
|
1433
|
+
* =========
|
|
1434
|
+
*
|
|
1435
|
+
* UPLO (input) CHARACTER*1
|
|
1436
|
+
* = 'U': Upper triangle of A is stored;
|
|
1437
|
+
* = 'L': Lower triangle of A is stored.
|
|
1438
|
+
*
|
|
1439
|
+
* N (input) INTEGER
|
|
1440
|
+
* The order of the matrix A. N >= 0.
|
|
1441
|
+
*
|
|
1442
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
1443
|
+
* On entry, the upper or lower triangle of the symmetric matrix
|
|
1444
|
+
* A, packed columnwise in a linear array. The j-th column of A
|
|
1445
|
+
* is stored in the array AP as follows:
|
|
1446
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
|
|
1447
|
+
* if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.
|
|
1448
|
+
*
|
|
1449
|
+
* On exit, the block diagonal matrix D and the multipliers used
|
|
1450
|
+
* to obtain the factor U or L, stored as a packed triangular
|
|
1451
|
+
* matrix overwriting A (see below for further details).
|
|
1452
|
+
*
|
|
1453
|
+
* IPIV (output) INTEGER array, dimension (N)
|
|
1454
|
+
* Details of the interchanges and the block structure of D.
|
|
1455
|
+
* If IPIV(k) > 0, then rows and columns k and IPIV(k) were
|
|
1456
|
+
* interchanged and D(k,k) is a 1-by-1 diagonal block.
|
|
1457
|
+
* If UPLO = 'U' and IPIV(k) = IPIV(k-1) < 0, then rows and
|
|
1458
|
+
* columns k-1 and -IPIV(k) were interchanged and D(k-1:k,k-1:k)
|
|
1459
|
+
* is a 2-by-2 diagonal block. If UPLO = 'L' and IPIV(k) =
|
|
1460
|
+
* IPIV(k+1) < 0, then rows and columns k+1 and -IPIV(k) were
|
|
1461
|
+
* interchanged and D(k:k+1,k:k+1) is a 2-by-2 diagonal block.
|
|
1462
|
+
*
|
|
1463
|
+
* INFO (output) INTEGER
|
|
1464
|
+
* = 0: successful exit
|
|
1465
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1466
|
+
* > 0: if INFO = i, D(i,i) is exactly zero. The factorization
|
|
1467
|
+
* has been completed, but the block diagonal matrix D is
|
|
1468
|
+
* exactly singular, and division by zero will occur if it
|
|
1469
|
+
* is used to solve a system of equations.
|
|
1470
|
+
*
|
|
1471
|
+
|
|
1472
|
+
* Further Details
|
|
1473
|
+
* ===============
|
|
1474
|
+
*
|
|
1475
|
+
* 5-96 - Based on modifications by J. Lewis, Boeing Computer Services
|
|
1476
|
+
* Company
|
|
1477
|
+
*
|
|
1478
|
+
* If UPLO = 'U', then A = U*D*U', where
|
|
1479
|
+
* U = P(n)*U(n)* ... *P(k)U(k)* ...,
|
|
1480
|
+
* i.e., U is a product of terms P(k)*U(k), where k decreases from n to
|
|
1481
|
+
* 1 in steps of 1 or 2, and D is a block diagonal matrix with 1-by-1
|
|
1482
|
+
* and 2-by-2 diagonal blocks D(k). P(k) is a permutation matrix as
|
|
1483
|
+
* defined by IPIV(k), and U(k) is a unit upper triangular matrix, such
|
|
1484
|
+
* that if the diagonal block D(k) is of order s (s = 1 or 2), then
|
|
1485
|
+
*
|
|
1486
|
+
* ( I v 0 ) k-s
|
|
1487
|
+
* U(k) = ( 0 I 0 ) s
|
|
1488
|
+
* ( 0 0 I ) n-k
|
|
1489
|
+
* k-s s n-k
|
|
1490
|
+
*
|
|
1491
|
+
* If s = 1, D(k) overwrites A(k,k), and v overwrites A(1:k-1,k).
|
|
1492
|
+
* If s = 2, the upper triangle of D(k) overwrites A(k-1,k-1), A(k-1,k),
|
|
1493
|
+
* and A(k,k), and v overwrites A(1:k-2,k-1:k).
|
|
1494
|
+
*
|
|
1495
|
+
* If UPLO = 'L', then A = L*D*L', where
|
|
1496
|
+
* L = P(1)*L(1)* ... *P(k)*L(k)* ...,
|
|
1497
|
+
* i.e., L is a product of terms P(k)*L(k), where k increases from 1 to
|
|
1498
|
+
* n in steps of 1 or 2, and D is a block diagonal matrix with 1-by-1
|
|
1499
|
+
* and 2-by-2 diagonal blocks D(k). P(k) is a permutation matrix as
|
|
1500
|
+
* defined by IPIV(k), and L(k) is a unit lower triangular matrix, such
|
|
1501
|
+
* that if the diagonal block D(k) is of order s (s = 1 or 2), then
|
|
1502
|
+
*
|
|
1503
|
+
* ( I 0 0 ) k-1
|
|
1504
|
+
* L(k) = ( 0 I 0 ) s
|
|
1505
|
+
* ( 0 v I ) n-k-s+1
|
|
1506
|
+
* k-1 s n-k-s+1
|
|
1507
|
+
*
|
|
1508
|
+
* If s = 1, D(k) overwrites A(k,k), and v overwrites A(k+1:n,k).
|
|
1509
|
+
* If s = 2, the lower triangle of D(k) overwrites A(k,k), A(k+1,k),
|
|
1510
|
+
* and A(k+1,k+1), and v overwrites A(k+2:n,k:k+1).
|
|
1511
|
+
*
|
|
1512
|
+
* =====================================================================
|
|
1513
|
+
*
|
|
1514
|
+
|
|
1515
|
+
|
|
1516
|
+
</PRE>
|
|
1517
|
+
<A HREF="#top">go to the page top</A>
|
|
1518
|
+
|
|
1519
|
+
<A NAME="ssptri"></A>
|
|
1520
|
+
<H2>ssptri</H2>
|
|
1521
|
+
<PRE>
|
|
1522
|
+
USAGE:
|
|
1523
|
+
info, ap = NumRu::Lapack.ssptri( uplo, ap, ipiv, [:usage => usage, :help => help])
|
|
1524
|
+
|
|
1525
|
+
|
|
1526
|
+
FORTRAN MANUAL
|
|
1527
|
+
SUBROUTINE SSPTRI( UPLO, N, AP, IPIV, WORK, INFO )
|
|
1528
|
+
|
|
1529
|
+
* Purpose
|
|
1530
|
+
* =======
|
|
1531
|
+
*
|
|
1532
|
+
* SSPTRI computes the inverse of a real symmetric indefinite matrix
|
|
1533
|
+
* A in packed storage using the factorization A = U*D*U**T or
|
|
1534
|
+
* A = L*D*L**T computed by SSPTRF.
|
|
1535
|
+
*
|
|
1536
|
+
|
|
1537
|
+
* Arguments
|
|
1538
|
+
* =========
|
|
1539
|
+
*
|
|
1540
|
+
* UPLO (input) CHARACTER*1
|
|
1541
|
+
* Specifies whether the details of the factorization are stored
|
|
1542
|
+
* as an upper or lower triangular matrix.
|
|
1543
|
+
* = 'U': Upper triangular, form is A = U*D*U**T;
|
|
1544
|
+
* = 'L': Lower triangular, form is A = L*D*L**T.
|
|
1545
|
+
*
|
|
1546
|
+
* N (input) INTEGER
|
|
1547
|
+
* The order of the matrix A. N >= 0.
|
|
1548
|
+
*
|
|
1549
|
+
* AP (input/output) REAL array, dimension (N*(N+1)/2)
|
|
1550
|
+
* On entry, the block diagonal matrix D and the multipliers
|
|
1551
|
+
* used to obtain the factor U or L as computed by SSPTRF,
|
|
1552
|
+
* stored as a packed triangular matrix.
|
|
1553
|
+
*
|
|
1554
|
+
* On exit, if INFO = 0, the (symmetric) inverse of the original
|
|
1555
|
+
* matrix, stored as a packed triangular matrix. The j-th column
|
|
1556
|
+
* of inv(A) is stored in the array AP as follows:
|
|
1557
|
+
* if UPLO = 'U', AP(i + (j-1)*j/2) = inv(A)(i,j) for 1<=i<=j;
|
|
1558
|
+
* if UPLO = 'L',
|
|
1559
|
+
* AP(i + (j-1)*(2n-j)/2) = inv(A)(i,j) for j<=i<=n.
|
|
1560
|
+
*
|
|
1561
|
+
* IPIV (input) INTEGER array, dimension (N)
|
|
1562
|
+
* Details of the interchanges and the block structure of D
|
|
1563
|
+
* as determined by SSPTRF.
|
|
1564
|
+
*
|
|
1565
|
+
* WORK (workspace) REAL array, dimension (N)
|
|
1566
|
+
*
|
|
1567
|
+
* INFO (output) INTEGER
|
|
1568
|
+
* = 0: successful exit
|
|
1569
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1570
|
+
* > 0: if INFO = i, D(i,i) = 0; the matrix is singular and its
|
|
1571
|
+
* inverse could not be computed.
|
|
1572
|
+
*
|
|
1573
|
+
|
|
1574
|
+
* =====================================================================
|
|
1575
|
+
*
|
|
1576
|
+
|
|
1577
|
+
|
|
1578
|
+
</PRE>
|
|
1579
|
+
<A HREF="#top">go to the page top</A>
|
|
1580
|
+
|
|
1581
|
+
<A NAME="ssptrs"></A>
|
|
1582
|
+
<H2>ssptrs</H2>
|
|
1583
|
+
<PRE>
|
|
1584
|
+
USAGE:
|
|
1585
|
+
info, b = NumRu::Lapack.ssptrs( uplo, ap, ipiv, b, [:usage => usage, :help => help])
|
|
1586
|
+
|
|
1587
|
+
|
|
1588
|
+
FORTRAN MANUAL
|
|
1589
|
+
SUBROUTINE SSPTRS( UPLO, N, NRHS, AP, IPIV, B, LDB, INFO )
|
|
1590
|
+
|
|
1591
|
+
* Purpose
|
|
1592
|
+
* =======
|
|
1593
|
+
*
|
|
1594
|
+
* SSPTRS solves a system of linear equations A*X = B with a real
|
|
1595
|
+
* symmetric matrix A stored in packed format using the factorization
|
|
1596
|
+
* A = U*D*U**T or A = L*D*L**T computed by SSPTRF.
|
|
1597
|
+
*
|
|
1598
|
+
|
|
1599
|
+
* Arguments
|
|
1600
|
+
* =========
|
|
1601
|
+
*
|
|
1602
|
+
* UPLO (input) CHARACTER*1
|
|
1603
|
+
* Specifies whether the details of the factorization are stored
|
|
1604
|
+
* as an upper or lower triangular matrix.
|
|
1605
|
+
* = 'U': Upper triangular, form is A = U*D*U**T;
|
|
1606
|
+
* = 'L': Lower triangular, form is A = L*D*L**T.
|
|
1607
|
+
*
|
|
1608
|
+
* N (input) INTEGER
|
|
1609
|
+
* The order of the matrix A. N >= 0.
|
|
1610
|
+
*
|
|
1611
|
+
* NRHS (input) INTEGER
|
|
1612
|
+
* The number of right hand sides, i.e., the number of columns
|
|
1613
|
+
* of the matrix B. NRHS >= 0.
|
|
1614
|
+
*
|
|
1615
|
+
* AP (input) REAL array, dimension (N*(N+1)/2)
|
|
1616
|
+
* The block diagonal matrix D and the multipliers used to
|
|
1617
|
+
* obtain the factor U or L as computed by SSPTRF, stored as a
|
|
1618
|
+
* packed triangular matrix.
|
|
1619
|
+
*
|
|
1620
|
+
* IPIV (input) INTEGER array, dimension (N)
|
|
1621
|
+
* Details of the interchanges and the block structure of D
|
|
1622
|
+
* as determined by SSPTRF.
|
|
1623
|
+
*
|
|
1624
|
+
* B (input/output) REAL array, dimension (LDB,NRHS)
|
|
1625
|
+
* On entry, the right hand side matrix B.
|
|
1626
|
+
* On exit, the solution matrix X.
|
|
1627
|
+
*
|
|
1628
|
+
* LDB (input) INTEGER
|
|
1629
|
+
* The leading dimension of the array B. LDB >= max(1,N).
|
|
1630
|
+
*
|
|
1631
|
+
* INFO (output) INTEGER
|
|
1632
|
+
* = 0: successful exit
|
|
1633
|
+
* < 0: if INFO = -i, the i-th argument had an illegal value
|
|
1634
|
+
*
|
|
1635
|
+
|
|
1636
|
+
* =====================================================================
|
|
1637
|
+
*
|
|
1638
|
+
|
|
1639
|
+
|
|
1640
|
+
</PRE>
|
|
1641
|
+
<A HREF="#top">go to the page top</A>
|
|
1642
|
+
|
|
1643
|
+
<HR />
|
|
1644
|
+
<A HREF="s.html">back to matrix types</A><BR>
|
|
1645
|
+
<A HREF="s.html">back to data types</A>
|
|
1646
|
+
</BODY>
|
|
1647
|
+
</HTML>
|