ruby-lapack 1.7 → 1.7.1

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