OSchip
/
OpenBLAS
Not watched
1
0
Fork 0
Code
Releases 66 Wiki evaluate Activity
Issues 0 Pull Requests 0 Datasets Model Cloudbrain HPC
You can not select more than 25 topics Topics must start with a chinese character,a letter or number, can include dashes ('-') and can be up to 35 characters long.
Branch: develop
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from 'develop'
${ noResults }
OpenBLAS/kernel/alpha/iamax.S

441 lines
8.6 kB
Raw Permalink Blame History 

  1. /*********************************************************************/

  2. /* Copyright 2009, 2010 The University of Texas at Austin.           */

  3. /* All rights reserved.                                              */

  4. /*                                                                   */

  5. /* Redistribution and use in source and binary forms, with or        */

  6. /* without modification, are permitted provided that the following   */

  7. /* conditions are met:                                               */

  8. /*                                                                   */

  9. /*   1. Redistributions of source code must retain the above         */

  10. /*      copyright notice, this list of conditions and the following  */

  11. /*      disclaimer.                                                  */

  12. /*                                                                   */

  13. /*   2. Redistributions in binary form must reproduce the above      */

  14. /*      copyright notice, this list of conditions and the following  */

  15. /*      disclaimer in the documentation and/or other materials       */

  16. /*      provided with the distribution.                              */

  17. /*                                                                   */

  18. /*    THIS  SOFTWARE IS PROVIDED  BY THE  UNIVERSITY OF  TEXAS AT    */

  19. /*    AUSTIN  ``AS IS''  AND ANY  EXPRESS OR  IMPLIED WARRANTIES,    */

  20. /*    INCLUDING, BUT  NOT LIMITED  TO, THE IMPLIED  WARRANTIES OF    */

  21. /*    MERCHANTABILITY  AND FITNESS FOR  A PARTICULAR  PURPOSE ARE    */

  22. /*    DISCLAIMED.  IN  NO EVENT SHALL THE UNIVERSITY  OF TEXAS AT    */

  23. /*    AUSTIN OR CONTRIBUTORS BE  LIABLE FOR ANY DIRECT, INDIRECT,    */

  24. /*    INCIDENTAL,  SPECIAL, EXEMPLARY,  OR  CONSEQUENTIAL DAMAGES    */

  25. /*    (INCLUDING, BUT  NOT LIMITED TO,  PROCUREMENT OF SUBSTITUTE    */

  26. /*    GOODS  OR  SERVICES; LOSS  OF  USE,  DATA,  OR PROFITS;  OR    */

  27. /*    BUSINESS INTERRUPTION) HOWEVER CAUSED  AND ON ANY THEORY OF    */

  28. /*    LIABILITY, WHETHER  IN CONTRACT, STRICT  LIABILITY, OR TORT    */

  29. /*    (INCLUDING NEGLIGENCE OR OTHERWISE)  ARISING IN ANY WAY OUT    */

  30. /*    OF  THE  USE OF  THIS  SOFTWARE,  EVEN  IF ADVISED  OF  THE    */

  31. /*    POSSIBILITY OF SUCH DAMAGE.                                    */

  32. /*                                                                   */

  33. /* The views and conclusions contained in the software and           */

  34. /* documentation are those of the authors and should not be          */

  35. /* interpreted as representing official policies, either expressed   */

  36. /* or implied, of The University of Texas at Austin.                 */

  37. /*********************************************************************/

  38. 

  39. #define ASSEMBLER

  40. #include "common.h"

  41. 

  42. 

  43. #define N	$16

  44. #define X	$17

  45. #define INCX	$18

  46. #define XX	$19

  47. 

  48. #ifndef USE_MIN

  49. #define CMPLT(a, b) cmptlt a, b

  50. #else

  51. #define CMPLT(a, b) cmptlt b, a

  52. #endif

  53. 

  54. #define STACKSIZE 6 * 8

  55. 

  56. 	PROLOGUE

  57. 	PROFCODE

  58. 	.frame	$sp, STACKSIZE, $26, 0

  59. 

  60. #ifdef F_INTERFACE

  61. 	ldl	N,     0(N)		# n

  62. 	ldl	INCX,  0(INCX)		# incx

  63. #endif

  64. 	lda	$sp, -STACKSIZE($sp)

  65. 	mov	X, XX

  66. 	.align 4

  67. 

  68. 	stt	$f2,   0($sp)

  69. 	fclr	$f16

  70. 	cmplt	$31, N,    $2

  71. 	unop

  72. 

  73. 	stt	$f3,   8($sp)

  74. 	fclr	$f17

  75. 	cmplt	$31, INCX, $3

  76. 	unop

  77. 

  78. 	stt	$f4,  16($sp)

  79. 	fclr	$f18

  80. 	SXADDQ	INCX, $31, INCX

  81. 	unop

  82. 

  83. 	stt	$f5,  24($sp)

  84. 	fclr	$f19

  85. 	and	$2,  $3,  $2

  86. 	clr	$0

  87. 

  88. 	stt	$f6,  32($sp)

  89. 	fclr	$f0

  90. 	sra	N, 3, $1

  91. 	beq	$2,  $End		# if (n <= 0) or (incx <= 0) return

  92. 	.align 4

  93. 

  94. 	LD	$f20,  0 * SIZE(X)

  95. 	unop

  96. 	fabs	$f20, $f0

  97. 	ble	$1,  $L15

  98. 	.align 4

  99. 

  100. 	fabs	$f20, $f1

  101. 	unop

  102. 	addq	X, INCX, X

  103. 	unop

  104. 

  105. 	LD	$f21,  0 * SIZE(X)

  106. 	fabs	$f20, $f2

  107. 	addq	X, INCX, X

  108. 	unop

  109. 

  110. 	LD	$f22,  0 * SIZE(X)

  111. 	fabs	$f20, $f3

  112. 	addq	X, INCX, X

  113. 	unop

  114. 

  115. 	LD	$f23,  0 * SIZE(X)

  116. 	fabs	$f20, $f4

  117. 	addq	X, INCX, X

  118. 	unop

  119. 

  120. 	LD	$f24,  0 * SIZE(X)

  121. 	addq	X, INCX, X

  122. 	fabs	$f20, $f5

  123. 	unop

  124. 

  125. 	LD	$f25,  0 * SIZE(X)

  126. 	fabs	$f20, $f6

  127. 	addq	X, INCX, X

  128. 	unop

  129. 

  130. 	LD	$f26,  0 * SIZE(X)

  131. 	fabs	$f20, $f28

  132. 	addq	X, INCX, X

  133. 	lda	$1,  -1($1)

  134. 

  135. 	LD	$f27,  0 * SIZE(X)

  136. 	unop

  137. 	addq	X, INCX, X

  138. 	ble	$1, $L13

  139. 	.align 4

  140. 

  141. $L12:

  142. 	fcmovne	$f16, $f12, $f4

  143. 	unop

  144. 	fabs	$f20, $f29

  145. 	ldl	$31, 56 * SIZE(X)

  146. 

  147. 	fcmovne	$f17, $f13, $f5

  148. 	LD	$f20,  0 * SIZE(X)

  149. 	fabs	$f21, $f30

  150. 	addq	X, INCX, X

  151. 

  152. 	fcmovne	$f18, $f14, $f6

  153. 	LD	$f21,  0 * SIZE(X)

  154. 	fabs	$f22, $f10

  155. 	addq	X, INCX, X

  156. 

  157. 	fcmovne	$f19, $f15, $f28

  158. 	LD	$f22,  0 * SIZE(X)

  159. 	fabs	$f23, $f11

  160. 	addq	X, INCX, X

  161. 

  162. 	fabs	$f24, $f12

  163. 	LD	$f23,  0 * SIZE(X)

  164. 	CMPLT($f0,  $f29),  $f16

  165. 	addq	X, INCX, X

  166. 

  167. 	fabs	$f25, $f13

  168. 	LD	$f24,  0 * SIZE(X)

  169. 	CMPLT($f1,  $f30),  $f17

  170. 	addq	X, INCX, X

  171. 

  172. 	fabs	$f26, $f14

  173. 	LD	$f25,  0 * SIZE(X)

  174. 	CMPLT($f2,  $f10), $f18

  175. 	addq	X, INCX, X

  176. 

  177. 	fabs	$f27, $f15

  178. 	LD	$f26,  0 * SIZE(X)

  179. 	CMPLT($f3,  $f11), $f19

  180. 	addq	X, INCX, X

  181. 

  182. 	fcmovne	$f16, $f29,  $f0

  183. 	LD	$f27,  0 * SIZE(X)

  184. 	CMPLT($f4,  $f12), $f16

  185. 	addq	X, INCX, X

  186. 

  187. 	fcmovne	$f17, $f30,  $f1

  188. 	unop

  189. 	CMPLT($f5,  $f13), $f17

  190. 	lda	$1,   -1($1)		# i --

  191. 

  192. 	fcmovne	$f18, $f10, $f2

  193. 	unop

  194. 	CMPLT($f6,  $f14), $f18

  195. 	unop

  196. 

  197. 	fcmovne	$f19, $f11, $f3

  198. 	unop

  199. 	CMPLT($f28,  $f15), $f19

  200. 	bgt	$1,$L12

  201. 	.align 4

  202. 

  203. $L13:

  204. 	fcmovne	$f16, $f12, $f4

  205. 	fabs	$f20, $f29

  206. 	fcmovne	$f17, $f13, $f5

  207. 	fabs	$f21, $f30

  208. 

  209. 	fcmovne	$f18, $f14, $f6

  210. 	fabs	$f22, $f10

  211. 	fcmovne	$f19, $f15, $f28

  212. 	fabs	$f23, $f11

  213. 

  214. 	fabs	$f24, $f12

  215. 	CMPLT($f0,  $f29),  $f16

  216. 	fabs	$f25, $f13

  217. 	CMPLT($f1,  $f30),  $f17

  218. 

  219. 	fabs	$f26, $f14

  220. 	CMPLT($f2,  $f10), $f18

  221. 	fabs	$f27, $f15

  222. 	CMPLT($f3,  $f11), $f19

  223. 

  224. 	fcmovne	$f16, $f29,  $f0

  225. 	CMPLT($f4,  $f12), $f16

  226. 	fcmovne	$f17, $f30,  $f1

  227. 	CMPLT($f5,  $f13), $f17

  228. 

  229. 	fcmovne	$f18, $f10, $f2

  230. 	CMPLT($f6,  $f14), $f18

  231. 	fcmovne	$f19, $f11, $f3

  232. 	CMPLT($f28,  $f15), $f19

  233. 

  234. 	fcmovne	$f16, $f12, $f4

  235. 	CMPLT($f0,  $f1), $f16

  236. 	fcmovne	$f17, $f13, $f5

  237. 	CMPLT($f2,  $f3), $f17

  238. 

  239. 	fcmovne	$f18, $f14, $f6

  240. 	CMPLT($f4,  $f5), $f18

  241. 	fcmovne	$f19, $f15, $f28

  242. 	CMPLT($f6,  $f28), $f19

  243. 

  244. 	fcmovne	$f16, $f1, $f0

  245. 	fcmovne	$f17, $f3, $f2

  246. 	fcmovne	$f18, $f5, $f4

  247. 	fcmovne	$f19, $f28, $f6

  248. 

  249. 	CMPLT($f0,  $f2), $f16

  250. 	CMPLT($f4,  $f6), $f17

  251. 

  252. 	fcmovne	$f16, $f2, $f0

  253. 	fcmovne	$f17, $f6, $f4

  254. 

  255. 	CMPLT($f0,  $f4), $f16

  256. 	fcmovne	$f16, $f4, $f0

  257. 	.align 4

  258. 

  259. $L15:

  260. 	and	N, 7, $1

  261. 	unop

  262. 	unop

  263. 	ble	$1,  $L20

  264. 	.align 4

  265. 

  266. $L16:

  267. 	LD	$f20,  0 * SIZE(X)

  268. 	addq	X, INCX, X

  269. 

  270. 	fabs	$f20, $f29

  271. 	CMPLT($f0,  $f29), $f16

  272. 	fcmovne	$f16, $f29, $f0

  273. 

  274. 	lda	$1,   -1($1)		# i --

  275. 	bgt	$1, $L16

  276. 	.align 4

  277. 

  278. $L20:

  279. 	sra	N, 3, $1

  280. 	ble	$1,  $L40

  281. 	.align 4

  282. 

  283. 	LD	$f10,  0 * SIZE(XX)

  284. 	addq	XX, INCX, XX

  285. 	LD	$f11,  0 * SIZE(XX)

  286. 	addq	XX, INCX, XX

  287. 

  288. 	LD	$f12,  0 * SIZE(XX)

  289. 	addq	XX, INCX, XX

  290. 	LD	$f13,  0 * SIZE(XX)

  291. 	addq	XX, INCX, XX

  292. 

  293. 	LD	$f14,  0 * SIZE(XX)

  294. 	addq	XX, INCX, XX

  295. 	LD	$f15,  0 * SIZE(XX)

  296. 	addq	XX, INCX, XX

  297. 

  298. 	LD	$f16,  0 * SIZE(XX)

  299. 	addq	XX, INCX, XX

  300. 	LD	$f17,  0 * SIZE(XX)

  301. 	addq	XX, INCX, XX

  302. 

  303. 	fabs	$f10, $f18

  304. 	fabs	$f11, $f19

  305. 	fabs	$f12, $f20

  306. 	fabs	$f13, $f21

  307. 

  308. 	lda	$1,  -1($1)

  309. 	ble	$1, $L23

  310. 	.align 4

  311. 

  312. $L22:

  313. 	LD	$f10,  0 * SIZE(XX)

  314. 	fabs	$f14, $f22

  315. 	addq	XX, INCX, XX

  316. 	cmpteq	$f0, $f18, $f2

  317. 

  318. 	LD	$f11,  0 * SIZE(XX)

  319. 	fabs	$f15, $f23

  320. 	addq	XX, INCX, XX

  321. 	cmpteq	$f0, $f19, $f3

  322. 

  323. 	LD	$f12,  0 * SIZE(XX)

  324. 	fabs	$f16, $f24

  325. 	addq	XX, INCX, XX

  326. 	cmpteq	$f0, $f20, $f4

  327. 

  328. 	LD	$f13,  0 * SIZE(XX)

  329. 	fabs	$f17, $f25

  330. 	addq	XX, INCX, XX

  331. 	cmpteq	$f0, $f21, $f5

  332. 

  333. 	LD	$f14,  0 * SIZE(XX)

  334. 	lda	$1,   -1($1)		# i --

  335. 	cmpteq	$f0, $f22, $f26

  336. 	addq	XX, INCX, XX

  337. 

  338. 	lda	$0,    1($0)

  339. 	fbne	$f2,  $End

  340. 

  341. 	LD	$f15,  0 * SIZE(XX)

  342. 	cmpteq	$f0, $f23, $f27

  343. 	lda	$0,    1($0)

  344. 	fbne	$f3,  $End

  345. 

  346. 	addq	XX, INCX, XX

  347. 	cmpteq	$f0, $f24, $f28

  348. 	lda	$0,    1($0)

  349. 	fbne	$f4,  $End

  350. 

  351. 	LD	$f16,  0 * SIZE(XX)

  352. 	cmpteq	$f0, $f25, $f29

  353. 	lda	$0,    1($0)

  354. 	fbne	$f5,  $End

  355. 

  356. 	addq	XX, INCX, XX

  357. 	lda	$0,    1($0)

  358. 	fabs	$f10, $f18

  359. 	fbne	$f26, $End

  360. 

  361. 	LD	$f17,  0 * SIZE(XX)

  362. 	lda	$0,    1($0)

  363. 	fabs	$f11, $f19

  364. 	fbne	$f27, $End

  365. 

  366. 	addq	XX, INCX, XX

  367. 	lda	$0,    1($0)

  368. 	fabs	$f12, $f20

  369. 	fbne	$f28, $End

  370. 

  371. 	lda	$0,    1($0)

  372. 	fabs	$f13, $f21

  373. 	fbne	$f29, $End

  374. 	bgt	$1,  $L22

  375. 	.align 4

  376. 

  377. $L23:

  378. 	fabs	$f14, $f22

  379. 	cmpteq	$f0, $f18, $f2

  380. 	fabs	$f15, $f23

  381. 	cmpteq	$f0, $f19, $f3

  382. 

  383. 	fabs	$f16, $f24

  384. 	cmpteq	$f0, $f20, $f4

  385. 	fabs	$f17, $f25

  386. 	cmpteq	$f0, $f21, $f5

  387. 

  388. 	cmpteq	$f0, $f22, $f26

  389. 	lda	$0,    1($0)

  390. 	unop

  391. 	fbne	$f2,  $End

  392. 

  393. 	cmpteq	$f0, $f23, $f27

  394. 	lda	$0,    1($0)

  395. 	unop

  396. 	fbne	$f3,  $End

  397. 

  398. 	cmpteq	$f0, $f24, $f28

  399. 	lda	$0,    1($0)

  400. 	unop

  401. 	fbne	$f4,  $End

  402. 

  403. 	cmpteq	$f0, $f25, $f29

  404. 	lda	$0,    1($0)

  405. 	unop

  406. 	fbne	$f5,  $End

  407. 

  408. 	lda	$0,    1($0)

  409. 	fbne	$f26, $End

  410. 	lda	$0,    1($0)

  411. 	fbne	$f27, $End

  412. 	lda	$0,    1($0)

  413. 	fbne	$f28, $End

  414. 	lda	$0,    1($0)

  415. 	fbne	$f29, $End

  416. 	.align 4

  417. 

  418. $L40:

  419. 	LD	$f20,  0 * SIZE(XX)

  420. 	addq	XX, INCX, XX

  421. 

  422. 	fabs	$f20, $f25

  423. 	cmpteq	$f0, $f25, $f29

  424. 

  425. 	lda	$0,    1($0)

  426. 	fbne	$f29,  $End

  427. 	br	$31, $L40

  428. 	.align 4

  429. 

  430. $End:

  431. 	ldt	$f2,   0($sp)

  432. 	ldt	$f3,   8($sp)

  433. 	ldt	$f4,  16($sp)

  434. 	ldt	$f5,  24($sp)

  435. 

  436. 	ldt	$f6,  32($sp)

  437. 	lda	$sp,  STACKSIZE($sp)

  438. 	ret

  439. 

  440. 	EPILOGUE