vector.S 9.9 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473
  1. #include <asm/processor.h>
  2. #include <asm/ppc_asm.h>
  3. #include <asm/reg.h>
  4. #include <asm/asm-offsets.h>
  5. #include <asm/cputable.h>
  6. #include <asm/thread_info.h>
  7. #include <asm/page.h>
  8. #include <asm/ptrace.h>
  9. #ifdef CONFIG_PPC_TRANSACTIONAL_MEM
  10. /* void do_load_up_transact_altivec(struct thread_struct *thread)
  11. *
  12. * This is similar to load_up_altivec but for the transactional version of the
  13. * vector regs. It doesn't mess with the task MSR or valid flags.
  14. * Furthermore, VEC laziness is not supported with TM currently.
  15. */
  16. _GLOBAL(do_load_up_transact_altivec)
  17. mfmsr r6
  18. oris r5,r6,MSR_VEC@h
  19. MTMSRD(r5)
  20. isync
  21. li r4,1
  22. stw r4,THREAD_USED_VR(r3)
  23. li r10,THREAD_TRANSACT_VRSTATE+VRSTATE_VSCR
  24. lvx vr0,r10,r3
  25. mtvscr vr0
  26. addi r10,r3,THREAD_TRANSACT_VRSTATE
  27. REST_32VRS(0,r4,r10)
  28. /* Disable VEC again. */
  29. MTMSRD(r6)
  30. isync
  31. blr
  32. #endif
  33. /*
  34. * Load state from memory into VMX registers including VSCR.
  35. * Assumes the caller has enabled VMX in the MSR.
  36. */
  37. _GLOBAL(load_vr_state)
  38. li r4,VRSTATE_VSCR
  39. lvx vr0,r4,r3
  40. mtvscr vr0
  41. REST_32VRS(0,r4,r3)
  42. blr
  43. /*
  44. * Store VMX state into memory, including VSCR.
  45. * Assumes the caller has enabled VMX in the MSR.
  46. */
  47. _GLOBAL(store_vr_state)
  48. SAVE_32VRS(0, r4, r3)
  49. mfvscr vr0
  50. li r4, VRSTATE_VSCR
  51. stvx vr0, r4, r3
  52. blr
  53. /*
  54. * Disable VMX for the task which had it previously,
  55. * and save its vector registers in its thread_struct.
  56. * Enables the VMX for use in the kernel on return.
  57. * On SMP we know the VMX is free, since we give it up every
  58. * switch (ie, no lazy save of the vector registers).
  59. */
  60. _GLOBAL(load_up_altivec)
  61. mfmsr r5 /* grab the current MSR */
  62. oris r5,r5,MSR_VEC@h
  63. MTMSRD(r5) /* enable use of AltiVec now */
  64. isync
  65. /*
  66. * For SMP, we don't do lazy VMX switching because it just gets too
  67. * horrendously complex, especially when a task switches from one CPU
  68. * to another. Instead we call giveup_altvec in switch_to.
  69. * VRSAVE isn't dealt with here, that is done in the normal context
  70. * switch code. Note that we could rely on vrsave value to eventually
  71. * avoid saving all of the VREGs here...
  72. */
  73. #ifndef CONFIG_SMP
  74. LOAD_REG_ADDRBASE(r3, last_task_used_altivec)
  75. toreal(r3)
  76. PPC_LL r4,ADDROFF(last_task_used_altivec)(r3)
  77. PPC_LCMPI 0,r4,0
  78. beq 1f
  79. /* Save VMX state to last_task_used_altivec's THREAD struct */
  80. toreal(r4)
  81. addi r4,r4,THREAD
  82. addi r7,r4,THREAD_VRSTATE
  83. SAVE_32VRS(0,r5,r7)
  84. mfvscr vr0
  85. li r10,VRSTATE_VSCR
  86. stvx vr0,r10,r7
  87. /* Disable VMX for last_task_used_altivec */
  88. PPC_LL r5,PT_REGS(r4)
  89. toreal(r5)
  90. PPC_LL r4,_MSR-STACK_FRAME_OVERHEAD(r5)
  91. lis r10,MSR_VEC@h
  92. andc r4,r4,r10
  93. PPC_STL r4,_MSR-STACK_FRAME_OVERHEAD(r5)
  94. 1:
  95. #endif /* CONFIG_SMP */
  96. /* Hack: if we get an altivec unavailable trap with VRSAVE
  97. * set to all zeros, we assume this is a broken application
  98. * that fails to set it properly, and thus we switch it to
  99. * all 1's
  100. */
  101. mfspr r4,SPRN_VRSAVE
  102. cmpwi 0,r4,0
  103. bne+ 1f
  104. li r4,-1
  105. mtspr SPRN_VRSAVE,r4
  106. 1:
  107. /* enable use of VMX after return */
  108. #ifdef CONFIG_PPC32
  109. mfspr r5,SPRN_SPRG_THREAD /* current task's THREAD (phys) */
  110. oris r9,r9,MSR_VEC@h
  111. #else
  112. ld r4,PACACURRENT(r13)
  113. addi r5,r4,THREAD /* Get THREAD */
  114. oris r12,r12,MSR_VEC@h
  115. std r12,_MSR(r1)
  116. #endif
  117. addi r7,r5,THREAD_VRSTATE
  118. li r4,1
  119. li r10,VRSTATE_VSCR
  120. stw r4,THREAD_USED_VR(r5)
  121. lvx vr0,r10,r7
  122. mtvscr vr0
  123. REST_32VRS(0,r4,r7)
  124. #ifndef CONFIG_SMP
  125. /* Update last_task_used_altivec to 'current' */
  126. subi r4,r5,THREAD /* Back to 'current' */
  127. fromreal(r4)
  128. PPC_STL r4,ADDROFF(last_task_used_altivec)(r3)
  129. #endif /* CONFIG_SMP */
  130. /* restore registers and return */
  131. blr
  132. _GLOBAL(giveup_altivec_notask)
  133. mfmsr r3
  134. andis. r4,r3,MSR_VEC@h
  135. bnelr /* Already enabled? */
  136. oris r3,r3,MSR_VEC@h
  137. SYNC
  138. MTMSRD(r3) /* enable use of VMX now */
  139. isync
  140. blr
  141. /*
  142. * giveup_altivec(tsk)
  143. * Disable VMX for the task given as the argument,
  144. * and save the vector registers in its thread_struct.
  145. * Enables the VMX for use in the kernel on return.
  146. */
  147. _GLOBAL(giveup_altivec)
  148. mfmsr r5
  149. oris r5,r5,MSR_VEC@h
  150. SYNC
  151. MTMSRD(r5) /* enable use of VMX now */
  152. isync
  153. PPC_LCMPI 0,r3,0
  154. beqlr /* if no previous owner, done */
  155. addi r3,r3,THREAD /* want THREAD of task */
  156. PPC_LL r7,THREAD_VRSAVEAREA(r3)
  157. PPC_LL r5,PT_REGS(r3)
  158. PPC_LCMPI 0,r7,0
  159. bne 2f
  160. addi r7,r3,THREAD_VRSTATE
  161. 2: PPC_LCMPI 0,r5,0
  162. SAVE_32VRS(0,r4,r7)
  163. mfvscr vr0
  164. li r4,VRSTATE_VSCR
  165. stvx vr0,r4,r7
  166. beq 1f
  167. PPC_LL r4,_MSR-STACK_FRAME_OVERHEAD(r5)
  168. #ifdef CONFIG_VSX
  169. BEGIN_FTR_SECTION
  170. lis r3,(MSR_VEC|MSR_VSX)@h
  171. FTR_SECTION_ELSE
  172. lis r3,MSR_VEC@h
  173. ALT_FTR_SECTION_END_IFSET(CPU_FTR_VSX)
  174. #else
  175. lis r3,MSR_VEC@h
  176. #endif
  177. andc r4,r4,r3 /* disable FP for previous task */
  178. PPC_STL r4,_MSR-STACK_FRAME_OVERHEAD(r5)
  179. 1:
  180. #ifndef CONFIG_SMP
  181. li r5,0
  182. LOAD_REG_ADDRBASE(r4,last_task_used_altivec)
  183. PPC_STL r5,ADDROFF(last_task_used_altivec)(r4)
  184. #endif /* CONFIG_SMP */
  185. blr
  186. #ifdef CONFIG_VSX
  187. #ifdef CONFIG_PPC32
  188. #error This asm code isn't ready for 32-bit kernels
  189. #endif
  190. /*
  191. * load_up_vsx(unused, unused, tsk)
  192. * Disable VSX for the task which had it previously,
  193. * and save its vector registers in its thread_struct.
  194. * Reuse the fp and vsx saves, but first check to see if they have
  195. * been saved already.
  196. */
  197. _GLOBAL(load_up_vsx)
  198. /* Load FP and VSX registers if they haven't been done yet */
  199. andi. r5,r12,MSR_FP
  200. beql+ load_up_fpu /* skip if already loaded */
  201. andis. r5,r12,MSR_VEC@h
  202. beql+ load_up_altivec /* skip if already loaded */
  203. #ifndef CONFIG_SMP
  204. ld r3,last_task_used_vsx@got(r2)
  205. ld r4,0(r3)
  206. cmpdi 0,r4,0
  207. beq 1f
  208. /* Disable VSX for last_task_used_vsx */
  209. addi r4,r4,THREAD
  210. ld r5,PT_REGS(r4)
  211. ld r4,_MSR-STACK_FRAME_OVERHEAD(r5)
  212. lis r6,MSR_VSX@h
  213. andc r6,r4,r6
  214. std r6,_MSR-STACK_FRAME_OVERHEAD(r5)
  215. 1:
  216. #endif /* CONFIG_SMP */
  217. ld r4,PACACURRENT(r13)
  218. addi r4,r4,THREAD /* Get THREAD */
  219. li r6,1
  220. stw r6,THREAD_USED_VSR(r4) /* ... also set thread used vsr */
  221. /* enable use of VSX after return */
  222. oris r12,r12,MSR_VSX@h
  223. std r12,_MSR(r1)
  224. #ifndef CONFIG_SMP
  225. /* Update last_task_used_vsx to 'current' */
  226. ld r4,PACACURRENT(r13)
  227. std r4,0(r3)
  228. #endif /* CONFIG_SMP */
  229. b fast_exception_return
  230. /*
  231. * __giveup_vsx(tsk)
  232. * Disable VSX for the task given as the argument.
  233. * Does NOT save vsx registers.
  234. * Enables the VSX for use in the kernel on return.
  235. */
  236. _GLOBAL(__giveup_vsx)
  237. mfmsr r5
  238. oris r5,r5,MSR_VSX@h
  239. mtmsrd r5 /* enable use of VSX now */
  240. isync
  241. cmpdi 0,r3,0
  242. beqlr- /* if no previous owner, done */
  243. addi r3,r3,THREAD /* want THREAD of task */
  244. ld r5,PT_REGS(r3)
  245. cmpdi 0,r5,0
  246. beq 1f
  247. ld r4,_MSR-STACK_FRAME_OVERHEAD(r5)
  248. lis r3,MSR_VSX@h
  249. andc r4,r4,r3 /* disable VSX for previous task */
  250. std r4,_MSR-STACK_FRAME_OVERHEAD(r5)
  251. 1:
  252. #ifndef CONFIG_SMP
  253. li r5,0
  254. ld r4,last_task_used_vsx@got(r2)
  255. std r5,0(r4)
  256. #endif /* CONFIG_SMP */
  257. blr
  258. #endif /* CONFIG_VSX */
  259. /*
  260. * The routines below are in assembler so we can closely control the
  261. * usage of floating-point registers. These routines must be called
  262. * with preempt disabled.
  263. */
  264. #ifdef CONFIG_PPC32
  265. .data
  266. fpzero:
  267. .long 0
  268. fpone:
  269. .long 0x3f800000 /* 1.0 in single-precision FP */
  270. fphalf:
  271. .long 0x3f000000 /* 0.5 in single-precision FP */
  272. #define LDCONST(fr, name) \
  273. lis r11,name@ha; \
  274. lfs fr,name@l(r11)
  275. #else
  276. .section ".toc","aw"
  277. fpzero:
  278. .tc FD_0_0[TC],0
  279. fpone:
  280. .tc FD_3ff00000_0[TC],0x3ff0000000000000 /* 1.0 */
  281. fphalf:
  282. .tc FD_3fe00000_0[TC],0x3fe0000000000000 /* 0.5 */
  283. #define LDCONST(fr, name) \
  284. lfd fr,name@toc(r2)
  285. #endif
  286. .text
  287. /*
  288. * Internal routine to enable floating point and set FPSCR to 0.
  289. * Don't call it from C; it doesn't use the normal calling convention.
  290. */
  291. fpenable:
  292. #ifdef CONFIG_PPC32
  293. stwu r1,-64(r1)
  294. #else
  295. stdu r1,-64(r1)
  296. #endif
  297. mfmsr r10
  298. ori r11,r10,MSR_FP
  299. mtmsr r11
  300. isync
  301. stfd fr0,24(r1)
  302. stfd fr1,16(r1)
  303. stfd fr31,8(r1)
  304. LDCONST(fr1, fpzero)
  305. mffs fr31
  306. MTFSF_L(fr1)
  307. blr
  308. fpdisable:
  309. mtlr r12
  310. MTFSF_L(fr31)
  311. lfd fr31,8(r1)
  312. lfd fr1,16(r1)
  313. lfd fr0,24(r1)
  314. mtmsr r10
  315. isync
  316. addi r1,r1,64
  317. blr
  318. /*
  319. * Vector add, floating point.
  320. */
  321. _GLOBAL(vaddfp)
  322. mflr r12
  323. bl fpenable
  324. li r0,4
  325. mtctr r0
  326. li r6,0
  327. 1: lfsx fr0,r4,r6
  328. lfsx fr1,r5,r6
  329. fadds fr0,fr0,fr1
  330. stfsx fr0,r3,r6
  331. addi r6,r6,4
  332. bdnz 1b
  333. b fpdisable
  334. /*
  335. * Vector subtract, floating point.
  336. */
  337. _GLOBAL(vsubfp)
  338. mflr r12
  339. bl fpenable
  340. li r0,4
  341. mtctr r0
  342. li r6,0
  343. 1: lfsx fr0,r4,r6
  344. lfsx fr1,r5,r6
  345. fsubs fr0,fr0,fr1
  346. stfsx fr0,r3,r6
  347. addi r6,r6,4
  348. bdnz 1b
  349. b fpdisable
  350. /*
  351. * Vector multiply and add, floating point.
  352. */
  353. _GLOBAL(vmaddfp)
  354. mflr r12
  355. bl fpenable
  356. stfd fr2,32(r1)
  357. li r0,4
  358. mtctr r0
  359. li r7,0
  360. 1: lfsx fr0,r4,r7
  361. lfsx fr1,r5,r7
  362. lfsx fr2,r6,r7
  363. fmadds fr0,fr0,fr2,fr1
  364. stfsx fr0,r3,r7
  365. addi r7,r7,4
  366. bdnz 1b
  367. lfd fr2,32(r1)
  368. b fpdisable
  369. /*
  370. * Vector negative multiply and subtract, floating point.
  371. */
  372. _GLOBAL(vnmsubfp)
  373. mflr r12
  374. bl fpenable
  375. stfd fr2,32(r1)
  376. li r0,4
  377. mtctr r0
  378. li r7,0
  379. 1: lfsx fr0,r4,r7
  380. lfsx fr1,r5,r7
  381. lfsx fr2,r6,r7
  382. fnmsubs fr0,fr0,fr2,fr1
  383. stfsx fr0,r3,r7
  384. addi r7,r7,4
  385. bdnz 1b
  386. lfd fr2,32(r1)
  387. b fpdisable
  388. /*
  389. * Vector reciprocal estimate. We just compute 1.0/x.
  390. * r3 -> destination, r4 -> source.
  391. */
  392. _GLOBAL(vrefp)
  393. mflr r12
  394. bl fpenable
  395. li r0,4
  396. LDCONST(fr1, fpone)
  397. mtctr r0
  398. li r6,0
  399. 1: lfsx fr0,r4,r6
  400. fdivs fr0,fr1,fr0
  401. stfsx fr0,r3,r6
  402. addi r6,r6,4
  403. bdnz 1b
  404. b fpdisable
  405. /*
  406. * Vector reciprocal square-root estimate, floating point.
  407. * We use the frsqrte instruction for the initial estimate followed
  408. * by 2 iterations of Newton-Raphson to get sufficient accuracy.
  409. * r3 -> destination, r4 -> source.
  410. */
  411. _GLOBAL(vrsqrtefp)
  412. mflr r12
  413. bl fpenable
  414. stfd fr2,32(r1)
  415. stfd fr3,40(r1)
  416. stfd fr4,48(r1)
  417. stfd fr5,56(r1)
  418. li r0,4
  419. LDCONST(fr4, fpone)
  420. LDCONST(fr5, fphalf)
  421. mtctr r0
  422. li r6,0
  423. 1: lfsx fr0,r4,r6
  424. frsqrte fr1,fr0 /* r = frsqrte(s) */
  425. fmuls fr3,fr1,fr0 /* r * s */
  426. fmuls fr2,fr1,fr5 /* r * 0.5 */
  427. fnmsubs fr3,fr1,fr3,fr4 /* 1 - s * r * r */
  428. fmadds fr1,fr2,fr3,fr1 /* r = r + 0.5 * r * (1 - s * r * r) */
  429. fmuls fr3,fr1,fr0 /* r * s */
  430. fmuls fr2,fr1,fr5 /* r * 0.5 */
  431. fnmsubs fr3,fr1,fr3,fr4 /* 1 - s * r * r */
  432. fmadds fr1,fr2,fr3,fr1 /* r = r + 0.5 * r * (1 - s * r * r) */
  433. stfsx fr1,r3,r6
  434. addi r6,r6,4
  435. bdnz 1b
  436. lfd fr5,56(r1)
  437. lfd fr4,48(r1)
  438. lfd fr3,40(r1)
  439. lfd fr2,32(r1)
  440. b fpdisable