ruby-lapack 1.6 → 1.7

Sign up to get free protection for your applications and to get access to all the features.
Files changed (1637) hide show
  1. checksums.yaml +7 -0
  2. data/Rakefile +1 -1
  3. data/dev/make_csrc.rb +8 -2
  4. data/ext/extconf.rb +15 -11
  5. data/ext/rb_lapack.h +7 -1
  6. metadata +48 -1700
  7. data/ext/cbbcsd.c +0 -283
  8. data/ext/cbdsqr.c +0 -182
  9. data/ext/cgbbrd.c +0 -157
  10. data/ext/cgbcon.c +0 -98
  11. data/ext/cgbequ.c +0 -98
  12. data/ext/cgbequb.c +0 -96
  13. data/ext/cgbrfs.c +0 -161
  14. data/ext/cgbrfsx.c +0 -249
  15. data/ext/cgbsv.c +0 -115
  16. data/ext/cgbsvx.c +0 -286
  17. data/ext/cgbsvxx.c +0 -289
  18. data/ext/cgbtf2.c +0 -93
  19. data/ext/cgbtrf.c +0 -93
  20. data/ext/cgbtrs.c +0 -111
  21. data/ext/cgebak.c +0 -101
  22. data/ext/cgebal.c +0 -91
  23. data/ext/cgebd2.c +0 -112
  24. data/ext/cgebrd.c +0 -127
  25. data/ext/cgecon.c +0 -78
  26. data/ext/cgeequ.c +0 -88
  27. data/ext/cgeequb.c +0 -88
  28. data/ext/cgees.c +0 -142
  29. data/ext/cgeesx.c +0 -152
  30. data/ext/cgeev.c +0 -132
  31. data/ext/cgeevx.c +0 -173
  32. data/ext/cgegs.c +0 -166
  33. data/ext/cgegv.c +0 -171
  34. data/ext/cgehd2.c +0 -92
  35. data/ext/cgehrd.c +0 -107
  36. data/ext/cgelq2.c +0 -86
  37. data/ext/cgelqf.c +0 -103
  38. data/ext/cgels.c +0 -137
  39. data/ext/cgelsd.c +0 -154
  40. data/ext/cgelss.c +0 -151
  41. data/ext/cgelsx.c +0 -139
  42. data/ext/cgelsy.c +0 -166
  43. data/ext/cgeql2.c +0 -88
  44. data/ext/cgeqlf.c +0 -103
  45. data/ext/cgeqp3.c +0 -129
  46. data/ext/cgeqpf.c +0 -114
  47. data/ext/cgeqr2.c +0 -88
  48. data/ext/cgeqr2p.c +0 -88
  49. data/ext/cgeqrf.c +0 -103
  50. data/ext/cgeqrfp.c +0 -103
  51. data/ext/cgerfs.c +0 -153
  52. data/ext/cgerfsx.c +0 -219
  53. data/ext/cgerq2.c +0 -86
  54. data/ext/cgerqf.c +0 -103
  55. data/ext/cgesc2.c +0 -108
  56. data/ext/cgesdd.c +0 -135
  57. data/ext/cgesv.c +0 -107
  58. data/ext/cgesvd.c +0 -146
  59. data/ext/cgesvx.c +0 -278
  60. data/ext/cgesvxx.c +0 -281
  61. data/ext/cgetc2.c +0 -89
  62. data/ext/cgetf2.c +0 -85
  63. data/ext/cgetrf.c +0 -85
  64. data/ext/cgetri.c +0 -103
  65. data/ext/cgetrs.c +0 -103
  66. data/ext/cggbak.c +0 -113
  67. data/ext/cggbal.c +0 -128
  68. data/ext/cgges.c +0 -192
  69. data/ext/cggesx.c +0 -230
  70. data/ext/cggev.c +0 -171
  71. data/ext/cggevx.c +0 -226
  72. data/ext/cggglm.c +0 -156
  73. data/ext/cgghrd.c +0 -167
  74. data/ext/cgglse.c +0 -171
  75. data/ext/cggqrf.c +0 -137
  76. data/ext/cggrqf.c +0 -141
  77. data/ext/cggsvd.c +0 -184
  78. data/ext/cggsvp.c +0 -174
  79. data/ext/cgtcon.c +0 -121
  80. data/ext/cgtrfs.c +0 -209
  81. data/ext/cgtsv.c +0 -142
  82. data/ext/cgtsvx.c +0 -256
  83. data/ext/cgttrf.c +0 -132
  84. data/ext/cgttrs.c +0 -137
  85. data/ext/cgtts2.c +0 -134
  86. data/ext/chbev.c +0 -110
  87. data/ext/chbevd.c +0 -158
  88. data/ext/chbevx.c +0 -160
  89. data/ext/chbgst.c +0 -120
  90. data/ext/chbgv.c +0 -140
  91. data/ext/chbgvd.c +0 -188
  92. data/ext/chbgvx.c +0 -189
  93. data/ext/chbtrd.c +0 -130
  94. data/ext/checon.c +0 -87
  95. data/ext/cheequb.c +0 -82
  96. data/ext/cheev.c +0 -110
  97. data/ext/cheevd.c +0 -143
  98. data/ext/cheevr.c +0 -190
  99. data/ext/cheevx.c +0 -160
  100. data/ext/chegs2.c +0 -95
  101. data/ext/chegst.c +0 -95
  102. data/ext/chegv.c +0 -140
  103. data/ext/chegvd.c +0 -173
  104. data/ext/chegvx.c +0 -190
  105. data/ext/cherfs.c +0 -153
  106. data/ext/cherfsx.c +0 -218
  107. data/ext/chesv.c +0 -123
  108. data/ext/chesvx.c +0 -183
  109. data/ext/chesvxx.c +0 -258
  110. data/ext/chetd2.c +0 -101
  111. data/ext/chetf2.c +0 -85
  112. data/ext/chetrd.c +0 -113
  113. data/ext/chetrf.c +0 -97
  114. data/ext/chetri.c +0 -92
  115. data/ext/chetrs.c +0 -103
  116. data/ext/chetrs2.c +0 -106
  117. data/ext/chfrk.c +0 -109
  118. data/ext/chgeqz.c +0 -208
  119. data/ext/chla_transtype.c +0 -51
  120. data/ext/chpcon.c +0 -85
  121. data/ext/chpev.c +0 -105
  122. data/ext/chpevd.c +0 -153
  123. data/ext/chpevx.c +0 -144
  124. data/ext/chpgst.c +0 -94
  125. data/ext/chpgv.c +0 -132
  126. data/ext/chpgvd.c +0 -170
  127. data/ext/chpgvx.c +0 -170
  128. data/ext/chprfs.c +0 -149
  129. data/ext/chpsv.c +0 -110
  130. data/ext/chpsvx.c +0 -163
  131. data/ext/chptrd.c +0 -100
  132. data/ext/chptrf.c +0 -84
  133. data/ext/chptri.c +0 -89
  134. data/ext/chptrs.c +0 -101
  135. data/ext/chsein.c +0 -185
  136. data/ext/chseqr.c +0 -145
  137. data/ext/cla_gbamv.c +0 -127
  138. data/ext/cla_gbrcond_c.c +0 -142
  139. data/ext/cla_gbrcond_x.c +0 -138
  140. data/ext/cla_gbrfsx_extended.c +0 -295
  141. data/ext/cla_gbrpvgrw.c +0 -87
  142. data/ext/cla_geamv.c +0 -117
  143. data/ext/cla_gercond_c.c +0 -134
  144. data/ext/cla_gercond_x.c +0 -130
  145. data/ext/cla_gerfsx_extended.c +0 -281
  146. data/ext/cla_heamv.c +0 -116
  147. data/ext/cla_hercond_c.c +0 -134
  148. data/ext/cla_hercond_x.c +0 -130
  149. data/ext/cla_herfsx_extended.c +0 -283
  150. data/ext/cla_herpvgrw.c +0 -107
  151. data/ext/cla_lin_berr.c +0 -84
  152. data/ext/cla_porcond_c.c +0 -122
  153. data/ext/cla_porcond_x.c +0 -118
  154. data/ext/cla_porfsx_extended.c +0 -271
  155. data/ext/cla_porpvgrw.c +0 -95
  156. data/ext/cla_rpvgrw.c +0 -79
  157. data/ext/cla_syamv.c +0 -115
  158. data/ext/cla_syrcond_c.c +0 -134
  159. data/ext/cla_syrcond_x.c +0 -130
  160. data/ext/cla_syrfsx_extended.c +0 -283
  161. data/ext/cla_syrpvgrw.c +0 -107
  162. data/ext/cla_wwaddw.c +0 -102
  163. data/ext/clabrd.c +0 -132
  164. data/ext/clacgv.c +0 -75
  165. data/ext/clacn2.c +0 -103
  166. data/ext/clacon.c +0 -80
  167. data/ext/clacp2.c +0 -77
  168. data/ext/clacpy.c +0 -77
  169. data/ext/clacrm.c +0 -90
  170. data/ext/clacrt.c +0 -108
  171. data/ext/cladiv.c +0 -57
  172. data/ext/claed0.c +0 -134
  173. data/ext/claed7.c +0 -247
  174. data/ext/claed8.c +0 -198
  175. data/ext/claein.c +0 -113
  176. data/ext/claesy.c +0 -74
  177. data/ext/claev2.c +0 -71
  178. data/ext/clag2z.c +0 -76
  179. data/ext/clags2.c +0 -92
  180. data/ext/clagtm.c +0 -132
  181. data/ext/clahef.c +0 -97
  182. data/ext/clahqr.c +0 -135
  183. data/ext/clahr2.c +0 -112
  184. data/ext/clahrd.c +0 -112
  185. data/ext/claic1.c +0 -90
  186. data/ext/clals0.c +0 -201
  187. data/ext/clalsa.c +0 -270
  188. data/ext/clalsd.c +0 -145
  189. data/ext/clangb.c +0 -76
  190. data/ext/clange.c +0 -74
  191. data/ext/clangt.c +0 -87
  192. data/ext/clanhb.c +0 -78
  193. data/ext/clanhe.c +0 -72
  194. data/ext/clanhf.c +0 -80
  195. data/ext/clanhp.c +0 -74
  196. data/ext/clanhs.c +0 -70
  197. data/ext/clanht.c +0 -75
  198. data/ext/clansb.c +0 -78
  199. data/ext/clansp.c +0 -74
  200. data/ext/clansy.c +0 -72
  201. data/ext/clantb.c +0 -80
  202. data/ext/clantp.c +0 -80
  203. data/ext/clantr.c +0 -82
  204. data/ext/clapll.c +0 -105
  205. data/ext/clapmr.c +0 -97
  206. data/ext/clapmt.c +0 -101
  207. data/ext/claqgb.c +0 -117
  208. data/ext/claqge.c +0 -109
  209. data/ext/claqhb.c +0 -97
  210. data/ext/claqhe.c +0 -97
  211. data/ext/claqhp.c +0 -94
  212. data/ext/claqp2.c +0 -158
  213. data/ext/claqps.c +0 -208
  214. data/ext/claqr0.c +0 -145
  215. data/ext/claqr1.c +0 -76
  216. data/ext/claqr2.c +0 -174
  217. data/ext/claqr3.c +0 -174
  218. data/ext/claqr4.c +0 -145
  219. data/ext/claqr5.c +0 -179
  220. data/ext/claqsb.c +0 -101
  221. data/ext/claqsp.c +0 -94
  222. data/ext/claqsy.c +0 -97
  223. data/ext/clar1v.c +0 -173
  224. data/ext/clar2v.c +0 -149
  225. data/ext/clarcm.c +0 -86
  226. data/ext/clarf.c +0 -102
  227. data/ext/clarfb.c +0 -123
  228. data/ext/clarfg.c +0 -84
  229. data/ext/clarfgp.c +0 -84
  230. data/ext/clarft.c +0 -105
  231. data/ext/clarfx.c +0 -94
  232. data/ext/clargv.c +0 -114
  233. data/ext/clarnv.c +0 -83
  234. data/ext/clarrv.c +0 -271
  235. data/ext/clarscl2.c +0 -82
  236. data/ext/clartg.c +0 -63
  237. data/ext/clartv.c +0 -130
  238. data/ext/clarz.c +0 -106
  239. data/ext/clarzb.c +0 -127
  240. data/ext/clarzt.c +0 -105
  241. data/ext/clascl.c +0 -97
  242. data/ext/clascl2.c +0 -82
  243. data/ext/claset.c +0 -88
  244. data/ext/clasr.c +0 -110
  245. data/ext/classq.c +0 -70
  246. data/ext/claswp.c +0 -94
  247. data/ext/clasyf.c +0 -97
  248. data/ext/clatbs.c +0 -130
  249. data/ext/clatdf.c +0 -119
  250. data/ext/clatps.c +0 -124
  251. data/ext/clatrd.c +0 -105
  252. data/ext/clatrs.c +0 -126
  253. data/ext/clatrz.c +0 -87
  254. data/ext/clatzm.c +0 -132
  255. data/ext/clauu2.c +0 -77
  256. data/ext/clauum.c +0 -77
  257. data/ext/cpbcon.c +0 -82
  258. data/ext/cpbequ.c +0 -83
  259. data/ext/cpbrfs.c +0 -145
  260. data/ext/cpbstf.c +0 -81
  261. data/ext/cpbsv.c +0 -107
  262. data/ext/cpbsvx.c +0 -201
  263. data/ext/cpbtf2.c +0 -81
  264. data/ext/cpbtrf.c +0 -81
  265. data/ext/cpbtrs.c +0 -95
  266. data/ext/cpftrf.c +0 -82
  267. data/ext/cpftri.c +0 -82
  268. data/ext/cpftrs.c +0 -97
  269. data/ext/cpocon.c +0 -78
  270. data/ext/cpoequ.c +0 -75
  271. data/ext/cpoequb.c +0 -75
  272. data/ext/cporfs.c +0 -141
  273. data/ext/cporfsx.c +0 -206
  274. data/ext/cposv.c +0 -103
  275. data/ext/cposvx.c +0 -197
  276. data/ext/cposvxx.c +0 -235
  277. data/ext/cpotf2.c +0 -77
  278. data/ext/cpotrf.c +0 -77
  279. data/ext/cpotri.c +0 -77
  280. data/ext/cpotrs.c +0 -91
  281. data/ext/cppcon.c +0 -78
  282. data/ext/cppequ.c +0 -79
  283. data/ext/cpprfs.c +0 -139
  284. data/ext/cppsv.c +0 -104
  285. data/ext/cppsvx.c +0 -191
  286. data/ext/cpptrf.c +0 -78
  287. data/ext/cpptri.c +0 -78
  288. data/ext/cpptrs.c +0 -93
  289. data/ext/cpstf2.c +0 -95
  290. data/ext/cpstrf.c +0 -95
  291. data/ext/cptcon.c +0 -81
  292. data/ext/cpteqr.c +0 -126
  293. data/ext/cptrfs.c +0 -161
  294. data/ext/cptsv.c +0 -119
  295. data/ext/cptsvx.c +0 -171
  296. data/ext/cpttrf.c +0 -93
  297. data/ext/cpttrs.c +0 -101
  298. data/ext/cptts2.c +0 -98
  299. data/ext/crot.c +0 -107
  300. data/ext/cspcon.c +0 -85
  301. data/ext/cspmv.c +0 -115
  302. data/ext/cspr.c +0 -96
  303. data/ext/csprfs.c +0 -149
  304. data/ext/cspsv.c +0 -110
  305. data/ext/cspsvx.c +0 -163
  306. data/ext/csptrf.c +0 -84
  307. data/ext/csptri.c +0 -89
  308. data/ext/csptrs.c +0 -101
  309. data/ext/csrscl.c +0 -79
  310. data/ext/cstedc.c +0 -177
  311. data/ext/cstegr.c +0 -188
  312. data/ext/cstein.c +0 -134
  313. data/ext/cstemr.c +0 -193
  314. data/ext/csteqr.c +0 -126
  315. data/ext/csycon.c +0 -87
  316. data/ext/csyconv.c +0 -84
  317. data/ext/csyequb.c +0 -82
  318. data/ext/csymv.c +0 -115
  319. data/ext/csyr.c +0 -95
  320. data/ext/csyrfs.c +0 -153
  321. data/ext/csyrfsx.c +0 -218
  322. data/ext/csysv.c +0 -129
  323. data/ext/csysvx.c +0 -183
  324. data/ext/csysvxx.c +0 -258
  325. data/ext/csyswapr.c +0 -82
  326. data/ext/csytf2.c +0 -85
  327. data/ext/csytrf.c +0 -97
  328. data/ext/csytri.c +0 -92
  329. data/ext/csytri2.c +0 -108
  330. data/ext/csytri2x.c +0 -96
  331. data/ext/csytrs.c +0 -103
  332. data/ext/csytrs2.c +0 -106
  333. data/ext/ctbcon.c +0 -86
  334. data/ext/ctbrfs.c +0 -127
  335. data/ext/ctbtrs.c +0 -103
  336. data/ext/ctfsm.c +0 -111
  337. data/ext/ctftri.c +0 -86
  338. data/ext/ctfttp.c +0 -79
  339. data/ext/ctfttr.c +0 -80
  340. data/ext/ctgevc.c +0 -156
  341. data/ext/ctgex2.c +0 -171
  342. data/ext/ctgexc.c +0 -172
  343. data/ext/ctgsen.c +0 -244
  344. data/ext/ctgsja.c +0 -227
  345. data/ext/ctgsna.c +0 -164
  346. data/ext/ctgsy2.c +0 -176
  347. data/ext/ctgsyl.c +0 -190
  348. data/ext/ctpcon.c +0 -82
  349. data/ext/ctprfs.c +0 -123
  350. data/ext/ctptri.c +0 -82
  351. data/ext/ctptrs.c +0 -101
  352. data/ext/ctpttf.c +0 -79
  353. data/ext/ctpttr.c +0 -76
  354. data/ext/ctrcon.c +0 -82
  355. data/ext/ctrevc.c +0 -154
  356. data/ext/ctrexc.c +0 -111
  357. data/ext/ctrrfs.c +0 -123
  358. data/ext/ctrsen.c +0 -154
  359. data/ext/ctrsna.c +0 -137
  360. data/ext/ctrsyl.c +0 -116
  361. data/ext/ctrti2.c +0 -81
  362. data/ext/ctrtri.c +0 -81
  363. data/ext/ctrtrs.c +0 -99
  364. data/ext/ctrttf.c +0 -77
  365. data/ext/ctrttp.c +0 -73
  366. data/ext/ctzrqf.c +0 -83
  367. data/ext/ctzrzf.c +0 -101
  368. data/ext/cunbdb.c +0 -232
  369. data/ext/cuncsd.c +0 -204
  370. data/ext/cung2l.c +0 -92
  371. data/ext/cung2r.c +0 -92
  372. data/ext/cungbr.c +0 -115
  373. data/ext/cunghr.c +0 -111
  374. data/ext/cungl2.c +0 -90
  375. data/ext/cunglq.c +0 -107
  376. data/ext/cungql.c +0 -107
  377. data/ext/cungqr.c +0 -107
  378. data/ext/cungr2.c +0 -90
  379. data/ext/cungrq.c +0 -107
  380. data/ext/cungtr.c +0 -107
  381. data/ext/cunm2l.c +0 -114
  382. data/ext/cunm2r.c +0 -114
  383. data/ext/cunmbr.c +0 -139
  384. data/ext/cunmhr.c +0 -133
  385. data/ext/cunml2.c +0 -110
  386. data/ext/cunmlq.c +0 -125
  387. data/ext/cunmql.c +0 -129
  388. data/ext/cunmqr.c +0 -129
  389. data/ext/cunmr2.c +0 -110
  390. data/ext/cunmr3.c +0 -114
  391. data/ext/cunmrq.c +0 -125
  392. data/ext/cunmrz.c +0 -129
  393. data/ext/cunmtr.c +0 -129
  394. data/ext/cupgtr.c +0 -91
  395. data/ext/cupmtr.c +0 -116
  396. data/ext/dbbcsd.c +0 -287
  397. data/ext/dbdsdc.c +0 -151
  398. data/ext/dbdsqr.c +0 -182
  399. data/ext/ddisna.c +0 -75
  400. data/ext/dgbbrd.c +0 -154
  401. data/ext/dgbcon.c +0 -98
  402. data/ext/dgbequ.c +0 -98
  403. data/ext/dgbequb.c +0 -96
  404. data/ext/dgbrfs.c +0 -161
  405. data/ext/dgbrfsx.c +0 -249
  406. data/ext/dgbsv.c +0 -115
  407. data/ext/dgbsvx.c +0 -286
  408. data/ext/dgbsvxx.c +0 -289
  409. data/ext/dgbtf2.c +0 -93
  410. data/ext/dgbtrf.c +0 -93
  411. data/ext/dgbtrs.c +0 -111
  412. data/ext/dgebak.c +0 -101
  413. data/ext/dgebal.c +0 -91
  414. data/ext/dgebd2.c +0 -112
  415. data/ext/dgebrd.c +0 -127
  416. data/ext/dgecon.c +0 -78
  417. data/ext/dgeequ.c +0 -88
  418. data/ext/dgeequb.c +0 -88
  419. data/ext/dgees.c +0 -148
  420. data/ext/dgeesx.c +0 -170
  421. data/ext/dgeev.c +0 -137
  422. data/ext/dgeevx.c +0 -181
  423. data/ext/dgegs.c +0 -171
  424. data/ext/dgegv.c +0 -171
  425. data/ext/dgehd2.c +0 -92
  426. data/ext/dgehrd.c +0 -107
  427. data/ext/dgejsv.c +0 -159
  428. data/ext/dgelq2.c +0 -86
  429. data/ext/dgelqf.c +0 -103
  430. data/ext/dgels.c +0 -137
  431. data/ext/dgelsd.c +0 -149
  432. data/ext/dgelss.c +0 -148
  433. data/ext/dgelsx.c +0 -136
  434. data/ext/dgelsy.c +0 -163
  435. data/ext/dgeql2.c +0 -88
  436. data/ext/dgeqlf.c +0 -103
  437. data/ext/dgeqp3.c +0 -126
  438. data/ext/dgeqpf.c +0 -111
  439. data/ext/dgeqr2.c +0 -88
  440. data/ext/dgeqr2p.c +0 -88
  441. data/ext/dgeqrf.c +0 -103
  442. data/ext/dgeqrfp.c +0 -103
  443. data/ext/dgerfs.c +0 -153
  444. data/ext/dgerfsx.c +0 -219
  445. data/ext/dgerq2.c +0 -86
  446. data/ext/dgerqf.c +0 -103
  447. data/ext/dgesc2.c +0 -108
  448. data/ext/dgesdd.c +0 -132
  449. data/ext/dgesv.c +0 -107
  450. data/ext/dgesvd.c +0 -143
  451. data/ext/dgesvj.c +0 -156
  452. data/ext/dgesvx.c +0 -278
  453. data/ext/dgesvxx.c +0 -281
  454. data/ext/dgetc2.c +0 -89
  455. data/ext/dgetf2.c +0 -85
  456. data/ext/dgetrf.c +0 -85
  457. data/ext/dgetri.c +0 -103
  458. data/ext/dgetrs.c +0 -103
  459. data/ext/dggbak.c +0 -113
  460. data/ext/dggbal.c +0 -128
  461. data/ext/dgges.c +0 -198
  462. data/ext/dggesx.c +0 -231
  463. data/ext/dggev.c +0 -171
  464. data/ext/dggevx.c +0 -229
  465. data/ext/dggglm.c +0 -156
  466. data/ext/dgghrd.c +0 -167
  467. data/ext/dgglse.c +0 -171
  468. data/ext/dggqrf.c +0 -137
  469. data/ext/dggrqf.c +0 -141
  470. data/ext/dggsvd.c +0 -181
  471. data/ext/dggsvp.c +0 -171
  472. data/ext/dgsvj0.c +0 -182
  473. data/ext/dgsvj1.c +0 -186
  474. data/ext/dgtcon.c +0 -124
  475. data/ext/dgtrfs.c +0 -209
  476. data/ext/dgtsv.c +0 -142
  477. data/ext/dgtsvx.c +0 -256
  478. data/ext/dgttrf.c +0 -132
  479. data/ext/dgttrs.c +0 -137
  480. data/ext/dgtts2.c +0 -134
  481. data/ext/dhgeqz.c +0 -213
  482. data/ext/dhsein.c +0 -205
  483. data/ext/dhseqr.c +0 -153
  484. data/ext/disnan.c +0 -51
  485. data/ext/dla_gbamv.c +0 -129
  486. data/ext/dla_gbrcond.c +0 -142
  487. data/ext/dla_gbrfsx_extended.c +0 -293
  488. data/ext/dla_gbrpvgrw.c +0 -87
  489. data/ext/dla_geamv.c +0 -119
  490. data/ext/dla_gercond.c +0 -134
  491. data/ext/dla_gerfsx_extended.c +0 -281
  492. data/ext/dla_lin_berr.c +0 -84
  493. data/ext/dla_porcond.c +0 -122
  494. data/ext/dla_porfsx_extended.c +0 -271
  495. data/ext/dla_porpvgrw.c +0 -95
  496. data/ext/dla_rpvgrw.c +0 -79
  497. data/ext/dla_syamv.c +0 -113
  498. data/ext/dla_syrcond.c +0 -134
  499. data/ext/dla_syrfsx_extended.c +0 -283
  500. data/ext/dla_syrpvgrw.c +0 -107
  501. data/ext/dla_wwaddw.c +0 -102
  502. data/ext/dlabad.c +0 -54
  503. data/ext/dlabrd.c +0 -132
  504. data/ext/dlacn2.c +0 -106
  505. data/ext/dlacon.c +0 -83
  506. data/ext/dlacpy.c +0 -77
  507. data/ext/dladiv.c +0 -66
  508. data/ext/dlae2.c +0 -62
  509. data/ext/dlaebz.c +0 -218
  510. data/ext/dlaed0.c +0 -127
  511. data/ext/dlaed1.c +0 -133
  512. data/ext/dlaed2.c +0 -189
  513. data/ext/dlaed3.c +0 -161
  514. data/ext/dlaed4.c +0 -90
  515. data/ext/dlaed5.c +0 -87
  516. data/ext/dlaed6.c +0 -90
  517. data/ext/dlaed7.c +0 -248
  518. data/ext/dlaed8.c +0 -206
  519. data/ext/dlaed9.c +0 -111
  520. data/ext/dlaeda.c +0 -160
  521. data/ext/dlaein.c +0 -143
  522. data/ext/dlaev2.c +0 -68
  523. data/ext/dlaexc.c +0 -118
  524. data/ext/dlag2.c +0 -91
  525. data/ext/dlag2s.c +0 -76
  526. data/ext/dlags2.c +0 -90
  527. data/ext/dlagtf.c +0 -140
  528. data/ext/dlagtm.c +0 -132
  529. data/ext/dlagts.c +0 -139
  530. data/ext/dlagv2.c +0 -132
  531. data/ext/dlahqr.c +0 -143
  532. data/ext/dlahr2.c +0 -112
  533. data/ext/dlahrd.c +0 -112
  534. data/ext/dlaic1.c +0 -89
  535. data/ext/dlaln2.c +0 -120
  536. data/ext/dlals0.c +0 -201
  537. data/ext/dlalsa.c +0 -270
  538. data/ext/dlalsd.c +0 -142
  539. data/ext/dlamrg.c +0 -80
  540. data/ext/dlaneg.c +0 -83
  541. data/ext/dlangb.c +0 -78
  542. data/ext/dlange.c +0 -74
  543. data/ext/dlangt.c +0 -87
  544. data/ext/dlanhs.c +0 -70
  545. data/ext/dlansb.c +0 -78
  546. data/ext/dlansf.c +0 -80
  547. data/ext/dlansp.c +0 -76
  548. data/ext/dlanst.c +0 -75
  549. data/ext/dlansy.c +0 -74
  550. data/ext/dlantb.c +0 -82
  551. data/ext/dlantp.c +0 -80
  552. data/ext/dlantr.c +0 -82
  553. data/ext/dlanv2.c +0 -82
  554. data/ext/dlapll.c +0 -105
  555. data/ext/dlapmr.c +0 -97
  556. data/ext/dlapmt.c +0 -101
  557. data/ext/dlapy2.c +0 -55
  558. data/ext/dlapy3.c +0 -59
  559. data/ext/dlaqgb.c +0 -117
  560. data/ext/dlaqge.c +0 -109
  561. data/ext/dlaqp2.c +0 -158
  562. data/ext/dlaqps.c +0 -208
  563. data/ext/dlaqr0.c +0 -153
  564. data/ext/dlaqr1.c +0 -82
  565. data/ext/dlaqr2.c +0 -182
  566. data/ext/dlaqr3.c +0 -182
  567. data/ext/dlaqr4.c +0 -153
  568. data/ext/dlaqr5.c +0 -200
  569. data/ext/dlaqsb.c +0 -101
  570. data/ext/dlaqsp.c +0 -94
  571. data/ext/dlaqsy.c +0 -97
  572. data/ext/dlaqtr.c +0 -114
  573. data/ext/dlar1v.c +0 -173
  574. data/ext/dlar2v.c +0 -149
  575. data/ext/dlarf.c +0 -101
  576. data/ext/dlarfb.c +0 -123
  577. data/ext/dlarfg.c +0 -83
  578. data/ext/dlarfgp.c +0 -83
  579. data/ext/dlarft.c +0 -105
  580. data/ext/dlarfx.c +0 -93
  581. data/ext/dlargv.c +0 -114
  582. data/ext/dlarnv.c +0 -83
  583. data/ext/dlarra.c +0 -124
  584. data/ext/dlarrb.c +0 -178
  585. data/ext/dlarrc.c +0 -96
  586. data/ext/dlarrd.c +0 -190
  587. data/ext/dlarre.c +0 -221
  588. data/ext/dlarrf.c +0 -176
  589. data/ext/dlarrj.c +0 -147
  590. data/ext/dlarrk.c +0 -97
  591. data/ext/dlarrr.c +0 -82
  592. data/ext/dlarrv.c +0 -271
  593. data/ext/dlarscl2.c +0 -82
  594. data/ext/dlartg.c +0 -61
  595. data/ext/dlartgp.c +0 -61
  596. data/ext/dlartgs.c +0 -62
  597. data/ext/dlartv.c +0 -130
  598. data/ext/dlaruv.c +0 -79
  599. data/ext/dlarz.c +0 -105
  600. data/ext/dlarzb.c +0 -127
  601. data/ext/dlarzt.c +0 -105
  602. data/ext/dlas2.c +0 -62
  603. data/ext/dlascl.c +0 -97
  604. data/ext/dlascl2.c +0 -82
  605. data/ext/dlasd0.c +0 -120
  606. data/ext/dlasd1.c +0 -162
  607. data/ext/dlasd2.c +0 -228
  608. data/ext/dlasd3.c +0 -202
  609. data/ext/dlasd4.c +0 -93
  610. data/ext/dlasd5.c +0 -90
  611. data/ext/dlasd6.c +0 -236
  612. data/ext/dlasd7.c +0 -225
  613. data/ext/dlasd8.c +0 -173
  614. data/ext/dlasda.c +0 -221
  615. data/ext/dlasdq.c +0 -186
  616. data/ext/dlasdt.c +0 -82
  617. data/ext/dlaset.c +0 -86
  618. data/ext/dlasq1.c +0 -96
  619. data/ext/dlasq2.c +0 -74
  620. data/ext/dlasq3.c +0 -138
  621. data/ext/dlasq4.c +0 -107
  622. data/ext/dlasq5.c +0 -94
  623. data/ext/dlasq6.c +0 -86
  624. data/ext/dlasr.c +0 -110
  625. data/ext/dlasrt.c +0 -74
  626. data/ext/dlassq.c +0 -70
  627. data/ext/dlasv2.c +0 -74
  628. data/ext/dlaswp.c +0 -94
  629. data/ext/dlasy2.c +0 -126
  630. data/ext/dlasyf.c +0 -97
  631. data/ext/dlat2s.c +0 -76
  632. data/ext/dlatbs.c +0 -130
  633. data/ext/dlatdf.c +0 -119
  634. data/ext/dlatps.c +0 -124
  635. data/ext/dlatrd.c +0 -105
  636. data/ext/dlatrs.c +0 -126
  637. data/ext/dlatrz.c +0 -87
  638. data/ext/dlatzm.c +0 -131
  639. data/ext/dlauu2.c +0 -77
  640. data/ext/dlauum.c +0 -77
  641. data/ext/dopgtr.c +0 -91
  642. data/ext/dopmtr.c +0 -116
  643. data/ext/dorbdb.c +0 -232
  644. data/ext/dorcsd.c +0 -197
  645. data/ext/dorg2l.c +0 -92
  646. data/ext/dorg2r.c +0 -92
  647. data/ext/dorgbr.c +0 -115
  648. data/ext/dorghr.c +0 -111
  649. data/ext/dorgl2.c +0 -90
  650. data/ext/dorglq.c +0 -107
  651. data/ext/dorgql.c +0 -107
  652. data/ext/dorgqr.c +0 -107
  653. data/ext/dorgr2.c +0 -90
  654. data/ext/dorgrq.c +0 -107
  655. data/ext/dorgtr.c +0 -107
  656. data/ext/dorm2l.c +0 -114
  657. data/ext/dorm2r.c +0 -114
  658. data/ext/dormbr.c +0 -139
  659. data/ext/dormhr.c +0 -133
  660. data/ext/dorml2.c +0 -110
  661. data/ext/dormlq.c +0 -125
  662. data/ext/dormql.c +0 -129
  663. data/ext/dormqr.c +0 -129
  664. data/ext/dormr2.c +0 -110
  665. data/ext/dormr3.c +0 -114
  666. data/ext/dormrq.c +0 -125
  667. data/ext/dormrz.c +0 -129
  668. data/ext/dormtr.c +0 -129
  669. data/ext/dpbcon.c +0 -82
  670. data/ext/dpbequ.c +0 -83
  671. data/ext/dpbrfs.c +0 -145
  672. data/ext/dpbstf.c +0 -81
  673. data/ext/dpbsv.c +0 -107
  674. data/ext/dpbsvx.c +0 -201
  675. data/ext/dpbtf2.c +0 -81
  676. data/ext/dpbtrf.c +0 -81
  677. data/ext/dpbtrs.c +0 -95
  678. data/ext/dpftrf.c +0 -82
  679. data/ext/dpftri.c +0 -82
  680. data/ext/dpftrs.c +0 -97
  681. data/ext/dpocon.c +0 -78
  682. data/ext/dpoequ.c +0 -75
  683. data/ext/dpoequb.c +0 -75
  684. data/ext/dporfs.c +0 -141
  685. data/ext/dporfsx.c +0 -206
  686. data/ext/dposv.c +0 -103
  687. data/ext/dposvx.c +0 -197
  688. data/ext/dposvxx.c +0 -235
  689. data/ext/dpotf2.c +0 -77
  690. data/ext/dpotrf.c +0 -77
  691. data/ext/dpotri.c +0 -77
  692. data/ext/dpotrs.c +0 -91
  693. data/ext/dppcon.c +0 -78
  694. data/ext/dppequ.c +0 -79
  695. data/ext/dpprfs.c +0 -139
  696. data/ext/dppsv.c +0 -104
  697. data/ext/dppsvx.c +0 -191
  698. data/ext/dpptrf.c +0 -78
  699. data/ext/dpptri.c +0 -78
  700. data/ext/dpptrs.c +0 -93
  701. data/ext/dpstf2.c +0 -95
  702. data/ext/dpstrf.c +0 -95
  703. data/ext/dptcon.c +0 -81
  704. data/ext/dpteqr.c +0 -126
  705. data/ext/dptrfs.c +0 -154
  706. data/ext/dptsv.c +0 -119
  707. data/ext/dptsvx.c +0 -168
  708. data/ext/dpttrf.c +0 -93
  709. data/ext/dpttrs.c +0 -97
  710. data/ext/dptts2.c +0 -94
  711. data/ext/drscl.c +0 -79
  712. data/ext/dsbev.c +0 -107
  713. data/ext/dsbevd.c +0 -140
  714. data/ext/dsbevx.c +0 -157
  715. data/ext/dsbgst.c +0 -117
  716. data/ext/dsbgv.c +0 -137
  717. data/ext/dsbgvd.c +0 -170
  718. data/ext/dsbgvx.c +0 -197
  719. data/ext/dsbtrd.c +0 -130
  720. data/ext/dsfrk.c +0 -109
  721. data/ext/dsgesv.c +0 -115
  722. data/ext/dspcon.c +0 -88
  723. data/ext/dspev.c +0 -102
  724. data/ext/dspevd.c +0 -135
  725. data/ext/dspevx.c +0 -141
  726. data/ext/dspgst.c +0 -94
  727. data/ext/dspgv.c +0 -129
  728. data/ext/dspgvd.c +0 -162
  729. data/ext/dspgvx.c +0 -168
  730. data/ext/dsposv.c +0 -111
  731. data/ext/dsprfs.c +0 -149
  732. data/ext/dspsv.c +0 -110
  733. data/ext/dspsvx.c +0 -163
  734. data/ext/dsptrd.c +0 -100
  735. data/ext/dsptrf.c +0 -84
  736. data/ext/dsptri.c +0 -89
  737. data/ext/dsptrs.c +0 -101
  738. data/ext/dstebz.c +0 -135
  739. data/ext/dstedc.c +0 -159
  740. data/ext/dstegr.c +0 -188
  741. data/ext/dstein.c +0 -134
  742. data/ext/dstemr.c +0 -193
  743. data/ext/dsteqr.c +0 -126
  744. data/ext/dsterf.c +0 -93
  745. data/ext/dstev.c +0 -111
  746. data/ext/dstevd.c +0 -144
  747. data/ext/dstevr.c +0 -188
  748. data/ext/dstevx.c +0 -158
  749. data/ext/dsycon.c +0 -90
  750. data/ext/dsyconv.c +0 -84
  751. data/ext/dsyequb.c +0 -82
  752. data/ext/dsyev.c +0 -107
  753. data/ext/dsyevd.c +0 -125
  754. data/ext/dsyevr.c +0 -172
  755. data/ext/dsyevx.c +0 -157
  756. data/ext/dsygs2.c +0 -95
  757. data/ext/dsygst.c +0 -95
  758. data/ext/dsygv.c +0 -137
  759. data/ext/dsygvd.c +0 -155
  760. data/ext/dsygvx.c +0 -187
  761. data/ext/dsyrfs.c +0 -153
  762. data/ext/dsyrfsx.c +0 -218
  763. data/ext/dsysv.c +0 -129
  764. data/ext/dsysvx.c +0 -183
  765. data/ext/dsysvxx.c +0 -258
  766. data/ext/dsyswapr.c +0 -82
  767. data/ext/dsytd2.c +0 -101
  768. data/ext/dsytf2.c +0 -85
  769. data/ext/dsytrd.c +0 -113
  770. data/ext/dsytrf.c +0 -97
  771. data/ext/dsytri.c +0 -92
  772. data/ext/dsytri2.c +0 -108
  773. data/ext/dsytri2x.c +0 -96
  774. data/ext/dsytrs.c +0 -103
  775. data/ext/dsytrs2.c +0 -106
  776. data/ext/dtbcon.c +0 -86
  777. data/ext/dtbrfs.c +0 -127
  778. data/ext/dtbtrs.c +0 -103
  779. data/ext/dtfsm.c +0 -110
  780. data/ext/dtftri.c +0 -86
  781. data/ext/dtfttp.c +0 -79
  782. data/ext/dtfttr.c +0 -80
  783. data/ext/dtgevc.c +0 -153
  784. data/ext/dtgex2.c +0 -180
  785. data/ext/dtgexc.c +0 -187
  786. data/ext/dtgsen.c +0 -252
  787. data/ext/dtgsja.c +0 -227
  788. data/ext/dtgsna.c +0 -164
  789. data/ext/dtgsy2.c +0 -182
  790. data/ext/dtgsyl.c +0 -190
  791. data/ext/dtpcon.c +0 -82
  792. data/ext/dtprfs.c +0 -123
  793. data/ext/dtptri.c +0 -82
  794. data/ext/dtptrs.c +0 -101
  795. data/ext/dtpttf.c +0 -79
  796. data/ext/dtpttr.c +0 -76
  797. data/ext/dtrcon.c +0 -82
  798. data/ext/dtrevc.c +0 -150
  799. data/ext/dtrexc.c +0 -116
  800. data/ext/dtrrfs.c +0 -123
  801. data/ext/dtrsen.c +0 -169
  802. data/ext/dtrsna.c +0 -137
  803. data/ext/dtrsyl.c +0 -116
  804. data/ext/dtrti2.c +0 -81
  805. data/ext/dtrtri.c +0 -81
  806. data/ext/dtrtrs.c +0 -99
  807. data/ext/dtrttf.c +0 -77
  808. data/ext/dtrttp.c +0 -73
  809. data/ext/dtzrqf.c +0 -83
  810. data/ext/dtzrzf.c +0 -101
  811. data/ext/dzsum1.c +0 -63
  812. data/ext/icmax1.c +0 -63
  813. data/ext/ieeeck.c +0 -59
  814. data/ext/ilaclc.c +0 -65
  815. data/ext/ilaclr.c +0 -65
  816. data/ext/iladiag.c +0 -51
  817. data/ext/iladlc.c +0 -65
  818. data/ext/iladlr.c +0 -65
  819. data/ext/ilaenv.c +0 -75
  820. data/ext/ilaprec.c +0 -51
  821. data/ext/ilaslc.c +0 -65
  822. data/ext/ilaslr.c +0 -65
  823. data/ext/ilatrans.c +0 -51
  824. data/ext/ilauplo.c +0 -51
  825. data/ext/ilaver.c +0 -53
  826. data/ext/ilazlc.c +0 -65
  827. data/ext/ilazlr.c +0 -65
  828. data/ext/iparmq.c +0 -75
  829. data/ext/izmax1.c +0 -63
  830. data/ext/lsamen.c +0 -59
  831. data/ext/rb_lapack.c +0 -3279
  832. data/ext/sbbcsd.c +0 -287
  833. data/ext/sbdsdc.c +0 -157
  834. data/ext/sbdsqr.c +0 -182
  835. data/ext/scsum1.c +0 -63
  836. data/ext/sdisna.c +0 -75
  837. data/ext/sgbbrd.c +0 -154
  838. data/ext/sgbcon.c +0 -98
  839. data/ext/sgbequ.c +0 -98
  840. data/ext/sgbequb.c +0 -96
  841. data/ext/sgbrfs.c +0 -161
  842. data/ext/sgbrfsx.c +0 -249
  843. data/ext/sgbsv.c +0 -115
  844. data/ext/sgbsvx.c +0 -286
  845. data/ext/sgbsvxx.c +0 -289
  846. data/ext/sgbtf2.c +0 -93
  847. data/ext/sgbtrf.c +0 -93
  848. data/ext/sgbtrs.c +0 -111
  849. data/ext/sgebak.c +0 -101
  850. data/ext/sgebal.c +0 -91
  851. data/ext/sgebd2.c +0 -112
  852. data/ext/sgebrd.c +0 -127
  853. data/ext/sgecon.c +0 -78
  854. data/ext/sgeequ.c +0 -88
  855. data/ext/sgeequb.c +0 -88
  856. data/ext/sgees.c +0 -148
  857. data/ext/sgeesx.c +0 -170
  858. data/ext/sgeev.c +0 -137
  859. data/ext/sgeevx.c +0 -181
  860. data/ext/sgegs.c +0 -171
  861. data/ext/sgegv.c +0 -171
  862. data/ext/sgehd2.c +0 -92
  863. data/ext/sgehrd.c +0 -107
  864. data/ext/sgejsv.c +0 -159
  865. data/ext/sgelq2.c +0 -86
  866. data/ext/sgelqf.c +0 -103
  867. data/ext/sgels.c +0 -137
  868. data/ext/sgelsd.c +0 -149
  869. data/ext/sgelss.c +0 -148
  870. data/ext/sgelsx.c +0 -136
  871. data/ext/sgelsy.c +0 -163
  872. data/ext/sgeql2.c +0 -88
  873. data/ext/sgeqlf.c +0 -103
  874. data/ext/sgeqp3.c +0 -126
  875. data/ext/sgeqpf.c +0 -111
  876. data/ext/sgeqr2.c +0 -88
  877. data/ext/sgeqr2p.c +0 -88
  878. data/ext/sgeqrf.c +0 -103
  879. data/ext/sgeqrfp.c +0 -103
  880. data/ext/sgerfs.c +0 -153
  881. data/ext/sgerfsx.c +0 -219
  882. data/ext/sgerq2.c +0 -86
  883. data/ext/sgerqf.c +0 -103
  884. data/ext/sgesc2.c +0 -108
  885. data/ext/sgesdd.c +0 -132
  886. data/ext/sgesv.c +0 -107
  887. data/ext/sgesvd.c +0 -143
  888. data/ext/sgesvj.c +0 -156
  889. data/ext/sgesvx.c +0 -278
  890. data/ext/sgesvxx.c +0 -281
  891. data/ext/sgetc2.c +0 -89
  892. data/ext/sgetf2.c +0 -85
  893. data/ext/sgetrf.c +0 -85
  894. data/ext/sgetri.c +0 -103
  895. data/ext/sgetrs.c +0 -103
  896. data/ext/sggbak.c +0 -113
  897. data/ext/sggbal.c +0 -128
  898. data/ext/sgges.c +0 -198
  899. data/ext/sggesx.c +0 -231
  900. data/ext/sggev.c +0 -171
  901. data/ext/sggevx.c +0 -229
  902. data/ext/sggglm.c +0 -156
  903. data/ext/sgghrd.c +0 -167
  904. data/ext/sgglse.c +0 -171
  905. data/ext/sggqrf.c +0 -137
  906. data/ext/sggrqf.c +0 -141
  907. data/ext/sggsvd.c +0 -181
  908. data/ext/sggsvp.c +0 -171
  909. data/ext/sgsvj0.c +0 -182
  910. data/ext/sgsvj1.c +0 -186
  911. data/ext/sgtcon.c +0 -124
  912. data/ext/sgtrfs.c +0 -209
  913. data/ext/sgtsv.c +0 -142
  914. data/ext/sgtsvx.c +0 -256
  915. data/ext/sgttrf.c +0 -132
  916. data/ext/sgttrs.c +0 -137
  917. data/ext/sgtts2.c +0 -134
  918. data/ext/shgeqz.c +0 -213
  919. data/ext/shsein.c +0 -205
  920. data/ext/shseqr.c +0 -153
  921. data/ext/sisnan.c +0 -51
  922. data/ext/sla_gbamv.c +0 -129
  923. data/ext/sla_gbrcond.c +0 -142
  924. data/ext/sla_gbrfsx_extended.c +0 -291
  925. data/ext/sla_gbrpvgrw.c +0 -87
  926. data/ext/sla_geamv.c +0 -119
  927. data/ext/sla_gercond.c +0 -134
  928. data/ext/sla_gerfsx_extended.c +0 -283
  929. data/ext/sla_lin_berr.c +0 -84
  930. data/ext/sla_porcond.c +0 -122
  931. data/ext/sla_porfsx_extended.c +0 -271
  932. data/ext/sla_porpvgrw.c +0 -95
  933. data/ext/sla_rpvgrw.c +0 -79
  934. data/ext/sla_syamv.c +0 -116
  935. data/ext/sla_syrcond.c +0 -134
  936. data/ext/sla_syrfsx_extended.c +0 -283
  937. data/ext/sla_syrpvgrw.c +0 -107
  938. data/ext/sla_wwaddw.c +0 -102
  939. data/ext/slabad.c +0 -54
  940. data/ext/slabrd.c +0 -132
  941. data/ext/slacn2.c +0 -106
  942. data/ext/slacon.c +0 -83
  943. data/ext/slacpy.c +0 -77
  944. data/ext/sladiv.c +0 -66
  945. data/ext/slae2.c +0 -62
  946. data/ext/slaebz.c +0 -218
  947. data/ext/slaed0.c +0 -127
  948. data/ext/slaed1.c +0 -133
  949. data/ext/slaed2.c +0 -189
  950. data/ext/slaed3.c +0 -161
  951. data/ext/slaed4.c +0 -90
  952. data/ext/slaed5.c +0 -87
  953. data/ext/slaed6.c +0 -90
  954. data/ext/slaed7.c +0 -248
  955. data/ext/slaed8.c +0 -206
  956. data/ext/slaed9.c +0 -111
  957. data/ext/slaeda.c +0 -160
  958. data/ext/slaein.c +0 -143
  959. data/ext/slaev2.c +0 -68
  960. data/ext/slaexc.c +0 -118
  961. data/ext/slag2.c +0 -91
  962. data/ext/slag2d.c +0 -76
  963. data/ext/slags2.c +0 -90
  964. data/ext/slagtf.c +0 -140
  965. data/ext/slagtm.c +0 -132
  966. data/ext/slagts.c +0 -139
  967. data/ext/slagv2.c +0 -132
  968. data/ext/slahqr.c +0 -143
  969. data/ext/slahr2.c +0 -112
  970. data/ext/slahrd.c +0 -114
  971. data/ext/slaic1.c +0 -89
  972. data/ext/slaln2.c +0 -120
  973. data/ext/slals0.c +0 -201
  974. data/ext/slalsa.c +0 -270
  975. data/ext/slalsd.c +0 -142
  976. data/ext/slamrg.c +0 -80
  977. data/ext/slaneg.c +0 -83
  978. data/ext/slangb.c +0 -78
  979. data/ext/slange.c +0 -74
  980. data/ext/slangt.c +0 -87
  981. data/ext/slanhs.c +0 -70
  982. data/ext/slansb.c +0 -78
  983. data/ext/slansf.c +0 -78
  984. data/ext/slansp.c +0 -76
  985. data/ext/slanst.c +0 -75
  986. data/ext/slansy.c +0 -74
  987. data/ext/slantb.c +0 -82
  988. data/ext/slantp.c +0 -80
  989. data/ext/slantr.c +0 -82
  990. data/ext/slanv2.c +0 -82
  991. data/ext/slapll.c +0 -105
  992. data/ext/slapmr.c +0 -97
  993. data/ext/slapmt.c +0 -101
  994. data/ext/slapy2.c +0 -55
  995. data/ext/slapy3.c +0 -59
  996. data/ext/slaqgb.c +0 -117
  997. data/ext/slaqge.c +0 -109
  998. data/ext/slaqp2.c +0 -158
  999. data/ext/slaqps.c +0 -208
  1000. data/ext/slaqr0.c +0 -153
  1001. data/ext/slaqr1.c +0 -82
  1002. data/ext/slaqr2.c +0 -182
  1003. data/ext/slaqr3.c +0 -182
  1004. data/ext/slaqr4.c +0 -153
  1005. data/ext/slaqr5.c +0 -200
  1006. data/ext/slaqsb.c +0 -101
  1007. data/ext/slaqsp.c +0 -94
  1008. data/ext/slaqsy.c +0 -97
  1009. data/ext/slaqtr.c +0 -114
  1010. data/ext/slar1v.c +0 -173
  1011. data/ext/slar2v.c +0 -149
  1012. data/ext/slarf.c +0 -101
  1013. data/ext/slarfb.c +0 -123
  1014. data/ext/slarfg.c +0 -83
  1015. data/ext/slarfgp.c +0 -83
  1016. data/ext/slarft.c +0 -105
  1017. data/ext/slarfx.c +0 -93
  1018. data/ext/slargv.c +0 -114
  1019. data/ext/slarnv.c +0 -83
  1020. data/ext/slarra.c +0 -124
  1021. data/ext/slarrb.c +0 -178
  1022. data/ext/slarrc.c +0 -96
  1023. data/ext/slarrd.c +0 -190
  1024. data/ext/slarre.c +0 -221
  1025. data/ext/slarrf.c +0 -176
  1026. data/ext/slarrj.c +0 -147
  1027. data/ext/slarrk.c +0 -97
  1028. data/ext/slarrr.c +0 -82
  1029. data/ext/slarrv.c +0 -271
  1030. data/ext/slarscl2.c +0 -82
  1031. data/ext/slartg.c +0 -61
  1032. data/ext/slartgp.c +0 -61
  1033. data/ext/slartgs.c +0 -62
  1034. data/ext/slartv.c +0 -130
  1035. data/ext/slaruv.c +0 -79
  1036. data/ext/slarz.c +0 -105
  1037. data/ext/slarzb.c +0 -127
  1038. data/ext/slarzt.c +0 -105
  1039. data/ext/slas2.c +0 -62
  1040. data/ext/slascl.c +0 -97
  1041. data/ext/slascl2.c +0 -82
  1042. data/ext/slasd0.c +0 -120
  1043. data/ext/slasd1.c +0 -160
  1044. data/ext/slasd2.c +0 -228
  1045. data/ext/slasd3.c +0 -212
  1046. data/ext/slasd4.c +0 -93
  1047. data/ext/slasd5.c +0 -90
  1048. data/ext/slasd6.c +0 -236
  1049. data/ext/slasd7.c +0 -225
  1050. data/ext/slasd8.c +0 -173
  1051. data/ext/slasda.c +0 -221
  1052. data/ext/slasdq.c +0 -186
  1053. data/ext/slasdt.c +0 -82
  1054. data/ext/slaset.c +0 -86
  1055. data/ext/slasq1.c +0 -96
  1056. data/ext/slasq2.c +0 -74
  1057. data/ext/slasq3.c +0 -138
  1058. data/ext/slasq4.c +0 -107
  1059. data/ext/slasq5.c +0 -94
  1060. data/ext/slasq6.c +0 -86
  1061. data/ext/slasr.c +0 -110
  1062. data/ext/slasrt.c +0 -74
  1063. data/ext/slassq.c +0 -70
  1064. data/ext/slasv2.c +0 -74
  1065. data/ext/slaswp.c +0 -94
  1066. data/ext/slasy2.c +0 -126
  1067. data/ext/slasyf.c +0 -97
  1068. data/ext/slatbs.c +0 -130
  1069. data/ext/slatdf.c +0 -119
  1070. data/ext/slatps.c +0 -124
  1071. data/ext/slatrd.c +0 -105
  1072. data/ext/slatrs.c +0 -126
  1073. data/ext/slatrz.c +0 -87
  1074. data/ext/slatzm.c +0 -131
  1075. data/ext/slauu2.c +0 -77
  1076. data/ext/slauum.c +0 -77
  1077. data/ext/sopgtr.c +0 -91
  1078. data/ext/sopmtr.c +0 -116
  1079. data/ext/sorbdb.c +0 -232
  1080. data/ext/sorcsd.c +0 -197
  1081. data/ext/sorg2l.c +0 -92
  1082. data/ext/sorg2r.c +0 -92
  1083. data/ext/sorgbr.c +0 -115
  1084. data/ext/sorghr.c +0 -111
  1085. data/ext/sorgl2.c +0 -90
  1086. data/ext/sorglq.c +0 -107
  1087. data/ext/sorgql.c +0 -107
  1088. data/ext/sorgqr.c +0 -107
  1089. data/ext/sorgr2.c +0 -90
  1090. data/ext/sorgrq.c +0 -107
  1091. data/ext/sorgtr.c +0 -107
  1092. data/ext/sorm2l.c +0 -114
  1093. data/ext/sorm2r.c +0 -114
  1094. data/ext/sormbr.c +0 -139
  1095. data/ext/sormhr.c +0 -133
  1096. data/ext/sorml2.c +0 -110
  1097. data/ext/sormlq.c +0 -125
  1098. data/ext/sormql.c +0 -129
  1099. data/ext/sormqr.c +0 -129
  1100. data/ext/sormr2.c +0 -110
  1101. data/ext/sormr3.c +0 -114
  1102. data/ext/sormrq.c +0 -125
  1103. data/ext/sormrz.c +0 -129
  1104. data/ext/sormtr.c +0 -129
  1105. data/ext/spbcon.c +0 -82
  1106. data/ext/spbequ.c +0 -83
  1107. data/ext/spbrfs.c +0 -145
  1108. data/ext/spbstf.c +0 -81
  1109. data/ext/spbsv.c +0 -107
  1110. data/ext/spbsvx.c +0 -201
  1111. data/ext/spbtf2.c +0 -81
  1112. data/ext/spbtrf.c +0 -81
  1113. data/ext/spbtrs.c +0 -95
  1114. data/ext/spftrf.c +0 -82
  1115. data/ext/spftri.c +0 -82
  1116. data/ext/spftrs.c +0 -97
  1117. data/ext/spocon.c +0 -78
  1118. data/ext/spoequ.c +0 -75
  1119. data/ext/spoequb.c +0 -75
  1120. data/ext/sporfs.c +0 -141
  1121. data/ext/sporfsx.c +0 -206
  1122. data/ext/sposv.c +0 -103
  1123. data/ext/sposvx.c +0 -197
  1124. data/ext/sposvxx.c +0 -235
  1125. data/ext/spotf2.c +0 -77
  1126. data/ext/spotrf.c +0 -77
  1127. data/ext/spotri.c +0 -77
  1128. data/ext/spotrs.c +0 -91
  1129. data/ext/sppcon.c +0 -78
  1130. data/ext/sppequ.c +0 -79
  1131. data/ext/spprfs.c +0 -139
  1132. data/ext/sppsv.c +0 -104
  1133. data/ext/sppsvx.c +0 -191
  1134. data/ext/spptrf.c +0 -78
  1135. data/ext/spptri.c +0 -78
  1136. data/ext/spptrs.c +0 -93
  1137. data/ext/spstf2.c +0 -95
  1138. data/ext/spstrf.c +0 -95
  1139. data/ext/sptcon.c +0 -81
  1140. data/ext/spteqr.c +0 -126
  1141. data/ext/sptrfs.c +0 -154
  1142. data/ext/sptsv.c +0 -119
  1143. data/ext/sptsvx.c +0 -168
  1144. data/ext/spttrf.c +0 -93
  1145. data/ext/spttrs.c +0 -97
  1146. data/ext/sptts2.c +0 -94
  1147. data/ext/srscl.c +0 -79
  1148. data/ext/ssbev.c +0 -107
  1149. data/ext/ssbevd.c +0 -140
  1150. data/ext/ssbevx.c +0 -157
  1151. data/ext/ssbgst.c +0 -117
  1152. data/ext/ssbgv.c +0 -137
  1153. data/ext/ssbgvd.c +0 -170
  1154. data/ext/ssbgvx.c +0 -197
  1155. data/ext/ssbtrd.c +0 -130
  1156. data/ext/ssfrk.c +0 -109
  1157. data/ext/sspcon.c +0 -88
  1158. data/ext/sspev.c +0 -102
  1159. data/ext/sspevd.c +0 -135
  1160. data/ext/sspevx.c +0 -141
  1161. data/ext/sspgst.c +0 -94
  1162. data/ext/sspgv.c +0 -129
  1163. data/ext/sspgvd.c +0 -162
  1164. data/ext/sspgvx.c +0 -168
  1165. data/ext/ssprfs.c +0 -149
  1166. data/ext/sspsv.c +0 -110
  1167. data/ext/sspsvx.c +0 -163
  1168. data/ext/ssptrd.c +0 -100
  1169. data/ext/ssptrf.c +0 -84
  1170. data/ext/ssptri.c +0 -89
  1171. data/ext/ssptrs.c +0 -101
  1172. data/ext/sstebz.c +0 -135
  1173. data/ext/sstedc.c +0 -159
  1174. data/ext/sstegr.c +0 -188
  1175. data/ext/sstein.c +0 -134
  1176. data/ext/sstemr.c +0 -193
  1177. data/ext/ssteqr.c +0 -126
  1178. data/ext/ssterf.c +0 -93
  1179. data/ext/sstev.c +0 -111
  1180. data/ext/sstevd.c +0 -144
  1181. data/ext/sstevr.c +0 -188
  1182. data/ext/sstevx.c +0 -158
  1183. data/ext/ssycon.c +0 -90
  1184. data/ext/ssyconv.c +0 -84
  1185. data/ext/ssyequb.c +0 -82
  1186. data/ext/ssyev.c +0 -107
  1187. data/ext/ssyevd.c +0 -125
  1188. data/ext/ssyevr.c +0 -172
  1189. data/ext/ssyevx.c +0 -157
  1190. data/ext/ssygs2.c +0 -95
  1191. data/ext/ssygst.c +0 -95
  1192. data/ext/ssygv.c +0 -137
  1193. data/ext/ssygvd.c +0 -155
  1194. data/ext/ssygvx.c +0 -191
  1195. data/ext/ssyrfs.c +0 -153
  1196. data/ext/ssyrfsx.c +0 -218
  1197. data/ext/ssysv.c +0 -129
  1198. data/ext/ssysvx.c +0 -183
  1199. data/ext/ssysvxx.c +0 -258
  1200. data/ext/ssyswapr.c +0 -82
  1201. data/ext/ssytd2.c +0 -101
  1202. data/ext/ssytf2.c +0 -85
  1203. data/ext/ssytrd.c +0 -113
  1204. data/ext/ssytrf.c +0 -97
  1205. data/ext/ssytri.c +0 -92
  1206. data/ext/ssytri2.c +0 -127
  1207. data/ext/ssytri2x.c +0 -96
  1208. data/ext/ssytrs.c +0 -103
  1209. data/ext/ssytrs2.c +0 -106
  1210. data/ext/stbcon.c +0 -86
  1211. data/ext/stbrfs.c +0 -127
  1212. data/ext/stbtrs.c +0 -103
  1213. data/ext/stfsm.c +0 -112
  1214. data/ext/stftri.c +0 -86
  1215. data/ext/stfttp.c +0 -79
  1216. data/ext/stfttr.c +0 -80
  1217. data/ext/stgevc.c +0 -153
  1218. data/ext/stgex2.c +0 -184
  1219. data/ext/stgexc.c +0 -191
  1220. data/ext/stgsen.c +0 -252
  1221. data/ext/stgsja.c +0 -227
  1222. data/ext/stgsna.c +0 -164
  1223. data/ext/stgsy2.c +0 -182
  1224. data/ext/stgsyl.c +0 -190
  1225. data/ext/stpcon.c +0 -82
  1226. data/ext/stprfs.c +0 -123
  1227. data/ext/stptri.c +0 -82
  1228. data/ext/stptrs.c +0 -101
  1229. data/ext/stpttf.c +0 -79
  1230. data/ext/stpttr.c +0 -76
  1231. data/ext/strcon.c +0 -82
  1232. data/ext/strevc.c +0 -150
  1233. data/ext/strexc.c +0 -116
  1234. data/ext/strrfs.c +0 -123
  1235. data/ext/strsen.c +0 -169
  1236. data/ext/strsna.c +0 -137
  1237. data/ext/strsyl.c +0 -116
  1238. data/ext/strti2.c +0 -81
  1239. data/ext/strtri.c +0 -81
  1240. data/ext/strtrs.c +0 -99
  1241. data/ext/strttf.c +0 -77
  1242. data/ext/strttp.c +0 -73
  1243. data/ext/stzrqf.c +0 -83
  1244. data/ext/stzrzf.c +0 -101
  1245. data/ext/xerbla.c +0 -52
  1246. data/ext/xerbla_array.c +0 -53
  1247. data/ext/zbbcsd.c +0 -283
  1248. data/ext/zbdsqr.c +0 -182
  1249. data/ext/zcgesv.c +0 -118
  1250. data/ext/zcposv.c +0 -114
  1251. data/ext/zdrscl.c +0 -79
  1252. data/ext/zgbbrd.c +0 -157
  1253. data/ext/zgbcon.c +0 -98
  1254. data/ext/zgbequ.c +0 -98
  1255. data/ext/zgbequb.c +0 -96
  1256. data/ext/zgbrfs.c +0 -161
  1257. data/ext/zgbrfsx.c +0 -249
  1258. data/ext/zgbsv.c +0 -115
  1259. data/ext/zgbsvx.c +0 -286
  1260. data/ext/zgbsvxx.c +0 -289
  1261. data/ext/zgbtf2.c +0 -93
  1262. data/ext/zgbtrf.c +0 -93
  1263. data/ext/zgbtrs.c +0 -111
  1264. data/ext/zgebak.c +0 -101
  1265. data/ext/zgebal.c +0 -91
  1266. data/ext/zgebd2.c +0 -112
  1267. data/ext/zgebrd.c +0 -127
  1268. data/ext/zgecon.c +0 -78
  1269. data/ext/zgeequ.c +0 -88
  1270. data/ext/zgeequb.c +0 -88
  1271. data/ext/zgees.c +0 -142
  1272. data/ext/zgeesx.c +0 -152
  1273. data/ext/zgeev.c +0 -132
  1274. data/ext/zgeevx.c +0 -173
  1275. data/ext/zgegs.c +0 -166
  1276. data/ext/zgegv.c +0 -171
  1277. data/ext/zgehd2.c +0 -92
  1278. data/ext/zgehrd.c +0 -107
  1279. data/ext/zgelq2.c +0 -86
  1280. data/ext/zgelqf.c +0 -103
  1281. data/ext/zgels.c +0 -137
  1282. data/ext/zgelsd.c +0 -154
  1283. data/ext/zgelss.c +0 -151
  1284. data/ext/zgelsx.c +0 -139
  1285. data/ext/zgelsy.c +0 -166
  1286. data/ext/zgeql2.c +0 -88
  1287. data/ext/zgeqlf.c +0 -103
  1288. data/ext/zgeqp3.c +0 -129
  1289. data/ext/zgeqpf.c +0 -114
  1290. data/ext/zgeqr2.c +0 -88
  1291. data/ext/zgeqr2p.c +0 -88
  1292. data/ext/zgeqrf.c +0 -103
  1293. data/ext/zgeqrfp.c +0 -103
  1294. data/ext/zgerfs.c +0 -153
  1295. data/ext/zgerfsx.c +0 -219
  1296. data/ext/zgerq2.c +0 -86
  1297. data/ext/zgerqf.c +0 -103
  1298. data/ext/zgesc2.c +0 -108
  1299. data/ext/zgesdd.c +0 -135
  1300. data/ext/zgesv.c +0 -107
  1301. data/ext/zgesvd.c +0 -146
  1302. data/ext/zgesvx.c +0 -278
  1303. data/ext/zgesvxx.c +0 -281
  1304. data/ext/zgetc2.c +0 -89
  1305. data/ext/zgetf2.c +0 -85
  1306. data/ext/zgetrf.c +0 -85
  1307. data/ext/zgetri.c +0 -103
  1308. data/ext/zgetrs.c +0 -103
  1309. data/ext/zggbak.c +0 -113
  1310. data/ext/zggbal.c +0 -128
  1311. data/ext/zgges.c +0 -192
  1312. data/ext/zggesx.c +0 -230
  1313. data/ext/zggev.c +0 -171
  1314. data/ext/zggevx.c +0 -226
  1315. data/ext/zggglm.c +0 -156
  1316. data/ext/zgghrd.c +0 -167
  1317. data/ext/zgglse.c +0 -171
  1318. data/ext/zggqrf.c +0 -137
  1319. data/ext/zggrqf.c +0 -141
  1320. data/ext/zggsvd.c +0 -184
  1321. data/ext/zggsvp.c +0 -174
  1322. data/ext/zgtcon.c +0 -121
  1323. data/ext/zgtrfs.c +0 -209
  1324. data/ext/zgtsv.c +0 -142
  1325. data/ext/zgtsvx.c +0 -256
  1326. data/ext/zgttrf.c +0 -132
  1327. data/ext/zgttrs.c +0 -137
  1328. data/ext/zgtts2.c +0 -134
  1329. data/ext/zhbev.c +0 -110
  1330. data/ext/zhbevd.c +0 -158
  1331. data/ext/zhbevx.c +0 -160
  1332. data/ext/zhbgst.c +0 -120
  1333. data/ext/zhbgv.c +0 -140
  1334. data/ext/zhbgvd.c +0 -188
  1335. data/ext/zhbgvx.c +0 -189
  1336. data/ext/zhbtrd.c +0 -130
  1337. data/ext/zhecon.c +0 -87
  1338. data/ext/zheequb.c +0 -82
  1339. data/ext/zheev.c +0 -110
  1340. data/ext/zheevd.c +0 -143
  1341. data/ext/zheevr.c +0 -190
  1342. data/ext/zheevx.c +0 -160
  1343. data/ext/zhegs2.c +0 -95
  1344. data/ext/zhegst.c +0 -95
  1345. data/ext/zhegv.c +0 -140
  1346. data/ext/zhegvd.c +0 -173
  1347. data/ext/zhegvx.c +0 -190
  1348. data/ext/zherfs.c +0 -153
  1349. data/ext/zherfsx.c +0 -218
  1350. data/ext/zhesv.c +0 -123
  1351. data/ext/zhesvx.c +0 -183
  1352. data/ext/zhesvxx.c +0 -258
  1353. data/ext/zhetd2.c +0 -101
  1354. data/ext/zhetf2.c +0 -85
  1355. data/ext/zhetrd.c +0 -113
  1356. data/ext/zhetrf.c +0 -97
  1357. data/ext/zhetri.c +0 -92
  1358. data/ext/zhetrs.c +0 -103
  1359. data/ext/zhetrs2.c +0 -106
  1360. data/ext/zhfrk.c +0 -109
  1361. data/ext/zhgeqz.c +0 -208
  1362. data/ext/zhpcon.c +0 -85
  1363. data/ext/zhpev.c +0 -105
  1364. data/ext/zhpevd.c +0 -153
  1365. data/ext/zhpevx.c +0 -144
  1366. data/ext/zhpgst.c +0 -94
  1367. data/ext/zhpgv.c +0 -132
  1368. data/ext/zhpgvd.c +0 -170
  1369. data/ext/zhpgvx.c +0 -170
  1370. data/ext/zhprfs.c +0 -149
  1371. data/ext/zhpsv.c +0 -110
  1372. data/ext/zhpsvx.c +0 -163
  1373. data/ext/zhptrd.c +0 -100
  1374. data/ext/zhptrf.c +0 -84
  1375. data/ext/zhptri.c +0 -89
  1376. data/ext/zhptrs.c +0 -101
  1377. data/ext/zhsein.c +0 -185
  1378. data/ext/zhseqr.c +0 -145
  1379. data/ext/zla_gbamv.c +0 -127
  1380. data/ext/zla_gbrcond_c.c +0 -142
  1381. data/ext/zla_gbrcond_x.c +0 -138
  1382. data/ext/zla_gbrfsx_extended.c +0 -295
  1383. data/ext/zla_gbrpvgrw.c +0 -87
  1384. data/ext/zla_geamv.c +0 -119
  1385. data/ext/zla_gercond_c.c +0 -134
  1386. data/ext/zla_gercond_x.c +0 -130
  1387. data/ext/zla_gerfsx_extended.c +0 -281
  1388. data/ext/zla_heamv.c +0 -116
  1389. data/ext/zla_hercond_c.c +0 -134
  1390. data/ext/zla_hercond_x.c +0 -130
  1391. data/ext/zla_herfsx_extended.c +0 -283
  1392. data/ext/zla_herpvgrw.c +0 -107
  1393. data/ext/zla_lin_berr.c +0 -84
  1394. data/ext/zla_porcond_c.c +0 -122
  1395. data/ext/zla_porcond_x.c +0 -118
  1396. data/ext/zla_porfsx_extended.c +0 -271
  1397. data/ext/zla_porpvgrw.c +0 -95
  1398. data/ext/zla_rpvgrw.c +0 -79
  1399. data/ext/zla_syamv.c +0 -116
  1400. data/ext/zla_syrcond_c.c +0 -134
  1401. data/ext/zla_syrcond_x.c +0 -130
  1402. data/ext/zla_syrfsx_extended.c +0 -283
  1403. data/ext/zla_syrpvgrw.c +0 -107
  1404. data/ext/zla_wwaddw.c +0 -102
  1405. data/ext/zlabrd.c +0 -132
  1406. data/ext/zlacgv.c +0 -75
  1407. data/ext/zlacn2.c +0 -103
  1408. data/ext/zlacon.c +0 -80
  1409. data/ext/zlacp2.c +0 -77
  1410. data/ext/zlacpy.c +0 -77
  1411. data/ext/zlacrm.c +0 -90
  1412. data/ext/zlacrt.c +0 -108
  1413. data/ext/zladiv.c +0 -57
  1414. data/ext/zlaed0.c +0 -134
  1415. data/ext/zlaed7.c +0 -247
  1416. data/ext/zlaed8.c +0 -198
  1417. data/ext/zlaein.c +0 -113
  1418. data/ext/zlaesy.c +0 -74
  1419. data/ext/zlaev2.c +0 -71
  1420. data/ext/zlag2c.c +0 -76
  1421. data/ext/zlags2.c +0 -92
  1422. data/ext/zlagtm.c +0 -132
  1423. data/ext/zlahef.c +0 -97
  1424. data/ext/zlahqr.c +0 -135
  1425. data/ext/zlahr2.c +0 -112
  1426. data/ext/zlahrd.c +0 -112
  1427. data/ext/zlaic1.c +0 -90
  1428. data/ext/zlals0.c +0 -201
  1429. data/ext/zlalsa.c +0 -270
  1430. data/ext/zlalsd.c +0 -145
  1431. data/ext/zlangb.c +0 -76
  1432. data/ext/zlange.c +0 -74
  1433. data/ext/zlangt.c +0 -87
  1434. data/ext/zlanhb.c +0 -78
  1435. data/ext/zlanhe.c +0 -74
  1436. data/ext/zlanhf.c +0 -80
  1437. data/ext/zlanhp.c +0 -76
  1438. data/ext/zlanhs.c +0 -70
  1439. data/ext/zlanht.c +0 -75
  1440. data/ext/zlansb.c +0 -78
  1441. data/ext/zlansp.c +0 -76
  1442. data/ext/zlansy.c +0 -74
  1443. data/ext/zlantb.c +0 -82
  1444. data/ext/zlantp.c +0 -80
  1445. data/ext/zlantr.c +0 -82
  1446. data/ext/zlapll.c +0 -105
  1447. data/ext/zlapmr.c +0 -97
  1448. data/ext/zlapmt.c +0 -101
  1449. data/ext/zlaqgb.c +0 -117
  1450. data/ext/zlaqge.c +0 -109
  1451. data/ext/zlaqhb.c +0 -97
  1452. data/ext/zlaqhe.c +0 -97
  1453. data/ext/zlaqhp.c +0 -94
  1454. data/ext/zlaqp2.c +0 -158
  1455. data/ext/zlaqps.c +0 -208
  1456. data/ext/zlaqr0.c +0 -153
  1457. data/ext/zlaqr1.c +0 -76
  1458. data/ext/zlaqr2.c +0 -174
  1459. data/ext/zlaqr3.c +0 -174
  1460. data/ext/zlaqr4.c +0 -147
  1461. data/ext/zlaqr5.c +0 -179
  1462. data/ext/zlaqsb.c +0 -101
  1463. data/ext/zlaqsp.c +0 -94
  1464. data/ext/zlaqsy.c +0 -97
  1465. data/ext/zlar1v.c +0 -173
  1466. data/ext/zlar2v.c +0 -149
  1467. data/ext/zlarcm.c +0 -86
  1468. data/ext/zlarf.c +0 -102
  1469. data/ext/zlarfb.c +0 -123
  1470. data/ext/zlarfg.c +0 -84
  1471. data/ext/zlarfgp.c +0 -84
  1472. data/ext/zlarft.c +0 -105
  1473. data/ext/zlarfx.c +0 -94
  1474. data/ext/zlargv.c +0 -114
  1475. data/ext/zlarnv.c +0 -83
  1476. data/ext/zlarrv.c +0 -271
  1477. data/ext/zlarscl2.c +0 -82
  1478. data/ext/zlartg.c +0 -63
  1479. data/ext/zlartv.c +0 -130
  1480. data/ext/zlarz.c +0 -106
  1481. data/ext/zlarzb.c +0 -127
  1482. data/ext/zlarzt.c +0 -105
  1483. data/ext/zlascl.c +0 -97
  1484. data/ext/zlascl2.c +0 -82
  1485. data/ext/zlaset.c +0 -88
  1486. data/ext/zlasr.c +0 -110
  1487. data/ext/zlassq.c +0 -70
  1488. data/ext/zlaswp.c +0 -94
  1489. data/ext/zlasyf.c +0 -97
  1490. data/ext/zlat2c.c +0 -76
  1491. data/ext/zlatbs.c +0 -130
  1492. data/ext/zlatdf.c +0 -119
  1493. data/ext/zlatps.c +0 -124
  1494. data/ext/zlatrd.c +0 -105
  1495. data/ext/zlatrs.c +0 -126
  1496. data/ext/zlatrz.c +0 -87
  1497. data/ext/zlatzm.c +0 -132
  1498. data/ext/zlauu2.c +0 -77
  1499. data/ext/zlauum.c +0 -77
  1500. data/ext/zpbcon.c +0 -82
  1501. data/ext/zpbequ.c +0 -83
  1502. data/ext/zpbrfs.c +0 -145
  1503. data/ext/zpbstf.c +0 -81
  1504. data/ext/zpbsv.c +0 -107
  1505. data/ext/zpbsvx.c +0 -201
  1506. data/ext/zpbtf2.c +0 -81
  1507. data/ext/zpbtrf.c +0 -81
  1508. data/ext/zpbtrs.c +0 -95
  1509. data/ext/zpftrf.c +0 -82
  1510. data/ext/zpftri.c +0 -82
  1511. data/ext/zpftrs.c +0 -97
  1512. data/ext/zpocon.c +0 -78
  1513. data/ext/zpoequ.c +0 -75
  1514. data/ext/zpoequb.c +0 -75
  1515. data/ext/zporfs.c +0 -141
  1516. data/ext/zporfsx.c +0 -206
  1517. data/ext/zposv.c +0 -103
  1518. data/ext/zposvx.c +0 -197
  1519. data/ext/zposvxx.c +0 -235
  1520. data/ext/zpotf2.c +0 -77
  1521. data/ext/zpotrf.c +0 -77
  1522. data/ext/zpotri.c +0 -77
  1523. data/ext/zpotrs.c +0 -91
  1524. data/ext/zppcon.c +0 -78
  1525. data/ext/zppequ.c +0 -79
  1526. data/ext/zpprfs.c +0 -139
  1527. data/ext/zppsv.c +0 -104
  1528. data/ext/zppsvx.c +0 -191
  1529. data/ext/zpptrf.c +0 -78
  1530. data/ext/zpptri.c +0 -78
  1531. data/ext/zpptrs.c +0 -93
  1532. data/ext/zpstf2.c +0 -95
  1533. data/ext/zpstrf.c +0 -95
  1534. data/ext/zptcon.c +0 -81
  1535. data/ext/zpteqr.c +0 -126
  1536. data/ext/zptrfs.c +0 -161
  1537. data/ext/zptsv.c +0 -123
  1538. data/ext/zptsvx.c +0 -171
  1539. data/ext/zpttrf.c +0 -93
  1540. data/ext/zpttrs.c +0 -101
  1541. data/ext/zptts2.c +0 -98
  1542. data/ext/zrot.c +0 -107
  1543. data/ext/zspcon.c +0 -85
  1544. data/ext/zspmv.c +0 -117
  1545. data/ext/zspr.c +0 -96
  1546. data/ext/zsprfs.c +0 -149
  1547. data/ext/zspsv.c +0 -110
  1548. data/ext/zspsvx.c +0 -163
  1549. data/ext/zsptrf.c +0 -84
  1550. data/ext/zsptri.c +0 -89
  1551. data/ext/zsptrs.c +0 -101
  1552. data/ext/zstedc.c +0 -177
  1553. data/ext/zstegr.c +0 -188
  1554. data/ext/zstein.c +0 -134
  1555. data/ext/zstemr.c +0 -193
  1556. data/ext/zsteqr.c +0 -126
  1557. data/ext/zsycon.c +0 -87
  1558. data/ext/zsyconv.c +0 -84
  1559. data/ext/zsyequb.c +0 -82
  1560. data/ext/zsymv.c +0 -115
  1561. data/ext/zsyr.c +0 -95
  1562. data/ext/zsyrfs.c +0 -153
  1563. data/ext/zsyrfsx.c +0 -218
  1564. data/ext/zsysv.c +0 -129
  1565. data/ext/zsysvx.c +0 -183
  1566. data/ext/zsysvxx.c +0 -258
  1567. data/ext/zsyswapr.c +0 -82
  1568. data/ext/zsytf2.c +0 -85
  1569. data/ext/zsytrf.c +0 -97
  1570. data/ext/zsytri.c +0 -92
  1571. data/ext/zsytri2.c +0 -104
  1572. data/ext/zsytri2x.c +0 -96
  1573. data/ext/zsytrs.c +0 -103
  1574. data/ext/zsytrs2.c +0 -106
  1575. data/ext/ztbcon.c +0 -86
  1576. data/ext/ztbrfs.c +0 -127
  1577. data/ext/ztbtrs.c +0 -103
  1578. data/ext/ztfsm.c +0 -111
  1579. data/ext/ztftri.c +0 -86
  1580. data/ext/ztfttp.c +0 -79
  1581. data/ext/ztfttr.c +0 -80
  1582. data/ext/ztgevc.c +0 -156
  1583. data/ext/ztgex2.c +0 -171
  1584. data/ext/ztgexc.c +0 -172
  1585. data/ext/ztgsen.c +0 -244
  1586. data/ext/ztgsja.c +0 -227
  1587. data/ext/ztgsna.c +0 -164
  1588. data/ext/ztgsy2.c +0 -176
  1589. data/ext/ztgsyl.c +0 -190
  1590. data/ext/ztpcon.c +0 -82
  1591. data/ext/ztprfs.c +0 -123
  1592. data/ext/ztptri.c +0 -82
  1593. data/ext/ztptrs.c +0 -101
  1594. data/ext/ztpttf.c +0 -79
  1595. data/ext/ztpttr.c +0 -76
  1596. data/ext/ztrcon.c +0 -82
  1597. data/ext/ztrevc.c +0 -154
  1598. data/ext/ztrexc.c +0 -111
  1599. data/ext/ztrrfs.c +0 -123
  1600. data/ext/ztrsen.c +0 -154
  1601. data/ext/ztrsna.c +0 -137
  1602. data/ext/ztrsyl.c +0 -116
  1603. data/ext/ztrti2.c +0 -81
  1604. data/ext/ztrtri.c +0 -81
  1605. data/ext/ztrtrs.c +0 -99
  1606. data/ext/ztrttf.c +0 -77
  1607. data/ext/ztrttp.c +0 -73
  1608. data/ext/ztzrqf.c +0 -83
  1609. data/ext/ztzrzf.c +0 -101
  1610. data/ext/zunbdb.c +0 -232
  1611. data/ext/zuncsd.c +0 -204
  1612. data/ext/zung2l.c +0 -92
  1613. data/ext/zung2r.c +0 -92
  1614. data/ext/zungbr.c +0 -115
  1615. data/ext/zunghr.c +0 -111
  1616. data/ext/zungl2.c +0 -90
  1617. data/ext/zunglq.c +0 -107
  1618. data/ext/zungql.c +0 -107
  1619. data/ext/zungqr.c +0 -107
  1620. data/ext/zungr2.c +0 -90
  1621. data/ext/zungrq.c +0 -107
  1622. data/ext/zungtr.c +0 -107
  1623. data/ext/zunm2l.c +0 -114
  1624. data/ext/zunm2r.c +0 -114
  1625. data/ext/zunmbr.c +0 -139
  1626. data/ext/zunmhr.c +0 -133
  1627. data/ext/zunml2.c +0 -110
  1628. data/ext/zunmlq.c +0 -125
  1629. data/ext/zunmql.c +0 -129
  1630. data/ext/zunmqr.c +0 -129
  1631. data/ext/zunmr2.c +0 -110
  1632. data/ext/zunmr3.c +0 -114
  1633. data/ext/zunmrq.c +0 -125
  1634. data/ext/zunmrz.c +0 -129
  1635. data/ext/zunmtr.c +0 -129
  1636. data/ext/zupgtr.c +0 -91
  1637. data/ext/zupmtr.c +0 -116
