ViewVC Help
View File | Revision Log | Show Annotations | Revision Graph | Root Listing
root/cebix/BasiliskII/src/uae_cpu/fpp.cpp
Revision: 1.1.1.1 (vendor branch)
Committed: 1999-10-03T14:16:26Z (25 years, 1 month ago) by cebix
Branch: cebix
CVS Tags: release-0_7-2, snapshot-21101999, start
Changes since 1.1: +0 -0 lines
Log Message:
Imported sources

File Contents

# User Rev Content
1 cebix 1.1 /*
2     * UAE - The Un*x Amiga Emulator
3     *
4     * MC68881 emulation
5     *
6     * Copyright 1996 Herman ten Brugge
7     */
8    
9     #include <math.h>
10     #include <stdio.h>
11    
12     #include "sysdeps.h"
13    
14     #include "memory.h"
15     #include "readcpu.h"
16     #include "newcpu.h"
17    
18     #if 1
19    
20     #define DEBUG_FPP 0
21    
22     /* single : S 8*E 23*F */
23     /* double : S 11*E 52*F */
24     /* extended : S 15*E 64*F */
25     /* E = 0 & F = 0 -> 0 */
26     /* E = MAX & F = 0 -> Infin */
27     /* E = MAX & F # 0 -> NotANumber */
28     /* E = biased by 127 (single) ,1023 (double) ,16383 (extended) */
29    
30     static __inline__ double to_single (uae_u32 value)
31     {
32     double frac;
33    
34     if ((value & 0x7fffffff) == 0)
35     return (0.0);
36     frac = (double) ((value & 0x7fffff) | 0x800000) / 8388608.0;
37     if (value & 0x80000000)
38     frac = -frac;
39     return (ldexp (frac, ((value >> 23) & 0xff) - 127));
40     }
41    
42     static __inline__ uae_u32 from_single (double src)
43     {
44     int expon;
45     uae_u32 tmp;
46     double frac;
47    
48     if (src == 0.0)
49     return 0;
50     if (src < 0) {
51     tmp = 0x80000000;
52     src = -src;
53     } else {
54     tmp = 0;
55     }
56     frac = frexp (src, &expon);
57     frac += 0.5 / 16777216.0;
58     if (frac >= 1.0) {
59     frac /= 2.0;
60     expon++;
61     }
62     return (tmp | (((expon + 127 - 1) & 0xff) << 23) |
63     (((int) (frac * 16777216.0)) & 0x7fffff));
64     }
65    
66     static __inline__ double to_exten(uae_u32 wrd1, uae_u32 wrd2, uae_u32 wrd3)
67     {
68     double frac;
69    
70     if ((wrd1 & 0x7fff0000) == 0 && wrd2 == 0 && wrd3 == 0)
71     return 0.0;
72     frac = (double) wrd2 / 2147483648.0 +
73     (double) wrd3 / 9223372036854775808.0;
74     if (wrd1 & 0x80000000)
75     frac = -frac;
76     return ldexp (frac, ((wrd1 >> 16) & 0x7fff) - 16383);
77     }
78    
79     static __inline__ void from_exten(double src, uae_u32 * wrd1, uae_u32 * wrd2, uae_u32 * wrd3)
80     {
81     int expon;
82     double frac;
83    
84     if (src == 0.0) {
85     *wrd1 = 0;
86     *wrd2 = 0;
87     *wrd3 = 0;
88     return;
89     }
90     if (src < 0) {
91     *wrd1 = 0x80000000;
92     src = -src;
93     } else {
94     *wrd1 = 0;
95     }
96     frac = frexp (src, &expon);
97     frac += 0.5 / 18446744073709551616.0;
98     if (frac >= 1.0) {
99     frac /= 2.0;
100     expon++;
101     }
102     *wrd1 |= (((expon + 16383 - 1) & 0x7fff) << 16);
103     *wrd2 = (uae_u32) (frac * 4294967296.0);
104     *wrd3 = (uae_u32) (frac * 18446744073709551616.0 - *wrd2 * 4294967296.0);
105     }
106    
107     static __inline__ double to_double(uae_u32 wrd1, uae_u32 wrd2)
108     {
109     double frac;
110    
111     if ((wrd1 & 0x7fffffff) == 0 && wrd2 == 0)
112     return 0.0;
113     frac = (double) ((wrd1 & 0xfffff) | 0x100000) / 1048576.0 +
114     (double) wrd2 / 4503599627370496.0;
115     if (wrd1 & 0x80000000)
116     frac = -frac;
117     return ldexp (frac, ((wrd1 >> 20) & 0x7ff) - 1023);
118     }
119    
120     static __inline__ void from_double(double src, uae_u32 * wrd1, uae_u32 * wrd2)
121     {
122     int expon;
123     int tmp;
124     double frac;
125    
126     if (src == 0.0) {
127     *wrd1 = 0;
128     *wrd2 = 0;
129     return;
130     }
131     if (src < 0) {
132     *wrd1 = 0x80000000;
133     src = -src;
134     } else {
135     *wrd1 = 0;
136     }
137     frac = frexp (src, &expon);
138     frac += 0.5 / 9007199254740992.0;
139     if (frac >= 1.0) {
140     frac /= 2.0;
141     expon++;
142     }
143     tmp = (uae_u32) (frac * 2097152.0);
144     *wrd1 |= (((expon + 1023 - 1) & 0x7ff) << 20) | (tmp & 0xfffff);
145     *wrd2 = (uae_u32) (frac * 9007199254740992.0 - tmp * 4294967296.0);
146     }
147    
148     static __inline__ double to_pack(uae_u32 wrd1, uae_u32 wrd2, uae_u32 wrd3)
149     {
150     double d;
151     char *cp;
152     char str[100];
153    
154     cp = str;
155     if (wrd1 & 0x80000000)
156     *cp++ = '-';
157     *cp++ = (wrd1 & 0xf) + '0';
158     *cp++ = '.';
159     *cp++ = ((wrd2 >> 28) & 0xf) + '0';
160     *cp++ = ((wrd2 >> 24) & 0xf) + '0';
161     *cp++ = ((wrd2 >> 20) & 0xf) + '0';
162     *cp++ = ((wrd2 >> 16) & 0xf) + '0';
163     *cp++ = ((wrd2 >> 12) & 0xf) + '0';
164     *cp++ = ((wrd2 >> 8) & 0xf) + '0';
165     *cp++ = ((wrd2 >> 4) & 0xf) + '0';
166     *cp++ = ((wrd2 >> 0) & 0xf) + '0';
167     *cp++ = ((wrd3 >> 28) & 0xf) + '0';
168     *cp++ = ((wrd3 >> 24) & 0xf) + '0';
169     *cp++ = ((wrd3 >> 20) & 0xf) + '0';
170     *cp++ = ((wrd3 >> 16) & 0xf) + '0';
171     *cp++ = ((wrd3 >> 12) & 0xf) + '0';
172     *cp++ = ((wrd3 >> 8) & 0xf) + '0';
173     *cp++ = ((wrd3 >> 4) & 0xf) + '0';
174     *cp++ = ((wrd3 >> 0) & 0xf) + '0';
175     *cp++ = 'E';
176     if (wrd1 & 0x40000000)
177     *cp++ = '-';
178     *cp++ = ((wrd1 >> 24) & 0xf) + '0';
179     *cp++ = ((wrd1 >> 20) & 0xf) + '0';
180     *cp++ = ((wrd1 >> 16) & 0xf) + '0';
181     *cp = 0;
182     sscanf(str, "%le", &d);
183     return d;
184     }
185    
186     static __inline__ void from_pack(double src, uae_u32 * wrd1, uae_u32 * wrd2, uae_u32 * wrd3)
187     {
188     int i;
189     int t;
190     char *cp;
191     char str[100];
192    
193     sprintf(str, "%.16e", src);
194     cp = str;
195     *wrd1 = *wrd2 = *wrd3 = 0;
196     if (*cp == '-') {
197     cp++;
198     *wrd1 = 0x80000000;
199     }
200     if (*cp == '+')
201     cp++;
202     *wrd1 |= (*cp++ - '0');
203     if (*cp == '.')
204     cp++;
205     for (i = 0; i < 8; i++) {
206     *wrd2 <<= 4;
207     if (*cp >= '0' && *cp <= '9')
208     *wrd2 |= *cp++ - '0';
209     }
210     for (i = 0; i < 8; i++) {
211     *wrd3 <<= 4;
212     if (*cp >= '0' && *cp <= '9')
213     *wrd3 |= *cp++ - '0';
214     }
215     if (*cp == 'e' || *cp == 'E') {
216     cp++;
217     if (*cp == '-') {
218     cp++;
219     *wrd1 |= 0x40000000;
220     }
221     if (*cp == '+')
222     cp++;
223     t = 0;
224     for (i = 0; i < 3; i++) {
225     if (*cp >= '0' && *cp <= '9')
226     t = (t << 4) | (*cp++ - '0');
227     }
228     *wrd1 |= t << 16;
229     }
230     }
231    
232     static __inline__ int get_fp_value (uae_u32 opcode, uae_u16 extra, double *src)
233     {
234     uaecptr tmppc;
235     uae_u16 tmp;
236     int size;
237     int mode;
238     int reg;
239     uae_u32 ad = 0;
240     static int sz1[8] =
241     {4, 4, 12, 12, 2, 8, 1, 0};
242     static int sz2[8] =
243     {4, 4, 12, 12, 2, 8, 2, 0};
244    
245     if ((extra & 0x4000) == 0) {
246     *src = regs.fp[(extra >> 10) & 7];
247     return 1;
248     }
249     mode = (opcode >> 3) & 7;
250     reg = opcode & 7;
251     size = (extra >> 10) & 7;
252     switch (mode) {
253     case 0:
254     switch (size) {
255     case 6:
256     *src = (double) (uae_s8) m68k_dreg (regs, reg);
257     break;
258     case 4:
259     *src = (double) (uae_s16) m68k_dreg (regs, reg);
260     break;
261     case 0:
262     *src = (double) (uae_s32) m68k_dreg (regs, reg);
263     break;
264     case 1:
265     *src = to_single(m68k_dreg (regs, reg));
266     break;
267     default:
268     return 0;
269     }
270     return 1;
271     case 1:
272     return 0;
273     case 2:
274     ad = m68k_areg (regs, reg);
275     break;
276     case 3:
277     ad = m68k_areg (regs, reg);
278     m68k_areg (regs, reg) += reg == 7 ? sz2[size] : sz1[size];
279     break;
280     case 4:
281     m68k_areg (regs, reg) -= reg == 7 ? sz2[size] : sz1[size];
282     ad = m68k_areg (regs, reg);
283     break;
284     case 5:
285     ad = m68k_areg (regs, reg) + (uae_s32) (uae_s16) next_iword();
286     break;
287     case 6:
288     ad = get_disp_ea_020 (m68k_areg (regs, reg), next_iword());
289     break;
290     case 7:
291     switch (reg) {
292     case 0:
293     ad = (uae_s32) (uae_s16) next_iword();
294     break;
295     case 1:
296     ad = next_ilong();
297     break;
298     case 2:
299     ad = m68k_getpc ();
300     ad += (uae_s32) (uae_s16) next_iword();
301     break;
302     case 3:
303     tmppc = m68k_getpc ();
304     tmp = next_iword();
305     ad = get_disp_ea_020 (tmppc, tmp);
306     break;
307     case 4:
308     ad = m68k_getpc ();
309     m68k_setpc (ad + sz2[size]);
310     break;
311     default:
312     return 0;
313     }
314     }
315     switch (size) {
316     case 0:
317     *src = (double) (uae_s32) get_long (ad);
318     break;
319     case 1:
320     *src = to_single(get_long (ad));
321     break;
322     case 2:{
323     uae_u32 wrd1, wrd2, wrd3;
324     wrd1 = get_long (ad);
325     ad += 4;
326     wrd2 = get_long (ad);
327     ad += 4;
328     wrd3 = get_long (ad);
329     *src = to_exten(wrd1, wrd2, wrd3);
330     }
331     break;
332     case 3:{
333     uae_u32 wrd1, wrd2, wrd3;
334     wrd1 = get_long (ad);
335     ad += 4;
336     wrd2 = get_long (ad);
337     ad += 4;
338     wrd3 = get_long (ad);
339     *src = to_pack(wrd1, wrd2, wrd3);
340     }
341     break;
342     case 4:
343     *src = (double) (uae_s16) get_word(ad);
344     break;
345     case 5:{
346     uae_u32 wrd1, wrd2;
347     wrd1 = get_long (ad);
348     ad += 4;
349     wrd2 = get_long (ad);
350     *src = to_double(wrd1, wrd2);
351     }
352     break;
353     case 6:
354     *src = (double) (uae_s8) get_byte(ad);
355     break;
356     default:
357     return 0;
358     }
359     return 1;
360     }
361    
362     static __inline__ int put_fp_value (double value, uae_u32 opcode, uae_u16 extra)
363     {
364     uae_u16 tmp;
365     uaecptr tmppc;
366     int size;
367     int mode;
368     int reg;
369     uae_u32 ad;
370     static int sz1[8] =
371     {4, 4, 12, 12, 2, 8, 1, 0};
372     static int sz2[8] =
373     {4, 4, 12, 12, 2, 8, 2, 0};
374    
375     if ((extra & 0x4000) == 0) {
376     regs.fp[(extra >> 10) & 7] = value;
377     return 1;
378     }
379     mode = (opcode >> 3) & 7;
380     reg = opcode & 7;
381     size = (extra >> 10) & 7;
382     ad = 0xffffffff;
383     switch (mode) {
384     case 0:
385     switch (size) {
386     case 6:
387     m68k_dreg (regs, reg) = (((int) value & 0xff)
388     | (m68k_dreg (regs, reg) & ~0xff));
389     break;
390     case 4:
391     m68k_dreg (regs, reg) = (((int) value & 0xffff)
392     | (m68k_dreg (regs, reg) & ~0xffff));
393     break;
394     case 0:
395     m68k_dreg (regs, reg) = (int) value;
396     break;
397     case 1:
398     m68k_dreg (regs, reg) = from_single(value);
399     break;
400     default:
401     return 0;
402     }
403     return 1;
404     case 1:
405     return 0;
406     case 2:
407     ad = m68k_areg (regs, reg);
408     break;
409     case 3:
410     ad = m68k_areg (regs, reg);
411     m68k_areg (regs, reg) += reg == 7 ? sz2[size] : sz1[size];
412     break;
413     case 4:
414     m68k_areg (regs, reg) -= reg == 7 ? sz2[size] : sz1[size];
415     ad = m68k_areg (regs, reg);
416     break;
417     case 5:
418     ad = m68k_areg (regs, reg) + (uae_s32) (uae_s16) next_iword();
419     break;
420     case 6:
421     ad = get_disp_ea_020 (m68k_areg (regs, reg), next_iword());
422     break;
423     case 7:
424     switch (reg) {
425     case 0:
426     ad = (uae_s32) (uae_s16) next_iword();
427     break;
428     case 1:
429     ad = next_ilong();
430     break;
431     case 2:
432     ad = m68k_getpc ();
433     ad += (uae_s32) (uae_s16) next_iword();
434     break;
435     case 3:
436     tmppc = m68k_getpc ();
437     tmp = next_iword();
438     ad = get_disp_ea_020 (tmppc, tmp);
439     break;
440     case 4:
441     ad = m68k_getpc ();
442     m68k_setpc (ad + sz2[size]);
443     break;
444     default:
445     return 0;
446     }
447     }
448     switch (size) {
449     case 0:
450     put_long (ad, (uae_s32) value);
451     break;
452     case 1:
453     put_long (ad, from_single(value));
454     break;
455     case 2:
456     {
457     uae_u32 wrd1, wrd2, wrd3;
458     from_exten(value, &wrd1, &wrd2, &wrd3);
459     put_long (ad, wrd1);
460     ad += 4;
461     put_long (ad, wrd2);
462     ad += 4;
463     put_long (ad, wrd3);
464     }
465     break;
466     case 3:
467     {
468     uae_u32 wrd1, wrd2, wrd3;
469     from_pack(value, &wrd1, &wrd2, &wrd3);
470     put_long (ad, wrd1);
471     ad += 4;
472     put_long (ad, wrd2);
473     ad += 4;
474     put_long (ad, wrd3);
475     }
476     break;
477     case 4:
478     put_word(ad, (uae_s16) value);
479     break;
480     case 5:{
481     uae_u32 wrd1, wrd2;
482     from_double(value, &wrd1, &wrd2);
483     put_long (ad, wrd1);
484     ad += 4;
485     put_long (ad, wrd2);
486     }
487     break;
488     case 6:
489     put_byte(ad, (uae_s8) value);
490     break;
491     default:
492     return 0;
493     }
494     return 1;
495     }
496    
497     static __inline__ int get_fp_ad(uae_u32 opcode, uae_u32 * ad)
498     {
499     uae_u16 tmp;
500     uaecptr tmppc;
501     int mode;
502     int reg;
503    
504     mode = (opcode >> 3) & 7;
505     reg = opcode & 7;
506     switch (mode) {
507     case 0:
508     case 1:
509     return 0;
510     case 2:
511     *ad = m68k_areg (regs, reg);
512     break;
513     case 3:
514     *ad = m68k_areg (regs, reg);
515     break;
516     case 4:
517     *ad = m68k_areg (regs, reg);
518     break;
519     case 5:
520     *ad = m68k_areg (regs, reg) + (uae_s32) (uae_s16) next_iword();
521     break;
522     case 6:
523     *ad = get_disp_ea_020 (m68k_areg (regs, reg), next_iword());
524     break;
525     case 7:
526     switch (reg) {
527     case 0:
528     *ad = (uae_s32) (uae_s16) next_iword();
529     break;
530     case 1:
531     *ad = next_ilong();
532     break;
533     case 2:
534     *ad = m68k_getpc ();
535     *ad += (uae_s32) (uae_s16) next_iword();
536     break;
537     case 3:
538     tmppc = m68k_getpc ();
539     tmp = next_iword();
540     *ad = get_disp_ea_020 (tmppc, tmp);
541     break;
542     default:
543     return 0;
544     }
545     }
546     return 1;
547     }
548    
549     static __inline__ int fpp_cond(uae_u32 opcode, int contition)
550     {
551     int N = (regs.fpsr & 0x8000000) != 0;
552     int Z = (regs.fpsr & 0x4000000) != 0;
553     /* int I = (regs.fpsr & 0x2000000) != 0; */
554     int NotANumber = (regs.fpsr & 0x1000000) != 0;
555    
556     switch (contition) {
557     case 0x00:
558     return 0;
559     case 0x01:
560     return Z;
561     case 0x02:
562     return !(NotANumber || Z || N);
563     case 0x03:
564     return Z || !(NotANumber || N);
565     case 0x04:
566     return N && !(NotANumber || Z);
567     case 0x05:
568     return Z || (N && !NotANumber);
569     case 0x06:
570     return !(NotANumber || Z);
571     case 0x07:
572     return !NotANumber;
573     case 0x08:
574     return NotANumber;
575     case 0x09:
576     return NotANumber || Z;
577     case 0x0a:
578     return NotANumber || !(N || Z);
579     case 0x0b:
580     return NotANumber || Z || !N;
581     case 0x0c:
582     return NotANumber || (N && !Z);
583     case 0x0d:
584     return NotANumber || Z || N;
585     case 0x0e:
586     return !Z;
587     case 0x0f:
588     return 1;
589     case 0x10:
590     return 0;
591     case 0x11:
592     return Z;
593     case 0x12:
594     return !(NotANumber || Z || N);
595     case 0x13:
596     return Z || !(NotANumber || N);
597     case 0x14:
598     return N && !(NotANumber || Z);
599     case 0x15:
600     return Z || (N && !NotANumber);
601     case 0x16:
602     return !(NotANumber || Z);
603     case 0x17:
604     return !NotANumber;
605     case 0x18:
606     return NotANumber;
607     case 0x19:
608     return NotANumber || Z;
609     case 0x1a:
610     return NotANumber || !(N || Z);
611     case 0x1b:
612     return NotANumber || Z || !N;
613     case 0x1c:
614     return NotANumber || (Z && N);
615     case 0x1d:
616     return NotANumber || Z || N;
617     case 0x1e:
618     return !Z;
619     case 0x1f:
620     return 1;
621     }
622     return -1;
623     }
624    
625     void fdbcc_opp(uae_u32 opcode, uae_u16 extra)
626     {
627     uaecptr pc = (uae_u32) m68k_getpc ();
628     uae_s32 disp = (uae_s32) (uae_s16) next_iword();
629     int cc;
630    
631     #if DEBUG_FPP
632     printf("fdbcc_opp at %08lx\n", m68k_getpc ());
633     fflush(stdout);
634     #endif
635     cc = fpp_cond(opcode, extra & 0x3f);
636     if (cc == -1) {
637     m68k_setpc (pc - 4);
638     op_illg (opcode);
639     } else if (!cc) {
640     int reg = opcode & 0x7;
641    
642     m68k_dreg (regs, reg) = ((m68k_dreg (regs, reg) & ~0xffff)
643     | ((m68k_dreg (regs, reg) - 1) & 0xffff));
644     if ((m68k_dreg (regs, reg) & 0xffff) == 0xffff)
645     m68k_setpc (pc + disp);
646     }
647     }
648    
649     void fscc_opp(uae_u32 opcode, uae_u16 extra)
650     {
651     uae_u32 ad;
652     int cc;
653    
654     #if DEBUG_FPP
655     printf("fscc_opp at %08lx\n", m68k_getpc ());
656     fflush(stdout);
657     #endif
658     cc = fpp_cond(opcode, extra & 0x3f);
659     if (cc == -1) {
660     m68k_setpc (m68k_getpc () - 4);
661     op_illg (opcode);
662     } else if ((opcode & 0x38) == 0) {
663     m68k_dreg (regs, opcode & 7) = (m68k_dreg (regs, opcode & 7) & ~0xff) |
664     (cc ? 0xff : 0x00);
665     } else {
666     if (get_fp_ad(opcode, &ad) == 0) {
667     m68k_setpc (m68k_getpc () - 4);
668     op_illg (opcode);
669     } else
670     put_byte(ad, cc ? 0xff : 0x00);
671     }
672     }
673    
674     void ftrapcc_opp(uae_u32 opcode, uaecptr oldpc)
675     {
676     int cc;
677    
678     #if DEBUG_FPP
679     printf("ftrapcc_opp at %08lx\n", m68k_getpc ());
680     fflush(stdout);
681     #endif
682     cc = fpp_cond(opcode, opcode & 0x3f);
683     if (cc == -1) {
684     m68k_setpc (oldpc);
685     op_illg (opcode);
686     }
687     if (cc)
688     Exception(7, oldpc - 2);
689     }
690    
691     void fbcc_opp(uae_u32 opcode, uaecptr pc, uae_u32 extra)
692     {
693     int cc;
694    
695     #if DEBUG_FPP
696     printf("fbcc_opp at %08lx\n", m68k_getpc ());
697     fflush(stdout);
698     #endif
699     cc = fpp_cond(opcode, opcode & 0x3f);
700     if (cc == -1) {
701     m68k_setpc (pc);
702     op_illg (opcode);
703     } else if (cc) {
704     if ((opcode & 0x40) == 0)
705     extra = (uae_s32) (uae_s16) extra;
706     m68k_setpc (pc + extra);
707     }
708     }
709    
710     void fsave_opp(uae_u32 opcode)
711     {
712     uae_u32 ad;
713     int incr = (opcode & 0x38) == 0x20 ? -1 : 1;
714     int i;
715    
716     #if DEBUG_FPP
717     printf("fsave_opp at %08lx\n", m68k_getpc ());
718     fflush(stdout);
719     #endif
720     if (get_fp_ad(opcode, &ad) == 0) {
721     m68k_setpc (m68k_getpc () - 2);
722     op_illg (opcode);
723     return;
724     }
725     if (incr < 0) {
726     ad -= 4;
727     put_long (ad, 0x70000000);
728     for (i = 0; i < 5; i++) {
729     ad -= 4;
730     put_long (ad, 0x00000000);
731     }
732     ad -= 4;
733     put_long (ad, 0x1f180000);
734     } else {
735     put_long (ad, 0x1f180000);
736     ad += 4;
737     for (i = 0; i < 5; i++) {
738     put_long (ad, 0x00000000);
739     ad += 4;
740     }
741     put_long (ad, 0x70000000);
742     ad += 4;
743     }
744     if ((opcode & 0x38) == 0x18)
745     m68k_areg (regs, opcode & 7) = ad;
746     if ((opcode & 0x38) == 0x20)
747     m68k_areg (regs, opcode & 7) = ad;
748     }
749    
750     void frestore_opp(uae_u32 opcode)
751     {
752     uae_u32 ad;
753     uae_u32 d;
754     int incr = (opcode & 0x38) == 0x20 ? -1 : 1;
755    
756     #if DEBUG_FPP
757     printf("frestore_opp at %08lx\n", m68k_getpc ());
758     fflush(stdout);
759     #endif
760     if (get_fp_ad(opcode, &ad) == 0) {
761     m68k_setpc (m68k_getpc () - 2);
762     op_illg (opcode);
763     return;
764     }
765     if (incr < 0) {
766     ad -= 4;
767     d = get_long (ad);
768     if ((d & 0xff000000) != 0) {
769     if ((d & 0x00ff0000) == 0x00180000)
770     ad -= 6 * 4;
771     else if ((d & 0x00ff0000) == 0x00380000)
772     ad -= 14 * 4;
773     else if ((d & 0x00ff0000) == 0x00b40000)
774     ad -= 45 * 4;
775     }
776     } else {
777     d = get_long (ad);
778     ad += 4;
779     if ((d & 0xff000000) != 0) {
780     if ((d & 0x00ff0000) == 0x00180000)
781     ad += 6 * 4;
782     else if ((d & 0x00ff0000) == 0x00380000)
783     ad += 14 * 4;
784     else if ((d & 0x00ff0000) == 0x00b40000)
785     ad += 45 * 4;
786     }
787     }
788     if ((opcode & 0x38) == 0x18)
789     m68k_areg (regs, opcode & 7) = ad;
790     if ((opcode & 0x38) == 0x20)
791     m68k_areg (regs, opcode & 7) = ad;
792     }
793    
794     void fpp_opp(uae_u32 opcode, uae_u16 extra)
795     {
796     int reg;
797     double src;
798    
799     #if DEBUG_FPP
800     printf("FPP %04lx %04x at %08lx\n", opcode & 0xffff, extra & 0xffff,
801     m68k_getpc () - 4);
802     fflush(stdout);
803     #endif
804     switch ((extra >> 13) & 0x7) {
805     case 3:
806     if (put_fp_value (regs.fp[(extra >> 7) & 7], opcode, extra) == 0) {
807     m68k_setpc (m68k_getpc () - 4);
808     op_illg (opcode);
809     }
810     return;
811     case 4:
812     case 5:
813     if ((opcode & 0x38) == 0) {
814     if (extra & 0x2000) {
815     if (extra & 0x1000)
816     m68k_dreg (regs, opcode & 7) = regs.fpcr;
817     if (extra & 0x0800)
818     m68k_dreg (regs, opcode & 7) = regs.fpsr;
819     if (extra & 0x0400)
820     m68k_dreg (regs, opcode & 7) = regs.fpiar;
821     } else {
822     if (extra & 0x1000)
823     regs.fpcr = m68k_dreg (regs, opcode & 7);
824     if (extra & 0x0800)
825     regs.fpsr = m68k_dreg (regs, opcode & 7);
826     if (extra & 0x0400)
827     regs.fpiar = m68k_dreg (regs, opcode & 7);
828     }
829     } else if ((opcode & 0x38) == 1) {
830     if (extra & 0x2000) {
831     if (extra & 0x1000)
832     m68k_areg (regs, opcode & 7) = regs.fpcr;
833     if (extra & 0x0800)
834     m68k_areg (regs, opcode & 7) = regs.fpsr;
835     if (extra & 0x0400)
836     m68k_areg (regs, opcode & 7) = regs.fpiar;
837     } else {
838     if (extra & 0x1000)
839     regs.fpcr = m68k_areg (regs, opcode & 7);
840     if (extra & 0x0800)
841     regs.fpsr = m68k_areg (regs, opcode & 7);
842     if (extra & 0x0400)
843     regs.fpiar = m68k_areg (regs, opcode & 7);
844     }
845     } else if ((opcode & 0x3f) == 0x3c) {
846     if ((extra & 0x2000) == 0) {
847     if (extra & 0x1000)
848     regs.fpcr = next_ilong();
849     if (extra & 0x0800)
850     regs.fpsr = next_ilong();
851     if (extra & 0x0400)
852     regs.fpiar = next_ilong();
853     }
854     } else if (extra & 0x2000) {
855     /* FMOVEM FPP->memory */
856     uae_u32 ad;
857     int incr = 0;
858    
859     if (get_fp_ad(opcode, &ad) == 0) {
860     m68k_setpc (m68k_getpc () - 4);
861     op_illg (opcode);
862     return;
863     }
864     if ((opcode & 0x38) == 0x20) {
865     if (extra & 0x1000)
866     incr += 4;
867     if (extra & 0x0800)
868     incr += 4;
869     if (extra & 0x0400)
870     incr += 4;
871     }
872     ad -= incr;
873     if (extra & 0x1000) {
874     put_long (ad, regs.fpcr);
875     ad += 4;
876     }
877     if (extra & 0x0800) {
878     put_long (ad, regs.fpsr);
879     ad += 4;
880     }
881     if (extra & 0x0400) {
882     put_long (ad, regs.fpiar);
883     ad += 4;
884     }
885     ad -= incr;
886     if ((opcode & 0x38) == 0x18)
887     m68k_areg (regs, opcode & 7) = ad;
888     if ((opcode & 0x38) == 0x20)
889     m68k_areg (regs, opcode & 7) = ad;
890     } else {
891     /* FMOVEM memory->FPP */
892     uae_u32 ad;
893    
894     if (get_fp_ad(opcode, &ad) == 0) {
895     m68k_setpc (m68k_getpc () - 4);
896     op_illg (opcode);
897     return;
898     }
899     ad = (opcode & 0x38) == 0x20 ? ad - 12 : ad;
900     if (extra & 0x1000) {
901     regs.fpcr = get_long (ad);
902     ad += 4;
903     }
904     if (extra & 0x0800) {
905     regs.fpsr = get_long (ad);
906     ad += 4;
907     }
908     if (extra & 0x0400) {
909     regs.fpiar = get_long (ad);
910     ad += 4;
911     }
912     if ((opcode & 0x38) == 0x18)
913     m68k_areg (regs, opcode & 7) = ad;
914     if ((opcode & 0x38) == 0x20)
915     m68k_areg (regs, opcode & 7) = ad - 12;
916     }
917     return;
918     case 6:
919     case 7:
920     {
921     uae_u32 ad, list = 0;
922     int incr = 0;
923     if (extra & 0x2000) {
924     /* FMOVEM FPP->memory */
925     if (get_fp_ad(opcode, &ad) == 0) {
926     m68k_setpc (m68k_getpc () - 4);
927     op_illg (opcode);
928     return;
929     }
930     switch ((extra >> 11) & 3) {
931     case 0: /* static pred */
932     list = extra & 0xff;
933     incr = -1;
934     break;
935     case 1: /* dynamic pred */
936     list = m68k_dreg (regs, (extra >> 4) & 3) & 0xff;
937     incr = -1;
938     break;
939     case 2: /* static postinc */
940     list = extra & 0xff;
941     incr = 1;
942     break;
943     case 3: /* dynamic postinc */
944     list = m68k_dreg (regs, (extra >> 4) & 3) & 0xff;
945     incr = 1;
946     break;
947     }
948     while (list) {
949     uae_u32 wrd1, wrd2, wrd3;
950     if (incr < 0) {
951     from_exten(regs.fp[fpp_movem_index2[list]],
952     &wrd1, &wrd2, &wrd3);
953     ad -= 4;
954     put_long (ad, wrd3);
955     ad -= 4;
956     put_long (ad, wrd2);
957     ad -= 4;
958     put_long (ad, wrd1);
959     } else {
960     from_exten(regs.fp[fpp_movem_index1[list]],
961     &wrd1, &wrd2, &wrd3);
962     put_long (ad, wrd1);
963     ad += 4;
964     put_long (ad, wrd2);
965     ad += 4;
966     put_long (ad, wrd3);
967     ad += 4;
968     }
969     list = fpp_movem_next[list];
970     }
971     if ((opcode & 0x38) == 0x18)
972     m68k_areg (regs, opcode & 7) = ad;
973     if ((opcode & 0x38) == 0x20)
974     m68k_areg (regs, opcode & 7) = ad;
975     } else {
976     /* FMOVEM memory->FPP */
977     if (get_fp_ad(opcode, &ad) == 0) {
978     m68k_setpc (m68k_getpc () - 4);
979     op_illg (opcode);
980     return;
981     }
982     switch ((extra >> 11) & 3) {
983     case 0: /* static pred */
984     list = extra & 0xff;
985     incr = -1;
986     break;
987     case 1: /* dynamic pred */
988     list = m68k_dreg (regs, (extra >> 4) & 3) & 0xff;
989     incr = -1;
990     break;
991     case 2: /* static postinc */
992     list = extra & 0xff;
993     incr = 1;
994     break;
995     case 3: /* dynamic postinc */
996     list = m68k_dreg (regs, (extra >> 4) & 3) & 0xff;
997     incr = 1;
998     break;
999     }
1000     while (list) {
1001     uae_u32 wrd1, wrd2, wrd3;
1002     if (incr < 0) {
1003     ad -= 4;
1004     wrd3 = get_long (ad);
1005     ad -= 4;
1006     wrd2 = get_long (ad);
1007     ad -= 4;
1008     wrd1 = get_long (ad);
1009     regs.fp[fpp_movem_index2[list]] = to_exten (wrd1, wrd2, wrd3);
1010     } else {
1011     wrd1 = get_long (ad);
1012     ad += 4;
1013     wrd2 = get_long (ad);
1014     ad += 4;
1015     wrd3 = get_long (ad);
1016     ad += 4;
1017     regs.fp[fpp_movem_index1[list]] = to_exten (wrd1, wrd2, wrd3);
1018     }
1019     list = fpp_movem_next[list];
1020     }
1021     if ((opcode & 0x38) == 0x18)
1022     m68k_areg (regs, opcode & 7) = ad;
1023     if ((opcode & 0x38) == 0x20)
1024     m68k_areg (regs, opcode & 7) = ad;
1025     }
1026     }
1027     return;
1028     case 0:
1029     case 2:
1030     reg = (extra >> 7) & 7;
1031     if ((extra & 0xfc00) == 0x5c00) {
1032     switch (extra & 0x7f) {
1033     case 0x00:
1034     regs.fp[reg] = 4.0 * atan(1.0);
1035     break;
1036     case 0x0b:
1037     regs.fp[reg] = log10 (2.0);
1038     break;
1039     case 0x0c:
1040     regs.fp[reg] = exp (1.0);
1041     break;
1042     case 0x0d:
1043     regs.fp[reg] = log (exp (1.0)) / log (2.0);
1044     break;
1045     case 0x0e:
1046     regs.fp[reg] = log (exp (1.0)) / log (10.0);
1047     break;
1048     case 0x0f:
1049     regs.fp[reg] = 0.0;
1050     break;
1051     case 0x30:
1052     regs.fp[reg] = log (2.0);
1053     break;
1054     case 0x31:
1055     regs.fp[reg] = log (10.0);
1056     break;
1057     case 0x32:
1058     regs.fp[reg] = 1.0e0;
1059     break;
1060     case 0x33:
1061     regs.fp[reg] = 1.0e1;
1062     break;
1063     case 0x34:
1064     regs.fp[reg] = 1.0e2;
1065     break;
1066     case 0x35:
1067     regs.fp[reg] = 1.0e4;
1068     break;
1069     case 0x36:
1070     regs.fp[reg] = 1.0e8;
1071     break;
1072     case 0x37:
1073     regs.fp[reg] = 1.0e16;
1074     break;
1075     case 0x38:
1076     regs.fp[reg] = 1.0e32;
1077     break;
1078     case 0x39:
1079     regs.fp[reg] = 1.0e64;
1080     break;
1081     case 0x3a:
1082     regs.fp[reg] = 1.0e128;
1083     break;
1084     case 0x3b:
1085     regs.fp[reg] = 1.0e256;
1086     break;
1087     #if 0
1088     case 0x3c:
1089     regs.fp[reg] = 1.0e512;
1090     break;
1091     case 0x3d:
1092     regs.fp[reg] = 1.0e1024;
1093     break;
1094     case 0x3e:
1095     regs.fp[reg] = 1.0e2048;
1096     break;
1097     case 0x3f:
1098     regs.fp[reg] = 1.0e4096;
1099     break;
1100     #endif
1101     default:
1102     m68k_setpc (m68k_getpc () - 4);
1103     op_illg (opcode);
1104     break;
1105     }
1106     return;
1107     }
1108     if (get_fp_value (opcode, extra, &src) == 0) {
1109     m68k_setpc (m68k_getpc () - 4);
1110     op_illg (opcode);
1111     return;
1112     }
1113     switch (extra & 0x7f) {
1114     case 0x00: /* FMOVE */
1115     regs.fp[reg] = src;
1116     break;
1117     case 0x01: /* FINT */
1118     regs.fp[reg] = (int) (src + 0.5);
1119     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1120     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1121     break;
1122     case 0x02: /* FSINH */
1123     regs.fp[reg] = sinh (src);
1124     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1125     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1126     break;
1127     case 0x03: /* FINTRZ */
1128     regs.fp[reg] = (int) src;
1129     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1130     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1131     break;
1132     case 0x04: /* FSQRT */
1133     regs.fp[reg] = sqrt (src);
1134     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1135     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1136     break;
1137     case 0x06: /* FLOGNP1 */
1138     regs.fp[reg] = log (src + 1.0);
1139     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1140     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1141     break;
1142     case 0x08: /* FETOXM1 */
1143     regs.fp[reg] = exp (src) - 1.0;
1144     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1145     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1146     break;
1147     case 0x09: /* FTANH */
1148     regs.fp[reg] = tanh (src);
1149     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1150     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1151     break;
1152     case 0x0a: /* FATAN */
1153     regs.fp[reg] = atan (src);
1154     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1155     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1156     break;
1157     case 0x0c: /* FASIN */
1158     regs.fp[reg] = asin (src);
1159     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1160     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1161     break;
1162     case 0x0d: /* FATANH */
1163     #if 1 /* The BeBox doesn't have atanh, and it isn't in the HPUX libm either */
1164     regs.fp[reg] = log ((1 + src) / (1 - src)) / 2;
1165     #else
1166     regs.fp[reg] = atanh (src);
1167     #endif
1168     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1169     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1170     break;
1171     case 0x0e: /* FSIN */
1172     regs.fp[reg] = sin (src);
1173     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1174     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1175     break;
1176     case 0x0f: /* FTAN */
1177     regs.fp[reg] = tan (src);
1178     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1179     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1180     break;
1181     case 0x10: /* FETOX */
1182     regs.fp[reg] = exp (src);
1183     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1184     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1185     break;
1186     case 0x11: /* FTWOTOX */
1187     regs.fp[reg] = pow(2.0, src);
1188     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1189     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1190     break;
1191     case 0x12: /* FTENTOX */
1192     regs.fp[reg] = pow(10.0, src);
1193     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1194     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1195     break;
1196     case 0x14: /* FLOGN */
1197     regs.fp[reg] = log (src);
1198     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1199     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1200     break;
1201     case 0x15: /* FLOG10 */
1202     regs.fp[reg] = log10 (src);
1203     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1204     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1205     break;
1206     case 0x16: /* FLOG2 */
1207     regs.fp[reg] = log (src) / log (2.0);
1208     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1209     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1210     break;
1211     case 0x18: /* FABS */
1212     regs.fp[reg] = src < 0 ? -src : src;
1213     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1214     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1215     break;
1216     case 0x19: /* FCOSH */
1217     regs.fp[reg] = cosh(src);
1218     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1219     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1220     break;
1221     case 0x1a: /* FNEG */
1222     regs.fp[reg] = -src;
1223     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1224     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1225     break;
1226     case 0x1c: /* FACOS */
1227     regs.fp[reg] = acos(src);
1228     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1229     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1230     break;
1231     case 0x1d: /* FCOS */
1232     regs.fp[reg] = cos(src);
1233     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1234     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1235     break;
1236     case 0x1e: /* FGETEXP */
1237     {
1238     int expon;
1239     frexp (src, &expon);
1240     regs.fp[reg] = (double) (expon - 1);
1241     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1242     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1243     }
1244     break;
1245     case 0x1f: /* FGETMAN */
1246     {
1247     int expon;
1248     regs.fp[reg] = frexp (src, &expon) * 2.0;
1249     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1250     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1251     }
1252     break;
1253     case 0x20: /* FDIV */
1254     regs.fp[reg] /= src;
1255     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1256     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1257     break;
1258     case 0x21: /* FMOD */
1259     regs.fp[reg] = regs.fp[reg] -
1260     (double) ((int) (regs.fp[reg] / src)) * src;
1261     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1262     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1263     break;
1264     case 0x22: /* FADD */
1265     regs.fp[reg] += src;
1266     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1267     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1268     break;
1269     case 0x23: /* FMUL */
1270     regs.fp[reg] *= src;
1271     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1272     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1273     break;
1274     case 0x24: /* FSGLDIV */
1275     regs.fp[reg] /= src;
1276     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1277     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1278     break;
1279     case 0x25: /* FREM */
1280     regs.fp[reg] = regs.fp[reg] -
1281     (double) ((int) (regs.fp[reg] / src + 0.5)) * src;
1282     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1283     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1284     break;
1285     case 0x26: /* FSCALE */
1286     regs.fp[reg] *= exp (log (2.0) * src);
1287     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1288     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1289     break;
1290     case 0x27: /* FSGLMUL */
1291     regs.fp[reg] *= src;
1292     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1293     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1294     break;
1295     case 0x28: /* FSUB */
1296     regs.fp[reg] -= src;
1297     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1298     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1299     break;
1300     case 0x30: /* FSINCOS */
1301     case 0x31:
1302     case 0x32:
1303     case 0x33:
1304     case 0x34:
1305     case 0x35:
1306     case 0x36:
1307     case 0x37:
1308     regs.fp[reg] = sin (src);
1309     regs.fp[extra & 7] = cos(src);
1310     regs.fpsr = (regs.fp[reg] == 0 ? 0x4000000 : 0) |
1311     (regs.fp[reg] < 0 ? 0x8000000 : 0);
1312     break;
1313     case 0x38: /* FCMP */
1314     {
1315     double tmp = regs.fp[reg] - src;
1316     regs.fpsr = (tmp == 0 ? 0x4000000 : 0) |
1317     (tmp < 0 ? 0x8000000 : 0);
1318     }
1319     break;
1320     case 0x3a: /* FTST */
1321     regs.fpsr = (src == 0 ? 0x4000000 : 0) |
1322     (src < 0 ? 0x8000000 : 0);
1323     break;
1324     default:
1325     m68k_setpc (m68k_getpc () - 4);
1326     op_illg (opcode);
1327     break;
1328     }
1329     return;
1330     }
1331     m68k_setpc (m68k_getpc () - 4);
1332     op_illg (opcode);
1333     }
1334    
1335     #endif