Merge remote-tracking branch 's5p/for-next'
[linux-2.6/next.git] / arch / mips / math-emu / cp1emu.c
blobdbf2f93a50911b5914a295a82b3179c36562013f
1 /*
2 * cp1emu.c: a MIPS coprocessor 1 (fpu) instruction emulator
4 * MIPS floating point support
5 * Copyright (C) 1994-2000 Algorithmics Ltd.
7 * Kevin D. Kissell, kevink@mips.com and Carsten Langgaard, carstenl@mips.com
8 * Copyright (C) 2000 MIPS Technologies, Inc.
10 * This program is free software; you can distribute it and/or modify it
11 * under the terms of the GNU General Public License (Version 2) as
12 * published by the Free Software Foundation.
14 * This program is distributed in the hope it will be useful, but WITHOUT
15 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
16 * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
17 * for more details.
19 * You should have received a copy of the GNU General Public License along
20 * with this program; if not, write to the Free Software Foundation, Inc.,
21 * 59 Temple Place - Suite 330, Boston MA 02111-1307, USA.
23 * A complete emulator for MIPS coprocessor 1 instructions. This is
24 * required for #float(switch) or #float(trap), where it catches all
25 * COP1 instructions via the "CoProcessor Unusable" exception.
27 * More surprisingly it is also required for #float(ieee), to help out
28 * the hardware fpu at the boundaries of the IEEE-754 representation
29 * (denormalised values, infinities, underflow, etc). It is made
30 * quite nasty because emulation of some non-COP1 instructions is
31 * required, e.g. in branch delay slots.
33 * Note if you know that you won't have an fpu, then you'll get much
34 * better performance by compiling with -msoft-float!
36 #include <linux/sched.h>
37 #include <linux/module.h>
38 #include <linux/debugfs.h>
39 #include <linux/perf_event.h>
41 #include <asm/inst.h>
42 #include <asm/bootinfo.h>
43 #include <asm/processor.h>
44 #include <asm/ptrace.h>
45 #include <asm/signal.h>
46 #include <asm/mipsregs.h>
47 #include <asm/fpu_emulator.h>
48 #include <asm/uaccess.h>
49 #include <asm/branch.h>
51 #include "ieee754.h"
53 /* Strap kernel emulator for full MIPS IV emulation */
55 #ifdef __mips
56 #undef __mips
57 #endif
58 #define __mips 4
60 /* Function which emulates a floating point instruction. */
62 static int fpu_emu(struct pt_regs *, struct mips_fpu_struct *,
63 mips_instruction);
65 #if __mips >= 4 && __mips != 32
66 static int fpux_emu(struct pt_regs *,
67 struct mips_fpu_struct *, mips_instruction, void *__user *);
68 #endif
70 /* Further private data for which no space exists in mips_fpu_struct */
72 #ifdef CONFIG_DEBUG_FS
73 DEFINE_PER_CPU(struct mips_fpu_emulator_stats, fpuemustats);
74 #endif
76 /* Control registers */
78 #define FPCREG_RID 0 /* $0 = revision id */
79 #define FPCREG_CSR 31 /* $31 = csr */
81 /* Determine rounding mode from the RM bits of the FCSR */
82 #define modeindex(v) ((v) & FPU_CSR_RM)
84 /* Convert Mips rounding mode (0..3) to IEEE library modes. */
85 static const unsigned char ieee_rm[4] = {
86 [FPU_CSR_RN] = IEEE754_RN,
87 [FPU_CSR_RZ] = IEEE754_RZ,
88 [FPU_CSR_RU] = IEEE754_RU,
89 [FPU_CSR_RD] = IEEE754_RD,
91 /* Convert IEEE library modes to Mips rounding mode (0..3). */
92 static const unsigned char mips_rm[4] = {
93 [IEEE754_RN] = FPU_CSR_RN,
94 [IEEE754_RZ] = FPU_CSR_RZ,
95 [IEEE754_RD] = FPU_CSR_RD,
96 [IEEE754_RU] = FPU_CSR_RU,
99 #if __mips >= 4
100 /* convert condition code register number to csr bit */
101 static const unsigned int fpucondbit[8] = {
102 FPU_CSR_COND0,
103 FPU_CSR_COND1,
104 FPU_CSR_COND2,
105 FPU_CSR_COND3,
106 FPU_CSR_COND4,
107 FPU_CSR_COND5,
108 FPU_CSR_COND6,
109 FPU_CSR_COND7
111 #endif
115 * Redundant with logic already in kernel/branch.c,
116 * embedded in compute_return_epc. At some point,
117 * a single subroutine should be used across both
118 * modules.
120 static int isBranchInstr(mips_instruction * i)
122 switch (MIPSInst_OPCODE(*i)) {
123 case spec_op:
124 switch (MIPSInst_FUNC(*i)) {
125 case jalr_op:
126 case jr_op:
127 return 1;
129 break;
131 case bcond_op:
132 switch (MIPSInst_RT(*i)) {
133 case bltz_op:
134 case bgez_op:
135 case bltzl_op:
136 case bgezl_op:
137 case bltzal_op:
138 case bgezal_op:
139 case bltzall_op:
140 case bgezall_op:
141 return 1;
143 break;
145 case j_op:
146 case jal_op:
147 case jalx_op:
148 case beq_op:
149 case bne_op:
150 case blez_op:
151 case bgtz_op:
152 case beql_op:
153 case bnel_op:
154 case blezl_op:
155 case bgtzl_op:
156 return 1;
158 case cop0_op:
159 case cop1_op:
160 case cop2_op:
161 case cop1x_op:
162 if (MIPSInst_RS(*i) == bc_op)
163 return 1;
164 break;
167 return 0;
171 * In the Linux kernel, we support selection of FPR format on the
172 * basis of the Status.FR bit. If an FPU is not present, the FR bit
173 * is hardwired to zero, which would imply a 32-bit FPU even for
174 * 64-bit CPUs. For 64-bit kernels with no FPU we use TIF_32BIT_REGS
175 * as a proxy for the FR bit so that a 64-bit FPU is emulated. In any
176 * case, for a 32-bit kernel which uses the O32 MIPS ABI, only the
177 * even FPRs are used (Status.FR = 0).
179 static inline int cop1_64bit(struct pt_regs *xcp)
181 if (cpu_has_fpu)
182 return xcp->cp0_status & ST0_FR;
183 #ifdef CONFIG_64BIT
184 return !test_thread_flag(TIF_32BIT_REGS);
185 #else
186 return 0;
187 #endif
190 #define SIFROMREG(si, x) ((si) = cop1_64bit(xcp) || !(x & 1) ? \
191 (int)ctx->fpr[x] : (int)(ctx->fpr[x & ~1] >> 32))
193 #define SITOREG(si, x) (ctx->fpr[x & ~(cop1_64bit(xcp) == 0)] = \
194 cop1_64bit(xcp) || !(x & 1) ? \
195 ctx->fpr[x & ~1] >> 32 << 32 | (u32)(si) : \
196 ctx->fpr[x & ~1] << 32 >> 32 | (u64)(si) << 32)
198 #define DIFROMREG(di, x) ((di) = ctx->fpr[x & ~(cop1_64bit(xcp) == 0)])
199 #define DITOREG(di, x) (ctx->fpr[x & ~(cop1_64bit(xcp) == 0)] = (di))
201 #define SPFROMREG(sp, x) SIFROMREG((sp).bits, x)
202 #define SPTOREG(sp, x) SITOREG((sp).bits, x)
203 #define DPFROMREG(dp, x) DIFROMREG((dp).bits, x)
204 #define DPTOREG(dp, x) DITOREG((dp).bits, x)
207 * Emulate the single floating point instruction pointed at by EPC.
208 * Two instructions if the instruction is in a branch delay slot.
211 static int cop1Emulate(struct pt_regs *xcp, struct mips_fpu_struct *ctx,
212 void *__user *fault_addr)
214 mips_instruction ir;
215 unsigned long emulpc, contpc;
216 unsigned int cond;
218 if (!access_ok(VERIFY_READ, xcp->cp0_epc, sizeof(mips_instruction))) {
219 MIPS_FPU_EMU_INC_STATS(errors);
220 *fault_addr = (mips_instruction __user *)xcp->cp0_epc;
221 return SIGBUS;
223 if (__get_user(ir, (mips_instruction __user *) xcp->cp0_epc)) {
224 MIPS_FPU_EMU_INC_STATS(errors);
225 *fault_addr = (mips_instruction __user *)xcp->cp0_epc;
226 return SIGSEGV;
229 /* XXX NEC Vr54xx bug workaround */
230 if ((xcp->cp0_cause & CAUSEF_BD) && !isBranchInstr(&ir))
231 xcp->cp0_cause &= ~CAUSEF_BD;
233 if (xcp->cp0_cause & CAUSEF_BD) {
235 * The instruction to be emulated is in a branch delay slot
236 * which means that we have to emulate the branch instruction
237 * BEFORE we do the cop1 instruction.
239 * This branch could be a COP1 branch, but in that case we
240 * would have had a trap for that instruction, and would not
241 * come through this route.
243 * Linux MIPS branch emulator operates on context, updating the
244 * cp0_epc.
246 emulpc = xcp->cp0_epc + 4; /* Snapshot emulation target */
248 if (__compute_return_epc(xcp)) {
249 #ifdef CP1DBG
250 printk("failed to emulate branch at %p\n",
251 (void *) (xcp->cp0_epc));
252 #endif
253 return SIGILL;
255 if (!access_ok(VERIFY_READ, emulpc, sizeof(mips_instruction))) {
256 MIPS_FPU_EMU_INC_STATS(errors);
257 *fault_addr = (mips_instruction __user *)emulpc;
258 return SIGBUS;
260 if (__get_user(ir, (mips_instruction __user *) emulpc)) {
261 MIPS_FPU_EMU_INC_STATS(errors);
262 *fault_addr = (mips_instruction __user *)emulpc;
263 return SIGSEGV;
265 /* __compute_return_epc() will have updated cp0_epc */
266 contpc = xcp->cp0_epc;
267 /* In order not to confuse ptrace() et al, tweak context */
268 xcp->cp0_epc = emulpc - 4;
269 } else {
270 emulpc = xcp->cp0_epc;
271 contpc = xcp->cp0_epc + 4;
274 emul:
275 perf_sw_event(PERF_COUNT_SW_EMULATION_FAULTS, 1, xcp, 0);
276 MIPS_FPU_EMU_INC_STATS(emulated);
277 switch (MIPSInst_OPCODE(ir)) {
278 case ldc1_op:{
279 u64 __user *va = (u64 __user *) (xcp->regs[MIPSInst_RS(ir)] +
280 MIPSInst_SIMM(ir));
281 u64 val;
283 MIPS_FPU_EMU_INC_STATS(loads);
285 if (!access_ok(VERIFY_READ, va, sizeof(u64))) {
286 MIPS_FPU_EMU_INC_STATS(errors);
287 *fault_addr = va;
288 return SIGBUS;
290 if (__get_user(val, va)) {
291 MIPS_FPU_EMU_INC_STATS(errors);
292 *fault_addr = va;
293 return SIGSEGV;
295 DITOREG(val, MIPSInst_RT(ir));
296 break;
299 case sdc1_op:{
300 u64 __user *va = (u64 __user *) (xcp->regs[MIPSInst_RS(ir)] +
301 MIPSInst_SIMM(ir));
302 u64 val;
304 MIPS_FPU_EMU_INC_STATS(stores);
305 DIFROMREG(val, MIPSInst_RT(ir));
306 if (!access_ok(VERIFY_WRITE, va, sizeof(u64))) {
307 MIPS_FPU_EMU_INC_STATS(errors);
308 *fault_addr = va;
309 return SIGBUS;
311 if (__put_user(val, va)) {
312 MIPS_FPU_EMU_INC_STATS(errors);
313 *fault_addr = va;
314 return SIGSEGV;
316 break;
319 case lwc1_op:{
320 u32 __user *va = (u32 __user *) (xcp->regs[MIPSInst_RS(ir)] +
321 MIPSInst_SIMM(ir));
322 u32 val;
324 MIPS_FPU_EMU_INC_STATS(loads);
325 if (!access_ok(VERIFY_READ, va, sizeof(u32))) {
326 MIPS_FPU_EMU_INC_STATS(errors);
327 *fault_addr = va;
328 return SIGBUS;
330 if (__get_user(val, va)) {
331 MIPS_FPU_EMU_INC_STATS(errors);
332 *fault_addr = va;
333 return SIGSEGV;
335 SITOREG(val, MIPSInst_RT(ir));
336 break;
339 case swc1_op:{
340 u32 __user *va = (u32 __user *) (xcp->regs[MIPSInst_RS(ir)] +
341 MIPSInst_SIMM(ir));
342 u32 val;
344 MIPS_FPU_EMU_INC_STATS(stores);
345 SIFROMREG(val, MIPSInst_RT(ir));
346 if (!access_ok(VERIFY_WRITE, va, sizeof(u32))) {
347 MIPS_FPU_EMU_INC_STATS(errors);
348 *fault_addr = va;
349 return SIGBUS;
351 if (__put_user(val, va)) {
352 MIPS_FPU_EMU_INC_STATS(errors);
353 *fault_addr = va;
354 return SIGSEGV;
356 break;
359 case cop1_op:
360 switch (MIPSInst_RS(ir)) {
362 #if defined(__mips64)
363 case dmfc_op:
364 /* copregister fs -> gpr[rt] */
365 if (MIPSInst_RT(ir) != 0) {
366 DIFROMREG(xcp->regs[MIPSInst_RT(ir)],
367 MIPSInst_RD(ir));
369 break;
371 case dmtc_op:
372 /* copregister fs <- rt */
373 DITOREG(xcp->regs[MIPSInst_RT(ir)], MIPSInst_RD(ir));
374 break;
375 #endif
377 case mfc_op:
378 /* copregister rd -> gpr[rt] */
379 if (MIPSInst_RT(ir) != 0) {
380 SIFROMREG(xcp->regs[MIPSInst_RT(ir)],
381 MIPSInst_RD(ir));
383 break;
385 case mtc_op:
386 /* copregister rd <- rt */
387 SITOREG(xcp->regs[MIPSInst_RT(ir)], MIPSInst_RD(ir));
388 break;
390 case cfc_op:{
391 /* cop control register rd -> gpr[rt] */
392 u32 value;
394 if (MIPSInst_RD(ir) == FPCREG_CSR) {
395 value = ctx->fcr31;
396 value = (value & ~FPU_CSR_RM) |
397 mips_rm[modeindex(value)];
398 #ifdef CSRTRACE
399 printk("%p gpr[%d]<-csr=%08x\n",
400 (void *) (xcp->cp0_epc),
401 MIPSInst_RT(ir), value);
402 #endif
404 else if (MIPSInst_RD(ir) == FPCREG_RID)
405 value = 0;
406 else
407 value = 0;
408 if (MIPSInst_RT(ir))
409 xcp->regs[MIPSInst_RT(ir)] = value;
410 break;
413 case ctc_op:{
414 /* copregister rd <- rt */
415 u32 value;
417 if (MIPSInst_RT(ir) == 0)
418 value = 0;
419 else
420 value = xcp->regs[MIPSInst_RT(ir)];
422 /* we only have one writable control reg
424 if (MIPSInst_RD(ir) == FPCREG_CSR) {
425 #ifdef CSRTRACE
426 printk("%p gpr[%d]->csr=%08x\n",
427 (void *) (xcp->cp0_epc),
428 MIPSInst_RT(ir), value);
429 #endif
432 * Don't write reserved bits,
433 * and convert to ieee library modes
435 ctx->fcr31 = (value &
436 ~(FPU_CSR_RSVD | FPU_CSR_RM)) |
437 ieee_rm[modeindex(value)];
439 if ((ctx->fcr31 >> 5) & ctx->fcr31 & FPU_CSR_ALL_E) {
440 return SIGFPE;
442 break;
445 case bc_op:{
446 int likely = 0;
448 if (xcp->cp0_cause & CAUSEF_BD)
449 return SIGILL;
451 #if __mips >= 4
452 cond = ctx->fcr31 & fpucondbit[MIPSInst_RT(ir) >> 2];
453 #else
454 cond = ctx->fcr31 & FPU_CSR_COND;
455 #endif
456 switch (MIPSInst_RT(ir) & 3) {
457 case bcfl_op:
458 likely = 1;
459 case bcf_op:
460 cond = !cond;
461 break;
462 case bctl_op:
463 likely = 1;
464 case bct_op:
465 break;
466 default:
467 /* thats an illegal instruction */
468 return SIGILL;
471 xcp->cp0_cause |= CAUSEF_BD;
472 if (cond) {
473 /* branch taken: emulate dslot
474 * instruction
476 xcp->cp0_epc += 4;
477 contpc = (xcp->cp0_epc +
478 (MIPSInst_SIMM(ir) << 2));
480 if (!access_ok(VERIFY_READ, xcp->cp0_epc,
481 sizeof(mips_instruction))) {
482 MIPS_FPU_EMU_INC_STATS(errors);
483 *fault_addr = (mips_instruction __user *)xcp->cp0_epc;
484 return SIGBUS;
486 if (__get_user(ir,
487 (mips_instruction __user *) xcp->cp0_epc)) {
488 MIPS_FPU_EMU_INC_STATS(errors);
489 *fault_addr = (mips_instruction __user *)xcp->cp0_epc;
490 return SIGSEGV;
493 switch (MIPSInst_OPCODE(ir)) {
494 case lwc1_op:
495 case swc1_op:
496 #if (__mips >= 2 || defined(__mips64))
497 case ldc1_op:
498 case sdc1_op:
499 #endif
500 case cop1_op:
501 #if __mips >= 4 && __mips != 32
502 case cop1x_op:
503 #endif
504 /* its one of ours */
505 goto emul;
506 #if __mips >= 4
507 case spec_op:
508 if (MIPSInst_FUNC(ir) == movc_op)
509 goto emul;
510 break;
511 #endif
515 * Single step the non-cp1
516 * instruction in the dslot
518 return mips_dsemul(xcp, ir, contpc);
520 else {
521 /* branch not taken */
522 if (likely) {
524 * branch likely nullifies
525 * dslot if not taken
527 xcp->cp0_epc += 4;
528 contpc += 4;
530 * else continue & execute
531 * dslot as normal insn
535 break;
538 default:
539 if (!(MIPSInst_RS(ir) & 0x10))
540 return SIGILL;
542 int sig;
544 /* a real fpu computation instruction */
545 if ((sig = fpu_emu(xcp, ctx, ir)))
546 return sig;
549 break;
551 #if __mips >= 4 && __mips != 32
552 case cop1x_op:{
553 int sig = fpux_emu(xcp, ctx, ir, fault_addr);
554 if (sig)
555 return sig;
556 break;
558 #endif
560 #if __mips >= 4
561 case spec_op:
562 if (MIPSInst_FUNC(ir) != movc_op)
563 return SIGILL;
564 cond = fpucondbit[MIPSInst_RT(ir) >> 2];
565 if (((ctx->fcr31 & cond) != 0) == ((MIPSInst_RT(ir) & 1) != 0))
566 xcp->regs[MIPSInst_RD(ir)] =
567 xcp->regs[MIPSInst_RS(ir)];
568 break;
569 #endif
571 default:
572 return SIGILL;
575 /* we did it !! */
576 xcp->cp0_epc = contpc;
577 xcp->cp0_cause &= ~CAUSEF_BD;
579 return 0;
583 * Conversion table from MIPS compare ops 48-63
584 * cond = ieee754dp_cmp(x,y,IEEE754_UN,sig);
586 static const unsigned char cmptab[8] = {
587 0, /* cmp_0 (sig) cmp_sf */
588 IEEE754_CUN, /* cmp_un (sig) cmp_ngle */
589 IEEE754_CEQ, /* cmp_eq (sig) cmp_seq */
590 IEEE754_CEQ | IEEE754_CUN, /* cmp_ueq (sig) cmp_ngl */
591 IEEE754_CLT, /* cmp_olt (sig) cmp_lt */
592 IEEE754_CLT | IEEE754_CUN, /* cmp_ult (sig) cmp_nge */
593 IEEE754_CLT | IEEE754_CEQ, /* cmp_ole (sig) cmp_le */
594 IEEE754_CLT | IEEE754_CEQ | IEEE754_CUN, /* cmp_ule (sig) cmp_ngt */
598 #if __mips >= 4 && __mips != 32
601 * Additional MIPS4 instructions
604 #define DEF3OP(name, p, f1, f2, f3) \
605 static ieee754##p fpemu_##p##_##name(ieee754##p r, ieee754##p s, \
606 ieee754##p t) \
608 struct _ieee754_csr ieee754_csr_save; \
609 s = f1(s, t); \
610 ieee754_csr_save = ieee754_csr; \
611 s = f2(s, r); \
612 ieee754_csr_save.cx |= ieee754_csr.cx; \
613 ieee754_csr_save.sx |= ieee754_csr.sx; \
614 s = f3(s); \
615 ieee754_csr.cx |= ieee754_csr_save.cx; \
616 ieee754_csr.sx |= ieee754_csr_save.sx; \
617 return s; \
620 static ieee754dp fpemu_dp_recip(ieee754dp d)
622 return ieee754dp_div(ieee754dp_one(0), d);
625 static ieee754dp fpemu_dp_rsqrt(ieee754dp d)
627 return ieee754dp_div(ieee754dp_one(0), ieee754dp_sqrt(d));
630 static ieee754sp fpemu_sp_recip(ieee754sp s)
632 return ieee754sp_div(ieee754sp_one(0), s);
635 static ieee754sp fpemu_sp_rsqrt(ieee754sp s)
637 return ieee754sp_div(ieee754sp_one(0), ieee754sp_sqrt(s));
640 DEF3OP(madd, sp, ieee754sp_mul, ieee754sp_add, );
641 DEF3OP(msub, sp, ieee754sp_mul, ieee754sp_sub, );
642 DEF3OP(nmadd, sp, ieee754sp_mul, ieee754sp_add, ieee754sp_neg);
643 DEF3OP(nmsub, sp, ieee754sp_mul, ieee754sp_sub, ieee754sp_neg);
644 DEF3OP(madd, dp, ieee754dp_mul, ieee754dp_add, );
645 DEF3OP(msub, dp, ieee754dp_mul, ieee754dp_sub, );
646 DEF3OP(nmadd, dp, ieee754dp_mul, ieee754dp_add, ieee754dp_neg);
647 DEF3OP(nmsub, dp, ieee754dp_mul, ieee754dp_sub, ieee754dp_neg);
649 static int fpux_emu(struct pt_regs *xcp, struct mips_fpu_struct *ctx,
650 mips_instruction ir, void *__user *fault_addr)
652 unsigned rcsr = 0; /* resulting csr */
654 MIPS_FPU_EMU_INC_STATS(cp1xops);
656 switch (MIPSInst_FMA_FFMT(ir)) {
657 case s_fmt:{ /* 0 */
659 ieee754sp(*handler) (ieee754sp, ieee754sp, ieee754sp);
660 ieee754sp fd, fr, fs, ft;
661 u32 __user *va;
662 u32 val;
664 switch (MIPSInst_FUNC(ir)) {
665 case lwxc1_op:
666 va = (void __user *) (xcp->regs[MIPSInst_FR(ir)] +
667 xcp->regs[MIPSInst_FT(ir)]);
669 MIPS_FPU_EMU_INC_STATS(loads);
670 if (!access_ok(VERIFY_READ, va, sizeof(u32))) {
671 MIPS_FPU_EMU_INC_STATS(errors);
672 *fault_addr = va;
673 return SIGBUS;
675 if (__get_user(val, va)) {
676 MIPS_FPU_EMU_INC_STATS(errors);
677 *fault_addr = va;
678 return SIGSEGV;
680 SITOREG(val, MIPSInst_FD(ir));
681 break;
683 case swxc1_op:
684 va = (void __user *) (xcp->regs[MIPSInst_FR(ir)] +
685 xcp->regs[MIPSInst_FT(ir)]);
687 MIPS_FPU_EMU_INC_STATS(stores);
689 SIFROMREG(val, MIPSInst_FS(ir));
690 if (!access_ok(VERIFY_WRITE, va, sizeof(u32))) {
691 MIPS_FPU_EMU_INC_STATS(errors);
692 *fault_addr = va;
693 return SIGBUS;
695 if (put_user(val, va)) {
696 MIPS_FPU_EMU_INC_STATS(errors);
697 *fault_addr = va;
698 return SIGSEGV;
700 break;
702 case madd_s_op:
703 handler = fpemu_sp_madd;
704 goto scoptop;
705 case msub_s_op:
706 handler = fpemu_sp_msub;
707 goto scoptop;
708 case nmadd_s_op:
709 handler = fpemu_sp_nmadd;
710 goto scoptop;
711 case nmsub_s_op:
712 handler = fpemu_sp_nmsub;
713 goto scoptop;
715 scoptop:
716 SPFROMREG(fr, MIPSInst_FR(ir));
717 SPFROMREG(fs, MIPSInst_FS(ir));
718 SPFROMREG(ft, MIPSInst_FT(ir));
719 fd = (*handler) (fr, fs, ft);
720 SPTOREG(fd, MIPSInst_FD(ir));
722 copcsr:
723 if (ieee754_cxtest(IEEE754_INEXACT))
724 rcsr |= FPU_CSR_INE_X | FPU_CSR_INE_S;
725 if (ieee754_cxtest(IEEE754_UNDERFLOW))
726 rcsr |= FPU_CSR_UDF_X | FPU_CSR_UDF_S;
727 if (ieee754_cxtest(IEEE754_OVERFLOW))
728 rcsr |= FPU_CSR_OVF_X | FPU_CSR_OVF_S;
729 if (ieee754_cxtest(IEEE754_INVALID_OPERATION))
730 rcsr |= FPU_CSR_INV_X | FPU_CSR_INV_S;
732 ctx->fcr31 = (ctx->fcr31 & ~FPU_CSR_ALL_X) | rcsr;
733 if ((ctx->fcr31 >> 5) & ctx->fcr31 & FPU_CSR_ALL_E) {
734 /*printk ("SIGFPE: fpu csr = %08x\n",
735 ctx->fcr31); */
736 return SIGFPE;
739 break;
741 default:
742 return SIGILL;
744 break;
747 case d_fmt:{ /* 1 */
748 ieee754dp(*handler) (ieee754dp, ieee754dp, ieee754dp);
749 ieee754dp fd, fr, fs, ft;
750 u64 __user *va;
751 u64 val;
753 switch (MIPSInst_FUNC(ir)) {
754 case ldxc1_op:
755 va = (void __user *) (xcp->regs[MIPSInst_FR(ir)] +
756 xcp->regs[MIPSInst_FT(ir)]);
758 MIPS_FPU_EMU_INC_STATS(loads);
759 if (!access_ok(VERIFY_READ, va, sizeof(u64))) {
760 MIPS_FPU_EMU_INC_STATS(errors);
761 *fault_addr = va;
762 return SIGBUS;
764 if (__get_user(val, va)) {
765 MIPS_FPU_EMU_INC_STATS(errors);
766 *fault_addr = va;
767 return SIGSEGV;
769 DITOREG(val, MIPSInst_FD(ir));
770 break;
772 case sdxc1_op:
773 va = (void __user *) (xcp->regs[MIPSInst_FR(ir)] +
774 xcp->regs[MIPSInst_FT(ir)]);
776 MIPS_FPU_EMU_INC_STATS(stores);
777 DIFROMREG(val, MIPSInst_FS(ir));
778 if (!access_ok(VERIFY_WRITE, va, sizeof(u64))) {
779 MIPS_FPU_EMU_INC_STATS(errors);
780 *fault_addr = va;
781 return SIGBUS;
783 if (__put_user(val, va)) {
784 MIPS_FPU_EMU_INC_STATS(errors);
785 *fault_addr = va;
786 return SIGSEGV;
788 break;
790 case madd_d_op:
791 handler = fpemu_dp_madd;
792 goto dcoptop;
793 case msub_d_op:
794 handler = fpemu_dp_msub;
795 goto dcoptop;
796 case nmadd_d_op:
797 handler = fpemu_dp_nmadd;
798 goto dcoptop;
799 case nmsub_d_op:
800 handler = fpemu_dp_nmsub;
801 goto dcoptop;
803 dcoptop:
804 DPFROMREG(fr, MIPSInst_FR(ir));
805 DPFROMREG(fs, MIPSInst_FS(ir));
806 DPFROMREG(ft, MIPSInst_FT(ir));
807 fd = (*handler) (fr, fs, ft);
808 DPTOREG(fd, MIPSInst_FD(ir));
809 goto copcsr;
811 default:
812 return SIGILL;
814 break;
817 case 0x7: /* 7 */
818 if (MIPSInst_FUNC(ir) != pfetch_op) {
819 return SIGILL;
821 /* ignore prefx operation */
822 break;
824 default:
825 return SIGILL;
828 return 0;
830 #endif
835 * Emulate a single COP1 arithmetic instruction.
837 static int fpu_emu(struct pt_regs *xcp, struct mips_fpu_struct *ctx,
838 mips_instruction ir)
840 int rfmt; /* resulting format */
841 unsigned rcsr = 0; /* resulting csr */
842 unsigned cond;
843 union {
844 ieee754dp d;
845 ieee754sp s;
846 int w;
847 #ifdef __mips64
848 s64 l;
849 #endif
850 } rv; /* resulting value */
852 MIPS_FPU_EMU_INC_STATS(cp1ops);
853 switch (rfmt = (MIPSInst_FFMT(ir) & 0xf)) {
854 case s_fmt:{ /* 0 */
855 union {
856 ieee754sp(*b) (ieee754sp, ieee754sp);
857 ieee754sp(*u) (ieee754sp);
858 } handler;
860 switch (MIPSInst_FUNC(ir)) {
861 /* binary ops */
862 case fadd_op:
863 handler.b = ieee754sp_add;
864 goto scopbop;
865 case fsub_op:
866 handler.b = ieee754sp_sub;
867 goto scopbop;
868 case fmul_op:
869 handler.b = ieee754sp_mul;
870 goto scopbop;
871 case fdiv_op:
872 handler.b = ieee754sp_div;
873 goto scopbop;
875 /* unary ops */
876 #if __mips >= 2 || defined(__mips64)
877 case fsqrt_op:
878 handler.u = ieee754sp_sqrt;
879 goto scopuop;
880 #endif
881 #if __mips >= 4 && __mips != 32
882 case frsqrt_op:
883 handler.u = fpemu_sp_rsqrt;
884 goto scopuop;
885 case frecip_op:
886 handler.u = fpemu_sp_recip;
887 goto scopuop;
888 #endif
889 #if __mips >= 4
890 case fmovc_op:
891 cond = fpucondbit[MIPSInst_FT(ir) >> 2];
892 if (((ctx->fcr31 & cond) != 0) !=
893 ((MIPSInst_FT(ir) & 1) != 0))
894 return 0;
895 SPFROMREG(rv.s, MIPSInst_FS(ir));
896 break;
897 case fmovz_op:
898 if (xcp->regs[MIPSInst_FT(ir)] != 0)
899 return 0;
900 SPFROMREG(rv.s, MIPSInst_FS(ir));
901 break;
902 case fmovn_op:
903 if (xcp->regs[MIPSInst_FT(ir)] == 0)
904 return 0;
905 SPFROMREG(rv.s, MIPSInst_FS(ir));
906 break;
907 #endif
908 case fabs_op:
909 handler.u = ieee754sp_abs;
910 goto scopuop;
911 case fneg_op:
912 handler.u = ieee754sp_neg;
913 goto scopuop;
914 case fmov_op:
915 /* an easy one */
916 SPFROMREG(rv.s, MIPSInst_FS(ir));
917 goto copcsr;
919 /* binary op on handler */
920 scopbop:
922 ieee754sp fs, ft;
924 SPFROMREG(fs, MIPSInst_FS(ir));
925 SPFROMREG(ft, MIPSInst_FT(ir));
927 rv.s = (*handler.b) (fs, ft);
928 goto copcsr;
930 scopuop:
932 ieee754sp fs;
934 SPFROMREG(fs, MIPSInst_FS(ir));
935 rv.s = (*handler.u) (fs);
936 goto copcsr;
938 copcsr:
939 if (ieee754_cxtest(IEEE754_INEXACT))
940 rcsr |= FPU_CSR_INE_X | FPU_CSR_INE_S;
941 if (ieee754_cxtest(IEEE754_UNDERFLOW))
942 rcsr |= FPU_CSR_UDF_X | FPU_CSR_UDF_S;
943 if (ieee754_cxtest(IEEE754_OVERFLOW))
944 rcsr |= FPU_CSR_OVF_X | FPU_CSR_OVF_S;
945 if (ieee754_cxtest(IEEE754_ZERO_DIVIDE))
946 rcsr |= FPU_CSR_DIV_X | FPU_CSR_DIV_S;
947 if (ieee754_cxtest(IEEE754_INVALID_OPERATION))
948 rcsr |= FPU_CSR_INV_X | FPU_CSR_INV_S;
949 break;
951 /* unary conv ops */
952 case fcvts_op:
953 return SIGILL; /* not defined */
954 case fcvtd_op:{
955 ieee754sp fs;
957 SPFROMREG(fs, MIPSInst_FS(ir));
958 rv.d = ieee754dp_fsp(fs);
959 rfmt = d_fmt;
960 goto copcsr;
962 case fcvtw_op:{
963 ieee754sp fs;
965 SPFROMREG(fs, MIPSInst_FS(ir));
966 rv.w = ieee754sp_tint(fs);
967 rfmt = w_fmt;
968 goto copcsr;
971 #if __mips >= 2 || defined(__mips64)
972 case fround_op:
973 case ftrunc_op:
974 case fceil_op:
975 case ffloor_op:{
976 unsigned int oldrm = ieee754_csr.rm;
977 ieee754sp fs;
979 SPFROMREG(fs, MIPSInst_FS(ir));
980 ieee754_csr.rm = ieee_rm[modeindex(MIPSInst_FUNC(ir))];
981 rv.w = ieee754sp_tint(fs);
982 ieee754_csr.rm = oldrm;
983 rfmt = w_fmt;
984 goto copcsr;
986 #endif /* __mips >= 2 */
988 #if defined(__mips64)
989 case fcvtl_op:{
990 ieee754sp fs;
992 SPFROMREG(fs, MIPSInst_FS(ir));
993 rv.l = ieee754sp_tlong(fs);
994 rfmt = l_fmt;
995 goto copcsr;
998 case froundl_op:
999 case ftruncl_op:
1000 case fceill_op:
1001 case ffloorl_op:{
1002 unsigned int oldrm = ieee754_csr.rm;
1003 ieee754sp fs;
1005 SPFROMREG(fs, MIPSInst_FS(ir));
1006 ieee754_csr.rm = ieee_rm[modeindex(MIPSInst_FUNC(ir))];
1007 rv.l = ieee754sp_tlong(fs);
1008 ieee754_csr.rm = oldrm;
1009 rfmt = l_fmt;
1010 goto copcsr;
1012 #endif /* defined(__mips64) */
1014 default:
1015 if (MIPSInst_FUNC(ir) >= fcmp_op) {
1016 unsigned cmpop = MIPSInst_FUNC(ir) - fcmp_op;
1017 ieee754sp fs, ft;
1019 SPFROMREG(fs, MIPSInst_FS(ir));
1020 SPFROMREG(ft, MIPSInst_FT(ir));
1021 rv.w = ieee754sp_cmp(fs, ft,
1022 cmptab[cmpop & 0x7], cmpop & 0x8);
1023 rfmt = -1;
1024 if ((cmpop & 0x8) && ieee754_cxtest
1025 (IEEE754_INVALID_OPERATION))
1026 rcsr = FPU_CSR_INV_X | FPU_CSR_INV_S;
1027 else
1028 goto copcsr;
1031 else {
1032 return SIGILL;
1034 break;
1036 break;
1039 case d_fmt:{
1040 union {
1041 ieee754dp(*b) (ieee754dp, ieee754dp);
1042 ieee754dp(*u) (ieee754dp);
1043 } handler;
1045 switch (MIPSInst_FUNC(ir)) {
1046 /* binary ops */
1047 case fadd_op:
1048 handler.b = ieee754dp_add;
1049 goto dcopbop;
1050 case fsub_op:
1051 handler.b = ieee754dp_sub;
1052 goto dcopbop;
1053 case fmul_op:
1054 handler.b = ieee754dp_mul;
1055 goto dcopbop;
1056 case fdiv_op:
1057 handler.b = ieee754dp_div;
1058 goto dcopbop;
1060 /* unary ops */
1061 #if __mips >= 2 || defined(__mips64)
1062 case fsqrt_op:
1063 handler.u = ieee754dp_sqrt;
1064 goto dcopuop;
1065 #endif
1066 #if __mips >= 4 && __mips != 32
1067 case frsqrt_op:
1068 handler.u = fpemu_dp_rsqrt;
1069 goto dcopuop;
1070 case frecip_op:
1071 handler.u = fpemu_dp_recip;
1072 goto dcopuop;
1073 #endif
1074 #if __mips >= 4
1075 case fmovc_op:
1076 cond = fpucondbit[MIPSInst_FT(ir) >> 2];
1077 if (((ctx->fcr31 & cond) != 0) !=
1078 ((MIPSInst_FT(ir) & 1) != 0))
1079 return 0;
1080 DPFROMREG(rv.d, MIPSInst_FS(ir));
1081 break;
1082 case fmovz_op:
1083 if (xcp->regs[MIPSInst_FT(ir)] != 0)
1084 return 0;
1085 DPFROMREG(rv.d, MIPSInst_FS(ir));
1086 break;
1087 case fmovn_op:
1088 if (xcp->regs[MIPSInst_FT(ir)] == 0)
1089 return 0;
1090 DPFROMREG(rv.d, MIPSInst_FS(ir));
1091 break;
1092 #endif
1093 case fabs_op:
1094 handler.u = ieee754dp_abs;
1095 goto dcopuop;
1097 case fneg_op:
1098 handler.u = ieee754dp_neg;
1099 goto dcopuop;
1101 case fmov_op:
1102 /* an easy one */
1103 DPFROMREG(rv.d, MIPSInst_FS(ir));
1104 goto copcsr;
1106 /* binary op on handler */
1107 dcopbop:{
1108 ieee754dp fs, ft;
1110 DPFROMREG(fs, MIPSInst_FS(ir));
1111 DPFROMREG(ft, MIPSInst_FT(ir));
1113 rv.d = (*handler.b) (fs, ft);
1114 goto copcsr;
1116 dcopuop:{
1117 ieee754dp fs;
1119 DPFROMREG(fs, MIPSInst_FS(ir));
1120 rv.d = (*handler.u) (fs);
1121 goto copcsr;
1124 /* unary conv ops */
1125 case fcvts_op:{
1126 ieee754dp fs;
1128 DPFROMREG(fs, MIPSInst_FS(ir));
1129 rv.s = ieee754sp_fdp(fs);
1130 rfmt = s_fmt;
1131 goto copcsr;
1133 case fcvtd_op:
1134 return SIGILL; /* not defined */
1136 case fcvtw_op:{
1137 ieee754dp fs;
1139 DPFROMREG(fs, MIPSInst_FS(ir));
1140 rv.w = ieee754dp_tint(fs); /* wrong */
1141 rfmt = w_fmt;
1142 goto copcsr;
1145 #if __mips >= 2 || defined(__mips64)
1146 case fround_op:
1147 case ftrunc_op:
1148 case fceil_op:
1149 case ffloor_op:{
1150 unsigned int oldrm = ieee754_csr.rm;
1151 ieee754dp fs;
1153 DPFROMREG(fs, MIPSInst_FS(ir));
1154 ieee754_csr.rm = ieee_rm[modeindex(MIPSInst_FUNC(ir))];
1155 rv.w = ieee754dp_tint(fs);
1156 ieee754_csr.rm = oldrm;
1157 rfmt = w_fmt;
1158 goto copcsr;
1160 #endif
1162 #if defined(__mips64)
1163 case fcvtl_op:{
1164 ieee754dp fs;
1166 DPFROMREG(fs, MIPSInst_FS(ir));
1167 rv.l = ieee754dp_tlong(fs);
1168 rfmt = l_fmt;
1169 goto copcsr;
1172 case froundl_op:
1173 case ftruncl_op:
1174 case fceill_op:
1175 case ffloorl_op:{
1176 unsigned int oldrm = ieee754_csr.rm;
1177 ieee754dp fs;
1179 DPFROMREG(fs, MIPSInst_FS(ir));
1180 ieee754_csr.rm = ieee_rm[modeindex(MIPSInst_FUNC(ir))];
1181 rv.l = ieee754dp_tlong(fs);
1182 ieee754_csr.rm = oldrm;
1183 rfmt = l_fmt;
1184 goto copcsr;
1186 #endif /* __mips >= 3 */
1188 default:
1189 if (MIPSInst_FUNC(ir) >= fcmp_op) {
1190 unsigned cmpop = MIPSInst_FUNC(ir) - fcmp_op;
1191 ieee754dp fs, ft;
1193 DPFROMREG(fs, MIPSInst_FS(ir));
1194 DPFROMREG(ft, MIPSInst_FT(ir));
1195 rv.w = ieee754dp_cmp(fs, ft,
1196 cmptab[cmpop & 0x7], cmpop & 0x8);
1197 rfmt = -1;
1198 if ((cmpop & 0x8)
1200 ieee754_cxtest
1201 (IEEE754_INVALID_OPERATION))
1202 rcsr = FPU_CSR_INV_X | FPU_CSR_INV_S;
1203 else
1204 goto copcsr;
1207 else {
1208 return SIGILL;
1210 break;
1212 break;
1215 case w_fmt:{
1216 ieee754sp fs;
1218 switch (MIPSInst_FUNC(ir)) {
1219 case fcvts_op:
1220 /* convert word to single precision real */
1221 SPFROMREG(fs, MIPSInst_FS(ir));
1222 rv.s = ieee754sp_fint(fs.bits);
1223 rfmt = s_fmt;
1224 goto copcsr;
1225 case fcvtd_op:
1226 /* convert word to double precision real */
1227 SPFROMREG(fs, MIPSInst_FS(ir));
1228 rv.d = ieee754dp_fint(fs.bits);
1229 rfmt = d_fmt;
1230 goto copcsr;
1231 default:
1232 return SIGILL;
1234 break;
1237 #if defined(__mips64)
1238 case l_fmt:{
1239 switch (MIPSInst_FUNC(ir)) {
1240 case fcvts_op:
1241 /* convert long to single precision real */
1242 rv.s = ieee754sp_flong(ctx->fpr[MIPSInst_FS(ir)]);
1243 rfmt = s_fmt;
1244 goto copcsr;
1245 case fcvtd_op:
1246 /* convert long to double precision real */
1247 rv.d = ieee754dp_flong(ctx->fpr[MIPSInst_FS(ir)]);
1248 rfmt = d_fmt;
1249 goto copcsr;
1250 default:
1251 return SIGILL;
1253 break;
1255 #endif
1257 default:
1258 return SIGILL;
1262 * Update the fpu CSR register for this operation.
1263 * If an exception is required, generate a tidy SIGFPE exception,
1264 * without updating the result register.
1265 * Note: cause exception bits do not accumulate, they are rewritten
1266 * for each op; only the flag/sticky bits accumulate.
1268 ctx->fcr31 = (ctx->fcr31 & ~FPU_CSR_ALL_X) | rcsr;
1269 if ((ctx->fcr31 >> 5) & ctx->fcr31 & FPU_CSR_ALL_E) {
1270 /*printk ("SIGFPE: fpu csr = %08x\n",ctx->fcr31); */
1271 return SIGFPE;
1275 * Now we can safely write the result back to the register file.
1277 switch (rfmt) {
1278 case -1:{
1279 #if __mips >= 4
1280 cond = fpucondbit[MIPSInst_FD(ir) >> 2];
1281 #else
1282 cond = FPU_CSR_COND;
1283 #endif
1284 if (rv.w)
1285 ctx->fcr31 |= cond;
1286 else
1287 ctx->fcr31 &= ~cond;
1288 break;
1290 case d_fmt:
1291 DPTOREG(rv.d, MIPSInst_FD(ir));
1292 break;
1293 case s_fmt:
1294 SPTOREG(rv.s, MIPSInst_FD(ir));
1295 break;
1296 case w_fmt:
1297 SITOREG(rv.w, MIPSInst_FD(ir));
1298 break;
1299 #if defined(__mips64)
1300 case l_fmt:
1301 DITOREG(rv.l, MIPSInst_FD(ir));
1302 break;
1303 #endif
1304 default:
1305 return SIGILL;
1308 return 0;
1311 int fpu_emulator_cop1Handler(struct pt_regs *xcp, struct mips_fpu_struct *ctx,
1312 int has_fpu, void *__user *fault_addr)
1314 unsigned long oldepc, prevepc;
1315 mips_instruction insn;
1316 int sig = 0;
1318 oldepc = xcp->cp0_epc;
1319 do {
1320 prevepc = xcp->cp0_epc;
1322 if (!access_ok(VERIFY_READ, xcp->cp0_epc, sizeof(mips_instruction))) {
1323 MIPS_FPU_EMU_INC_STATS(errors);
1324 *fault_addr = (mips_instruction __user *)xcp->cp0_epc;
1325 return SIGBUS;
1327 if (__get_user(insn, (mips_instruction __user *) xcp->cp0_epc)) {
1328 MIPS_FPU_EMU_INC_STATS(errors);
1329 *fault_addr = (mips_instruction __user *)xcp->cp0_epc;
1330 return SIGSEGV;
1332 if (insn == 0)
1333 xcp->cp0_epc += 4; /* skip nops */
1334 else {
1336 * The 'ieee754_csr' is an alias of
1337 * ctx->fcr31. No need to copy ctx->fcr31 to
1338 * ieee754_csr. But ieee754_csr.rm is ieee
1339 * library modes. (not mips rounding mode)
1341 /* convert to ieee library modes */
1342 ieee754_csr.rm = ieee_rm[ieee754_csr.rm];
1343 sig = cop1Emulate(xcp, ctx, fault_addr);
1344 /* revert to mips rounding mode */
1345 ieee754_csr.rm = mips_rm[ieee754_csr.rm];
1348 if (has_fpu)
1349 break;
1350 if (sig)
1351 break;
1353 cond_resched();
1354 } while (xcp->cp0_epc > prevepc);
1356 /* SIGILL indicates a non-fpu instruction */
1357 if (sig == SIGILL && xcp->cp0_epc != oldepc)
1358 /* but if epc has advanced, then ignore it */
1359 sig = 0;
1361 return sig;
1364 #ifdef CONFIG_DEBUG_FS
1366 static int fpuemu_stat_get(void *data, u64 *val)
1368 int cpu;
1369 unsigned long sum = 0;
1370 for_each_online_cpu(cpu) {
1371 struct mips_fpu_emulator_stats *ps;
1372 local_t *pv;
1373 ps = &per_cpu(fpuemustats, cpu);
1374 pv = (void *)ps + (unsigned long)data;
1375 sum += local_read(pv);
1377 *val = sum;
1378 return 0;
1380 DEFINE_SIMPLE_ATTRIBUTE(fops_fpuemu_stat, fpuemu_stat_get, NULL, "%llu\n");
1382 extern struct dentry *mips_debugfs_dir;
1383 static int __init debugfs_fpuemu(void)
1385 struct dentry *d, *dir;
1387 if (!mips_debugfs_dir)
1388 return -ENODEV;
1389 dir = debugfs_create_dir("fpuemustats", mips_debugfs_dir);
1390 if (!dir)
1391 return -ENOMEM;
1393 #define FPU_STAT_CREATE(M) \
1394 do { \
1395 d = debugfs_create_file(#M , S_IRUGO, dir, \
1396 (void *)offsetof(struct mips_fpu_emulator_stats, M), \
1397 &fops_fpuemu_stat); \
1398 if (!d) \
1399 return -ENOMEM; \
1400 } while (0)
1402 FPU_STAT_CREATE(emulated);
1403 FPU_STAT_CREATE(loads);
1404 FPU_STAT_CREATE(stores);
1405 FPU_STAT_CREATE(cp1ops);
1406 FPU_STAT_CREATE(cp1xops);
1407 FPU_STAT_CREATE(errors);
1409 return 0;
1411 __initcall(debugfs_fpuemu);
1412 #endif