Merge pull request #542 from gameblabla/mdec_fix
[pcsx_rearmed.git] / libpcsxcore / gte.c
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/>                           *
5  *                                                                         *
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.                                   *
10  *                                                                         *
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.                          *
15  *                                                                         *
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  ***************************************************************************/
21
22 /*
23 * GTE functions.
24 */
25
26 #include "gte.h"
27 #include "psxmem.h"
28
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)
44
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; */
48
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])
111
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])
164
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)
174
175 #define gteop (psxRegs.code & 0x1ffffff)
176
177 #ifndef FLAGLESS
178
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;
184         }
185         return n_value;
186 }
187
188 static inline s32 LIM_(psxCP2Regs *regs, s32 value, s32 max, s32 min, u32 flag) {
189         s32 ret = value;
190         if (value > max) {
191                 gteFLAG |= flag;
192                 ret = max;
193         } else if (value < min) {
194                 gteFLAG |= flag;
195                 ret = min;
196         }
197         return ret;
198 }
199
200 static inline u32 limE_(psxCP2Regs *regs, u32 result) {
201         if (result > 0x1ffff) {
202                 gteFLAG |= (1 << 31) | (1 << 17);
203                 return 0x1ffff;
204         }
205         return result;
206 }
207
208 #else
209
210 #define BOUNDS_(regs, a, ...) (a)
211
212 static inline s32 LIM_(psxCP2Regs *regs, s32 value, s32 max, s32 min, u32 flag_unused) {
213         s32 ret = value;
214         if (value > max)
215                 ret = max;
216         else if (value < min)
217                 ret = min;
218         return ret;
219 }
220
221 static inline u32 limE_(psxCP2Regs *regs, u32 result) {
222         if (result > 0x1ffff)
223                 return 0x1ffff;
224         return result;
225 }
226
227 #endif
228
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) \
234         limE_(regs,result)
235
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))
246
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))
251
252 #ifndef __arm__
253 #define A1U A1
254 #define A2U A2
255 #define A3U A3
256 #else
257 /* these are unlikely to be hit and usually waste cycles, don't want them on ARM */
258 #define A1U(x) (x)
259 #define A2U(x) (x)
260 #define A3U(x) (x)
261 #endif
262
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) {
266         if (n < d * 2) {
267                 return ((u32)n << 16) / d;
268         }
269         return 0xffffffff;
270 }
271 #else
272 #include "gte_divider.h"
273 #endif // GTE_USE_NATIVE_DIVIDE
274
275 #ifndef FLAGLESS
276
277 u32 MFC2(int reg) {
278         psxCP2Regs *regs = &psxRegs.CP2;
279         switch (reg) {
280                 case 1:
281                 case 3:
282                 case 5:
283                 case 8:
284                 case 9:
285                 case 10:
286                 case 11:
287                         psxRegs.CP2D.r[reg] = (s32)psxRegs.CP2D.p[reg].sw.l;
288                         break;
289
290                 case 7:
291                 case 16:
292                 case 17:
293                 case 18:
294                 case 19:
295                         psxRegs.CP2D.r[reg] = (u32)psxRegs.CP2D.p[reg].w.l;
296                         break;
297
298                 case 15:
299                         psxRegs.CP2D.r[reg] = gteSXY2;
300                         break;
301
302                 case 28:
303                 case 29:
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);
307                         break;
308         }
309         return psxRegs.CP2D.r[reg];
310 }
311
312 void MTC2(u32 value, int reg) {
313         psxCP2Regs *regs = &psxRegs.CP2;
314         switch (reg) {
315                 case 15:
316                         gteSXY0 = gteSXY1;
317                         gteSXY1 = gteSXY2;
318                         gteSXY2 = value;
319                         gteSXYP = value;
320                         break;
321
322                 case 28:
323                         gteIRGB = value;
324                         gteIR1 = (value & 0x1f) << 7;
325                         gteIR2 = (value & 0x3e0) << 2;
326                         gteIR3 = (value & 0x7c00) >> 3;
327                         break;
328
329                 case 30:
330                         {
331                                 int a;
332                                 gteLZCS = value;
333
334                                 a = gteLZCS;
335                                 if (a > 0) {
336                                         int i;
337                                         for (i = 31; (a & (1 << i)) == 0 && i >= 0; i--);
338                                         gteLZCR = 31 - i;
339                                 } else if (a < 0) {
340                                         int i;
341                                         a ^= 0xffffffff;
342                                         for (i = 31; (a & (1 << i)) == 0 && i >= 0; i--);
343                                         gteLZCR = 31 - i;
344                                 } else {
345                                         gteLZCR = 32;
346                                 }
347                         }
348                         break;
349
350                 case 31:
351                         return;
352
353                 default:
354                         psxRegs.CP2D.r[reg] = value;
355         }
356 }
357
358 void CTC2(u32 value, int reg) {
359         switch (reg) {
360                 case 4:
361                 case 12:
362                 case 20:
363                 case 26:
364                 case 27:
365                 case 29:
366                 case 30:
367                         value = (s32)(s16)value;
368                         break;
369
370                 case 31:
371                         value = value & 0x7ffff000;
372                         if (value & 0x7f87e000) value |= 0x80000000;
373                         break;
374         }
375
376         psxRegs.CP2C.r[reg] = value;
377 }
378
379 void gteMFC2() {
380         if (!_Rt_) return;
381         psxRegs.GPR.r[_Rt_] = MFC2(_Rd_);
382 }
383
384 void gteCFC2() {
385         if (!_Rt_) return;
386         psxRegs.GPR.r[_Rt_] = psxRegs.CP2C.r[_Rd_];
387 }
388
389 void gteMTC2() {
390         MTC2(psxRegs.GPR.r[_Rt_], _Rd_);
391 }
392
393 void gteCTC2() {
394         CTC2(psxRegs.GPR.r[_Rt_], _Rd_);
395 }
396
397 #define _oB_ (psxRegs.GPR.r[_Rs_] + _Imm_)
398
399 void gteLWC2() {
400         MTC2(psxMemRead32(_oB_), _Rt_);
401 }
402
403 void gteSWC2() {
404         psxMemWrite32(_oB_, MFC2(_Rt_));
405 }
406
407 #endif // FLAGLESS
408
409 #if 0
410 #define DIVIDE DIVIDE_
411 static u32 DIVIDE_(s16 n, u16 d) {
412         if (n >= 0 && n < d * 2) {
413                 s32 n_ = n;
414                 return ((n_ << 16) + d / 2) / d;
415                 //return (u32)((float)(n_ << 16) / (float)d + (float)0.5);
416         }
417         return 0xffffffff;
418 }
419 #endif
420
421 void gteRTPS(psxCP2Regs *regs) {
422         int quotient;
423         s64 tmp;
424
425 #ifdef GTE_LOG
426         GTE_LOG("GTE RTPS\n");
427 #endif
428         gteFLAG = 0;
429
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);
436         gteSZ0 = gteSZ1;
437         gteSZ1 = gteSZ2;
438         gteSZ2 = gteSZ3;
439         gteSZ3 = limD(gteMAC3);
440         quotient = limE(DIVIDE(gteH, gteSZ3));
441         gteSXY0 = gteSXY1;
442         gteSXY1 = gteSXY2;
443         gteSX2 = limG1(F((s64)gteOFX + ((s64)gteIR1 * quotient)) >> 16);
444         gteSY2 = limG2(F((s64)gteOFY + ((s64)gteIR2 * quotient)) >> 16);
445
446         tmp = (s64)gteDQB + ((s64)gteDQA * quotient);
447         gteMAC0 = F(tmp);
448         gteIR0 = limH(tmp >> 12);
449 }
450
451 void gteRTPT(psxCP2Regs *regs) {
452         int quotient;
453         int v;
454         s32 vx, vy, vz;
455         s64 tmp;
456
457 #ifdef GTE_LOG
458         GTE_LOG("GTE RTPT\n");
459 #endif
460         gteFLAG = 0;
461
462         gteSZ0 = gteSZ3;
463         for (v = 0; v < 3; v++) {
464                 vx = VX(v);
465                 vy = VY(v);
466                 vz = VZ(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);
477         }
478
479         tmp = (s64)gteDQB + ((s64)gteDQA * quotient);
480         gteMAC0 = F(tmp);
481         gteIR0 = limH(tmp >> 12);
482 }
483
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);
490         s32 vx = VX(v);
491         s32 vy = VY(v);
492         s32 vz = VZ(v);
493
494 #ifdef GTE_LOG
495         GTE_LOG("GTE MVMVA\n");
496 #endif
497         gteFLAG = 0;
498
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);
502
503         gteIR1 = limB1(gteMAC1, lm);
504         gteIR2 = limB2(gteMAC2, lm);
505         gteIR3 = limB3(gteMAC3, lm);
506 }
507
508 void gteNCLIP(psxCP2Regs *regs) {
509 #ifdef GTE_LOG
510         GTE_LOG("GTE NCLIP\n");
511 #endif
512         gteFLAG = 0;
513
514         gteMAC0 = F((s64)gteSX0 * (gteSY1 - gteSY2) +
515                                 gteSX1 * (gteSY2 - gteSY0) +
516                                 gteSX2 * (gteSY0 - gteSY1));
517 }
518
519 void gteAVSZ3(psxCP2Regs *regs) {
520 #ifdef GTE_LOG
521         GTE_LOG("GTE AVSZ3\n");
522 #endif
523         gteFLAG = 0;
524
525         gteMAC0 = F((s64)gteZSF3 * (gteSZ1 + gteSZ2 + gteSZ3));
526         gteOTZ = limD(gteMAC0 >> 12);
527 }
528
529 void gteAVSZ4(psxCP2Regs *regs) {
530 #ifdef GTE_LOG
531         GTE_LOG("GTE AVSZ4\n");
532 #endif
533         gteFLAG = 0;
534
535         gteMAC0 = F((s64)gteZSF4 * (gteSZ0 + gteSZ1 + gteSZ2 + gteSZ3));
536         gteOTZ = limD(gteMAC0 >> 12);
537 }
538
539 void gteSQR(psxCP2Regs *regs) {
540         int shift = 12 * GTE_SF(gteop);
541         int lm = GTE_LM(gteop);
542
543 #ifdef GTE_LOG
544         GTE_LOG("GTE SQR\n");
545 #endif
546         gteFLAG = 0;
547
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);
554 }
555
556 void gteNCCS(psxCP2Regs *regs) {
557 #ifdef GTE_LOG
558         GTE_LOG("GTE NCCS\n");
559 #endif
560         gteFLAG = 0;
561
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;
577         gteIR1 = gteMAC1;
578         gteIR2 = gteMAC2;
579         gteIR3 = gteMAC3;
580
581         gteRGB0 = gteRGB1;
582         gteRGB1 = gteRGB2;
583         gteCODE2 = gteCODE;
584         gteR2 = limC1(gteMAC1 >> 4);
585         gteG2 = limC2(gteMAC2 >> 4);
586         gteB2 = limC3(gteMAC3 >> 4);
587 }
588
589 void gteNCCT(psxCP2Regs *regs) {
590         int v;
591         s32 vx, vy, vz;
592
593 #ifdef GTE_LOG
594         GTE_LOG("GTE NCCT\n");
595 #endif
596         gteFLAG = 0;
597
598         for (v = 0; v < 3; v++) {
599                 vx = VX(v);
600                 vy = VY(v);
601                 vz = VZ(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;
617
618                 gteRGB0 = gteRGB1;
619                 gteRGB1 = gteRGB2;
620                 gteCODE2 = gteCODE;
621                 gteR2 = limC1(gteMAC1 >> 4);
622                 gteG2 = limC2(gteMAC2 >> 4);
623                 gteB2 = limC3(gteMAC3 >> 4);
624         }
625         gteIR1 = gteMAC1;
626         gteIR2 = gteMAC2;
627         gteIR3 = gteMAC3;
628 }
629
630 void gteNCDS(psxCP2Regs *regs) {
631 #ifdef GTE_LOG
632         GTE_LOG("GTE NCDS\n");
633 #endif
634         gteFLAG = 0;
635
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);
654
655         gteRGB0 = gteRGB1;
656         gteRGB1 = gteRGB2;
657         gteCODE2 = gteCODE;
658         gteR2 = limC1(gteMAC1 >> 4);
659         gteG2 = limC2(gteMAC2 >> 4);
660         gteB2 = limC3(gteMAC3 >> 4);
661 }
662
663 void gteNCDT(psxCP2Regs *regs) {
664         int v;
665         s32 vx, vy, vz;
666
667 #ifdef GTE_LOG
668         GTE_LOG("GTE NCDT\n");
669 #endif
670         gteFLAG = 0;
671
672         for (v = 0; v < 3; v++) {
673                 vx = VX(v);
674                 vy = VY(v);
675                 vz = VZ(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;
691
692                 gteRGB0 = gteRGB1;
693                 gteRGB1 = gteRGB2;
694                 gteCODE2 = gteCODE;
695                 gteR2 = limC1(gteMAC1 >> 4);
696                 gteG2 = limC2(gteMAC2 >> 4);
697                 gteB2 = limC3(gteMAC3 >> 4);
698         }
699         gteIR1 = limB1(gteMAC1, 1);
700         gteIR2 = limB2(gteMAC2, 1);
701         gteIR3 = limB3(gteMAC3, 1);
702 }
703
704 void gteOP(psxCP2Regs *regs) {
705         int shift = 12 * GTE_SF(gteop);
706         int lm = GTE_LM(gteop);
707
708 #ifdef GTE_LOG
709         GTE_LOG("GTE OP\n");
710 #endif
711         gteFLAG = 0;
712
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);
719 }
720
721 void gteDCPL(psxCP2Regs *regs) {
722         int lm = GTE_LM(gteop);
723
724         s32 RIR1 = ((s32)gteR * gteIR1) >> 8;
725         s32 GIR2 = ((s32)gteG * gteIR2) >> 8;
726         s32 BIR3 = ((s32)gteB * gteIR3) >> 8;
727
728 #ifdef GTE_LOG
729         GTE_LOG("GTE DCPL\n");
730 #endif
731         gteFLAG = 0;
732
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);
736
737         gteIR1 = limB1(gteMAC1, lm);
738         gteIR2 = limB2(gteMAC2, lm);
739         gteIR3 = limB3(gteMAC3, lm);
740
741         gteRGB0 = gteRGB1;
742         gteRGB1 = gteRGB2;
743         gteCODE2 = gteCODE;
744         gteR2 = limC1(gteMAC1 >> 4);
745         gteG2 = limC2(gteMAC2 >> 4);
746         gteB2 = limC3(gteMAC3 >> 4);
747 }
748
749 void gteGPF(psxCP2Regs *regs) {
750         int shift = 12 * GTE_SF(gteop);
751
752 #ifdef GTE_LOG
753         GTE_LOG("GTE GPF\n");
754 #endif
755         gteFLAG = 0;
756
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);
763
764         gteRGB0 = gteRGB1;
765         gteRGB1 = gteRGB2;
766         gteCODE2 = gteCODE;
767         gteR2 = limC1(gteMAC1 >> 4);
768         gteG2 = limC2(gteMAC2 >> 4);
769         gteB2 = limC3(gteMAC3 >> 4);
770 }
771
772 void gteGPL(psxCP2Regs *regs) {
773         int shift = 12 * GTE_SF(gteop);
774
775 #ifdef GTE_LOG
776         GTE_LOG("GTE GPL\n");
777 #endif
778         gteFLAG = 0;
779
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);
786
787         gteRGB0 = gteRGB1;
788         gteRGB1 = gteRGB2;
789         gteCODE2 = gteCODE;
790         gteR2 = limC1(gteMAC1 >> 4);
791         gteG2 = limC2(gteMAC2 >> 4);
792         gteB2 = limC3(gteMAC3 >> 4);
793 }
794
795 void gteDPCS(psxCP2Regs *regs) {
796         int shift = 12 * GTE_SF(gteop);
797
798 #ifdef GTE_LOG
799         GTE_LOG("GTE DPCS\n");
800 #endif
801         gteFLAG = 0;
802
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;
806
807         gteIR1 = limB1(gteMAC1, 0);
808         gteIR2 = limB2(gteMAC2, 0);
809         gteIR3 = limB3(gteMAC3, 0);
810         gteRGB0 = gteRGB1;
811         gteRGB1 = gteRGB2;
812         gteCODE2 = gteCODE;
813         gteR2 = limC1(gteMAC1 >> 4);
814         gteG2 = limC2(gteMAC2 >> 4);
815         gteB2 = limC3(gteMAC3 >> 4);
816 }
817
818 void gteDPCT(psxCP2Regs *regs) {
819         int v;
820
821 #ifdef GTE_LOG
822         GTE_LOG("GTE DPCT\n");
823 #endif
824         gteFLAG = 0;
825
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;
830
831                 gteRGB0 = gteRGB1;
832                 gteRGB1 = gteRGB2;
833                 gteCODE2 = gteCODE;
834                 gteR2 = limC1(gteMAC1 >> 4);
835                 gteG2 = limC2(gteMAC2 >> 4);
836                 gteB2 = limC3(gteMAC3 >> 4);
837         }
838         gteIR1 = limB1(gteMAC1, 0);
839         gteIR2 = limB2(gteMAC2, 0);
840         gteIR3 = limB3(gteMAC3, 0);
841 }
842
843 void gteNCS(psxCP2Regs *regs) {
844 #ifdef GTE_LOG
845         GTE_LOG("GTE NCS\n");
846 #endif
847         gteFLAG = 0;
848
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);
861
862         gteRGB0 = gteRGB1;
863         gteRGB1 = gteRGB2;
864         gteCODE2 = gteCODE;
865         gteR2 = limC1(gteMAC1 >> 4);
866         gteG2 = limC2(gteMAC2 >> 4);
867         gteB2 = limC3(gteMAC3 >> 4);
868 }
869
870 void gteNCT(psxCP2Regs *regs) {
871         int v;
872         s32 vx, vy, vz;
873
874 #ifdef GTE_LOG
875         GTE_LOG("GTE NCT\n");
876 #endif
877         gteFLAG = 0;
878
879         for (v = 0; v < 3; v++) {
880                 vx = VX(v);
881                 vy = VY(v);
882                 vz = VZ(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);
892                 gteRGB0 = gteRGB1;
893                 gteRGB1 = gteRGB2;
894                 gteCODE2 = gteCODE;
895                 gteR2 = limC1(gteMAC1 >> 4);
896                 gteG2 = limC2(gteMAC2 >> 4);
897                 gteB2 = limC3(gteMAC3 >> 4);
898         }
899         gteIR1 = limB1(gteMAC1, 1);
900         gteIR2 = limB2(gteMAC2, 1);
901         gteIR3 = limB3(gteMAC3, 1);
902 }
903
904 void gteCC(psxCP2Regs *regs) {
905 #ifdef GTE_LOG
906         GTE_LOG("GTE CC\n");
907 #endif
908         gteFLAG = 0;
909
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);
922
923         gteRGB0 = gteRGB1;
924         gteRGB1 = gteRGB2;
925         gteCODE2 = gteCODE;
926         gteR2 = limC1(gteMAC1 >> 4);
927         gteG2 = limC2(gteMAC2 >> 4);
928         gteB2 = limC3(gteMAC3 >> 4);
929 }
930
931 void gteINTPL(psxCP2Regs *regs) {
932         int shift = 12 * GTE_SF(gteop);
933         int lm = GTE_LM(gteop);
934
935 #ifdef GTE_LOG
936         GTE_LOG("GTE INTPL\n");
937 #endif
938         gteFLAG = 0;
939
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);
946         gteRGB0 = gteRGB1;
947         gteRGB1 = gteRGB2;
948         gteCODE2 = gteCODE;
949         gteR2 = limC1(gteMAC1 >> 4);
950         gteG2 = limC2(gteMAC2 >> 4);
951         gteB2 = limC3(gteMAC3 >> 4);
952 }
953
954 void gteCDP(psxCP2Regs *regs) {
955 #ifdef GTE_LOG
956         GTE_LOG("GTE CDP\n");
957 #endif
958         gteFLAG = 0;
959
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);
972
973         gteRGB0 = gteRGB1;
974         gteRGB1 = gteRGB2;
975         gteCODE2 = gteCODE;
976         gteR2 = limC1(gteMAC1 >> 4);
977         gteG2 = limC2(gteMAC2 >> 4);
978         gteB2 = limC3(gteMAC3 >> 4);
979 }
980
981 /* decomposed/parametrized versions for the recompiler */
982
983 #ifndef FLAGLESS
984
985 void gteSQR_part_noshift(psxCP2Regs *regs) {
986         gteFLAG = 0;
987
988         gteMAC1 = gteIR1 * gteIR1;
989         gteMAC2 = gteIR2 * gteIR2;
990         gteMAC3 = gteIR3 * gteIR3;
991 }
992
993 void gteSQR_part_shift(psxCP2Regs *regs) {
994         gteFLAG = 0;
995
996         gteMAC1 = (gteIR1 * gteIR1) >> 12;
997         gteMAC2 = (gteIR2 * gteIR2) >> 12;
998         gteMAC3 = (gteIR3 * gteIR3) >> 12;
999 }
1000
1001 void gteOP_part_noshift(psxCP2Regs *regs) {
1002         gteFLAG = 0;
1003
1004         gteMAC1 = (gteR22 * gteIR3) - (gteR33 * gteIR2);
1005         gteMAC2 = (gteR33 * gteIR1) - (gteR11 * gteIR3);
1006         gteMAC3 = (gteR11 * gteIR2) - (gteR22 * gteIR1);
1007 }
1008
1009 void gteOP_part_shift(psxCP2Regs *regs) {
1010         gteFLAG = 0;
1011
1012         gteMAC1 = ((gteR22 * gteIR3) - (gteR33 * gteIR2)) >> 12;
1013         gteMAC2 = ((gteR33 * gteIR1) - (gteR11 * gteIR3)) >> 12;
1014         gteMAC3 = ((gteR11 * gteIR2) - (gteR22 * gteIR1)) >> 12;
1015 }
1016
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;
1021
1022         gteFLAG = 0;
1023
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);
1027 }
1028
1029 void gteGPF_part_noshift(psxCP2Regs *regs) {
1030         gteFLAG = 0;
1031
1032         gteMAC1 = gteIR0 * gteIR1;
1033         gteMAC2 = gteIR0 * gteIR2;
1034         gteMAC3 = gteIR0 * gteIR3;
1035 }
1036
1037 void gteGPF_part_shift(psxCP2Regs *regs) {
1038         gteFLAG = 0;
1039
1040         gteMAC1 = (gteIR0 * gteIR1) >> 12;
1041         gteMAC2 = (gteIR0 * gteIR2) >> 12;
1042         gteMAC3 = (gteIR0 * gteIR3) >> 12;
1043 }
1044
1045 #endif // !FLAGLESS
1046
1047 void gteGPL_part_noshift(psxCP2Regs *regs) {
1048         gteFLAG = 0;
1049
1050         gteMAC1 = A1((s64)gteMAC1 + (gteIR0 * gteIR1));
1051         gteMAC2 = A2((s64)gteMAC2 + (gteIR0 * gteIR2));
1052         gteMAC3 = A3((s64)gteMAC3 + (gteIR0 * gteIR3));
1053 }
1054
1055 void gteGPL_part_shift(psxCP2Regs *regs) {
1056         gteFLAG = 0;
1057
1058         gteMAC1 = A1((s64)gteMAC1 + ((gteIR0 * gteIR1) >> 12));
1059         gteMAC2 = A2((s64)gteMAC2 + ((gteIR0 * gteIR2) >> 12));
1060         gteMAC3 = A3((s64)gteMAC3 + ((gteIR0 * gteIR3) >> 12));
1061 }
1062
1063 void gteDPCS_part_noshift(psxCP2Regs *regs) {
1064         int shift = 0;
1065
1066         gteFLAG = 0;
1067
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;
1071 }
1072
1073 void gteDPCS_part_shift(psxCP2Regs *regs) {
1074         int shift = 12;
1075
1076         gteFLAG = 0;
1077
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;
1081 }
1082
1083 void gteINTPL_part_noshift(psxCP2Regs *regs) {
1084         gteFLAG = 0;
1085
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)));
1089 }
1090
1091 void gteINTPL_part_shift(psxCP2Regs *regs) {
1092         gteFLAG = 0;
1093
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;
1097 }
1098
1099 void gteMACtoRGB(psxCP2Regs *regs) {
1100         gteRGB0 = gteRGB1;
1101         gteRGB1 = gteRGB2;
1102         gteCODE2 = gteCODE;
1103         gteR2 = limC1(gteMAC1 >> 4);
1104         gteG2 = limC2(gteMAC2 >> 4);
1105         gteB2 = limC3(gteMAC3 >> 4);
1106 }
1107