@@ -1,95 +0,0 @@
1
- #include "rb_lapack.h"
2
-
3
- extern real sla_porpvgrw_(char* uplo, integer* ncols, real* a, integer* lda, real* af, integer* ldaf, real* work);
4
-
5
-
6
- static VALUE
7
- rblapack_sla_porpvgrw(int argc, VALUE *argv, VALUE self){
8
- VALUE rblapack_uplo;
9
- char uplo;
10
- VALUE rblapack_ncols;
11
- integer ncols;
12
- VALUE rblapack_a;
13
- real *a;
14
- VALUE rblapack_af;
15
- real *af;
16
- VALUE rblapack_work;
17
- real *work;
18
- VALUE rblapack___out__;
19
- real __out__;
20
-
21
- integer lda;
22
- integer n;
23
- integer ldaf;
24
-
25
- VALUE rblapack_options;
26
- if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
27
- argc--;
28
- rblapack_options = argv[argc];
29
- if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
30
- printf("%s\n", "USAGE:\n __out__ = NumRu::Lapack.sla_porpvgrw( uplo, ncols, a, af, work, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n REAL FUNCTION SLA_PORPVGRW( UPLO, NCOLS, A, LDA, AF, LDAF, WORK )\n\n* Purpose\n* =======\n* \n* SLA_PORPVGRW computes the reciprocal pivot growth factor\n* norm(A)/norm(U). The \"max absolute element\" norm is used. If this is\n* much less than 1, the stability of the LU factorization of the\n* (equilibrated) matrix A could be poor. This also means that the\n* solution X, estimated condition numbers, and error bounds could be\n* unreliable.\n*\n\n* Arguments\n* =========\n*\n* UPLO (input) CHARACTER*1\n* = 'U': Upper triangle of A is stored;\n* = 'L': Lower triangle of A is stored.\n*\n* NCOLS (input) INTEGER\n* The number of columns of the matrix A. NCOLS >= 0.\n*\n* A (input) REAL array, dimension (LDA,N)\n* On entry, the N-by-N matrix A.\n*\n* LDA (input) INTEGER\n* The leading dimension of the array A. LDA >= max(1,N).\n*\n* AF (input) REAL array, dimension (LDAF,N)\n* The triangular factor U or L from the Cholesky factorization\n* A = U**T*U or A = L*L**T, as computed by SPOTRF.\n*\n* LDAF (input) INTEGER\n* The leading dimension of the array AF. LDAF >= max(1,N).\n*\n* WORK (input) REAL array, dimension (2*N)\n*\n\n* =====================================================================\n*\n* .. Local Scalars ..\n INTEGER I, J\n REAL AMAX, UMAX, RPVGRW\n LOGICAL UPPER\n* ..\n* .. Intrinsic Functions ..\n INTRINSIC ABS, MAX, MIN\n* ..\n* .. External Functions ..\n EXTERNAL LSAME, SLASET\n LOGICAL LSAME\n* ..\n\n");
31
- return Qnil;
32
- }
33
- if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
34
- printf("%s\n", "USAGE:\n __out__ = NumRu::Lapack.sla_porpvgrw( uplo, ncols, a, af, work, [:usage => usage, :help => help])\n");
35
- return Qnil;
36
- }
37
- } else
38
- rblapack_options = Qnil;
39
- if (argc != 5 && argc != 5)
40
- rb_raise(rb_eArgError,"wrong number of arguments (%d for 5)", argc);
41
- rblapack_uplo = argv[0];
42
- rblapack_ncols = argv[1];
43
- rblapack_a = argv[2];
44
- rblapack_af = argv[3];
45
- rblapack_work = argv[4];
46
- if (argc == 5) {
47
- } else if (rblapack_options != Qnil) {
48
- } else {
49
- }
50
-
51
- uplo = StringValueCStr(rblapack_uplo)[0];
52
- if (!NA_IsNArray(rblapack_a))
53
- rb_raise(rb_eArgError, "a (3th argument) must be NArray");
54
- if (NA_RANK(rblapack_a) != 2)
55
- rb_raise(rb_eArgError, "rank of a (3th argument) must be %d", 2);
56
- lda = NA_SHAPE0(rblapack_a);
57
- n = NA_SHAPE1(rblapack_a);
58
- if (NA_TYPE(rblapack_a) != NA_SFLOAT)
59
- rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
60
- a = NA_PTR_TYPE(rblapack_a, real*);
61
- ncols = NUM2INT(rblapack_ncols);
62
- if (!NA_IsNArray(rblapack_af))
63
- rb_raise(rb_eArgError, "af (4th argument) must be NArray");
64
- if (NA_RANK(rblapack_af) != 2)
65
- rb_raise(rb_eArgError, "rank of af (4th argument) must be %d", 2);
66
- ldaf = NA_SHAPE0(rblapack_af);
67
- if (NA_SHAPE1(rblapack_af) != n)
68
- rb_raise(rb_eRuntimeError, "shape 1 of af must be the same as shape 1 of a");
69
- if (NA_TYPE(rblapack_af) != NA_SFLOAT)
70
- rblapack_af = na_change_type(rblapack_af, NA_SFLOAT);
71
- af = NA_PTR_TYPE(rblapack_af, real*);
72
- if (!NA_IsNArray(rblapack_work))
73
- rb_raise(rb_eArgError, "work (5th argument) must be NArray");
74
- if (NA_RANK(rblapack_work) != 1)
75
- rb_raise(rb_eArgError, "rank of work (5th argument) must be %d", 1);
76
- if (NA_SHAPE0(rblapack_work) != (2*n))
77
- rb_raise(rb_eRuntimeError, "shape 0 of work must be %d", 2*n);
78
- if (NA_TYPE(rblapack_work) != NA_SFLOAT)
79
- rblapack_work = na_change_type(rblapack_work, NA_SFLOAT);
80
- work = NA_PTR_TYPE(rblapack_work, real*);
81
-
82
- __out__ = sla_porpvgrw_(&uplo, &ncols, a, &lda, af, &ldaf, work);
83
-
84
- rblapack___out__ = rb_float_new((double)__out__);
85
- return rblapack___out__;
86
- }
87
-
88
- void
89
- init_lapack_sla_porpvgrw(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
90
- sHelp = sH;
91
- sUsage = sU;
92
- rblapack_ZERO = zero;
93
-
94
- rb_define_module_function(mLapack, "sla_porpvgrw", rblapack_sla_porpvgrw, -1);
95
- }
@@ -1,79 +0,0 @@
1
- #include "rb_lapack.h"
2
-
3
- extern real sla_rpvgrw_(integer* n, integer* ncols, real* a, integer* lda, real* af, integer* ldaf);
4
-
5
-
6
- static VALUE
7
- rblapack_sla_rpvgrw(int argc, VALUE *argv, VALUE self){
8
- VALUE rblapack_ncols;
9
- integer ncols;
10
- VALUE rblapack_a;
11
- real *a;
12
- VALUE rblapack_af;
13
- real *af;
14
- VALUE rblapack___out__;
15
- real __out__;
16
-
17
- integer lda;
18
- integer n;
19
- integer ldaf;
20
-
21
- VALUE rblapack_options;
22
- if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
23
- argc--;
24
- rblapack_options = argv[argc];
25
- if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
26
- printf("%s\n", "USAGE:\n __out__ = NumRu::Lapack.sla_rpvgrw( ncols, a, af, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n REAL FUNCTION SLA_RPVGRW( N, NCOLS, A, LDA, AF, LDAF )\n\n* Purpose\n* =======\n*\n* SLA_RPVGRW computes the reciprocal pivot growth factor\n* norm(A)/norm(U). The \"max absolute element\" norm is used. If this is\n* much less than 1, the stability of the LU factorization of the\n* (equilibrated) matrix A could be poor. This also means that the\n* solution X, estimated condition numbers, and error bounds could be\n* unreliable.\n*\n\n* Arguments\n* =========\n*\n* N (input) INTEGER\n* The number of linear equations, i.e., the order of the\n* matrix A. N >= 0.\n*\n* NCOLS (input) INTEGER\n* The number of columns of the matrix A. NCOLS >= 0.\n*\n* A (input) REAL array, dimension (LDA,N)\n* On entry, the N-by-N matrix A.\n*\n* LDA (input) INTEGER\n* The leading dimension of the array A. LDA >= max(1,N).\n*\n* AF (input) REAL array, dimension (LDAF,N)\n* The factors L and U from the factorization\n* A = P*L*U as computed by SGETRF.\n*\n* LDAF (input) INTEGER\n* The leading dimension of the array AF. LDAF >= max(1,N).\n*\n\n* =====================================================================\n*\n* .. Local Scalars ..\n INTEGER I, J\n REAL AMAX, UMAX, RPVGRW\n* ..\n* .. Intrinsic Functions ..\n INTRINSIC ABS, MAX, MIN\n* ..\n\n");
27
- return Qnil;
28
- }
29
- if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
30
- printf("%s\n", "USAGE:\n __out__ = NumRu::Lapack.sla_rpvgrw( ncols, a, af, [:usage => usage, :help => help])\n");
31
- return Qnil;
32
- }
33
- } else
34
- rblapack_options = Qnil;
35
- if (argc != 3 && argc != 3)
36
- rb_raise(rb_eArgError,"wrong number of arguments (%d for 3)", argc);
37
- rblapack_ncols = argv[0];
38
- rblapack_a = argv[1];
39
- rblapack_af = argv[2];
40
- if (argc == 3) {
41
- } else if (rblapack_options != Qnil) {
42
- } else {
43
- }
44
-
45
- ncols = NUM2INT(rblapack_ncols);
46
- if (!NA_IsNArray(rblapack_af))
47
- rb_raise(rb_eArgError, "af (3th argument) must be NArray");
48
- if (NA_RANK(rblapack_af) != 2)
49
- rb_raise(rb_eArgError, "rank of af (3th argument) must be %d", 2);
50
- ldaf = NA_SHAPE0(rblapack_af);
51
- n = NA_SHAPE1(rblapack_af);
52
- if (NA_TYPE(rblapack_af) != NA_SFLOAT)
53
- rblapack_af = na_change_type(rblapack_af, NA_SFLOAT);
54
- af = NA_PTR_TYPE(rblapack_af, real*);
55
- if (!NA_IsNArray(rblapack_a))
56
- rb_raise(rb_eArgError, "a (2th argument) must be NArray");
57
- if (NA_RANK(rblapack_a) != 2)
58
- rb_raise(rb_eArgError, "rank of a (2th argument) must be %d", 2);
59
- lda = NA_SHAPE0(rblapack_a);
60
- if (NA_SHAPE1(rblapack_a) != n)
61
- rb_raise(rb_eRuntimeError, "shape 1 of a must be the same as shape 1 of af");
62
- if (NA_TYPE(rblapack_a) != NA_SFLOAT)
63
- rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
64
- a = NA_PTR_TYPE(rblapack_a, real*);
65
-
66
- __out__ = sla_rpvgrw_(&n, &ncols, a, &lda, af, &ldaf);
67
-
68
- rblapack___out__ = rb_float_new((double)__out__);
69
- return rblapack___out__;
70
- }
71
-
72
- void
73
- init_lapack_sla_rpvgrw(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
74
- sHelp = sH;
75
- sUsage = sU;
76
- rblapack_ZERO = zero;
77
-
78
- rb_define_module_function(mLapack, "sla_rpvgrw", rblapack_sla_rpvgrw, -1);
79
- }
@@ -1,116 +0,0 @@
1
- #include "rb_lapack.h"
2
-
3
- extern VOID sla_syamv_(integer* uplo, integer* n, real* alpha, real* a, integer* lda, real* x, integer* incx, real* beta, real* y, integer* incy);
4
-
5
-
6
- static VALUE
7
- rblapack_sla_syamv(int argc, VALUE *argv, VALUE self){
8
- VALUE rblapack_uplo;
9
- integer uplo;
10
- VALUE rblapack_alpha;
11
- real alpha;
12
- VALUE rblapack_a;
13
- real *a;
14
- VALUE rblapack_x;
15
- real *x;
16
- VALUE rblapack_incx;
17
- integer incx;
18
- VALUE rblapack_beta;
19
- real beta;
20
- VALUE rblapack_y;
21
- real *y;
22
- VALUE rblapack_incy;
23
- integer incy;
24
- VALUE rblapack_y_out__;
25
- real *y_out__;
26
-
27
- integer lda;
28
- integer n;
29
-
30
- VALUE rblapack_options;
31
- if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
32
- argc--;
33
- rblapack_options = argv[argc];
34
- if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
35
- printf("%s\n", "USAGE:\n y = NumRu::Lapack.sla_syamv( uplo, alpha, a, x, incx, beta, y, incy, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n SUBROUTINE SLA_SYAMV( UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY )\n\n* Purpose\n* =======\n*\n* SLA_SYAMV performs the matrix-vector operation\n*\n* y := alpha*abs(A)*abs(x) + beta*abs(y),\n*\n* where alpha and beta are scalars, x and y are vectors and A is an\n* n by n symmetric matrix.\n*\n* This function is primarily used in calculating error bounds.\n* To protect against underflow during evaluation, components in\n* the resulting vector are perturbed away from zero by (N+1)\n* times the underflow threshold. To prevent unnecessarily large\n* errors for block-structure embedded in general matrices,\n* \"symbolically\" zero components are not perturbed. A zero\n* entry is considered \"symbolic\" if all multiplications involved\n* in computing that entry have at least one zero multiplicand.\n*\n\n* Arguments\n* ==========\n*\n* UPLO (input) INTEGER\n* On entry, UPLO specifies whether the upper or lower\n* triangular part of the array A is to be referenced as\n* follows:\n*\n* UPLO = BLAS_UPPER Only the upper triangular part of A\n* is to be referenced.\n*\n* UPLO = BLAS_LOWER Only the lower triangular part of A\n* is to be referenced.\n*\n* Unchanged on exit.\n*\n* N (input) INTEGER\n* On entry, N specifies the number of columns of the matrix A.\n* N must be at least zero.\n* Unchanged on exit.\n*\n* ALPHA (input) REAL .\n* On entry, ALPHA specifies the scalar alpha.\n* Unchanged on exit.\n*\n* A - REAL array of DIMENSION ( LDA, n ).\n* Before entry, the leading m by n part of the array A must\n* contain the matrix of coefficients.\n* Unchanged on exit.\n*\n* LDA (input) INTEGER\n* On entry, LDA specifies the first dimension of A as declared\n* in the calling (sub) program. LDA must be at least\n* max( 1, n ).\n* Unchanged on exit.\n*\n* X (input) REAL array, dimension\n* ( 1 + ( n - 1 )*abs( INCX ) )\n* Before entry, the incremented array X must contain the\n* vector x.\n* Unchanged on exit.\n*\n* INCX (input) INTEGER\n* On entry, INCX specifies the increment for the elements of\n* X. INCX must not be zero.\n* Unchanged on exit.\n*\n* BETA (input) REAL .\n* On entry, BETA specifies the scalar beta. When BETA is\n* supplied as zero then Y need not be set on input.\n* Unchanged on exit.\n*\n* Y (input/output) REAL array, dimension\n* ( 1 + ( n - 1 )*abs( INCY ) )\n* Before entry with BETA non-zero, the incremented array Y\n* must contain the vector y. On exit, Y is overwritten by the\n* updated vector y.\n*\n* INCY (input) INTEGER\n* On entry, INCY specifies the increment for the elements of\n* Y. INCY must not be zero.\n* Unchanged on exit.\n*\n\n* Further Details\n* ===============\n*\n* Level 2 Blas routine.\n*\n* -- Written on 22-October-1986.\n* Jack Dongarra, Argonne National Lab.\n* Jeremy Du Croz, Nag Central Office.\n* Sven Hammarling, Nag Central Office.\n* Richard Hanson, Sandia National Labs.\n* -- Modified for the absolute-value product, April 2006\n* Jason Riedy, UC Berkeley\n*\n* =====================================================================\n*\n\n");
36
- return Qnil;
37
- }
38
- if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
39
- printf("%s\n", "USAGE:\n y = NumRu::Lapack.sla_syamv( uplo, alpha, a, x, incx, beta, y, incy, [:usage => usage, :help => help])\n");
40
- return Qnil;
41
- }
42
- } else
43
- rblapack_options = Qnil;
44
- if (argc != 8 && argc != 8)
45
- rb_raise(rb_eArgError,"wrong number of arguments (%d for 8)", argc);
46
- rblapack_uplo = argv[0];
47
- rblapack_alpha = argv[1];
48
- rblapack_a = argv[2];
49
- rblapack_x = argv[3];
50
- rblapack_incx = argv[4];
51
- rblapack_beta = argv[5];
52
- rblapack_y = argv[6];
53
- rblapack_incy = argv[7];
54
- if (argc == 8) {
55
- } else if (rblapack_options != Qnil) {
56
- } else {
57
- }
58
-
59
- uplo = NUM2INT(rblapack_uplo);
60
- if (!NA_IsNArray(rblapack_a))
61
- rb_raise(rb_eArgError, "a (3th argument) must be NArray");
62
- if (NA_RANK(rblapack_a) != 2)
63
- rb_raise(rb_eArgError, "rank of a (3th argument) must be %d", 2);
64
- lda = NA_SHAPE0(rblapack_a);
65
- if (lda != (MAX(1, n)))
66
- rb_raise(rb_eRuntimeError, "shape 0 of a must be %d", MAX(1, n));
67
- n = NA_SHAPE1(rblapack_a);
68
- if (NA_TYPE(rblapack_a) != NA_SFLOAT)
69
- rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
70
- a = NA_PTR_TYPE(rblapack_a, real*);
71
- incx = NUM2INT(rblapack_incx);
72
- incy = NUM2INT(rblapack_incy);
73
- alpha = (real)NUM2DBL(rblapack_alpha);
74
- beta = (real)NUM2DBL(rblapack_beta);
75
- lda = MAX(1, n);
76
- if (!NA_IsNArray(rblapack_x))
77
- rb_raise(rb_eArgError, "x (4th argument) must be NArray");
78
- if (NA_RANK(rblapack_x) != 1)
79
- rb_raise(rb_eArgError, "rank of x (4th argument) must be %d", 1);
80
- if (NA_SHAPE0(rblapack_x) != (1 + ( n - 1 )*abs( incx )))
81
- rb_raise(rb_eRuntimeError, "shape 0 of x must be %d", 1 + ( n - 1 )*abs( incx ));
82
- if (NA_TYPE(rblapack_x) != NA_SFLOAT)
83
- rblapack_x = na_change_type(rblapack_x, NA_SFLOAT);
84
- x = NA_PTR_TYPE(rblapack_x, real*);
85
- if (!NA_IsNArray(rblapack_y))
86
- rb_raise(rb_eArgError, "y (7th argument) must be NArray");
87
- if (NA_RANK(rblapack_y) != 1)
88
- rb_raise(rb_eArgError, "rank of y (7th argument) must be %d", 1);
89
- if (NA_SHAPE0(rblapack_y) != (1 + ( n - 1 )*abs( incy )))
90
- rb_raise(rb_eRuntimeError, "shape 0 of y must be %d", 1 + ( n - 1 )*abs( incy ));
91
- if (NA_TYPE(rblapack_y) != NA_SFLOAT)
92
- rblapack_y = na_change_type(rblapack_y, NA_SFLOAT);
93
- y = NA_PTR_TYPE(rblapack_y, real*);
94
- {
95
- int shape[1];
96
- shape[0] = 1 + ( n - 1 )*abs( incy );
97
- rblapack_y_out__ = na_make_object(NA_SFLOAT, 1, shape, cNArray);
98
- }
99
- y_out__ = NA_PTR_TYPE(rblapack_y_out__, real*);
100
- MEMCPY(y_out__, y, real, NA_TOTAL(rblapack_y));
101
- rblapack_y = rblapack_y_out__;
102
- y = y_out__;
103
-
104
- sla_syamv_(&uplo, &n, &alpha, a, &lda, x, &incx, &beta, y, &incy);
105
-
106
- return rblapack_y;
107
- }
108
-
109
- void
110
- init_lapack_sla_syamv(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
111
- sHelp = sH;
112
- sUsage = sU;
113
- rblapack_ZERO = zero;
114
-
115
- rb_define_module_function(mLapack, "sla_syamv", rblapack_sla_syamv, -1);
116
- }
@@ -1,134 +0,0 @@
1
- #include "rb_lapack.h"
2
-
3
- extern real sla_syrcond_(char* uplo, integer* n, real* a, integer* lda, real* af, integer* ldaf, integer* ipiv, integer* cmode, real* c, integer* info, real* work, integer* iwork);
4
-
5
-
6
- static VALUE
7
- rblapack_sla_syrcond(int argc, VALUE *argv, VALUE self){
8
- VALUE rblapack_uplo;
9
- char uplo;
10
- VALUE rblapack_a;
11
- real *a;
12
- VALUE rblapack_af;
13
- real *af;
14
- VALUE rblapack_ipiv;
15
- integer *ipiv;
16
- VALUE rblapack_cmode;
17
- integer cmode;
18
- VALUE rblapack_c;
19
- real *c;
20
- VALUE rblapack_work;
21
- real *work;
22
- VALUE rblapack_iwork;
23
- integer *iwork;
24
- VALUE rblapack_info;
25
- integer info;
26
- VALUE rblapack___out__;
27
- real __out__;
28
-
29
- integer lda;
30
- integer n;
31
- integer ldaf;
32
-
33
- VALUE rblapack_options;
34
- if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
35
- argc--;
36
- rblapack_options = argv[argc];
37
- if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
38
- printf("%s\n", "USAGE:\n info, __out__ = NumRu::Lapack.sla_syrcond( uplo, a, af, ipiv, cmode, c, work, iwork, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n REAL FUNCTION SLA_SYRCOND( UPLO, N, A, LDA, AF, LDAF, IPIV, CMODE, C, INFO, WORK, IWORK )\n\n* Purpose\n* =======\n*\n* SLA_SYRCOND estimates the Skeel condition number of op(A) * op2(C)\n* where op2 is determined by CMODE as follows\n* CMODE = 1 op2(C) = C\n* CMODE = 0 op2(C) = I\n* CMODE = -1 op2(C) = inv(C)\n* The Skeel condition number cond(A) = norminf( |inv(A)||A| )\n* is computed by computing scaling factors R such that\n* diag(R)*A*op2(C) is row equilibrated and computing the standard\n* infinity-norm condition number.\n*\n\n* Arguments\n* ==========\n*\n* UPLO (input) CHARACTER*1\n* = 'U': Upper triangle of A is stored;\n* = 'L': Lower triangle of A is stored.\n*\n* N (input) INTEGER\n* The number of linear equations, i.e., the order of the\n* matrix A. N >= 0.\n*\n* A (input) REAL array, dimension (LDA,N)\n* On entry, the N-by-N matrix A.\n*\n* LDA (input) INTEGER\n* The leading dimension of the array A. LDA >= max(1,N).\n*\n* AF (input) REAL array, dimension (LDAF,N)\n* The block diagonal matrix D and the multipliers used to\n* obtain the factor U or L as computed by SSYTRF.\n*\n* LDAF (input) INTEGER\n* The leading dimension of the array AF. LDAF >= max(1,N).\n*\n* IPIV (input) INTEGER array, dimension (N)\n* Details of the interchanges and the block structure of D\n* as determined by SSYTRF.\n*\n* CMODE (input) INTEGER\n* Determines op2(C) in the formula op(A) * op2(C) as follows:\n* CMODE = 1 op2(C) = C\n* CMODE = 0 op2(C) = I\n* CMODE = -1 op2(C) = inv(C)\n*\n* C (input) REAL array, dimension (N)\n* The vector C in the formula op(A) * op2(C).\n*\n* INFO (output) INTEGER\n* = 0: Successful exit.\n* i > 0: The ith argument is invalid.\n*\n* WORK (input) REAL array, dimension (3*N).\n* Workspace.\n*\n* IWORK (input) INTEGER array, dimension (N).\n* Workspace.\n*\n\n* =====================================================================\n*\n* .. Local Scalars ..\n CHARACTER NORMIN\n INTEGER KASE, I, J\n REAL AINVNM, SMLNUM, TMP\n LOGICAL UP\n* ..\n* .. Local Arrays ..\n INTEGER ISAVE( 3 )\n* ..\n* .. External Functions ..\n LOGICAL LSAME\n INTEGER ISAMAX\n REAL SLAMCH\n EXTERNAL LSAME, ISAMAX, SLAMCH\n* ..\n* .. External Subroutines ..\n EXTERNAL SLACN2, SLATRS, SRSCL, XERBLA, SSYTRS\n* ..\n* .. Intrinsic Functions ..\n INTRINSIC ABS, MAX\n* ..\n\n");
39
- return Qnil;
40
- }
41
- if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
42
- printf("%s\n", "USAGE:\n info, __out__ = NumRu::Lapack.sla_syrcond( uplo, a, af, ipiv, cmode, c, work, iwork, [:usage => usage, :help => help])\n");
43
- return Qnil;
44
- }
45
- } else
46
- rblapack_options = Qnil;
47
- if (argc != 8 && argc != 8)
48
- rb_raise(rb_eArgError,"wrong number of arguments (%d for 8)", argc);
49
- rblapack_uplo = argv[0];
50
- rblapack_a = argv[1];
51
- rblapack_af = argv[2];
52
- rblapack_ipiv = argv[3];
53
- rblapack_cmode = argv[4];
54
- rblapack_c = argv[5];
55
- rblapack_work = argv[6];
56
- rblapack_iwork = argv[7];
57
- if (argc == 8) {
58
- } else if (rblapack_options != Qnil) {
59
- } else {
60
- }
61
-
62
- uplo = StringValueCStr(rblapack_uplo)[0];
63
- if (!NA_IsNArray(rblapack_af))
64
- rb_raise(rb_eArgError, "af (3th argument) must be NArray");
65
- if (NA_RANK(rblapack_af) != 2)
66
- rb_raise(rb_eArgError, "rank of af (3th argument) must be %d", 2);
67
- ldaf = NA_SHAPE0(rblapack_af);
68
- n = NA_SHAPE1(rblapack_af);
69
- if (NA_TYPE(rblapack_af) != NA_SFLOAT)
70
- rblapack_af = na_change_type(rblapack_af, NA_SFLOAT);
71
- af = NA_PTR_TYPE(rblapack_af, real*);
72
- cmode = NUM2INT(rblapack_cmode);
73
- if (!NA_IsNArray(rblapack_iwork))
74
- rb_raise(rb_eArgError, "iwork (8th argument) must be NArray");
75
- if (NA_RANK(rblapack_iwork) != 1)
76
- rb_raise(rb_eArgError, "rank of iwork (8th argument) must be %d", 1);
77
- if (NA_SHAPE0(rblapack_iwork) != n)
78
- rb_raise(rb_eRuntimeError, "shape 0 of iwork must be the same as shape 1 of af");
79
- if (NA_TYPE(rblapack_iwork) != NA_LINT)
80
- rblapack_iwork = na_change_type(rblapack_iwork, NA_LINT);
81
- iwork = NA_PTR_TYPE(rblapack_iwork, integer*);
82
- if (!NA_IsNArray(rblapack_a))
83
- rb_raise(rb_eArgError, "a (2th argument) must be NArray");
84
- if (NA_RANK(rblapack_a) != 2)
85
- rb_raise(rb_eArgError, "rank of a (2th argument) must be %d", 2);
86
- lda = NA_SHAPE0(rblapack_a);
87
- if (NA_SHAPE1(rblapack_a) != n)
88
- rb_raise(rb_eRuntimeError, "shape 1 of a must be the same as shape 1 of af");
89
- if (NA_TYPE(rblapack_a) != NA_SFLOAT)
90
- rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
91
- a = NA_PTR_TYPE(rblapack_a, real*);
92
- if (!NA_IsNArray(rblapack_c))
93
- rb_raise(rb_eArgError, "c (6th argument) must be NArray");
94
- if (NA_RANK(rblapack_c) != 1)
95
- rb_raise(rb_eArgError, "rank of c (6th argument) must be %d", 1);
96
- if (NA_SHAPE0(rblapack_c) != n)
97
- rb_raise(rb_eRuntimeError, "shape 0 of c must be the same as shape 1 of af");
98
- if (NA_TYPE(rblapack_c) != NA_SFLOAT)
99
- rblapack_c = na_change_type(rblapack_c, NA_SFLOAT);
100
- c = NA_PTR_TYPE(rblapack_c, real*);
101
- if (!NA_IsNArray(rblapack_ipiv))
102
- rb_raise(rb_eArgError, "ipiv (4th argument) must be NArray");
103
- if (NA_RANK(rblapack_ipiv) != 1)
104
- rb_raise(rb_eArgError, "rank of ipiv (4th argument) must be %d", 1);
105
- if (NA_SHAPE0(rblapack_ipiv) != n)
106
- rb_raise(rb_eRuntimeError, "shape 0 of ipiv must be the same as shape 1 of af");
107
- if (NA_TYPE(rblapack_ipiv) != NA_LINT)
108
- rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
109
- ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
110
- if (!NA_IsNArray(rblapack_work))
111
- rb_raise(rb_eArgError, "work (7th argument) must be NArray");
112
- if (NA_RANK(rblapack_work) != 1)
113
- rb_raise(rb_eArgError, "rank of work (7th argument) must be %d", 1);
114
- if (NA_SHAPE0(rblapack_work) != (3*n))
115
- rb_raise(rb_eRuntimeError, "shape 0 of work must be %d", 3*n);
116
- if (NA_TYPE(rblapack_work) != NA_SFLOAT)
117
- rblapack_work = na_change_type(rblapack_work, NA_SFLOAT);
118
- work = NA_PTR_TYPE(rblapack_work, real*);
119
-
120
- __out__ = sla_syrcond_(&uplo, &n, a, &lda, af, &ldaf, ipiv, &cmode, c, &info, work, iwork);
121
-
122
- rblapack_info = INT2NUM(info);
123
- rblapack___out__ = rb_float_new((double)__out__);
124
- return rb_ary_new3(2, rblapack_info, rblapack___out__);
125
- }
126
-
127
- void
128
- init_lapack_sla_syrcond(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
129
- sHelp = sH;
130
- sUsage = sU;
131
- rblapack_ZERO = zero;
132
-
133
- rb_define_module_function(mLapack, "sla_syrcond", rblapack_sla_syrcond, -1);
134
- }
@@ -1,283 +0,0 @@
1
- #include "rb_lapack.h"
2
-
3
- extern VOID sla_syrfsx_extended_(integer* prec_type, char* uplo, integer* n, integer* nrhs, real* a, integer* lda, real* af, integer* ldaf, integer* ipiv, logical* colequ, real* c, real* b, integer* ldb, real* y, integer* ldy, real* berr_out, integer* n_norms, real* err_bnds_norm, real* err_bnds_comp, real* res, real* ayb, real* dy, real* y_tail, real* rcond, integer* ithresh, real* rthresh, real* dz_ub, logical* ignore_cwise, integer* info);
4
-
5
-
6
- static VALUE
7
- rblapack_sla_syrfsx_extended(int argc, VALUE *argv, VALUE self){
8
- VALUE rblapack_prec_type;
9
- integer prec_type;
10
- VALUE rblapack_uplo;
11
- char uplo;
12
- VALUE rblapack_a;
13
- real *a;
14
- VALUE rblapack_af;
15
- real *af;
16
- VALUE rblapack_ipiv;
17
- integer *ipiv;
18
- VALUE rblapack_colequ;
19
- logical colequ;
20
- VALUE rblapack_c;
21
- real *c;
22
- VALUE rblapack_b;
23
- real *b;
24
- VALUE rblapack_y;
25
- real *y;
26
- VALUE rblapack_n_norms;
27
- integer n_norms;
28
- VALUE rblapack_err_bnds_norm;
29
- real *err_bnds_norm;
30
- VALUE rblapack_err_bnds_comp;
31
- real *err_bnds_comp;
32
- VALUE rblapack_res;
33
- real *res;
34
- VALUE rblapack_ayb;
35
- real *ayb;
36
- VALUE rblapack_dy;
37
- real *dy;
38
- VALUE rblapack_y_tail;
39
- real *y_tail;
40
- VALUE rblapack_rcond;
41
- real rcond;
42
- VALUE rblapack_ithresh;
43
- integer ithresh;
44
- VALUE rblapack_rthresh;
45
- real rthresh;
46
- VALUE rblapack_dz_ub;
47
- real dz_ub;
48
- VALUE rblapack_ignore_cwise;
49
- logical ignore_cwise;
50
- VALUE rblapack_berr_out;
51
- real *berr_out;
52
- VALUE rblapack_info;
53
- integer info;
54
- VALUE rblapack_y_out__;
55
- real *y_out__;
56
- VALUE rblapack_err_bnds_norm_out__;
57
- real *err_bnds_norm_out__;
58
- VALUE rblapack_err_bnds_comp_out__;
59
- real *err_bnds_comp_out__;
60
-
61
- integer lda;
62
- integer n;
63
- integer ldaf;
64
- integer ldb;
65
- integer nrhs;
66
- integer ldy;
67
- integer n_err_bnds;
68
-
69
- VALUE rblapack_options;
70
- if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
71
- argc--;
72
- rblapack_options = argv[argc];
73
- if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
74
- printf("%s\n", "USAGE:\n berr_out, info, y, err_bnds_norm, err_bnds_comp = NumRu::Lapack.sla_syrfsx_extended( prec_type, uplo, a, af, ipiv, colequ, c, b, y, n_norms, err_bnds_norm, err_bnds_comp, res, ayb, dy, y_tail, rcond, ithresh, rthresh, dz_ub, ignore_cwise, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n SUBROUTINE SLA_SYRFSX_EXTENDED( PREC_TYPE, UPLO, N, NRHS, A, LDA, AF, LDAF, IPIV, COLEQU, C, B, LDB, Y, LDY, BERR_OUT, N_NORMS, ERR_BNDS_NORM, ERR_BNDS_COMP, RES, AYB, DY, Y_TAIL, RCOND, ITHRESH, RTHRESH, DZ_UB, IGNORE_CWISE, INFO )\n\n* Purpose\n* =======\n* \n* SLA_SYRFSX_EXTENDED improves the computed solution to a system of\n* linear equations by performing extra-precise iterative refinement\n* and provides error bounds and backward error estimates for the solution.\n* This subroutine is called by SSYRFSX to perform iterative refinement.\n* In addition to normwise error bound, the code provides maximum\n* componentwise error bound if possible. See comments for ERR_BNDS_NORM\n* and ERR_BNDS_COMP for details of the error bounds. Note that this\n* subroutine is only resonsible for setting the second fields of\n* ERR_BNDS_NORM and ERR_BNDS_COMP.\n*\n\n* Arguments\n* =========\n*\n* PREC_TYPE (input) INTEGER\n* Specifies the intermediate precision to be used in refinement.\n* The value is defined by ILAPREC(P) where P is a CHARACTER and\n* P = 'S': Single\n* = 'D': Double\n* = 'I': Indigenous\n* = 'X', 'E': Extra\n*\n* UPLO (input) CHARACTER*1\n* = 'U': Upper triangle of A is stored;\n* = 'L': Lower triangle of A is stored.\n*\n* N (input) INTEGER\n* The number of linear equations, i.e., the order of the\n* matrix A. N >= 0.\n*\n* NRHS (input) INTEGER\n* The number of right-hand-sides, i.e., the number of columns of the\n* matrix B.\n*\n* A (input) REAL array, dimension (LDA,N)\n* On entry, the N-by-N matrix A.\n*\n* LDA (input) INTEGER\n* The leading dimension of the array A. LDA >= max(1,N).\n*\n* AF (input) REAL array, dimension (LDAF,N)\n* The block diagonal matrix D and the multipliers used to\n* obtain the factor U or L as computed by SSYTRF.\n*\n* LDAF (input) INTEGER\n* The leading dimension of the array AF. LDAF >= max(1,N).\n*\n* IPIV (input) INTEGER array, dimension (N)\n* Details of the interchanges and the block structure of D\n* as determined by SSYTRF.\n*\n* COLEQU (input) LOGICAL\n* If .TRUE. then column equilibration was done to A before calling\n* this routine. This is needed to compute the solution and error\n* bounds correctly.\n*\n* C (input) REAL array, dimension (N)\n* The column scale factors for A. If COLEQU = .FALSE., C\n* is not accessed. If C is input, each element of C should be a power\n* of the radix to ensure a reliable solution and error estimates.\n* Scaling by powers of the radix does not cause rounding errors unless\n* the result underflows or overflows. Rounding errors during scaling\n* lead to refining with a matrix that is not equivalent to the\n* input matrix, producing error estimates that may not be\n* reliable.\n*\n* B (input) REAL array, dimension (LDB,NRHS)\n* The right-hand-side matrix B.\n*\n* LDB (input) INTEGER\n* The leading dimension of the array B. LDB >= max(1,N).\n*\n* Y (input/output) REAL array, dimension (LDY,NRHS)\n* On entry, the solution matrix X, as computed by SSYTRS.\n* On exit, the improved solution matrix Y.\n*\n* LDY (input) INTEGER\n* The leading dimension of the array Y. LDY >= max(1,N).\n*\n* BERR_OUT (output) REAL array, dimension (NRHS)\n* On exit, BERR_OUT(j) contains the componentwise relative backward\n* error for right-hand-side j from the formula\n* max(i) ( abs(RES(i)) / ( abs(op(A_s))*abs(Y) + abs(B_s) )(i) )\n* where abs(Z) is the componentwise absolute value of the matrix\n* or vector Z. This is computed by SLA_LIN_BERR.\n*\n* N_NORMS (input) INTEGER\n* Determines which error bounds to return (see ERR_BNDS_NORM\n* and ERR_BNDS_COMP).\n* If N_NORMS >= 1 return normwise error bounds.\n* If N_NORMS >= 2 return componentwise error bounds.\n*\n* ERR_BNDS_NORM (input/output) REAL array, dimension (NRHS, N_ERR_BNDS)\n* For each right-hand side, this array contains information about\n* various error bounds and condition numbers corresponding to the\n* normwise relative error, which is defined as follows:\n*\n* Normwise relative error in the ith solution vector:\n* max_j (abs(XTRUE(j,i) - X(j,i)))\n* ------------------------------\n* max_j abs(X(j,i))\n*\n* The array is indexed by the type of error information as described\n* below. There currently are up to three pieces of information\n* returned.\n*\n* The first index in ERR_BNDS_NORM(i,:) corresponds to the ith\n* right-hand side.\n*\n* The second index in ERR_BNDS_NORM(:,err) contains the following\n* three fields:\n* err = 1 \"Trust/don't trust\" boolean. Trust the answer if the\n* reciprocal condition number is less than the threshold\n* sqrt(n) * slamch('Epsilon').\n*\n* err = 2 \"Guaranteed\" error bound: The estimated forward error,\n* almost certainly within a factor of 10 of the true error\n* so long as the next entry is greater than the threshold\n* sqrt(n) * slamch('Epsilon'). This error bound should only\n* be trusted if the previous boolean is true.\n*\n* err = 3 Reciprocal condition number: Estimated normwise\n* reciprocal condition number. Compared with the threshold\n* sqrt(n) * slamch('Epsilon') to determine if the error\n* estimate is \"guaranteed\". These reciprocal condition\n* numbers are 1 / (norm(Z^{-1},inf) * norm(Z,inf)) for some\n* appropriately scaled matrix Z.\n* Let Z = S*A, where S scales each row by a power of the\n* radix so all absolute row sums of Z are approximately 1.\n*\n* This subroutine is only responsible for setting the second field\n* above.\n* See Lapack Working Note 165 for further details and extra\n* cautions.\n*\n* ERR_BNDS_COMP (input/output) REAL array, dimension (NRHS, N_ERR_BNDS)\n* For each right-hand side, this array contains information about\n* various error bounds and condition numbers corresponding to the\n* componentwise relative error, which is defined as follows:\n*\n* Componentwise relative error in the ith solution vector:\n* abs(XTRUE(j,i) - X(j,i))\n* max_j ----------------------\n* abs(X(j,i))\n*\n* The array is indexed by the right-hand side i (on which the\n* componentwise relative error depends), and the type of error\n* information as described below. There currently are up to three\n* pieces of information returned for each right-hand side. If\n* componentwise accuracy is not requested (PARAMS(3) = 0.0), then\n* ERR_BNDS_COMP is not accessed. If N_ERR_BNDS .LT. 3, then at most\n* the first (:,N_ERR_BNDS) entries are returned.\n*\n* The first index in ERR_BNDS_COMP(i,:) corresponds to the ith\n* right-hand side.\n*\n* The second index in ERR_BNDS_COMP(:,err) contains the following\n* three fields:\n* err = 1 \"Trust/don't trust\" boolean. Trust the answer if the\n* reciprocal condition number is less than the threshold\n* sqrt(n) * slamch('Epsilon').\n*\n* err = 2 \"Guaranteed\" error bound: The estimated forward error,\n* almost certainly within a factor of 10 of the true error\n* so long as the next entry is greater than the threshold\n* sqrt(n) * slamch('Epsilon'). This error bound should only\n* be trusted if the previous boolean is true.\n*\n* err = 3 Reciprocal condition number: Estimated componentwise\n* reciprocal condition number. Compared with the threshold\n* sqrt(n) * slamch('Epsilon') to determine if the error\n* estimate is \"guaranteed\". These reciprocal condition\n* numbers are 1 / (norm(Z^{-1},inf) * norm(Z,inf)) for some\n* appropriately scaled matrix Z.\n* Let Z = S*(A*diag(x)), where x is the solution for the\n* current right-hand side and S scales each row of\n* A*diag(x) by a power of the radix so all absolute row\n* sums of Z are approximately 1.\n*\n* This subroutine is only responsible for setting the second field\n* above.\n* See Lapack Working Note 165 for further details and extra\n* cautions.\n*\n* RES (input) REAL array, dimension (N)\n* Workspace to hold the intermediate residual.\n*\n* AYB (input) REAL array, dimension (N)\n* Workspace. This can be the same workspace passed for Y_TAIL.\n*\n* DY (input) REAL array, dimension (N)\n* Workspace to hold the intermediate solution.\n*\n* Y_TAIL (input) REAL array, dimension (N)\n* Workspace to hold the trailing bits of the intermediate solution.\n*\n* RCOND (input) REAL\n* Reciprocal scaled condition number. This is an estimate of the\n* reciprocal Skeel condition number of the matrix A after\n* equilibration (if done). If this is less than the machine\n* precision (in particular, if it is zero), the matrix is singular\n* to working precision. Note that the error may still be small even\n* if this number is very small and the matrix appears ill-\n* conditioned.\n*\n* ITHRESH (input) INTEGER\n* The maximum number of residual computations allowed for\n* refinement. The default is 10. For 'aggressive' set to 100 to\n* permit convergence using approximate factorizations or\n* factorizations other than LU. If the factorization uses a\n* technique other than Gaussian elimination, the guarantees in\n* ERR_BNDS_NORM and ERR_BNDS_COMP may no longer be trustworthy.\n*\n* RTHRESH (input) REAL\n* Determines when to stop refinement if the error estimate stops\n* decreasing. Refinement will stop when the next solution no longer\n* satisfies norm(dx_{i+1}) < RTHRESH * norm(dx_i) where norm(Z) is\n* the infinity norm of Z. RTHRESH satisfies 0 < RTHRESH <= 1. The\n* default value is 0.5. For 'aggressive' set to 0.9 to permit\n* convergence on extremely ill-conditioned matrices. See LAWN 165\n* for more details.\n*\n* DZ_UB (input) REAL\n* Determines when to start considering componentwise convergence.\n* Componentwise convergence is only considered after each component\n* of the solution Y is stable, which we definte as the relative\n* change in each component being less than DZ_UB. The default value\n* is 0.25, requiring the first bit to be stable. See LAWN 165 for\n* more details.\n*\n* IGNORE_CWISE (input) LOGICAL\n* If .TRUE. then ignore componentwise convergence. Default value\n* is .FALSE..\n*\n* INFO (output) INTEGER\n* = 0: Successful exit.\n* < 0: if INFO = -i, the ith argument to SSYTRS had an illegal\n* value\n*\n\n* =====================================================================\n*\n* .. Local Scalars ..\n INTEGER UPLO2, CNT, I, J, X_STATE, Z_STATE\n REAL YK, DYK, YMIN, NORMY, NORMX, NORMDX, DXRAT,\n $ DZRAT, PREVNORMDX, PREV_DZ_Z, DXRATMAX,\n $ DZRATMAX, DX_X, DZ_Z, FINAL_DX_X, FINAL_DZ_Z,\n $ EPS, HUGEVAL, INCR_THRESH\n LOGICAL INCR_PREC\n* ..\n\n");
75
- return Qnil;
76
- }
77
- if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
78
- printf("%s\n", "USAGE:\n berr_out, info, y, err_bnds_norm, err_bnds_comp = NumRu::Lapack.sla_syrfsx_extended( prec_type, uplo, a, af, ipiv, colequ, c, b, y, n_norms, err_bnds_norm, err_bnds_comp, res, ayb, dy, y_tail, rcond, ithresh, rthresh, dz_ub, ignore_cwise, [:usage => usage, :help => help])\n");
79
- return Qnil;
80
- }
81
- } else
82
- rblapack_options = Qnil;
83
- if (argc != 21 && argc != 21)
84
- rb_raise(rb_eArgError,"wrong number of arguments (%d for 21)", argc);
85
- rblapack_prec_type = argv[0];
86
- rblapack_uplo = argv[1];
87
- rblapack_a = argv[2];
88
- rblapack_af = argv[3];
89
- rblapack_ipiv = argv[4];
90
- rblapack_colequ = argv[5];
91
- rblapack_c = argv[6];
92
- rblapack_b = argv[7];
93
- rblapack_y = argv[8];
94
- rblapack_n_norms = argv[9];
95
- rblapack_err_bnds_norm = argv[10];
96
- rblapack_err_bnds_comp = argv[11];
97
- rblapack_res = argv[12];
98
- rblapack_ayb = argv[13];
99
- rblapack_dy = argv[14];
100
- rblapack_y_tail = argv[15];
101
- rblapack_rcond = argv[16];
102
- rblapack_ithresh = argv[17];
103
- rblapack_rthresh = argv[18];
104
- rblapack_dz_ub = argv[19];
105
- rblapack_ignore_cwise = argv[20];
106
- if (argc == 21) {
107
- } else if (rblapack_options != Qnil) {
108
- } else {
109
- }
110
-
111
- prec_type = NUM2INT(rblapack_prec_type);
112
- if (!NA_IsNArray(rblapack_a))
113
- rb_raise(rb_eArgError, "a (3th argument) must be NArray");
114
- if (NA_RANK(rblapack_a) != 2)
115
- rb_raise(rb_eArgError, "rank of a (3th argument) must be %d", 2);
116
- lda = NA_SHAPE0(rblapack_a);
117
- n = NA_SHAPE1(rblapack_a);
118
- if (NA_TYPE(rblapack_a) != NA_SFLOAT)
119
- rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
120
- a = NA_PTR_TYPE(rblapack_a, real*);
121
- if (!NA_IsNArray(rblapack_ipiv))
122
- rb_raise(rb_eArgError, "ipiv (5th argument) must be NArray");
123
- if (NA_RANK(rblapack_ipiv) != 1)
124
- rb_raise(rb_eArgError, "rank of ipiv (5th argument) must be %d", 1);
125
- if (NA_SHAPE0(rblapack_ipiv) != n)
126
- rb_raise(rb_eRuntimeError, "shape 0 of ipiv must be the same as shape 1 of a");
127
- if (NA_TYPE(rblapack_ipiv) != NA_LINT)
128
- rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
129
- ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
130
- if (!NA_IsNArray(rblapack_c))
131
- rb_raise(rb_eArgError, "c (7th argument) must be NArray");
132
- if (NA_RANK(rblapack_c) != 1)
133
- rb_raise(rb_eArgError, "rank of c (7th argument) must be %d", 1);
134
- if (NA_SHAPE0(rblapack_c) != n)
135
- rb_raise(rb_eRuntimeError, "shape 0 of c must be the same as shape 1 of a");
136
- if (NA_TYPE(rblapack_c) != NA_SFLOAT)
137
- rblapack_c = na_change_type(rblapack_c, NA_SFLOAT);
138
- c = NA_PTR_TYPE(rblapack_c, real*);
139
- if (!NA_IsNArray(rblapack_y))
140
- rb_raise(rb_eArgError, "y (9th argument) must be NArray");
141
- if (NA_RANK(rblapack_y) != 2)
142
- rb_raise(rb_eArgError, "rank of y (9th argument) must be %d", 2);
143
- ldy = NA_SHAPE0(rblapack_y);
144
- nrhs = NA_SHAPE1(rblapack_y);
145
- if (NA_TYPE(rblapack_y) != NA_SFLOAT)
146
- rblapack_y = na_change_type(rblapack_y, NA_SFLOAT);
147
- y = NA_PTR_TYPE(rblapack_y, real*);
148
- if (!NA_IsNArray(rblapack_err_bnds_norm))
149
- rb_raise(rb_eArgError, "err_bnds_norm (11th argument) must be NArray");
150
- if (NA_RANK(rblapack_err_bnds_norm) != 2)
151
- rb_raise(rb_eArgError, "rank of err_bnds_norm (11th argument) must be %d", 2);
152
- if (NA_SHAPE0(rblapack_err_bnds_norm) != nrhs)
153
- rb_raise(rb_eRuntimeError, "shape 0 of err_bnds_norm must be the same as shape 1 of y");
154
- n_err_bnds = NA_SHAPE1(rblapack_err_bnds_norm);
155
- if (NA_TYPE(rblapack_err_bnds_norm) != NA_SFLOAT)
156
- rblapack_err_bnds_norm = na_change_type(rblapack_err_bnds_norm, NA_SFLOAT);
157
- err_bnds_norm = NA_PTR_TYPE(rblapack_err_bnds_norm, real*);
158
- if (!NA_IsNArray(rblapack_res))
159
- rb_raise(rb_eArgError, "res (13th argument) must be NArray");
160
- if (NA_RANK(rblapack_res) != 1)
161
- rb_raise(rb_eArgError, "rank of res (13th argument) must be %d", 1);
162
- if (NA_SHAPE0(rblapack_res) != n)
163
- rb_raise(rb_eRuntimeError, "shape 0 of res must be the same as shape 1 of a");
164
- if (NA_TYPE(rblapack_res) != NA_SFLOAT)
165
- rblapack_res = na_change_type(rblapack_res, NA_SFLOAT);
166
- res = NA_PTR_TYPE(rblapack_res, real*);
167
- if (!NA_IsNArray(rblapack_dy))
168
- rb_raise(rb_eArgError, "dy (15th argument) must be NArray");
169
- if (NA_RANK(rblapack_dy) != 1)
170
- rb_raise(rb_eArgError, "rank of dy (15th argument) must be %d", 1);
171
- if (NA_SHAPE0(rblapack_dy) != n)
172
- rb_raise(rb_eRuntimeError, "shape 0 of dy must be the same as shape 1 of a");
173
- if (NA_TYPE(rblapack_dy) != NA_SFLOAT)
174
- rblapack_dy = na_change_type(rblapack_dy, NA_SFLOAT);
175
- dy = NA_PTR_TYPE(rblapack_dy, real*);
176
- rcond = (real)NUM2DBL(rblapack_rcond);
177
- rthresh = (real)NUM2DBL(rblapack_rthresh);
178
- ignore_cwise = (rblapack_ignore_cwise == Qtrue);
179
- uplo = StringValueCStr(rblapack_uplo)[0];
180
- colequ = (rblapack_colequ == Qtrue);
181
- n_norms = NUM2INT(rblapack_n_norms);
182
- if (!NA_IsNArray(rblapack_ayb))
183
- rb_raise(rb_eArgError, "ayb (14th argument) must be NArray");
184
- if (NA_RANK(rblapack_ayb) != 1)
185
- rb_raise(rb_eArgError, "rank of ayb (14th argument) must be %d", 1);
186
- if (NA_SHAPE0(rblapack_ayb) != n)
187
- rb_raise(rb_eRuntimeError, "shape 0 of ayb must be the same as shape 1 of a");
188
- if (NA_TYPE(rblapack_ayb) != NA_SFLOAT)
189
- rblapack_ayb = na_change_type(rblapack_ayb, NA_SFLOAT);
190
- ayb = NA_PTR_TYPE(rblapack_ayb, real*);
191
- ithresh = NUM2INT(rblapack_ithresh);
192
- if (!NA_IsNArray(rblapack_af))
193
- rb_raise(rb_eArgError, "af (4th argument) must be NArray");
194
- if (NA_RANK(rblapack_af) != 2)
195
- rb_raise(rb_eArgError, "rank of af (4th argument) must be %d", 2);
196
- ldaf = NA_SHAPE0(rblapack_af);
197
- if (NA_SHAPE1(rblapack_af) != n)
198
- rb_raise(rb_eRuntimeError, "shape 1 of af must be the same as shape 1 of a");
199
- if (NA_TYPE(rblapack_af) != NA_SFLOAT)
200
- rblapack_af = na_change_type(rblapack_af, NA_SFLOAT);
201
- af = NA_PTR_TYPE(rblapack_af, real*);
202
- if (!NA_IsNArray(rblapack_err_bnds_comp))
203
- rb_raise(rb_eArgError, "err_bnds_comp (12th argument) must be NArray");
204
- if (NA_RANK(rblapack_err_bnds_comp) != 2)
205
- rb_raise(rb_eArgError, "rank of err_bnds_comp (12th argument) must be %d", 2);
206
- if (NA_SHAPE0(rblapack_err_bnds_comp) != nrhs)
207
- rb_raise(rb_eRuntimeError, "shape 0 of err_bnds_comp must be the same as shape 1 of y");
208
- if (NA_SHAPE1(rblapack_err_bnds_comp) != n_err_bnds)
209
- rb_raise(rb_eRuntimeError, "shape 1 of err_bnds_comp must be the same as shape 1 of err_bnds_norm");
210
- if (NA_TYPE(rblapack_err_bnds_comp) != NA_SFLOAT)
211
- rblapack_err_bnds_comp = na_change_type(rblapack_err_bnds_comp, NA_SFLOAT);
212
- err_bnds_comp = NA_PTR_TYPE(rblapack_err_bnds_comp, real*);
213
- dz_ub = (real)NUM2DBL(rblapack_dz_ub);
214
- if (!NA_IsNArray(rblapack_b))
215
- rb_raise(rb_eArgError, "b (8th argument) must be NArray");
216
- if (NA_RANK(rblapack_b) != 2)
217
- rb_raise(rb_eArgError, "rank of b (8th argument) must be %d", 2);
218
- ldb = NA_SHAPE0(rblapack_b);
219
- if (NA_SHAPE1(rblapack_b) != nrhs)
220
- rb_raise(rb_eRuntimeError, "shape 1 of b must be the same as shape 1 of y");
221
- if (NA_TYPE(rblapack_b) != NA_SFLOAT)
222
- rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
223
- b = NA_PTR_TYPE(rblapack_b, real*);
224
- if (!NA_IsNArray(rblapack_y_tail))
225
- rb_raise(rb_eArgError, "y_tail (16th argument) must be NArray");
226
- if (NA_RANK(rblapack_y_tail) != 1)
227
- rb_raise(rb_eArgError, "rank of y_tail (16th argument) must be %d", 1);
228
- if (NA_SHAPE0(rblapack_y_tail) != n)
229
- rb_raise(rb_eRuntimeError, "shape 0 of y_tail must be the same as shape 1 of a");
230
- if (NA_TYPE(rblapack_y_tail) != NA_SFLOAT)
231
- rblapack_y_tail = na_change_type(rblapack_y_tail, NA_SFLOAT);
232
- y_tail = NA_PTR_TYPE(rblapack_y_tail, real*);
233
- {
234
- int shape[1];
235
- shape[0] = nrhs;
236
- rblapack_berr_out = na_make_object(NA_SFLOAT, 1, shape, cNArray);
237
- }
238
- berr_out = NA_PTR_TYPE(rblapack_berr_out, real*);
239
- {
240
- int shape[2];
241
- shape[0] = ldy;
242
- shape[1] = nrhs;
243
- rblapack_y_out__ = na_make_object(NA_SFLOAT, 2, shape, cNArray);
244
- }
245
- y_out__ = NA_PTR_TYPE(rblapack_y_out__, real*);
246
- MEMCPY(y_out__, y, real, NA_TOTAL(rblapack_y));
247
- rblapack_y = rblapack_y_out__;
248
- y = y_out__;
249
- {
250
- int shape[2];
251
- shape[0] = nrhs;
252
- shape[1] = n_err_bnds;
253
- rblapack_err_bnds_norm_out__ = na_make_object(NA_SFLOAT, 2, shape, cNArray);
254
- }
255
- err_bnds_norm_out__ = NA_PTR_TYPE(rblapack_err_bnds_norm_out__, real*);
256
- MEMCPY(err_bnds_norm_out__, err_bnds_norm, real, NA_TOTAL(rblapack_err_bnds_norm));
257
- rblapack_err_bnds_norm = rblapack_err_bnds_norm_out__;
258
- err_bnds_norm = err_bnds_norm_out__;
259
- {
260
- int shape[2];
261
- shape[0] = nrhs;
262
- shape[1] = n_err_bnds;
263
- rblapack_err_bnds_comp_out__ = na_make_object(NA_SFLOAT, 2, shape, cNArray);
264
- }
265
- err_bnds_comp_out__ = NA_PTR_TYPE(rblapack_err_bnds_comp_out__, real*);
266
- MEMCPY(err_bnds_comp_out__, err_bnds_comp, real, NA_TOTAL(rblapack_err_bnds_comp));
267
- rblapack_err_bnds_comp = rblapack_err_bnds_comp_out__;
268
- err_bnds_comp = err_bnds_comp_out__;
269
-
270
- sla_syrfsx_extended_(&prec_type, &uplo, &n, &nrhs, a, &lda, af, &ldaf, ipiv, &colequ, c, b, &ldb, y, &ldy, berr_out, &n_norms, err_bnds_norm, err_bnds_comp, res, ayb, dy, y_tail, &rcond, &ithresh, &rthresh, &dz_ub, &ignore_cwise, &info);
271
-
272
- rblapack_info = INT2NUM(info);
273
- return rb_ary_new3(5, rblapack_berr_out, rblapack_info, rblapack_y, rblapack_err_bnds_norm, rblapack_err_bnds_comp);
274
- }
275
-
276
- void
277
- init_lapack_sla_syrfsx_extended(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
278
- sHelp = sH;
279
- sUsage = sU;
280
- rblapack_ZERO = zero;
281
-
282
- rb_define_module_function(mLapack, "sla_syrfsx_extended", rblapack_sla_syrfsx_extended, -1);
283
- }