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 ? psxRegs.CP2D.p[n << 1].sw.l : psxRegs.CP2D.p[9].sw.l)
30 #define VY(n) (n < 3 ? psxRegs.CP2D.p[n << 1].sw.h : psxRegs.CP2D.p[10].sw.l)
31 #define VZ(n) (n < 3 ? psxRegs.CP2D.p[(n << 1) + 1].sw.l : psxRegs.CP2D.p[11].sw.l)
32 #define MX11(n) (n < 3 ? psxRegs.CP2C.p[(n << 3)].sw.l : 0)
33 #define MX12(n) (n < 3 ? psxRegs.CP2C.p[(n << 3)].sw.h : 0)
34 #define MX13(n) (n < 3 ? psxRegs.CP2C.p[(n << 3) + 1].sw.l : 0)
35 #define MX21(n) (n < 3 ? psxRegs.CP2C.p[(n << 3) + 1].sw.h : 0)
36 #define MX22(n) (n < 3 ? psxRegs.CP2C.p[(n << 3) + 2].sw.l : 0)
37 #define MX23(n) (n < 3 ? psxRegs.CP2C.p[(n << 3) + 2].sw.h : 0)
38 #define MX31(n) (n < 3 ? psxRegs.CP2C.p[(n << 3) + 3].sw.l : 0)
39 #define MX32(n) (n < 3 ? psxRegs.CP2C.p[(n << 3) + 3].sw.h : 0)
40 #define MX33(n) (n < 3 ? psxRegs.CP2C.p[(n << 3) + 4].sw.l : 0)
41 #define CV1(n) (n < 3 ? (s32)psxRegs.CP2C.r[(n << 3) + 5] : 0)
42 #define CV2(n) (n < 3 ? (s32)psxRegs.CP2C.r[(n << 3) + 6] : 0)
43 #define CV3(n) (n < 3 ? (s32)psxRegs.CP2C.r[(n << 3) + 7] : 0)
45 #define fSX(n) ((psxRegs.CP2D.p)[((n) + 12)].sw.l)
46 #define fSY(n) ((psxRegs.CP2D.p)[((n) + 12)].sw.h)
47 #define fSZ(n) ((psxRegs.CP2D.p)[((n) + 17)].w.l) /* (n == 0) => SZ1; */
49 #define gteVXY0 (psxRegs.CP2D.r[0])
50 #define gteVX0 (psxRegs.CP2D.p[0].sw.l)
51 #define gteVY0 (psxRegs.CP2D.p[0].sw.h)
52 #define gteVZ0 (psxRegs.CP2D.p[1].sw.l)
53 #define gteVXY1 (psxRegs.CP2D.r[2])
54 #define gteVX1 (psxRegs.CP2D.p[2].sw.l)
55 #define gteVY1 (psxRegs.CP2D.p[2].sw.h)
56 #define gteVZ1 (psxRegs.CP2D.p[3].sw.l)
57 #define gteVXY2 (psxRegs.CP2D.r[4])
58 #define gteVX2 (psxRegs.CP2D.p[4].sw.l)
59 #define gteVY2 (psxRegs.CP2D.p[4].sw.h)
60 #define gteVZ2 (psxRegs.CP2D.p[5].sw.l)
61 #define gteRGB (psxRegs.CP2D.r[6])
62 #define gteR (psxRegs.CP2D.p[6].b.l)
63 #define gteG (psxRegs.CP2D.p[6].b.h)
64 #define gteB (psxRegs.CP2D.p[6].b.h2)
65 #define gteCODE (psxRegs.CP2D.p[6].b.h3)
66 #define gteOTZ (psxRegs.CP2D.p[7].w.l)
67 #define gteIR0 (psxRegs.CP2D.p[8].sw.l)
68 #define gteIR1 (psxRegs.CP2D.p[9].sw.l)
69 #define gteIR2 (psxRegs.CP2D.p[10].sw.l)
70 #define gteIR3 (psxRegs.CP2D.p[11].sw.l)
71 #define gteSXY0 (psxRegs.CP2D.r[12])
72 #define gteSX0 (psxRegs.CP2D.p[12].sw.l)
73 #define gteSY0 (psxRegs.CP2D.p[12].sw.h)
74 #define gteSXY1 (psxRegs.CP2D.r[13])
75 #define gteSX1 (psxRegs.CP2D.p[13].sw.l)
76 #define gteSY1 (psxRegs.CP2D.p[13].sw.h)
77 #define gteSXY2 (psxRegs.CP2D.r[14])
78 #define gteSX2 (psxRegs.CP2D.p[14].sw.l)
79 #define gteSY2 (psxRegs.CP2D.p[14].sw.h)
80 #define gteSXYP (psxRegs.CP2D.r[15])
81 #define gteSXP (psxRegs.CP2D.p[15].sw.l)
82 #define gteSYP (psxRegs.CP2D.p[15].sw.h)
83 #define gteSZ0 (psxRegs.CP2D.p[16].w.l)
84 #define gteSZ1 (psxRegs.CP2D.p[17].w.l)
85 #define gteSZ2 (psxRegs.CP2D.p[18].w.l)
86 #define gteSZ3 (psxRegs.CP2D.p[19].w.l)
87 #define gteRGB0 (psxRegs.CP2D.r[20])
88 #define gteR0 (psxRegs.CP2D.p[20].b.l)
89 #define gteG0 (psxRegs.CP2D.p[20].b.h)
90 #define gteB0 (psxRegs.CP2D.p[20].b.h2)
91 #define gteCODE0 (psxRegs.CP2D.p[20].b.h3)
92 #define gteRGB1 (psxRegs.CP2D.r[21])
93 #define gteR1 (psxRegs.CP2D.p[21].b.l)
94 #define gteG1 (psxRegs.CP2D.p[21].b.h)
95 #define gteB1 (psxRegs.CP2D.p[21].b.h2)
96 #define gteCODE1 (psxRegs.CP2D.p[21].b.h3)
97 #define gteRGB2 (psxRegs.CP2D.r[22])
98 #define gteR2 (psxRegs.CP2D.p[22].b.l)
99 #define gteG2 (psxRegs.CP2D.p[22].b.h)
100 #define gteB2 (psxRegs.CP2D.p[22].b.h2)
101 #define gteCODE2 (psxRegs.CP2D.p[22].b.h3)
102 #define gteRES1 (psxRegs.CP2D.r[23])
103 #define gteMAC0 (((s32 *)psxRegs.CP2D.r)[24])
104 #define gteMAC1 (((s32 *)psxRegs.CP2D.r)[25])
105 #define gteMAC2 (((s32 *)psxRegs.CP2D.r)[26])
106 #define gteMAC3 (((s32 *)psxRegs.CP2D.r)[27])
107 #define gteIRGB (psxRegs.CP2D.r[28])
108 #define gteORGB (psxRegs.CP2D.r[29])
109 #define gteLZCS (psxRegs.CP2D.r[30])
110 #define gteLZCR (psxRegs.CP2D.r[31])
112 #define gteR11R12 (((s32 *)psxRegs.CP2C.r)[0])
113 #define gteR22R23 (((s32 *)psxRegs.CP2C.r)[2])
114 #define gteR11 (psxRegs.CP2C.p[0].sw.l)
115 #define gteR12 (psxRegs.CP2C.p[0].sw.h)
116 #define gteR13 (psxRegs.CP2C.p[1].sw.l)
117 #define gteR21 (psxRegs.CP2C.p[1].sw.h)
118 #define gteR22 (psxRegs.CP2C.p[2].sw.l)
119 #define gteR23 (psxRegs.CP2C.p[2].sw.h)
120 #define gteR31 (psxRegs.CP2C.p[3].sw.l)
121 #define gteR32 (psxRegs.CP2C.p[3].sw.h)
122 #define gteR33 (psxRegs.CP2C.p[4].sw.l)
123 #define gteTRX (((s32 *)psxRegs.CP2C.r)[5])
124 #define gteTRY (((s32 *)psxRegs.CP2C.r)[6])
125 #define gteTRZ (((s32 *)psxRegs.CP2C.r)[7])
126 #define gteL11 (psxRegs.CP2C.p[8].sw.l)
127 #define gteL12 (psxRegs.CP2C.p[8].sw.h)
128 #define gteL13 (psxRegs.CP2C.p[9].sw.l)
129 #define gteL21 (psxRegs.CP2C.p[9].sw.h)
130 #define gteL22 (psxRegs.CP2C.p[10].sw.l)
131 #define gteL23 (psxRegs.CP2C.p[10].sw.h)
132 #define gteL31 (psxRegs.CP2C.p[11].sw.l)
133 #define gteL32 (psxRegs.CP2C.p[11].sw.h)
134 #define gteL33 (psxRegs.CP2C.p[12].sw.l)
135 #define gteRBK (((s32 *)psxRegs.CP2C.r)[13])
136 #define gteGBK (((s32 *)psxRegs.CP2C.r)[14])
137 #define gteBBK (((s32 *)psxRegs.CP2C.r)[15])
138 #define gteLR1 (psxRegs.CP2C.p[16].sw.l)
139 #define gteLR2 (psxRegs.CP2C.p[16].sw.h)
140 #define gteLR3 (psxRegs.CP2C.p[17].sw.l)
141 #define gteLG1 (psxRegs.CP2C.p[17].sw.h)
142 #define gteLG2 (psxRegs.CP2C.p[18].sw.l)
143 #define gteLG3 (psxRegs.CP2C.p[18].sw.h)
144 #define gteLB1 (psxRegs.CP2C.p[19].sw.l)
145 #define gteLB2 (psxRegs.CP2C.p[19].sw.h)
146 #define gteLB3 (psxRegs.CP2C.p[20].sw.l)
147 #define gteRFC (((s32 *)psxRegs.CP2C.r)[21])
148 #define gteGFC (((s32 *)psxRegs.CP2C.r)[22])
149 #define gteBFC (((s32 *)psxRegs.CP2C.r)[23])
150 #define gteOFX (((s32 *)psxRegs.CP2C.r)[24])
151 #define gteOFY (((s32 *)psxRegs.CP2C.r)[25])
152 #define gteH (psxRegs.CP2C.p[26].sw.l)
153 #define gteDQA (psxRegs.CP2C.p[27].sw.l)
154 #define gteDQB (((s32 *)psxRegs.CP2C.r)[28])
155 #define gteZSF3 (psxRegs.CP2C.p[29].sw.l)
156 #define gteZSF4 (psxRegs.CP2C.p[30].sw.l)
157 #define gteFLAG (psxRegs.CP2C.r[31])
159 #define GTE_OP(op) ((op >> 20) & 31)
160 #define GTE_SF(op) ((op >> 19) & 1)
161 #define GTE_MX(op) ((op >> 17) & 3)
162 #define GTE_V(op) ((op >> 15) & 3)
163 #define GTE_CV(op) ((op >> 13) & 3)
164 #define GTE_CD(op) ((op >> 11) & 3) /* not used */
165 #define GTE_LM(op) ((op >> 10) & 1)
166 #define GTE_CT(op) ((op >> 6) & 15) /* not used */
167 #define GTE_FUNCT(op) (op & 63)
169 #define gteop (psxRegs.code & 0x1ffffff)
173 static inline s64 BOUNDS(s64 n_value, s64 n_max, int n_maxflag, s64 n_min, int n_minflag) {
174 if (n_value > n_max) {
175 gteFLAG |= n_maxflag;
176 } else if (n_value < n_min) {
177 gteFLAG |= n_minflag;
182 static inline s32 LIM(s32 value, s32 max, s32 min, u32 flag) {
187 } else if (value < min) {
194 static inline u32 limE(u32 result) {
195 if (result > 0x1ffff) {
196 gteFLAG |= (1 << 31) | (1 << 17);
204 #define BOUNDS(a, ...) (a)
206 static inline s32 LIM(s32 value, s32 max, s32 min, u32 flag_unused) {
210 else if (value < min)
215 static inline u32 limE(u32 result) {
216 if (result > 0x1ffff)
223 #define A1(a) BOUNDS((a), 0x7fffffff, (1 << 30), -(s64)0x80000000, (1 << 31) | (1 << 27))
224 #define A2(a) BOUNDS((a), 0x7fffffff, (1 << 29), -(s64)0x80000000, (1 << 31) | (1 << 26))
225 #define A3(a) BOUNDS((a), 0x7fffffff, (1 << 28), -(s64)0x80000000, (1 << 31) | (1 << 25))
226 #define limB1(a, l) LIM((a), 0x7fff, -0x8000 * !l, (1 << 31) | (1 << 24))
227 #define limB2(a, l) LIM((a), 0x7fff, -0x8000 * !l, (1 << 31) | (1 << 23))
228 #define limB3(a, l) LIM((a), 0x7fff, -0x8000 * !l, (1 << 22))
229 #define limC1(a) LIM((a), 0x00ff, 0x0000, (1 << 21))
230 #define limC2(a) LIM((a), 0x00ff, 0x0000, (1 << 20))
231 #define limC3(a) LIM((a), 0x00ff, 0x0000, (1 << 19))
232 #define limD(a) LIM((a), 0xffff, 0x0000, (1 << 31) | (1 << 18))
234 #define F(a) BOUNDS((a), 0x7fffffff, (1 << 31) | (1 << 16), -(s64)0x80000000, (1 << 31) | (1 << 15))
235 #define limG1(a) LIM((a), 0x3ff, -0x400, (1 << 31) | (1 << 14))
236 #define limG2(a) LIM((a), 0x3ff, -0x400, (1 << 31) | (1 << 13))
237 #define limH(a) LIM((a), 0x1000, 0x0000, (1 << 12))
239 #include "gte_divider.h"
243 static inline u32 MFC2(int reg) {
252 psxRegs.CP2D.r[reg] = (s32)psxRegs.CP2D.p[reg].sw.l;
260 psxRegs.CP2D.r[reg] = (u32)psxRegs.CP2D.p[reg].w.l;
264 psxRegs.CP2D.r[reg] = gteSXY2;
269 psxRegs.CP2D.r[reg] = LIM(gteIR1 >> 7, 0x1f, 0, 0) |
270 (LIM(gteIR2 >> 7, 0x1f, 0, 0) << 5) |
271 (LIM(gteIR3 >> 7, 0x1f, 0, 0) << 10);
274 return psxRegs.CP2D.r[reg];
277 static inline void MTC2(u32 value, int reg) {
288 gteIR1 = (value & 0x1f) << 7;
289 gteIR2 = (value & 0x3e0) << 2;
290 gteIR3 = (value & 0x7c00) >> 3;
301 for (i = 31; (a & (1 << i)) == 0 && i >= 0; i--);
306 for (i = 31; (a & (1 << i)) == 0 && i >= 0; i--);
318 psxRegs.CP2D.r[reg] = value;
322 static inline void CTC2(u32 value, int reg) {
331 value = (s32)(s16)value;
335 value = value & 0x7ffff000;
336 if (value & 0x7f87e000) value |= 0x80000000;
340 psxRegs.CP2C.r[reg] = value;
345 psxRegs.GPR.r[_Rt_] = MFC2(_Rd_);
350 psxRegs.GPR.r[_Rt_] = psxRegs.CP2C.r[_Rd_];
354 MTC2(psxRegs.GPR.r[_Rt_], _Rd_);
358 CTC2(psxRegs.GPR.r[_Rt_], _Rd_);
361 #define _oB_ (psxRegs.GPR.r[_Rs_] + _Imm_)
364 MTC2(psxMemRead32(_oB_), _Rt_);
368 psxMemWrite32(_oB_, MFC2(_Rt_));
374 #define DIVIDE DIVIDE_
375 static u32 DIVIDE_(s16 n, u16 d) {
376 if (n >= 0 && n < d * 2) {
378 return ((n_ << 16) + d / 2) / d;
379 //return (u32)((float)(n_ << 16) / (float)d + (float)0.5);
389 GTE_LOG("GTE RTPS\n");
393 gteMAC1 = A1((((s64)gteTRX << 12) + (gteR11 * gteVX0) + (gteR12 * gteVY0) + (gteR13 * gteVZ0)) >> 12);
394 gteMAC2 = A2((((s64)gteTRY << 12) + (gteR21 * gteVX0) + (gteR22 * gteVY0) + (gteR23 * gteVZ0)) >> 12);
395 gteMAC3 = A3((((s64)gteTRZ << 12) + (gteR31 * gteVX0) + (gteR32 * gteVY0) + (gteR33 * gteVZ0)) >> 12);
396 gteIR1 = limB1(gteMAC1, 0);
397 gteIR2 = limB2(gteMAC2, 0);
398 gteIR3 = limB3(gteMAC3, 0);
402 gteSZ3 = limD(gteMAC3);
403 quotient = limE(DIVIDE(gteH, gteSZ3));
406 gteSX2 = limG1(F((s64)gteOFX + ((s64)gteIR1 * quotient)) >> 16);
407 gteSY2 = limG2(F((s64)gteOFY + ((s64)gteIR2 * quotient)) >> 16);
409 gteMAC0 = F((s64)gteDQB + ((s64)gteDQA * quotient));
410 gteIR0 = limH(gteMAC0 >> 12);
419 GTE_LOG("GTE RTPT\n");
424 for (v = 0; v < 3; v++) {
428 gteMAC1 = A1((((s64)gteTRX << 12) + (gteR11 * vx) + (gteR12 * vy) + (gteR13 * vz)) >> 12);
429 gteMAC2 = A2((((s64)gteTRY << 12) + (gteR21 * vx) + (gteR22 * vy) + (gteR23 * vz)) >> 12);
430 gteMAC3 = A3((((s64)gteTRZ << 12) + (gteR31 * vx) + (gteR32 * vy) + (gteR33 * vz)) >> 12);
431 gteIR1 = limB1(gteMAC1, 0);
432 gteIR2 = limB2(gteMAC2, 0);
433 gteIR3 = limB3(gteMAC3, 0);
434 fSZ(v) = limD(gteMAC3);
435 quotient = limE(DIVIDE(gteH, fSZ(v)));
436 fSX(v) = limG1(F((s64)gteOFX + ((s64)gteIR1 * quotient)) >> 16);
437 fSY(v) = limG2(F((s64)gteOFY + ((s64)gteIR2 * quotient)) >> 16);
439 gteMAC0 = F((s64)gteDQB + ((s64)gteDQA * quotient));
440 gteIR0 = limH(gteMAC0 >> 12);
444 int shift = 12 * GTE_SF(gteop);
445 int mx = GTE_MX(gteop);
446 int v = GTE_V(gteop);
447 int cv = GTE_CV(gteop);
448 int lm = GTE_LM(gteop);
454 GTE_LOG("GTE MVMVA\n");
458 gteMAC1 = A1((((s64)CV1(cv) << 12) + (MX11(mx) * vx) + (MX12(mx) * vy) + (MX13(mx) * vz)) >> shift);
459 gteMAC2 = A2((((s64)CV2(cv) << 12) + (MX21(mx) * vx) + (MX22(mx) * vy) + (MX23(mx) * vz)) >> shift);
460 gteMAC3 = A3((((s64)CV3(cv) << 12) + (MX31(mx) * vx) + (MX32(mx) * vy) + (MX33(mx) * vz)) >> shift);
462 gteIR1 = limB1(gteMAC1, lm);
463 gteIR2 = limB2(gteMAC2, lm);
464 gteIR3 = limB3(gteMAC3, lm);
469 GTE_LOG("GTE NCLIP\n");
473 gteMAC0 = F((s64)gteSX0 * (gteSY1 - gteSY2) +
474 gteSX1 * (gteSY2 - gteSY0) +
475 gteSX2 * (gteSY0 - gteSY1));
480 GTE_LOG("GTE AVSZ3\n");
484 gteMAC0 = F((s64)(gteZSF3 * gteSZ1) + (gteZSF3 * gteSZ2) + (gteZSF3 * gteSZ3));
485 gteOTZ = limD(gteMAC0 >> 12);
490 GTE_LOG("GTE AVSZ4\n");
494 gteMAC0 = F((s64)(gteZSF4 * (gteSZ0 + gteSZ1 + gteSZ2 + gteSZ3)));
495 gteOTZ = limD(gteMAC0 >> 12);
499 int shift = 12 * GTE_SF(gteop);
500 int lm = GTE_LM(gteop);
503 GTE_LOG("GTE SQR\n");
507 gteMAC1 = A1((gteIR1 * gteIR1) >> shift);
508 gteMAC2 = A2((gteIR2 * gteIR2) >> shift);
509 gteMAC3 = A3((gteIR3 * gteIR3) >> shift);
510 gteIR1 = limB1(gteMAC1, lm);
511 gteIR2 = limB2(gteMAC2, lm);
512 gteIR3 = limB3(gteMAC3, lm);
517 GTE_LOG("GTE NCCS\n");
521 gteMAC1 = A1((((s64)gteL11 * gteVX0) + (gteL12 * gteVY0) + (gteL13 * gteVZ0)) >> 12);
522 gteMAC2 = A2((((s64)gteL21 * gteVX0) + (gteL22 * gteVY0) + (gteL23 * gteVZ0)) >> 12);
523 gteMAC3 = A3((((s64)gteL31 * gteVX0) + (gteL32 * gteVY0) + (gteL33 * gteVZ0)) >> 12);
524 gteIR1 = limB1(gteMAC1, 1);
525 gteIR2 = limB2(gteMAC2, 1);
526 gteIR3 = limB3(gteMAC3, 1);
527 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
528 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
529 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
530 gteIR1 = limB1(gteMAC1, 1);
531 gteIR2 = limB2(gteMAC2, 1);
532 gteIR3 = limB3(gteMAC3, 1);
533 gteMAC1 = A1(((s64)gteR * gteIR1) >> 8);
534 gteMAC2 = A2(((s64)gteG * gteIR2) >> 8);
535 gteMAC3 = A3(((s64)gteB * gteIR3) >> 8);
536 gteIR1 = limB1(gteMAC1, 1);
537 gteIR2 = limB2(gteMAC2, 1);
538 gteIR3 = limB3(gteMAC3, 1);
543 gteR2 = limC1(gteMAC1 >> 4);
544 gteG2 = limC2(gteMAC2 >> 4);
545 gteB2 = limC3(gteMAC3 >> 4);
553 GTE_LOG("GTE NCCT\n");
557 for (v = 0; v < 3; v++) {
561 gteMAC1 = A1((((s64)gteL11 * vx) + (gteL12 * vy) + (gteL13 * vz)) >> 12);
562 gteMAC2 = A2((((s64)gteL21 * vx) + (gteL22 * vy) + (gteL23 * vz)) >> 12);
563 gteMAC3 = A3((((s64)gteL31 * vx) + (gteL32 * vy) + (gteL33 * vz)) >> 12);
564 gteIR1 = limB1(gteMAC1, 1);
565 gteIR2 = limB2(gteMAC2, 1);
566 gteIR3 = limB3(gteMAC3, 1);
567 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
568 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
569 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
570 gteIR1 = limB1(gteMAC1, 1);
571 gteIR2 = limB2(gteMAC2, 1);
572 gteIR3 = limB3(gteMAC3, 1);
573 gteMAC1 = A1(((s64)gteR * gteIR1) >> 8);
574 gteMAC2 = A2(((s64)gteG * gteIR2) >> 8);
575 gteMAC3 = A3(((s64)gteB * gteIR3) >> 8);
580 gteR2 = limC1(gteMAC1 >> 4);
581 gteG2 = limC2(gteMAC2 >> 4);
582 gteB2 = limC3(gteMAC3 >> 4);
584 gteIR1 = limB1(gteMAC1, 1);
585 gteIR2 = limB2(gteMAC2, 1);
586 gteIR3 = limB3(gteMAC3, 1);
591 GTE_LOG("GTE NCDS\n");
595 gteMAC1 = A1((((s64)gteL11 * gteVX0) + (gteL12 * gteVY0) + (gteL13 * gteVZ0)) >> 12);
596 gteMAC2 = A2((((s64)gteL21 * gteVX0) + (gteL22 * gteVY0) + (gteL23 * gteVZ0)) >> 12);
597 gteMAC3 = A3((((s64)gteL31 * gteVX0) + (gteL32 * gteVY0) + (gteL33 * gteVZ0)) >> 12);
598 gteIR1 = limB1(gteMAC1, 1);
599 gteIR2 = limB2(gteMAC2, 1);
600 gteIR3 = limB3(gteMAC3, 1);
601 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
602 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
603 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
604 gteIR1 = limB1(gteMAC1, 1);
605 gteIR2 = limB2(gteMAC2, 1);
606 gteIR3 = limB3(gteMAC3, 1);
607 gteMAC1 = A1(((((s64)gteR << 4) * gteIR1) + (gteIR0 * limB1(gteRFC - ((gteR * gteIR1) >> 8), 0))) >> 12);
608 gteMAC2 = A2(((((s64)gteG << 4) * gteIR2) + (gteIR0 * limB2(gteGFC - ((gteG * gteIR2) >> 8), 0))) >> 12);
609 gteMAC3 = A3(((((s64)gteB << 4) * gteIR3) + (gteIR0 * limB3(gteBFC - ((gteB * gteIR3) >> 8), 0))) >> 12);
610 gteIR1 = limB1(gteMAC1, 1);
611 gteIR2 = limB2(gteMAC2, 1);
612 gteIR3 = limB3(gteMAC3, 1);
617 gteR2 = limC1(gteMAC1 >> 4);
618 gteG2 = limC2(gteMAC2 >> 4);
619 gteB2 = limC3(gteMAC3 >> 4);
627 GTE_LOG("GTE NCDT\n");
631 for (v = 0; v < 3; v++) {
635 gteMAC1 = A1((((s64)gteL11 * vx) + (gteL12 * vy) + (gteL13 * vz)) >> 12);
636 gteMAC2 = A2((((s64)gteL21 * vx) + (gteL22 * vy) + (gteL23 * vz)) >> 12);
637 gteMAC3 = A3((((s64)gteL31 * vx) + (gteL32 * vy) + (gteL33 * vz)) >> 12);
638 gteIR1 = limB1(gteMAC1, 1);
639 gteIR2 = limB2(gteMAC2, 1);
640 gteIR3 = limB3(gteMAC3, 1);
641 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
642 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
643 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
644 gteIR1 = limB1(gteMAC1, 1);
645 gteIR2 = limB2(gteMAC2, 1);
646 gteIR3 = limB3(gteMAC3, 1);
647 gteMAC1 = A1(((((s64)gteR << 4) * gteIR1) + (gteIR0 * limB1(gteRFC - ((gteR * gteIR1) >> 8), 0))) >> 12);
648 gteMAC2 = A2(((((s64)gteG << 4) * gteIR2) + (gteIR0 * limB2(gteGFC - ((gteG * gteIR2) >> 8), 0))) >> 12);
649 gteMAC3 = A3(((((s64)gteB << 4) * gteIR3) + (gteIR0 * limB3(gteBFC - ((gteB * gteIR3) >> 8), 0))) >> 12);
654 gteR2 = limC1(gteMAC1 >> 4);
655 gteG2 = limC2(gteMAC2 >> 4);
656 gteB2 = limC3(gteMAC3 >> 4);
658 gteIR1 = limB1(gteMAC1, 1);
659 gteIR2 = limB2(gteMAC2, 1);
660 gteIR3 = limB3(gteMAC3, 1);
664 int shift = 12 * GTE_SF(gteop);
665 int lm = GTE_LM(gteop);
672 gteMAC1 = A1(((s64)(gteR22 * gteIR3) - (gteR33 * gteIR2)) >> shift);
673 gteMAC2 = A2(((s64)(gteR33 * gteIR1) - (gteR11 * gteIR3)) >> shift);
674 gteMAC3 = A3(((s64)(gteR11 * gteIR2) - (gteR22 * gteIR1)) >> shift);
675 gteIR1 = limB1(gteMAC1, lm);
676 gteIR2 = limB2(gteMAC2, lm);
677 gteIR3 = limB3(gteMAC3, lm);
681 int lm = GTE_LM(gteop);
683 s64 RIR1 = ((s64)gteR * gteIR1) >> 8;
684 s64 GIR2 = ((s64)gteG * gteIR2) >> 8;
685 s64 BIR3 = ((s64)gteB * gteIR3) >> 8;
688 GTE_LOG("GTE DCPL\n");
692 gteMAC1 = A1(RIR1 + ((gteIR0 * limB1(gteRFC - RIR1, 0)) >> 12));
693 gteMAC2 = A2(GIR2 + ((gteIR0 * limB1(gteGFC - GIR2, 0)) >> 12));
694 gteMAC3 = A3(BIR3 + ((gteIR0 * limB1(gteBFC - BIR3, 0)) >> 12));
696 gteIR1 = limB1(gteMAC1, lm);
697 gteIR2 = limB2(gteMAC2, lm);
698 gteIR3 = limB3(gteMAC3, lm);
703 gteR2 = limC1(gteMAC1 >> 4);
704 gteG2 = limC2(gteMAC2 >> 4);
705 gteB2 = limC3(gteMAC3 >> 4);
709 int shift = 12 * GTE_SF(gteop);
712 GTE_LOG("GTE GPF\n");
716 gteMAC1 = A1(((s64)gteIR0 * gteIR1) >> shift);
717 gteMAC2 = A2(((s64)gteIR0 * gteIR2) >> shift);
718 gteMAC3 = A3(((s64)gteIR0 * gteIR3) >> shift);
719 gteIR1 = limB1(gteMAC1, 0);
720 gteIR2 = limB2(gteMAC2, 0);
721 gteIR3 = limB3(gteMAC3, 0);
726 gteR2 = limC1(gteMAC1 >> 4);
727 gteG2 = limC2(gteMAC2 >> 4);
728 gteB2 = limC3(gteMAC3 >> 4);
732 int shift = 12 * GTE_SF(gteop);
735 GTE_LOG("GTE GPL\n");
739 gteMAC1 = A1((((s64)gteMAC1 << shift) + (gteIR0 * gteIR1)) >> shift);
740 gteMAC2 = A2((((s64)gteMAC2 << shift) + (gteIR0 * gteIR2)) >> shift);
741 gteMAC3 = A3((((s64)gteMAC3 << shift) + (gteIR0 * gteIR3)) >> shift);
742 gteIR1 = limB1(gteMAC1, 0);
743 gteIR2 = limB2(gteMAC2, 0);
744 gteIR3 = limB3(gteMAC3, 0);
749 gteR2 = limC1(gteMAC1 >> 4);
750 gteG2 = limC2(gteMAC2 >> 4);
751 gteB2 = limC3(gteMAC3 >> 4);
755 int shift = 12 * GTE_SF(gteop);
758 GTE_LOG("GTE DPCS\n");
762 gteMAC1 = A1(((gteR << 16) + (gteIR0 * limB1(A1((s64)gteRFC - (gteR << 4)) << (12 - shift), 0))) >> 12);
763 gteMAC2 = A2(((gteG << 16) + (gteIR0 * limB2(A2((s64)gteGFC - (gteG << 4)) << (12 - shift), 0))) >> 12);
764 gteMAC3 = A3(((gteB << 16) + (gteIR0 * limB3(A3((s64)gteBFC - (gteB << 4)) << (12 - shift), 0))) >> 12);
766 gteIR1 = limB1(gteMAC1, 0);
767 gteIR2 = limB2(gteMAC2, 0);
768 gteIR3 = limB3(gteMAC3, 0);
772 gteR2 = limC1(gteMAC1 >> 4);
773 gteG2 = limC2(gteMAC2 >> 4);
774 gteB2 = limC3(gteMAC3 >> 4);
781 GTE_LOG("GTE DPCT\n");
785 for (v = 0; v < 3; v++) {
786 gteMAC1 = A1((((s64)gteR0 << 16) + ((s64)gteIR0 * (limB1(gteRFC - (gteR0 << 4), 0)))) >> 12);
787 gteMAC2 = A2((((s64)gteG0 << 16) + ((s64)gteIR0 * (limB1(gteGFC - (gteG0 << 4), 0)))) >> 12);
788 gteMAC3 = A3((((s64)gteB0 << 16) + ((s64)gteIR0 * (limB1(gteBFC - (gteB0 << 4), 0)))) >> 12);
793 gteR2 = limC1(gteMAC1 >> 4);
794 gteG2 = limC2(gteMAC2 >> 4);
795 gteB2 = limC3(gteMAC3 >> 4);
797 gteIR1 = limB1(gteMAC1, 0);
798 gteIR2 = limB2(gteMAC2, 0);
799 gteIR3 = limB3(gteMAC3, 0);
804 GTE_LOG("GTE NCS\n");
808 gteMAC1 = A1((((s64)gteL11 * gteVX0) + (gteL12 * gteVY0) + (gteL13 * gteVZ0)) >> 12);
809 gteMAC2 = A2((((s64)gteL21 * gteVX0) + (gteL22 * gteVY0) + (gteL23 * gteVZ0)) >> 12);
810 gteMAC3 = A3((((s64)gteL31 * gteVX0) + (gteL32 * gteVY0) + (gteL33 * gteVZ0)) >> 12);
811 gteIR1 = limB1(gteMAC1, 1);
812 gteIR2 = limB2(gteMAC2, 1);
813 gteIR3 = limB3(gteMAC3, 1);
814 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
815 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
816 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
817 gteIR1 = limB1(gteMAC1, 1);
818 gteIR2 = limB2(gteMAC2, 1);
819 gteIR3 = limB3(gteMAC3, 1);
824 gteR2 = limC1(gteMAC1 >> 4);
825 gteG2 = limC2(gteMAC2 >> 4);
826 gteB2 = limC3(gteMAC3 >> 4);
834 GTE_LOG("GTE NCT\n");
838 for (v = 0; v < 3; v++) {
842 gteMAC1 = A1((((s64)gteL11 * vx) + (gteL12 * vy) + (gteL13 * vz)) >> 12);
843 gteMAC2 = A2((((s64)gteL21 * vx) + (gteL22 * vy) + (gteL23 * vz)) >> 12);
844 gteMAC3 = A3((((s64)gteL31 * vx) + (gteL32 * vy) + (gteL33 * vz)) >> 12);
845 gteIR1 = limB1(gteMAC1, 1);
846 gteIR2 = limB2(gteMAC2, 1);
847 gteIR3 = limB3(gteMAC3, 1);
848 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
849 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
850 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
854 gteR2 = limC1(gteMAC1 >> 4);
855 gteG2 = limC2(gteMAC2 >> 4);
856 gteB2 = limC3(gteMAC3 >> 4);
858 gteIR1 = limB1(gteMAC1, 1);
859 gteIR2 = limB2(gteMAC2, 1);
860 gteIR3 = limB3(gteMAC3, 1);
869 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
870 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
871 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
872 gteIR1 = limB1(gteMAC1, 1);
873 gteIR2 = limB2(gteMAC2, 1);
874 gteIR3 = limB3(gteMAC3, 1);
875 gteMAC1 = A1(((s64)gteR * gteIR1) >> 8);
876 gteMAC2 = A2(((s64)gteG * gteIR2) >> 8);
877 gteMAC3 = A3(((s64)gteB * gteIR3) >> 8);
878 gteIR1 = limB1(gteMAC1, 1);
879 gteIR2 = limB2(gteMAC2, 1);
880 gteIR3 = limB3(gteMAC3, 1);
885 gteR2 = limC1(gteMAC1 >> 4);
886 gteG2 = limC2(gteMAC2 >> 4);
887 gteB2 = limC3(gteMAC3 >> 4);
891 int shift = 12 * GTE_SF(gteop);
892 int lm = GTE_LM(gteop);
895 GTE_LOG("GTE INTPL\n");
899 gteMAC1 = A1(((gteIR1 << 12) + (gteIR0 * limB1(((s64)gteRFC - gteIR1), 0))) >> shift);
900 gteMAC2 = A2(((gteIR2 << 12) + (gteIR0 * limB2(((s64)gteGFC - gteIR2), 0))) >> shift);
901 gteMAC3 = A3(((gteIR3 << 12) + (gteIR0 * limB3(((s64)gteBFC - gteIR3), 0))) >> shift);
902 gteIR1 = limB1(gteMAC1, lm);
903 gteIR2 = limB2(gteMAC2, lm);
904 gteIR3 = limB3(gteMAC3, lm);
908 gteR2 = limC1(gteMAC1 >> 4);
909 gteG2 = limC2(gteMAC2 >> 4);
910 gteB2 = limC3(gteMAC3 >> 4);
915 GTE_LOG("GTE CDP\n");
919 gteMAC1 = A1((((s64)gteRBK << 12) + (gteLR1 * gteIR1) + (gteLR2 * gteIR2) + (gteLR3 * gteIR3)) >> 12);
920 gteMAC2 = A2((((s64)gteGBK << 12) + (gteLG1 * gteIR1) + (gteLG2 * gteIR2) + (gteLG3 * gteIR3)) >> 12);
921 gteMAC3 = A3((((s64)gteBBK << 12) + (gteLB1 * gteIR1) + (gteLB2 * gteIR2) + (gteLB3 * gteIR3)) >> 12);
922 gteIR1 = limB1(gteMAC1, 1);
923 gteIR2 = limB2(gteMAC2, 1);
924 gteIR3 = limB3(gteMAC3, 1);
925 gteMAC1 = A1(((((s64)gteR << 4) * gteIR1) + (gteIR0 * limB1(gteRFC - ((gteR * gteIR1) >> 8), 0))) >> 12);
926 gteMAC2 = A2(((((s64)gteG << 4) * gteIR2) + (gteIR0 * limB2(gteGFC - ((gteG * gteIR2) >> 8), 0))) >> 12);
927 gteMAC3 = A3(((((s64)gteB << 4) * gteIR3) + (gteIR0 * limB3(gteBFC - ((gteB * gteIR3) >> 8), 0))) >> 12);
928 gteIR1 = limB1(gteMAC1, 1);
929 gteIR2 = limB2(gteMAC2, 1);
930 gteIR3 = limB3(gteMAC3, 1);
935 gteR2 = limC1(gteMAC1 >> 4);
936 gteG2 = limC2(gteMAC2 >> 4);
937 gteB2 = limC3(gteMAC3 >> 4);