1 /***************************************************************************
2 * PCSX-Revolution - PlayStation Emulator for Nintendo Wii *
3 * Copyright (C) 2009-2010 PCSX-Revolution Dev Team *
4 * <http://code.google.com/p/pcsx-revolution/> *
6 * This program is free software; you can redistribute it and/or modify *
7 * it under the terms of the GNU General Public License as published by *
8 * the Free Software Foundation; either version 2 of the License, or *
9 * (at your option) any later version. *
11 * This program is distributed in the hope that it will be useful, *
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
14 * GNU General Public License for more details. *
16 * You should have received a copy of the GNU General Public License *
17 * along with this program; if not, write to the *
18 * Free Software Foundation, Inc., *
19 * 51 Franklin Street, Fifth Floor, Boston, MA 02111-1307 USA. *
20 ***************************************************************************/
29 #define VX(n) (n < 3 ? regs->CP2D.p[n << 1].sw.l : regs->CP2D.p[9].sw.l)
30 #define VY(n) (n < 3 ? regs->CP2D.p[n << 1].sw.h : regs->CP2D.p[10].sw.l)
31 #define VZ(n) (n < 3 ? regs->CP2D.p[(n << 1) + 1].sw.l : regs->CP2D.p[11].sw.l)
32 #define MX11(n) (n < 3 ? regs->CP2C.p[(n << 3)].sw.l : 0)
33 #define MX12(n) (n < 3 ? regs->CP2C.p[(n << 3)].sw.h : 0)
34 #define MX13(n) (n < 3 ? regs->CP2C.p[(n << 3) + 1].sw.l : 0)
35 #define MX21(n) (n < 3 ? regs->CP2C.p[(n << 3) + 1].sw.h : 0)
36 #define MX22(n) (n < 3 ? regs->CP2C.p[(n << 3) + 2].sw.l : 0)
37 #define MX23(n) (n < 3 ? regs->CP2C.p[(n << 3) + 2].sw.h : 0)
38 #define MX31(n) (n < 3 ? regs->CP2C.p[(n << 3) + 3].sw.l : 0)
39 #define MX32(n) (n < 3 ? regs->CP2C.p[(n << 3) + 3].sw.h : 0)
40 #define MX33(n) (n < 3 ? regs->CP2C.p[(n << 3) + 4].sw.l : 0)
41 #define CV1(n) (n < 3 ? (s32)regs->CP2C.r[(n << 3) + 5] : 0)
42 #define CV2(n) (n < 3 ? (s32)regs->CP2C.r[(n << 3) + 6] : 0)
43 #define CV3(n) (n < 3 ? (s32)regs->CP2C.r[(n << 3) + 7] : 0)
45 #define fSX(n) ((regs->CP2D.p)[((n) + 12)].sw.l)
46 #define fSY(n) ((regs->CP2D.p)[((n) + 12)].sw.h)
47 #define fSZ(n) ((regs->CP2D.p)[((n) + 17)].w.l) /* (n == 0) => SZ1; */
49 #define gteVXY0 (regs->CP2D.r[0])
50 #define gteVX0 (regs->CP2D.p[0].sw.l)
51 #define gteVY0 (regs->CP2D.p[0].sw.h)
52 #define gteVZ0 (regs->CP2D.p[1].sw.l)
53 #define gteVXY1 (regs->CP2D.r[2])
54 #define gteVX1 (regs->CP2D.p[2].sw.l)
55 #define gteVY1 (regs->CP2D.p[2].sw.h)
56 #define gteVZ1 (regs->CP2D.p[3].sw.l)
57 #define gteVXY2 (regs->CP2D.r[4])
58 #define gteVX2 (regs->CP2D.p[4].sw.l)
59 #define gteVY2 (regs->CP2D.p[4].sw.h)
60 #define gteVZ2 (regs->CP2D.p[5].sw.l)
61 #define gteRGB (regs->CP2D.r[6])
62 #define gteR (regs->CP2D.p[6].b.l)
63 #define gteG (regs->CP2D.p[6].b.h)
64 #define gteB (regs->CP2D.p[6].b.h2)
65 #define gteCODE (regs->CP2D.p[6].b.h3)
66 #define gteOTZ (regs->CP2D.p[7].w.l)
67 #define gteIR0 (regs->CP2D.p[8].sw.l)
68 #define gteIR1 (regs->CP2D.p[9].sw.l)
69 #define gteIR2 (regs->CP2D.p[10].sw.l)
70 #define gteIR3 (regs->CP2D.p[11].sw.l)
71 #define gteSXY0 (regs->CP2D.r[12])
72 #define gteSX0 (regs->CP2D.p[12].sw.l)
73 #define gteSY0 (regs->CP2D.p[12].sw.h)
74 #define gteSXY1 (regs->CP2D.r[13])
75 #define gteSX1 (regs->CP2D.p[13].sw.l)
76 #define gteSY1 (regs->CP2D.p[13].sw.h)
77 #define gteSXY2 (regs->CP2D.r[14])
78 #define gteSX2 (regs->CP2D.p[14].sw.l)
79 #define gteSY2 (regs->CP2D.p[14].sw.h)
80 #define gteSXYP (regs->CP2D.r[15])
81 #define gteSXP (regs->CP2D.p[15].sw.l)
82 #define gteSYP (regs->CP2D.p[15].sw.h)
83 #define gteSZ0 (regs->CP2D.p[16].w.l)
84 #define gteSZ1 (regs->CP2D.p[17].w.l)
85 #define gteSZ2 (regs->CP2D.p[18].w.l)
86 #define gteSZ3 (regs->CP2D.p[19].w.l)
87 #define gteRGB0 (regs->CP2D.r[20])
88 #define gteR0 (regs->CP2D.p[20].b.l)
89 #define gteG0 (regs->CP2D.p[20].b.h)
90 #define gteB0 (regs->CP2D.p[20].b.h2)
91 #define gteCODE0 (regs->CP2D.p[20].b.h3)
92 #define gteRGB1 (regs->CP2D.r[21])
93 #define gteR1 (regs->CP2D.p[21].b.l)
94 #define gteG1 (regs->CP2D.p[21].b.h)
95 #define gteB1 (regs->CP2D.p[21].b.h2)
96 #define gteCODE1 (regs->CP2D.p[21].b.h3)
97 #define gteRGB2 (regs->CP2D.r[22])
98 #define gteR2 (regs->CP2D.p[22].b.l)
99 #define gteG2 (regs->CP2D.p[22].b.h)
100 #define gteB2 (regs->CP2D.p[22].b.h2)
101 #define gteCODE2 (regs->CP2D.p[22].b.h3)
102 #define gteRES1 (regs->CP2D.r[23])
103 #define gteMAC0 (((s32 *)regs->CP2D.r)[24])
104 #define gteMAC1 (((s32 *)regs->CP2D.r)[25])
105 #define gteMAC2 (((s32 *)regs->CP2D.r)[26])
106 #define gteMAC3 (((s32 *)regs->CP2D.r)[27])
107 #define gteIRGB (regs->CP2D.r[28])
108 #define gteORGB (regs->CP2D.r[29])
109 #define gteLZCS (regs->CP2D.r[30])
110 #define gteLZCR (regs->CP2D.r[31])
112 #define gteR11R12 (((s32 *)regs->CP2C.r)[0])
113 #define gteR22R23 (((s32 *)regs->CP2C.r)[2])
114 #define gteR11 (regs->CP2C.p[0].sw.l)
115 #define gteR12 (regs->CP2C.p[0].sw.h)
116 #define gteR13 (regs->CP2C.p[1].sw.l)
117 #define gteR21 (regs->CP2C.p[1].sw.h)
118 #define gteR22 (regs->CP2C.p[2].sw.l)
119 #define gteR23 (regs->CP2C.p[2].sw.h)
120 #define gteR31 (regs->CP2C.p[3].sw.l)
121 #define gteR32 (regs->CP2C.p[3].sw.h)
122 #define gteR33 (regs->CP2C.p[4].sw.l)
123 #define gteTRX (((s32 *)regs->CP2C.r)[5])
124 #define gteTRY (((s32 *)regs->CP2C.r)[6])
125 #define gteTRZ (((s32 *)regs->CP2C.r)[7])
126 #define gteL11 (regs->CP2C.p[8].sw.l)
127 #define gteL12 (regs->CP2C.p[8].sw.h)
128 #define gteL13 (regs->CP2C.p[9].sw.l)
129 #define gteL21 (regs->CP2C.p[9].sw.h)
130 #define gteL22 (regs->CP2C.p[10].sw.l)
131 #define gteL23 (regs->CP2C.p[10].sw.h)
132 #define gteL31 (regs->CP2C.p[11].sw.l)
133 #define gteL32 (regs->CP2C.p[11].sw.h)
134 #define gteL33 (regs->CP2C.p[12].sw.l)
135 #define gteRBK (((s32 *)regs->CP2C.r)[13])
136 #define gteGBK (((s32 *)regs->CP2C.r)[14])
137 #define gteBBK (((s32 *)regs->CP2C.r)[15])
138 #define gteLR1 (regs->CP2C.p[16].sw.l)
139 #define gteLR2 (regs->CP2C.p[16].sw.h)
140 #define gteLR3 (regs->CP2C.p[17].sw.l)
141 #define gteLG1 (regs->CP2C.p[17].sw.h)
142 #define gteLG2 (regs->CP2C.p[18].sw.l)
143 #define gteLG3 (regs->CP2C.p[18].sw.h)
144 #define gteLB1 (regs->CP2C.p[19].sw.l)
145 #define gteLB2 (regs->CP2C.p[19].sw.h)
146 #define gteLB3 (regs->CP2C.p[20].sw.l)
147 #define gteRFC (((s32 *)regs->CP2C.r)[21])
148 #define gteGFC (((s32 *)regs->CP2C.r)[22])
149 #define gteBFC (((s32 *)regs->CP2C.r)[23])
150 #define gteOFX (((s32 *)regs->CP2C.r)[24])
151 #define gteOFY (((s32 *)regs->CP2C.r)[25])
152 // senquack - gteH register is u16, not s16, and used in GTE that way.
153 // HOWEVER when read back by CPU using CFC2, it will be incorrectly
154 // sign-extended by bug in original hardware, according to Nocash docs
155 // GTE section 'Screen Offset and Distance'. The emulator does this
156 // sign extension when it is loaded to GTE by CTC2.
157 //#define gteH (regs->CP2C.p[26].sw.l)
158 #define gteH (regs->CP2C.p[26].w.l)
159 #define gteDQA (regs->CP2C.p[27].sw.l)
160 #define gteDQB (((s32 *)regs->CP2C.r)[28])
161 #define gteZSF3 (regs->CP2C.p[29].sw.l)
162 #define gteZSF4 (regs->CP2C.p[30].sw.l)
163 #define gteFLAG (regs->CP2C.r[31])
165 #define GTE_OP(op) ((op >> 20) & 31)
166 #define GTE_SF(op) ((op >> 19) & 1)
167 #define GTE_MX(op) ((op >> 17) & 3)
168 #define GTE_V(op) ((op >> 15) & 3)
169 #define GTE_CV(op) ((op >> 13) & 3)
170 #define GTE_CD(op) ((op >> 11) & 3) /* not used */
171 #define GTE_LM(op) ((op >> 10) & 1)
172 #define GTE_CT(op) ((op >> 6) & 15) /* not used */
173 #define GTE_FUNCT(op) (op & 63)
175 #define gteop (psxRegs.code & 0x1ffffff)
179 static inline s64 BOUNDS_(psxCP2Regs *regs, s64 n_value, s64 n_max, int n_maxflag, s64 n_min, int n_minflag) {
180 if (n_value > n_max) {
181 gteFLAG |= n_maxflag;
182 } else if (n_value < n_min) {
183 gteFLAG |= n_minflag;
188 static inline s32 LIM_(psxCP2Regs *regs, s32 value, s32 max, s32 min, u32 flag) {
193 } else if (value < min) {
200 static inline u32 limE_(psxCP2Regs *regs, u32 result) {
201 if (result > 0x1ffff) {
202 gteFLAG |= (1 << 31) | (1 << 17);
210 #define BOUNDS_(regs, a, ...) (a)
212 static inline s32 LIM_(psxCP2Regs *regs, s32 value, s32 max, s32 min, u32 flag_unused) {
216 else if (value < min)
221 static inline u32 limE_(psxCP2Regs *regs, u32 result) {
222 if (result > 0x1ffff)
229 #define BOUNDS(n_value,n_max,n_maxflag,n_min,n_minflag) \
230 BOUNDS_(regs,n_value,n_max,n_maxflag,n_min,n_minflag)
231 #define LIM(value,max,min,flag) \
232 LIM_(regs,value,max,min,flag)
233 #define limE(result) \
236 #define A1(a) BOUNDS((a), 0x7fffffff, (1 << 30), -(s64)0x80000000, (1 << 31) | (1 << 27))
237 #define A2(a) BOUNDS((a), 0x7fffffff, (1 << 29), -(s64)0x80000000, (1 << 31) | (1 << 26))
238 #define A3(a) BOUNDS((a), 0x7fffffff, (1 << 28), -(s64)0x80000000, (1 << 31) | (1 << 25))
239 #define limB1(a, l) LIM((a), 0x7fff, -0x8000 * !l, (1 << 31) | (1 << 24))
240 #define limB2(a, l) LIM((a), 0x7fff, -0x8000 * !l, (1 << 31) | (1 << 23))
241 #define limB3(a, l) LIM((a), 0x7fff, -0x8000 * !l, (1 << 22))
242 #define limC1(a) LIM((a), 0x00ff, 0x0000, (1 << 21))
243 #define limC2(a) LIM((a), 0x00ff, 0x0000, (1 << 20))
244 #define limC3(a) LIM((a), 0x00ff, 0x0000, (1 << 19))
245 #define limD(a) LIM((a), 0xffff, 0x0000, (1 << 31) | (1 << 18))
247 #define F(a) BOUNDS((a), 0x7fffffff, (1 << 31) | (1 << 16), -(s64)0x80000000, (1 << 31) | (1 << 15))
248 #define limG1(a) LIM((a), 0x3ff, -0x400, (1 << 31) | (1 << 14))
249 #define limG2(a) LIM((a), 0x3ff, -0x400, (1 << 31) | (1 << 13))
250 #define limH(a) LIM((a), 0x1000, 0x0000, (1 << 12))
257 /* these are unlikely to be hit and usually waste cycles, don't want them on ARM */
263 //senquack - n param should be unsigned (will be 'gteH' reg which is u16)
264 #ifdef GTE_USE_NATIVE_DIVIDE
265 INLINE u32 DIVIDE(u16 n, u16 d) {
267 return ((u32)n << 16) / d;
272 #include "gte_divider.h"
273 #endif // GTE_USE_NATIVE_DIVIDE
278 psxCP2Regs *regs = &psxRegs.CP2;
287 psxRegs.CP2D.r[reg] = (s32)psxRegs.CP2D.p[reg].sw.l;
295 psxRegs.CP2D.r[reg] = (u32)psxRegs.CP2D.p[reg].w.l;
299 psxRegs.CP2D.r[reg] = gteSXY2;
304 psxRegs.CP2D.r[reg] = LIM(gteIR1 >> 7, 0x1f, 0, 0) |
305 (LIM(gteIR2 >> 7, 0x1f, 0, 0) << 5) |
306 (LIM(gteIR3 >> 7, 0x1f, 0, 0) << 10);
309 return psxRegs.CP2D.r[reg];
312 void MTC2(u32 value, int reg) {
313 psxCP2Regs *regs = &psxRegs.CP2;
324 gteIR1 = (value & 0x1f) << 7;
325 gteIR2 = (value & 0x3e0) << 2;
326 gteIR3 = (value & 0x7c00) >> 3;
337 for (i = 31; (a & (1 << i)) == 0 && i >= 0; i--);
342 for (i = 31; (a & (1 << i)) == 0 && i >= 0; i--);
354 psxRegs.CP2D.r[reg] = value;
358 void CTC2(u32 value, int reg) {
367 value = (s32)(s16)value;
371 value = value & 0x7ffff000;
372 if (value & 0x7f87e000) value |= 0x80000000;
376 psxRegs.CP2C.r[reg] = value;
381 psxRegs.GPR.r[_Rt_] = MFC2(_Rd_);
386 psxRegs.GPR.r[_Rt_] = psxRegs.CP2C.r[_Rd_];
390 MTC2(psxRegs.GPR.r[_Rt_], _Rd_);
394 CTC2(psxRegs.GPR.r[_Rt_], _Rd_);
397 #define _oB_ (psxRegs.GPR.r[_Rs_] + _Imm_)
400 MTC2(psxMemRead32(_oB_), _Rt_);
404 psxMemWrite32(_oB_, MFC2(_Rt_));
410 #define DIVIDE DIVIDE_
411 static u32 DIVIDE_(s16 n, u16 d) {
412 if (n >= 0 && n < d * 2) {
414 return ((n_ << 16) + d / 2) / d;
415 //return (u32)((float)(n_ << 16) / (float)d + (float)0.5);
421 void gteRTPS(psxCP2Regs *regs) {
426 GTE_LOG("GTE RTPS\n");
430 gteMAC1 = A1((((s64)gteTRX << 12) + (gteR11 * gteVX0) + (gteR12 * gteVY0) + (gteR13 * gteVZ0)) >> 12);
431 gteMAC2 = A2((((s64)gteTRY << 12) + (gteR21 * gteVX0) + (gteR22 * gteVY0) + (gteR23 * gteVZ0)) >> 12);
432 gteMAC3 = A3((((s64)gteTRZ << 12) + (gteR31 * gteVX0) + (gteR32 * gteVY0) + (gteR33 * gteVZ0)) >> 12);
433 gteIR1 = limB1(gteMAC1, 0);
434 gteIR2 = limB2(gteMAC2, 0);
435 gteIR3 = limB3(gteMAC3, 0);
439 gteSZ3 = limD(gteMAC3);
440 quotient = limE(DIVIDE(gteH, gteSZ3));
443 gteSX2 = limG1(F((s64)gteOFX + ((s64)gteIR1 * quotient)) >> 16);
444 gteSY2 = limG2(F((s64)gteOFY + ((s64)gteIR2 * quotient)) >> 16);
446 tmp = (s64)gteDQB + ((s64)gteDQA * quotient);
448 gteIR0 = limH(tmp >> 12);
451 void gteRTPT(psxCP2Regs *regs) {
458 GTE_LOG("GTE RTPT\n");
463 for (v = 0; v < 3; v++) {
467 gteMAC1 = A1((((s64)gteTRX << 12) + (gteR11 * vx) + (gteR12 * vy) + (gteR13 * vz)) >> 12);
468 gteMAC2 = A2((((s64)gteTRY << 12) + (gteR21 * vx) + (gteR22 * vy) + (gteR23 * vz)) >> 12);
469 gteMAC3 = A3((((s64)gteTRZ << 12) + (gteR31 * vx) + (gteR32 * vy) + (gteR33 * vz)) >> 12);
470 gteIR1 = limB1(gteMAC1, 0);
471 gteIR2 = limB2(gteMAC2, 0);
472 gteIR3 = limB3(gteMAC3, 0);
473 fSZ(v) = limD(gteMAC3);
474 quotient = limE(DIVIDE(gteH, fSZ(v)));
475 fSX(v) = limG1(F((s64)gteOFX + ((s64)gteIR1 * quotient)) >> 16);
476 fSY(v) = limG2(F((s64)gteOFY + ((s64)gteIR2 * quotient)) >> 16);
479 tmp = (s64)gteDQB + ((s64)gteDQA * quotient);
481 gteIR0 = limH(tmp >> 12);
484 void gteMVMVA(psxCP2Regs *regs) {
485 int shift = 12 * GTE_SF(gteop);
486 int mx = GTE_MX(gteop);
487 int v = GTE_V(gteop);
488 int cv = GTE_CV(gteop);
489 int lm = GTE_LM(gteop);
495 GTE_LOG("GTE MVMVA\n");
499 gteMAC1 = A1((((s64)CV1(cv) << 12) + (MX11(mx) * vx) + (MX12(mx) * vy) + (MX13(mx) * vz)) >> shift);
500 gteMAC2 = A2((((s64)CV2(cv) << 12) + (MX21(mx) * vx) + (MX22(mx) * vy) + (MX23(mx) * vz)) >> shift);
501 gteMAC3 = A3((((s64)CV3(cv) << 12) + (MX31(mx) * vx) + (MX32(mx) * vy) + (MX33(mx) * vz)) >> shift);
503 gteIR1 = limB1(gteMAC1, lm);
504 gteIR2 = limB2(gteMAC2, lm);
505 gteIR3 = limB3(gteMAC3, lm);
508 void gteNCLIP(psxCP2Regs *regs) {
510 GTE_LOG("GTE NCLIP\n");
514 gteMAC0 = F((s64)gteSX0 * (gteSY1 - gteSY2) +
515 gteSX1 * (gteSY2 - gteSY0) +
516 gteSX2 * (gteSY0 - gteSY1));
519 void gteAVSZ3(psxCP2Regs *regs) {
521 GTE_LOG("GTE AVSZ3\n");
525 gteMAC0 = F((s64)gteZSF3 * (gteSZ1 + gteSZ2 + gteSZ3));
526 gteOTZ = limD(gteMAC0 >> 12);
529 void gteAVSZ4(psxCP2Regs *regs) {
531 GTE_LOG("GTE AVSZ4\n");
535 gteMAC0 = F((s64)gteZSF4 * (gteSZ0 + gteSZ1 + gteSZ2 + gteSZ3));
536 gteOTZ = limD(gteMAC0 >> 12);
539 void gteSQR(psxCP2Regs *regs) {
540 int shift = 12 * GTE_SF(gteop);
541 int lm = GTE_LM(gteop);
544 GTE_LOG("GTE SQR\n");
548 gteMAC1 = (gteIR1 * gteIR1) >> shift;
549 gteMAC2 = (gteIR2 * gteIR2) >> shift;
550 gteMAC3 = (gteIR3 * gteIR3) >> shift;
551 gteIR1 = limB1(gteMAC1, lm);
552 gteIR2 = limB2(gteMAC2, lm);
553 gteIR3 = limB3(gteMAC3, lm);
556 void gteNCCS(psxCP2Regs *regs) {
558 GTE_LOG("GTE NCCS\n");
562 gteMAC1 = ((s64)(gteL11 * gteVX0) + (gteL12 * gteVY0) + (gteL13 * gteVZ0)) >> 12;
563 gteMAC2 = ((s64)(gteL21 * gteVX0) + (gteL22 * gteVY0) + (gteL23 * gteVZ0)) >> 12;
564 gteMAC3 = ((s64)(gteL31 * gteVX0) + (gteL32 * gteVY0) + (gteL33 * gteVZ0)) >> 12;
565 gteIR1 = limB1(gteMAC1, 1);
566 gteIR2 = limB2(gteMAC2, 1);
567 gteIR3 = limB3(gteMAC3, 1);
568 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
569 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
570 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
571 gteIR1 = limB1(gteMAC1, 1);
572 gteIR2 = limB2(gteMAC2, 1);
573 gteIR3 = limB3(gteMAC3, 1);
574 gteMAC1 = ((s32)gteR * gteIR1) >> 8;
575 gteMAC2 = ((s32)gteG * gteIR2) >> 8;
576 gteMAC3 = ((s32)gteB * gteIR3) >> 8;
584 gteR2 = limC1(gteMAC1 >> 4);
585 gteG2 = limC2(gteMAC2 >> 4);
586 gteB2 = limC3(gteMAC3 >> 4);
589 void gteNCCT(psxCP2Regs *regs) {
594 GTE_LOG("GTE NCCT\n");
598 for (v = 0; v < 3; v++) {
602 gteMAC1 = ((s64)(gteL11 * vx) + (gteL12 * vy) + (gteL13 * vz)) >> 12;
603 gteMAC2 = ((s64)(gteL21 * vx) + (gteL22 * vy) + (gteL23 * vz)) >> 12;
604 gteMAC3 = ((s64)(gteL31 * vx) + (gteL32 * vy) + (gteL33 * vz)) >> 12;
605 gteIR1 = limB1(gteMAC1, 1);
606 gteIR2 = limB2(gteMAC2, 1);
607 gteIR3 = limB3(gteMAC3, 1);
608 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
609 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
610 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
611 gteIR1 = limB1(gteMAC1, 1);
612 gteIR2 = limB2(gteMAC2, 1);
613 gteIR3 = limB3(gteMAC3, 1);
614 gteMAC1 = ((s32)gteR * gteIR1) >> 8;
615 gteMAC2 = ((s32)gteG * gteIR2) >> 8;
616 gteMAC3 = ((s32)gteB * gteIR3) >> 8;
621 gteR2 = limC1(gteMAC1 >> 4);
622 gteG2 = limC2(gteMAC2 >> 4);
623 gteB2 = limC3(gteMAC3 >> 4);
630 void gteNCDS(psxCP2Regs *regs) {
632 GTE_LOG("GTE NCDS\n");
636 gteMAC1 = ((s64)(gteL11 * gteVX0) + (gteL12 * gteVY0) + (gteL13 * gteVZ0)) >> 12;
637 gteMAC2 = ((s64)(gteL21 * gteVX0) + (gteL22 * gteVY0) + (gteL23 * gteVZ0)) >> 12;
638 gteMAC3 = ((s64)(gteL31 * gteVX0) + (gteL32 * gteVY0) + (gteL33 * gteVZ0)) >> 12;
639 gteIR1 = limB1(gteMAC1, 1);
640 gteIR2 = limB2(gteMAC2, 1);
641 gteIR3 = limB3(gteMAC3, 1);
642 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
643 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
644 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
645 gteIR1 = limB1(gteMAC1, 1);
646 gteIR2 = limB2(gteMAC2, 1);
647 gteIR3 = limB3(gteMAC3, 1);
648 gteMAC1 = (((gteR << 4) * gteIR1) + (gteIR0 * limB1(A1U((s64)gteRFC - ((gteR * gteIR1) >> 8)), 0))) >> 12;
649 gteMAC2 = (((gteG << 4) * gteIR2) + (gteIR0 * limB2(A2U((s64)gteGFC - ((gteG * gteIR2) >> 8)), 0))) >> 12;
650 gteMAC3 = (((gteB << 4) * gteIR3) + (gteIR0 * limB3(A3U((s64)gteBFC - ((gteB * gteIR3) >> 8)), 0))) >> 12;
651 gteIR1 = limB1(gteMAC1, 1);
652 gteIR2 = limB2(gteMAC2, 1);
653 gteIR3 = limB3(gteMAC3, 1);
658 gteR2 = limC1(gteMAC1 >> 4);
659 gteG2 = limC2(gteMAC2 >> 4);
660 gteB2 = limC3(gteMAC3 >> 4);
663 void gteNCDT(psxCP2Regs *regs) {
668 GTE_LOG("GTE NCDT\n");
672 for (v = 0; v < 3; v++) {
676 gteMAC1 = ((s64)(gteL11 * vx) + (gteL12 * vy) + (gteL13 * vz)) >> 12;
677 gteMAC2 = ((s64)(gteL21 * vx) + (gteL22 * vy) + (gteL23 * vz)) >> 12;
678 gteMAC3 = ((s64)(gteL31 * vx) + (gteL32 * vy) + (gteL33 * vz)) >> 12;
679 gteIR1 = limB1(gteMAC1, 1);
680 gteIR2 = limB2(gteMAC2, 1);
681 gteIR3 = limB3(gteMAC3, 1);
682 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
683 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
684 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
685 gteIR1 = limB1(gteMAC1, 1);
686 gteIR2 = limB2(gteMAC2, 1);
687 gteIR3 = limB3(gteMAC3, 1);
688 gteMAC1 = (((gteR << 4) * gteIR1) + (gteIR0 * limB1(A1U((s64)gteRFC - ((gteR * gteIR1) >> 8)), 0))) >> 12;
689 gteMAC2 = (((gteG << 4) * gteIR2) + (gteIR0 * limB2(A2U((s64)gteGFC - ((gteG * gteIR2) >> 8)), 0))) >> 12;
690 gteMAC3 = (((gteB << 4) * gteIR3) + (gteIR0 * limB3(A3U((s64)gteBFC - ((gteB * gteIR3) >> 8)), 0))) >> 12;
695 gteR2 = limC1(gteMAC1 >> 4);
696 gteG2 = limC2(gteMAC2 >> 4);
697 gteB2 = limC3(gteMAC3 >> 4);
699 gteIR1 = limB1(gteMAC1, 1);
700 gteIR2 = limB2(gteMAC2, 1);
701 gteIR3 = limB3(gteMAC3, 1);
704 void gteOP(psxCP2Regs *regs) {
705 int shift = 12 * GTE_SF(gteop);
706 int lm = GTE_LM(gteop);
713 gteMAC1 = ((gteR22 * gteIR3) - (gteR33 * gteIR2)) >> shift;
714 gteMAC2 = ((gteR33 * gteIR1) - (gteR11 * gteIR3)) >> shift;
715 gteMAC3 = ((gteR11 * gteIR2) - (gteR22 * gteIR1)) >> shift;
716 gteIR1 = limB1(gteMAC1, lm);
717 gteIR2 = limB2(gteMAC2, lm);
718 gteIR3 = limB3(gteMAC3, lm);
721 void gteDCPL(psxCP2Regs *regs) {
722 int lm = GTE_LM(gteop);
724 s32 RIR1 = ((s32)gteR * gteIR1) >> 8;
725 s32 GIR2 = ((s32)gteG * gteIR2) >> 8;
726 s32 BIR3 = ((s32)gteB * gteIR3) >> 8;
729 GTE_LOG("GTE DCPL\n");
733 gteMAC1 = RIR1 + ((gteIR0 * limB1(A1U((s64)gteRFC - RIR1), 0)) >> 12);
734 gteMAC2 = GIR2 + ((gteIR0 * limB1(A2U((s64)gteGFC - GIR2), 0)) >> 12);
735 gteMAC3 = BIR3 + ((gteIR0 * limB1(A3U((s64)gteBFC - BIR3), 0)) >> 12);
737 gteIR1 = limB1(gteMAC1, lm);
738 gteIR2 = limB2(gteMAC2, lm);
739 gteIR3 = limB3(gteMAC3, lm);
744 gteR2 = limC1(gteMAC1 >> 4);
745 gteG2 = limC2(gteMAC2 >> 4);
746 gteB2 = limC3(gteMAC3 >> 4);
749 void gteGPF(psxCP2Regs *regs) {
750 int shift = 12 * GTE_SF(gteop);
753 GTE_LOG("GTE GPF\n");
757 gteMAC1 = (gteIR0 * gteIR1) >> shift;
758 gteMAC2 = (gteIR0 * gteIR2) >> shift;
759 gteMAC3 = (gteIR0 * gteIR3) >> shift;
760 gteIR1 = limB1(gteMAC1, 0);
761 gteIR2 = limB2(gteMAC2, 0);
762 gteIR3 = limB3(gteMAC3, 0);
767 gteR2 = limC1(gteMAC1 >> 4);
768 gteG2 = limC2(gteMAC2 >> 4);
769 gteB2 = limC3(gteMAC3 >> 4);
772 void gteGPL(psxCP2Regs *regs) {
773 int shift = 12 * GTE_SF(gteop);
776 GTE_LOG("GTE GPL\n");
780 gteMAC1 = A1((((s64)gteMAC1 << shift) + (gteIR0 * gteIR1)) >> shift);
781 gteMAC2 = A2((((s64)gteMAC2 << shift) + (gteIR0 * gteIR2)) >> shift);
782 gteMAC3 = A3((((s64)gteMAC3 << shift) + (gteIR0 * gteIR3)) >> shift);
783 gteIR1 = limB1(gteMAC1, 0);
784 gteIR2 = limB2(gteMAC2, 0);
785 gteIR3 = limB3(gteMAC3, 0);
790 gteR2 = limC1(gteMAC1 >> 4);
791 gteG2 = limC2(gteMAC2 >> 4);
792 gteB2 = limC3(gteMAC3 >> 4);
795 void gteDPCS(psxCP2Regs *regs) {
796 int shift = 12 * GTE_SF(gteop);
799 GTE_LOG("GTE DPCS\n");
803 gteMAC1 = ((gteR << 16) + (gteIR0 * limB1(A1U(((s64)gteRFC - (gteR << 4)) << (12 - shift)), 0))) >> 12;
804 gteMAC2 = ((gteG << 16) + (gteIR0 * limB2(A2U(((s64)gteGFC - (gteG << 4)) << (12 - shift)), 0))) >> 12;
805 gteMAC3 = ((gteB << 16) + (gteIR0 * limB3(A3U(((s64)gteBFC - (gteB << 4)) << (12 - shift)), 0))) >> 12;
807 gteIR1 = limB1(gteMAC1, 0);
808 gteIR2 = limB2(gteMAC2, 0);
809 gteIR3 = limB3(gteMAC3, 0);
813 gteR2 = limC1(gteMAC1 >> 4);
814 gteG2 = limC2(gteMAC2 >> 4);
815 gteB2 = limC3(gteMAC3 >> 4);
818 void gteDPCT(psxCP2Regs *regs) {
822 GTE_LOG("GTE DPCT\n");
826 for (v = 0; v < 3; v++) {
827 gteMAC1 = ((gteR0 << 16) + (gteIR0 * limB1(A1U((s64)gteRFC - (gteR0 << 4)), 0))) >> 12;
828 gteMAC2 = ((gteG0 << 16) + (gteIR0 * limB1(A2U((s64)gteGFC - (gteG0 << 4)), 0))) >> 12;
829 gteMAC3 = ((gteB0 << 16) + (gteIR0 * limB1(A3U((s64)gteBFC - (gteB0 << 4)), 0))) >> 12;
834 gteR2 = limC1(gteMAC1 >> 4);
835 gteG2 = limC2(gteMAC2 >> 4);
836 gteB2 = limC3(gteMAC3 >> 4);
838 gteIR1 = limB1(gteMAC1, 0);
839 gteIR2 = limB2(gteMAC2, 0);
840 gteIR3 = limB3(gteMAC3, 0);
843 void gteNCS(psxCP2Regs *regs) {
845 GTE_LOG("GTE NCS\n");
849 gteMAC1 = ((s64)(gteL11 * gteVX0) + (gteL12 * gteVY0) + (gteL13 * gteVZ0)) >> 12;
850 gteMAC2 = ((s64)(gteL21 * gteVX0) + (gteL22 * gteVY0) + (gteL23 * gteVZ0)) >> 12;
851 gteMAC3 = ((s64)(gteL31 * gteVX0) + (gteL32 * gteVY0) + (gteL33 * gteVZ0)) >> 12;
852 gteIR1 = limB1(gteMAC1, 1);
853 gteIR2 = limB2(gteMAC2, 1);
854 gteIR3 = limB3(gteMAC3, 1);
855 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
856 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
857 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
858 gteIR1 = limB1(gteMAC1, 1);
859 gteIR2 = limB2(gteMAC2, 1);
860 gteIR3 = limB3(gteMAC3, 1);
865 gteR2 = limC1(gteMAC1 >> 4);
866 gteG2 = limC2(gteMAC2 >> 4);
867 gteB2 = limC3(gteMAC3 >> 4);
870 void gteNCT(psxCP2Regs *regs) {
875 GTE_LOG("GTE NCT\n");
879 for (v = 0; v < 3; v++) {
883 gteMAC1 = ((s64)(gteL11 * vx) + (gteL12 * vy) + (gteL13 * vz)) >> 12;
884 gteMAC2 = ((s64)(gteL21 * vx) + (gteL22 * vy) + (gteL23 * vz)) >> 12;
885 gteMAC3 = ((s64)(gteL31 * vx) + (gteL32 * vy) + (gteL33 * vz)) >> 12;
886 gteIR1 = limB1(gteMAC1, 1);
887 gteIR2 = limB2(gteMAC2, 1);
888 gteIR3 = limB3(gteMAC3, 1);
889 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
890 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
891 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
895 gteR2 = limC1(gteMAC1 >> 4);
896 gteG2 = limC2(gteMAC2 >> 4);
897 gteB2 = limC3(gteMAC3 >> 4);
899 gteIR1 = limB1(gteMAC1, 1);
900 gteIR2 = limB2(gteMAC2, 1);
901 gteIR3 = limB3(gteMAC3, 1);
904 void gteCC(psxCP2Regs *regs) {
910 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
911 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
912 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
913 gteIR1 = limB1(gteMAC1, 1);
914 gteIR2 = limB2(gteMAC2, 1);
915 gteIR3 = limB3(gteMAC3, 1);
916 gteMAC1 = ((s32)gteR * gteIR1) >> 8;
917 gteMAC2 = ((s32)gteG * gteIR2) >> 8;
918 gteMAC3 = ((s32)gteB * gteIR3) >> 8;
919 gteIR1 = limB1(gteMAC1, 1);
920 gteIR2 = limB2(gteMAC2, 1);
921 gteIR3 = limB3(gteMAC3, 1);
926 gteR2 = limC1(gteMAC1 >> 4);
927 gteG2 = limC2(gteMAC2 >> 4);
928 gteB2 = limC3(gteMAC3 >> 4);
931 void gteINTPL(psxCP2Regs *regs) {
932 int shift = 12 * GTE_SF(gteop);
933 int lm = GTE_LM(gteop);
936 GTE_LOG("GTE INTPL\n");
940 gteMAC1 = ((gteIR1 << 12) + (gteIR0 * limB1(A1U((s64)gteRFC - gteIR1), 0))) >> shift;
941 gteMAC2 = ((gteIR2 << 12) + (gteIR0 * limB2(A2U((s64)gteGFC - gteIR2), 0))) >> shift;
942 gteMAC3 = ((gteIR3 << 12) + (gteIR0 * limB3(A3U((s64)gteBFC - gteIR3), 0))) >> shift;
943 gteIR1 = limB1(gteMAC1, lm);
944 gteIR2 = limB2(gteMAC2, lm);
945 gteIR3 = limB3(gteMAC3, lm);
949 gteR2 = limC1(gteMAC1 >> 4);
950 gteG2 = limC2(gteMAC2 >> 4);
951 gteB2 = limC3(gteMAC3 >> 4);
954 void gteCDP(psxCP2Regs *regs) {
956 GTE_LOG("GTE CDP\n");
960 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
961 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
962 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
963 gteIR1 = limB1(gteMAC1, 1);
964 gteIR2 = limB2(gteMAC2, 1);
965 gteIR3 = limB3(gteMAC3, 1);
966 gteMAC1 = (((gteR << 4) * gteIR1) + (gteIR0 * limB1(A1U((s64)gteRFC - ((gteR * gteIR1) >> 8)), 0))) >> 12;
967 gteMAC2 = (((gteG << 4) * gteIR2) + (gteIR0 * limB2(A2U((s64)gteGFC - ((gteG * gteIR2) >> 8)), 0))) >> 12;
968 gteMAC3 = (((gteB << 4) * gteIR3) + (gteIR0 * limB3(A3U((s64)gteBFC - ((gteB * gteIR3) >> 8)), 0))) >> 12;
969 gteIR1 = limB1(gteMAC1, 1);
970 gteIR2 = limB2(gteMAC2, 1);
971 gteIR3 = limB3(gteMAC3, 1);
976 gteR2 = limC1(gteMAC1 >> 4);
977 gteG2 = limC2(gteMAC2 >> 4);
978 gteB2 = limC3(gteMAC3 >> 4);
981 /* decomposed/parametrized versions for the recompiler */
985 void gteSQR_part_noshift(psxCP2Regs *regs) {
988 gteMAC1 = gteIR1 * gteIR1;
989 gteMAC2 = gteIR2 * gteIR2;
990 gteMAC3 = gteIR3 * gteIR3;
993 void gteSQR_part_shift(psxCP2Regs *regs) {
996 gteMAC1 = (gteIR1 * gteIR1) >> 12;
997 gteMAC2 = (gteIR2 * gteIR2) >> 12;
998 gteMAC3 = (gteIR3 * gteIR3) >> 12;
1001 void gteOP_part_noshift(psxCP2Regs *regs) {
1004 gteMAC1 = (gteR22 * gteIR3) - (gteR33 * gteIR2);
1005 gteMAC2 = (gteR33 * gteIR1) - (gteR11 * gteIR3);
1006 gteMAC3 = (gteR11 * gteIR2) - (gteR22 * gteIR1);
1009 void gteOP_part_shift(psxCP2Regs *regs) {
1012 gteMAC1 = ((gteR22 * gteIR3) - (gteR33 * gteIR2)) >> 12;
1013 gteMAC2 = ((gteR33 * gteIR1) - (gteR11 * gteIR3)) >> 12;
1014 gteMAC3 = ((gteR11 * gteIR2) - (gteR22 * gteIR1)) >> 12;
1017 void gteDCPL_part(psxCP2Regs *regs) {
1018 s32 RIR1 = ((s32)gteR * gteIR1) >> 8;
1019 s32 GIR2 = ((s32)gteG * gteIR2) >> 8;
1020 s32 BIR3 = ((s32)gteB * gteIR3) >> 8;
1024 gteMAC1 = RIR1 + ((gteIR0 * limB1(A1U((s64)gteRFC - RIR1), 0)) >> 12);
1025 gteMAC2 = GIR2 + ((gteIR0 * limB1(A2U((s64)gteGFC - GIR2), 0)) >> 12);
1026 gteMAC3 = BIR3 + ((gteIR0 * limB1(A3U((s64)gteBFC - BIR3), 0)) >> 12);
1029 void gteGPF_part_noshift(psxCP2Regs *regs) {
1032 gteMAC1 = gteIR0 * gteIR1;
1033 gteMAC2 = gteIR0 * gteIR2;
1034 gteMAC3 = gteIR0 * gteIR3;
1037 void gteGPF_part_shift(psxCP2Regs *regs) {
1040 gteMAC1 = (gteIR0 * gteIR1) >> 12;
1041 gteMAC2 = (gteIR0 * gteIR2) >> 12;
1042 gteMAC3 = (gteIR0 * gteIR3) >> 12;
1047 void gteGPL_part_noshift(psxCP2Regs *regs) {
1050 gteMAC1 = A1((s64)gteMAC1 + (gteIR0 * gteIR1));
1051 gteMAC2 = A2((s64)gteMAC2 + (gteIR0 * gteIR2));
1052 gteMAC3 = A3((s64)gteMAC3 + (gteIR0 * gteIR3));
1055 void gteGPL_part_shift(psxCP2Regs *regs) {
1058 gteMAC1 = A1((s64)gteMAC1 + ((gteIR0 * gteIR1) >> 12));
1059 gteMAC2 = A2((s64)gteMAC2 + ((gteIR0 * gteIR2) >> 12));
1060 gteMAC3 = A3((s64)gteMAC3 + ((gteIR0 * gteIR3) >> 12));
1063 void gteDPCS_part_noshift(psxCP2Regs *regs) {
1068 gteMAC1 = ((gteR << 16) + (gteIR0 * limB1(A1U((s64)gteRFC - (gteR << 4)) << (12 - shift), 0))) >> 12;
1069 gteMAC2 = ((gteG << 16) + (gteIR0 * limB2(A2U((s64)gteGFC - (gteG << 4)) << (12 - shift), 0))) >> 12;
1070 gteMAC3 = ((gteB << 16) + (gteIR0 * limB3(A3U((s64)gteBFC - (gteB << 4)) << (12 - shift), 0))) >> 12;
1073 void gteDPCS_part_shift(psxCP2Regs *regs) {
1078 gteMAC1 = ((gteR << 16) + (gteIR0 * limB1(A1U((s64)gteRFC - (gteR << 4)) << (12 - shift), 0))) >> 12;
1079 gteMAC2 = ((gteG << 16) + (gteIR0 * limB2(A2U((s64)gteGFC - (gteG << 4)) << (12 - shift), 0))) >> 12;
1080 gteMAC3 = ((gteB << 16) + (gteIR0 * limB3(A3U((s64)gteBFC - (gteB << 4)) << (12 - shift), 0))) >> 12;
1083 void gteINTPL_part_noshift(psxCP2Regs *regs) {
1086 gteMAC1 = ((gteIR1 << 12) + (gteIR0 * limB1(A1U((s64)gteRFC - gteIR1), 0)));
1087 gteMAC2 = ((gteIR2 << 12) + (gteIR0 * limB2(A2U((s64)gteGFC - gteIR2), 0)));
1088 gteMAC3 = ((gteIR3 << 12) + (gteIR0 * limB3(A3U((s64)gteBFC - gteIR3), 0)));
1091 void gteINTPL_part_shift(psxCP2Regs *regs) {
1094 gteMAC1 = ((gteIR1 << 12) + (gteIR0 * limB1(A1U((s64)gteRFC - gteIR1), 0))) >> 12;
1095 gteMAC2 = ((gteIR2 << 12) + (gteIR0 * limB2(A2U((s64)gteGFC - gteIR2), 0))) >> 12;
1096 gteMAC3 = ((gteIR3 << 12) + (gteIR0 * limB3(A3U((s64)gteBFC - gteIR3), 0))) >> 12;
1099 void gteMACtoRGB(psxCP2Regs *regs) {
1103 gteR2 = limC1(gteMAC1 >> 4);
1104 gteG2 = limC2(gteMAC2 >> 4);
1105 gteB2 = limC3(gteMAC3 >> 4);