op-i386.c

   1 /*
   2  *  i386 micro operations
   3  *
   4  *  Copyright (c) 2003 Fabrice Bellard
   5  *
   6  * This library is free software; you can redistribute it and/or
   7  * modify it under the terms of the GNU Lesser General Public
   8  * License as published by the Free Software Foundation; either
   9  * version 2 of the License, or (at your option) any later version.
  10  *
  11  * This library 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 GNU
  14  * Lesser General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU Lesser General Public
  17  * License along with this library; if not, write to the Free Software
  18  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  19  */
  20 #include "exec-i386.h"
  21
  22 /* NOTE: data are not static to force relocation generation by GCC */
  23
  24 uint8_t parity_table[256] = {
  25     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  26     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  27     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  28     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  29     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  30     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  31     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  32     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  33     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  34     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  35     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  36     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  37     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  38     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  39     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  40     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  41     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  42     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  43     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  44     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  45     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  46     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  47     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  48     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  49     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  50     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  51     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  52     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  53     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  54     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  55     CC_P, 0, 0, CC_P, 0, CC_P, CC_P, 0,
  56     0, CC_P, CC_P, 0, CC_P, 0, 0, CC_P,
  57 };
  58
  59 /* modulo 17 table */
  60 const uint8_t rclw_table[32] = {
  61     0, 1, 2, 3, 4, 5, 6, 7,
  62     8, 9,10,11,12,13,14,15,
  63    16, 0, 1, 2, 3, 4, 5, 6,
  64     7, 8, 9,10,11,12,13,14,
  65 };
  66
  67 /* modulo 9 table */
  68 const uint8_t rclb_table[32] = {
  69     0, 1, 2, 3, 4, 5, 6, 7,
  70     8, 0, 1, 2, 3, 4, 5, 6,
  71     7, 8, 0, 1, 2, 3, 4, 5,
  72     6, 7, 8, 0, 1, 2, 3, 4,
  73 };
  74
  75 #ifdef USE_X86LDOUBLE
  76 /* an array of Intel 80-bit FP constants, to be loaded via integer ops */
  77 typedef unsigned short f15ld[5];
  78 const f15ld f15rk[] =
  79 {
  80 /*0*/   {0x0000,0x0000,0x0000,0x0000,0x0000},
  81 /*1*/   {0x0000,0x0000,0x0000,0x8000,0x3fff},
  82 /*pi*/  {0xc235,0x2168,0xdaa2,0xc90f,0x4000},
  83 /*lg2*/ {0xf799,0xfbcf,0x9a84,0x9a20,0x3ffd},
  84 /*ln2*/ {0x79ac,0xd1cf,0x17f7,0xb172,0x3ffe},
  85 /*l2e*/ {0xf0bc,0x5c17,0x3b29,0xb8aa,0x3fff},
  86 /*l2t*/ {0x8afe,0xcd1b,0x784b,0xd49a,0x4000}
  87 };
  88 #else
  89 /* the same, 64-bit version */
  90 typedef unsigned short f15ld[4];
  91 const f15ld f15rk[] =
  92 {
  93 #ifndef WORDS_BIGENDIAN
  94 /*0*/   {0x0000,0x0000,0x0000,0x0000},
  95 /*1*/   {0x0000,0x0000,0x0000,0x3ff0},
  96 /*pi*/  {0x2d18,0x5444,0x21fb,0x4009},
  97 /*lg2*/ {0x79ff,0x509f,0x4413,0x3fd3},
  98 /*ln2*/ {0x39ef,0xfefa,0x2e42,0x3fe6},
  99 /*l2e*/ {0x82fe,0x652b,0x1547,0x3ff7},
 100 /*l2t*/ {0xa371,0x0979,0x934f,0x400a}
 101 #else
 102 /*0*/   {0x0000,0x0000,0x0000,0x0000},
 103 /*1*/   {0x3ff0,0x0000,0x0000,0x0000},
 104 /*pi*/  {0x4009,0x21fb,0x5444,0x2d18},
 105 /*lg2*/ {0x3fd3,0x4413,0x509f,0x79ff},
 106 /*ln2*/ {0x3fe6,0x2e42,0xfefa,0x39ef},
 107 /*l2e*/ {0x3ff7,0x1547,0x652b,0x82fe},
 108 /*l2t*/ {0x400a,0x934f,0x0979,0xa371}
 109 #endif
 110 };
 111 #endif
 112
 113 /* n must be a constant to be efficient */
 114 static inline int lshift(int x, int n)
 115 {
 116     if (n >= 0)
 117         return x << n;
 118     else
 119         return x >> (-n);
 120 }
 121
 122 /* we define the various pieces of code used by the JIT */
 123
 124 #define REG EAX
 125 #define REGNAME _EAX
 126 #include "opreg_template.h"
 127 #undef REG
 128 #undef REGNAME
 129
 130 #define REG ECX
 131 #define REGNAME _ECX
 132 #include "opreg_template.h"
 133 #undef REG
 134 #undef REGNAME
 135
 136 #define REG EDX
 137 #define REGNAME _EDX
 138 #include "opreg_template.h"
 139 #undef REG
 140 #undef REGNAME
 141
 142 #define REG EBX
 143 #define REGNAME _EBX
 144 #include "opreg_template.h"
 145 #undef REG
 146 #undef REGNAME
 147
 148 #define REG ESP
 149 #define REGNAME _ESP
 150 #include "opreg_template.h"
 151 #undef REG
 152 #undef REGNAME
 153
 154 #define REG EBP
 155 #define REGNAME _EBP
 156 #include "opreg_template.h"
 157 #undef REG
 158 #undef REGNAME
 159
 160 #define REG ESI
 161 #define REGNAME _ESI
 162 #include "opreg_template.h"
 163 #undef REG
 164 #undef REGNAME
 165
 166 #define REG EDI
 167 #define REGNAME _EDI
 168 #include "opreg_template.h"
 169 #undef REG
 170 #undef REGNAME
 171
 172 /* operations with flags */
 173
 174 void OPPROTO op_addl_T0_T1_cc(void)
 175 {
 176     CC_SRC = T0;
 177     T0 += T1;
 178     CC_DST = T0;
 179 }
 180
 181 void OPPROTO op_orl_T0_T1_cc(void)
 182 {
 183     T0 |= T1;
 184     CC_DST = T0;
 185 }
 186
 187 void OPPROTO op_andl_T0_T1_cc(void)
 188 {
 189     T0 &= T1;
 190     CC_DST = T0;
 191 }
 192
 193 void OPPROTO op_subl_T0_T1_cc(void)
 194 {
 195     CC_SRC = T0;
 196     T0 -= T1;
 197     CC_DST = T0;
 198 }
 199
 200 void OPPROTO op_xorl_T0_T1_cc(void)
 201 {
 202     T0 ^= T1;
 203     CC_DST = T0;
 204 }
 205
 206 void OPPROTO op_cmpl_T0_T1_cc(void)
 207 {
 208     CC_SRC = T0;
 209     CC_DST = T0 - T1;
 210 }
 211
 212 void OPPROTO op_negl_T0_cc(void)
 213 {
 214     CC_SRC = 0;
 215     T0 = -T0;
 216     CC_DST = T0;
 217 }
 218
 219 void OPPROTO op_incl_T0_cc(void)
 220 {
 221     CC_SRC = cc_table[CC_OP].compute_c();
 222     T0++;
 223     CC_DST = T0;
 224 }
 225
 226 void OPPROTO op_decl_T0_cc(void)
 227 {
 228     CC_SRC = cc_table[CC_OP].compute_c();
 229     T0--;
 230     CC_DST = T0;
 231 }
 232
 233 void OPPROTO op_testl_T0_T1_cc(void)
 234 {
 235     CC_DST = T0 & T1;
 236 }
 237
 238 /* operations without flags */
 239
 240 void OPPROTO op_addl_T0_T1(void)
 241 {
 242     T0 += T1;
 243 }
 244
 245 void OPPROTO op_orl_T0_T1(void)
 246 {
 247     T0 |= T1;
 248 }
 249
 250 void OPPROTO op_andl_T0_T1(void)
 251 {
 252     T0 &= T1;
 253 }
 254
 255 void OPPROTO op_subl_T0_T1(void)
 256 {
 257     T0 -= T1;
 258 }
 259
 260 void OPPROTO op_xorl_T0_T1(void)
 261 {
 262     T0 ^= T1;
 263 }
 264
 265 void OPPROTO op_negl_T0(void)
 266 {
 267     T0 = -T0;
 268 }
 269
 270 void OPPROTO op_incl_T0(void)
 271 {
 272     T0++;
 273 }
 274
 275 void OPPROTO op_decl_T0(void)
 276 {
 277     T0--;
 278 }
 279
 280 void OPPROTO op_notl_T0(void)
 281 {
 282     T0 = ~T0;
 283 }
 284
 285 void OPPROTO op_bswapl_T0(void)
 286 {
 287     T0 = bswap32(T0);
 288 }
 289
 290 /* multiply/divide */
 291 void OPPROTO op_mulb_AL_T0(void)
 292 {
 293     unsigned int res;
 294     res = (uint8_t)EAX * (uint8_t)T0;
 295     EAX = (EAX & 0xffff0000) | res;
 296     CC_SRC = (res & 0xff00);
 297 }
 298
 299 void OPPROTO op_imulb_AL_T0(void)
 300 {
 301     int res;
 302     res = (int8_t)EAX * (int8_t)T0;
 303     EAX = (EAX & 0xffff0000) | (res & 0xffff);
 304     CC_SRC = (res != (int8_t)res);
 305 }
 306
 307 void OPPROTO op_mulw_AX_T0(void)
 308 {
 309     unsigned int res;
 310     res = (uint16_t)EAX * (uint16_t)T0;
 311     EAX = (EAX & 0xffff0000) | (res & 0xffff);
 312     EDX = (EDX & 0xffff0000) | ((res >> 16) & 0xffff);
 313     CC_SRC = res >> 16;
 314 }
 315
 316 void OPPROTO op_imulw_AX_T0(void)
 317 {
 318     int res;
 319     res = (int16_t)EAX * (int16_t)T0;
 320     EAX = (EAX & 0xffff0000) | (res & 0xffff);
 321     EDX = (EDX & 0xffff0000) | ((res >> 16) & 0xffff);
 322     CC_SRC = (res != (int16_t)res);
 323 }
 324
 325 void OPPROTO op_mull_EAX_T0(void)
 326 {
 327     uint64_t res;
 328     res = (uint64_t)((uint32_t)EAX) * (uint64_t)((uint32_t)T0);
 329     EAX = res;
 330     EDX = res >> 32;
 331     CC_SRC = res >> 32;
 332 }
 333
 334 void OPPROTO op_imull_EAX_T0(void)
 335 {
 336     int64_t res;
 337     res = (int64_t)((int32_t)EAX) * (int64_t)((int32_t)T0);
 338     EAX = res;
 339     EDX = res >> 32;
 340     CC_SRC = (res != (int32_t)res);
 341 }
 342
 343 void OPPROTO op_imulw_T0_T1(void)
 344 {
 345     int res;
 346     res = (int16_t)T0 * (int16_t)T1;
 347     T0 = res;
 348     CC_SRC = (res != (int16_t)res);
 349 }
 350
 351 void OPPROTO op_imull_T0_T1(void)
 352 {
 353     int64_t res;
 354     res = (int64_t)((int32_t)T0) * (int64_t)((int32_t)T1);
 355     T0 = res;
 356     CC_SRC = (res != (int32_t)res);
 357 }
 358
 359 /* division, flags are undefined */
 360 /* XXX: add exceptions for overflow */
 361 void OPPROTO op_divb_AL_T0(void)
 362 {
 363     unsigned int num, den, q, r;
 364
 365     num = (EAX & 0xffff);
 366     den = (T0 & 0xff);
 367     if (den == 0)
 368         raise_exception(EXCP00_DIVZ);
 369     q = (num / den) & 0xff;
 370     r = (num % den) & 0xff;
 371     EAX = (EAX & 0xffff0000) | (r << 8) | q;
 372 }
 373
 374 void OPPROTO op_idivb_AL_T0(void)
 375 {
 376     int num, den, q, r;
 377
 378     num = (int16_t)EAX;
 379     den = (int8_t)T0;
 380     if (den == 0)
 381         raise_exception(EXCP00_DIVZ);
 382     q = (num / den) & 0xff;
 383     r = (num % den) & 0xff;
 384     EAX = (EAX & 0xffff0000) | (r << 8) | q;
 385 }
 386
 387 void OPPROTO op_divw_AX_T0(void)
 388 {
 389     unsigned int num, den, q, r;
 390
 391     num = (EAX & 0xffff) | ((EDX & 0xffff) << 16);
 392     den = (T0 & 0xffff);
 393     if (den == 0)
 394         raise_exception(EXCP00_DIVZ);
 395     q = (num / den) & 0xffff;
 396     r = (num % den) & 0xffff;
 397     EAX = (EAX & 0xffff0000) | q;
 398     EDX = (EDX & 0xffff0000) | r;
 399 }
 400
 401 void OPPROTO op_idivw_AX_T0(void)
 402 {
 403     int num, den, q, r;
 404
 405     num = (EAX & 0xffff) | ((EDX & 0xffff) << 16);
 406     den = (int16_t)T0;
 407     if (den == 0)
 408         raise_exception(EXCP00_DIVZ);
 409     q = (num / den) & 0xffff;
 410     r = (num % den) & 0xffff;
 411     EAX = (EAX & 0xffff0000) | q;
 412     EDX = (EDX & 0xffff0000) | r;
 413 }
 414
 415 void OPPROTO op_divl_EAX_T0(void)
 416 {
 417     unsigned int den, q, r;
 418     uint64_t num;
 419
 420     num = EAX | ((uint64_t)EDX << 32);
 421     den = T0;
 422     if (den == 0)
 423         raise_exception(EXCP00_DIVZ);
 424     q = (num / den);
 425     r = (num % den);
 426     EAX = q;
 427     EDX = r;
 428 }
 429
 430 void OPPROTO op_idivl_EAX_T0(void)
 431 {
 432     int den, q, r;
 433     int64_t num;
 434
 435     num = EAX | ((uint64_t)EDX << 32);
 436     den = T0;
 437     if (den == 0)
 438         raise_exception(EXCP00_DIVZ);
 439     q = (num / den);
 440     r = (num % den);
 441     EAX = q;
 442     EDX = r;
 443 }
 444
 445 /* constant load & misc op */
 446
 447 void OPPROTO op_movl_T0_im(void)
 448 {
 449     T0 = PARAM1;
 450 }
 451
 452 void OPPROTO op_addl_T0_im(void)
 453 {
 454     T0 += PARAM1;
 455 }
 456
 457 void OPPROTO op_andl_T0_ffff(void)
 458 {
 459     T0 = T0 & 0xffff;
 460 }
 461
 462 void OPPROTO op_movl_T0_T1(void)
 463 {
 464     T0 = T1;
 465 }
 466
 467 void OPPROTO op_movl_T1_im(void)
 468 {
 469     T1 = PARAM1;
 470 }
 471
 472 void OPPROTO op_addl_T1_im(void)
 473 {
 474     T1 += PARAM1;
 475 }
 476
 477 void OPPROTO op_movl_T1_A0(void)
 478 {
 479     T1 = A0;
 480 }
 481
 482 void OPPROTO op_movl_A0_im(void)
 483 {
 484     A0 = PARAM1;
 485 }
 486
 487 void OPPROTO op_addl_A0_im(void)
 488 {
 489     A0 += PARAM1;
 490 }
 491
 492 void OPPROTO op_addl_A0_AL(void)
 493 {
 494     A0 += (EAX & 0xff);
 495 }
 496
 497 void OPPROTO op_andl_A0_ffff(void)
 498 {
 499     A0 = A0 & 0xffff;
 500 }
 501
 502 /* memory access */
 503
 504 void OPPROTO op_ldub_T0_A0(void)
 505 {
 506     T0 = ldub((uint8_t *)A0);
 507 }
 508
 509 void OPPROTO op_ldsb_T0_A0(void)
 510 {
 511     T0 = ldsb((int8_t *)A0);
 512 }
 513
 514 void OPPROTO op_lduw_T0_A0(void)
 515 {
 516     T0 = lduw((uint8_t *)A0);
 517 }
 518
 519 void OPPROTO op_ldsw_T0_A0(void)
 520 {
 521     T0 = ldsw((int8_t *)A0);
 522 }
 523
 524 void OPPROTO op_ldl_T0_A0(void)
 525 {
 526     T0 = ldl((uint8_t *)A0);
 527 }
 528
 529 void OPPROTO op_ldub_T1_A0(void)
 530 {
 531     T1 = ldub((uint8_t *)A0);
 532 }
 533
 534 void OPPROTO op_ldsb_T1_A0(void)
 535 {
 536     T1 = ldsb((int8_t *)A0);
 537 }
 538
 539 void OPPROTO op_lduw_T1_A0(void)
 540 {
 541     T1 = lduw((uint8_t *)A0);
 542 }
 543
 544 void OPPROTO op_ldsw_T1_A0(void)
 545 {
 546     T1 = ldsw((int8_t *)A0);
 547 }
 548
 549 void OPPROTO op_ldl_T1_A0(void)
 550 {
 551     T1 = ldl((uint8_t *)A0);
 552 }
 553
 554 void OPPROTO op_stb_T0_A0(void)
 555 {
 556     stb((uint8_t *)A0, T0);
 557 }
 558
 559 void OPPROTO op_stw_T0_A0(void)
 560 {
 561     stw((uint8_t *)A0, T0);
 562 }
 563
 564 void OPPROTO op_stl_T0_A0(void)
 565 {
 566     stl((uint8_t *)A0, T0);
 567 }
 568
 569 /* used for bit operations */
 570
 571 void OPPROTO op_add_bitw_A0_T1(void)
 572 {
 573     A0 += ((int32_t)T1 >> 4) << 1;
 574 }
 575
 576 void OPPROTO op_add_bitl_A0_T1(void)
 577 {
 578     A0 += ((int32_t)T1 >> 5) << 2;
 579 }
 580
 581 /* indirect jump */
 582
 583 void OPPROTO op_jmp_T0(void)
 584 {
 585     EIP = T0;
 586 }
 587
 588 void OPPROTO op_jmp_im(void)
 589 {
 590     EIP = PARAM1;
 591 }
 592
 593 void OPPROTO op_int_im(void)
 594 {
 595     EIP = PARAM1;
 596     raise_exception(EXCP0D_GPF);
 597 }
 598
 599 void OPPROTO op_int3(void)
 600 {
 601     EIP = PARAM1;
 602     raise_exception(EXCP03_INT3);
 603 }
 604
 605 void OPPROTO op_into(void)
 606 {
 607     int eflags;
 608     eflags = cc_table[CC_OP].compute_all();
 609     if (eflags & CC_O) {
 610         EIP = PARAM1;
 611         raise_exception(EXCP04_INTO);
 612     } else {
 613         EIP = PARAM2;
 614     }
 615 }
 616
 617 /* string ops */
 618
 619 #define ldul ldl
 620
 621 #define SHIFT 0
 622 #include "ops_template.h"
 623 #undef SHIFT
 624
 625 #define SHIFT 1
 626 #include "ops_template.h"
 627 #undef SHIFT
 628
 629 #define SHIFT 2
 630 #include "ops_template.h"
 631 #undef SHIFT
 632
 633 /* sign extend */
 634
 635 void OPPROTO op_movsbl_T0_T0(void)
 636 {
 637     T0 = (int8_t)T0;
 638 }
 639
 640 void OPPROTO op_movzbl_T0_T0(void)
 641 {
 642     T0 = (uint8_t)T0;
 643 }
 644
 645 void OPPROTO op_movswl_T0_T0(void)
 646 {
 647     T0 = (int16_t)T0;
 648 }
 649
 650 void OPPROTO op_movzwl_T0_T0(void)
 651 {
 652     T0 = (uint16_t)T0;
 653 }
 654
 655 void OPPROTO op_movswl_EAX_AX(void)
 656 {
 657     EAX = (int16_t)EAX;
 658 }
 659
 660 void OPPROTO op_movsbw_AX_AL(void)
 661 {
 662     EAX = (EAX & 0xffff0000) | ((int8_t)EAX & 0xffff);
 663 }
 664
 665 void OPPROTO op_movslq_EDX_EAX(void)
 666 {
 667     EDX = (int32_t)EAX >> 31;
 668 }
 669
 670 void OPPROTO op_movswl_DX_AX(void)
 671 {
 672     EDX = (EDX & 0xffff0000) | (((int16_t)EAX >> 15) & 0xffff);
 673 }
 674
 675 /* push/pop */
 676
 677 void op_pushl_T0(void)
 678 {
 679     uint32_t offset;
 680     offset = ESP - 4;
 681     stl((void *)offset, T0);
 682     /* modify ESP after to handle exceptions correctly */
 683     ESP = offset;
 684 }
 685
 686 void op_pushw_T0(void)
 687 {
 688     uint32_t offset;
 689     offset = ESP - 2;
 690     stw((void *)offset, T0);
 691     /* modify ESP after to handle exceptions correctly */
 692     ESP = offset;
 693 }
 694
 695 void op_pushl_ss32_T0(void)
 696 {
 697     uint32_t offset;
 698     offset = ESP - 4;
 699     stl(env->seg_cache[R_SS].base + offset, T0);
 700     /* modify ESP after to handle exceptions correctly */
 701     ESP = offset;
 702 }
 703
 704 void op_pushw_ss32_T0(void)
 705 {
 706     uint32_t offset;
 707     offset = ESP - 2;
 708     stw(env->seg_cache[R_SS].base + offset, T0);
 709     /* modify ESP after to handle exceptions correctly */
 710     ESP = offset;
 711 }
 712
 713 void op_pushl_ss16_T0(void)
 714 {
 715     uint32_t offset;
 716     offset = (ESP - 4) & 0xffff;
 717     stl(env->seg_cache[R_SS].base + offset, T0);
 718     /* modify ESP after to handle exceptions correctly */
 719     ESP = (ESP & ~0xffff) | offset;
 720 }
 721
 722 void op_pushw_ss16_T0(void)
 723 {
 724     uint32_t offset;
 725     offset = (ESP - 2) & 0xffff;
 726     stw(env->seg_cache[R_SS].base + offset, T0);
 727     /* modify ESP after to handle exceptions correctly */
 728     ESP = (ESP & ~0xffff) | offset;
 729 }
 730
 731 /* NOTE: ESP update is done after */
 732 void op_popl_T0(void)
 733 {
 734     T0 = ldl((void *)ESP);
 735 }
 736
 737 void op_popw_T0(void)
 738 {
 739     T0 = lduw((void *)ESP);
 740 }
 741
 742 void op_popl_ss32_T0(void)
 743 {
 744     T0 = ldl(env->seg_cache[R_SS].base + ESP);
 745 }
 746
 747 void op_popw_ss32_T0(void)
 748 {
 749     T0 = lduw(env->seg_cache[R_SS].base + ESP);
 750 }
 751
 752 void op_popl_ss16_T0(void)
 753 {
 754     T0 = ldl(env->seg_cache[R_SS].base + (ESP & 0xffff));
 755 }
 756
 757 void op_popw_ss16_T0(void)
 758 {
 759     T0 = lduw(env->seg_cache[R_SS].base + (ESP & 0xffff));
 760 }
 761
 762 void op_addl_ESP_4(void)
 763 {
 764     ESP += 4;
 765 }
 766
 767 void op_addl_ESP_2(void)
 768 {
 769     ESP += 2;
 770 }
 771
 772 void op_addw_ESP_4(void)
 773 {
 774     ESP = (ESP & ~0xffff) | ((ESP + 4) & 0xffff);
 775 }
 776
 777 void op_addw_ESP_2(void)
 778 {
 779     ESP = (ESP & ~0xffff) | ((ESP + 2) & 0xffff);
 780 }
 781
 782 void op_addl_ESP_im(void)
 783 {
 784     ESP += PARAM1;
 785 }
 786
 787 void op_addw_ESP_im(void)
 788 {
 789     ESP = (ESP & ~0xffff) | ((ESP + PARAM1) & 0xffff);
 790 }
 791
 792 /* rdtsc */
 793 #ifndef __i386__
 794 uint64_t emu_time;
 795 #endif
 796 void op_rdtsc(void)
 797 {
 798     uint64_t val;
 799 #ifdef __i386__
 800     asm("rdtsc" : "=A" (val));
 801 #else
 802     /* better than nothing: the time increases */
 803     val = emu_time++;
 804 #endif
 805     EAX = val;
 806     EDX = val >> 32;
 807 }
 808
 809 /* bcd */
 810
 811 /* XXX: exception */
 812 void OPPROTO op_aam(void)
 813 {
 814     int base = PARAM1;
 815     int al, ah;
 816     al = EAX & 0xff;
 817     ah = al / base;
 818     al = al % base;
 819     EAX = (EAX & ~0xffff) | al | (ah << 8);
 820     CC_DST = al;
 821 }
 822
 823 void OPPROTO op_aad(void)
 824 {
 825     int base = PARAM1;
 826     int al, ah;
 827     al = EAX & 0xff;
 828     ah = (EAX >> 8) & 0xff;
 829     al = ((ah * base) + al) & 0xff;
 830     EAX = (EAX & ~0xffff) | al;
 831     CC_DST = al;
 832 }
 833
 834 void OPPROTO op_aaa(void)
 835 {
 836     int icarry;
 837     int al, ah, af;
 838     int eflags;
 839
 840     eflags = cc_table[CC_OP].compute_all();
 841     af = eflags & CC_A;
 842     al = EAX & 0xff;
 843     ah = (EAX >> 8) & 0xff;
 844
 845     icarry = (al > 0xf9);
 846     if (((al & 0x0f) > 9 ) || af) {
 847         al = (al + 6) & 0x0f;
 848         ah = (ah + 1 + icarry) & 0xff;
 849         eflags |= CC_C | CC_A;
 850     } else {
 851         eflags &= ~(CC_C | CC_A);
 852         al &= 0x0f;
 853     }
 854     EAX = (EAX & ~0xffff) | al | (ah << 8);
 855     CC_SRC = eflags;
 856 }
 857
 858 void OPPROTO op_aas(void)
 859 {
 860     int icarry;
 861     int al, ah, af;
 862     int eflags;
 863
 864     eflags = cc_table[CC_OP].compute_all();
 865     af = eflags & CC_A;
 866     al = EAX & 0xff;
 867     ah = (EAX >> 8) & 0xff;
 868
 869     icarry = (al < 6);
 870     if (((al & 0x0f) > 9 ) || af) {
 871         al = (al - 6) & 0x0f;
 872         ah = (ah - 1 - icarry) & 0xff;
 873         eflags |= CC_C | CC_A;
 874     } else {
 875         eflags &= ~(CC_C | CC_A);
 876         al &= 0x0f;
 877     }
 878     EAX = (EAX & ~0xffff) | al | (ah << 8);
 879     CC_SRC = eflags;
 880 }
 881
 882 void OPPROTO op_daa(void)
 883 {
 884     int al, af, cf;
 885     int eflags;
 886
 887     eflags = cc_table[CC_OP].compute_all();
 888     cf = eflags & CC_C;
 889     af = eflags & CC_A;
 890     al = EAX & 0xff;
 891
 892     eflags = 0;
 893     if (((al & 0x0f) > 9 ) || af) {
 894         al = (al + 6) & 0xff;
 895         eflags |= CC_A;
 896     }
 897     if ((al > 0x9f) || cf) {
 898         al = (al + 0x60) & 0xff;
 899         eflags |= CC_C;
 900     }
 901     EAX = (EAX & ~0xff) | al;
 902     /* well, speed is not an issue here, so we compute the flags by hand */
 903     eflags |= (al == 0) << 6; /* zf */
 904     eflags |= parity_table[al]; /* pf */
 905     eflags |= (al & 0x80); /* sf */
 906     CC_SRC = eflags;
 907 }
 908
 909 void OPPROTO op_das(void)
 910 {
 911     int al, al1, af, cf;
 912     int eflags;
 913
 914     eflags = cc_table[CC_OP].compute_all();
 915     cf = eflags & CC_C;
 916     af = eflags & CC_A;
 917     al = EAX & 0xff;
 918
 919     eflags = 0;
 920     al1 = al;
 921     if (((al & 0x0f) > 9 ) || af) {
 922         eflags |= CC_A;
 923         if (al < 6 || cf)
 924             eflags |= CC_C;
 925         al = (al - 6) & 0xff;
 926     }
 927     if ((al1 > 0x99) || cf) {
 928         al = (al - 0x60) & 0xff;
 929         eflags |= CC_C;
 930     }
 931     EAX = (EAX & ~0xff) | al;
 932     /* well, speed is not an issue here, so we compute the flags by hand */
 933     eflags |= (al == 0) << 6; /* zf */
 934     eflags |= parity_table[al]; /* pf */
 935     eflags |= (al & 0x80); /* sf */
 936     CC_SRC = eflags;
 937 }
 938
 939 /* segment handling */
 940
 941 void load_seg(int seg_reg, int selector)
 942 {
 943     SegmentCache *sc;
 944     SegmentDescriptorTable *dt;
 945     int index;
 946     uint32_t e1, e2;
 947     uint8_t *ptr;
 948
 949     env->segs[seg_reg] = selector;
 950     sc = &env->seg_cache[seg_reg];
 951     if (env->vm86) {
 952         sc->base = (void *)(selector << 4);
 953         sc->limit = 0xffff;
 954         sc->seg_32bit = 0;
 955     } else {
 956         if (selector & 0x4)
 957             dt = &env->ldt;
 958         else
 959             dt = &env->gdt;
 960         index = selector & ~7;
 961         if ((index + 7) > dt->limit)
 962             raise_exception(EXCP0D_GPF);
 963         ptr = dt->base + index;
 964         e1 = ldl(ptr);
 965         e2 = ldl(ptr + 4);
 966         sc->base = (void *)((e1 >> 16) | ((e2 & 0xff) << 16) | (e2 & 0xff000000));
 967         sc->limit = (e1 & 0xffff) | (e2 & 0x000f0000);
 968         if (e2 & (1 << 23))
 969             sc->limit = (sc->limit << 12) | 0xfff;
 970         sc->seg_32bit = (e2 >> 22) & 1;
 971 #if 0
 972         fprintf(logfile, "load_seg: sel=0x%04x base=0x%08lx limit=0x%08lx seg_32bit=%d\n",
 973                 selector, (unsigned long)sc->base, sc->limit, sc->seg_32bit);
 974 #endif
 975     }
 976 }
 977
 978 void OPPROTO op_movl_seg_T0(void)
 979 {
 980     load_seg(PARAM1, T0 & 0xffff);
 981 }
 982
 983 void OPPROTO op_movl_T0_seg(void)
 984 {
 985     T0 = env->segs[PARAM1];
 986 }
 987
 988 void OPPROTO op_addl_A0_seg(void)
 989 {
 990     A0 += *(unsigned long *)((char *)env + PARAM1);
 991 }
 992
 993 /* flags handling */
 994
 995 /* slow jumps cases (compute x86 flags) */
 996 void OPPROTO op_jo_cc(void)
 997 {
 998     int eflags;
 999     eflags = cc_table[CC_OP].compute_all();
1000     if (eflags & CC_O)
1001         EIP = PARAM1;
1002     else
1003         EIP = PARAM2;
1004     FORCE_RET();
1005 }
1006
1007 void OPPROTO op_jb_cc(void)
1008 {
1009     if (cc_table[CC_OP].compute_c())
1010         EIP = PARAM1;
1011     else
1012         EIP = PARAM2;
1013     FORCE_RET();
1014 }
1015
1016 void OPPROTO op_jz_cc(void)
1017 {
1018     int eflags;
1019     eflags = cc_table[CC_OP].compute_all();
1020     if (eflags & CC_Z)
1021         EIP = PARAM1;
1022     else
1023         EIP = PARAM2;
1024     FORCE_RET();
1025 }
1026
1027 void OPPROTO op_jbe_cc(void)
1028 {
1029     int eflags;
1030     eflags = cc_table[CC_OP].compute_all();
1031     if (eflags & (CC_Z | CC_C))
1032         EIP = PARAM1;
1033     else
1034         EIP = PARAM2;
1035     FORCE_RET();
1036 }
1037
1038 void OPPROTO op_js_cc(void)
1039 {
1040     int eflags;
1041     eflags = cc_table[CC_OP].compute_all();
1042     if (eflags & CC_S)
1043         EIP = PARAM1;
1044     else
1045         EIP = PARAM2;
1046     FORCE_RET();
1047 }
1048
1049 void OPPROTO op_jp_cc(void)
1050 {
1051     int eflags;
1052     eflags = cc_table[CC_OP].compute_all();
1053     if (eflags & CC_P)
1054         EIP = PARAM1;
1055     else
1056         EIP = PARAM2;
1057     FORCE_RET();
1058 }
1059
1060 void OPPROTO op_jl_cc(void)
1061 {
1062     int eflags;
1063     eflags = cc_table[CC_OP].compute_all();
1064     if ((eflags ^ (eflags >> 4)) & 0x80)
1065         EIP = PARAM1;
1066     else
1067         EIP = PARAM2;
1068     FORCE_RET();
1069 }
1070
1071 void OPPROTO op_jle_cc(void)
1072 {
1073     int eflags;
1074     eflags = cc_table[CC_OP].compute_all();
1075     if (((eflags ^ (eflags >> 4)) & 0x80) || (eflags & CC_Z))
1076         EIP = PARAM1;
1077     else
1078         EIP = PARAM2;
1079     FORCE_RET();
1080 }
1081
1082 /* slow set cases (compute x86 flags) */
1083 void OPPROTO op_seto_T0_cc(void)
1084 {
1085     int eflags;
1086     eflags = cc_table[CC_OP].compute_all();
1087     T0 = (eflags >> 11) & 1;
1088 }
1089
1090 void OPPROTO op_setb_T0_cc(void)
1091 {
1092     T0 = cc_table[CC_OP].compute_c();
1093 }
1094
1095 void OPPROTO op_setz_T0_cc(void)
1096 {
1097     int eflags;
1098     eflags = cc_table[CC_OP].compute_all();
1099     T0 = (eflags >> 6) & 1;
1100 }
1101
1102 void OPPROTO op_setbe_T0_cc(void)
1103 {
1104     int eflags;
1105     eflags = cc_table[CC_OP].compute_all();
1106     T0 = (eflags & (CC_Z | CC_C)) != 0;
1107 }
1108
1109 void OPPROTO op_sets_T0_cc(void)
1110 {
1111     int eflags;
1112     eflags = cc_table[CC_OP].compute_all();
1113     T0 = (eflags >> 7) & 1;
1114 }
1115
1116 void OPPROTO op_setp_T0_cc(void)
1117 {
1118     int eflags;
1119     eflags = cc_table[CC_OP].compute_all();
1120     T0 = (eflags >> 2) & 1;
1121 }
1122
1123 void OPPROTO op_setl_T0_cc(void)
1124 {
1125     int eflags;
1126     eflags = cc_table[CC_OP].compute_all();
1127     T0 = ((eflags ^ (eflags >> 4)) >> 7) & 1;
1128 }
1129
1130 void OPPROTO op_setle_T0_cc(void)
1131 {
1132     int eflags;
1133     eflags = cc_table[CC_OP].compute_all();
1134     T0 = (((eflags ^ (eflags >> 4)) & 0x80) || (eflags & CC_Z)) != 0;
1135 }
1136
1137 void OPPROTO op_xor_T0_1(void)
1138 {
1139     T0 ^= 1;
1140 }
1141
1142 void OPPROTO op_set_cc_op(void)
1143 {
1144     CC_OP = PARAM1;
1145 }
1146
1147 void OPPROTO op_movl_eflags_T0(void)
1148 {
1149     CC_SRC = T0;
1150     DF = 1 - (2 * ((T0 >> 10) & 1));
1151 }
1152
1153 /* XXX: compute only O flag */
1154 void OPPROTO op_movb_eflags_T0(void)
1155 {
1156     int of;
1157     of = cc_table[CC_OP].compute_all() & CC_O;
1158     CC_SRC = T0 | of;
1159 }
1160
1161 void OPPROTO op_movl_T0_eflags(void)
1162 {
1163     T0 = cc_table[CC_OP].compute_all();
1164     T0 |= (DF & DIRECTION_FLAG);
1165 }
1166
1167 void OPPROTO op_cld(void)
1168 {
1169     DF = 1;
1170 }
1171
1172 void OPPROTO op_std(void)
1173 {
1174     DF = -1;
1175 }
1176
1177 void OPPROTO op_clc(void)
1178 {
1179     int eflags;
1180     eflags = cc_table[CC_OP].compute_all();
1181     eflags &= ~CC_C;
1182     CC_SRC = eflags;
1183 }
1184
1185 void OPPROTO op_stc(void)
1186 {
1187     int eflags;
1188     eflags = cc_table[CC_OP].compute_all();
1189     eflags |= CC_C;
1190     CC_SRC = eflags;
1191 }
1192
1193 void OPPROTO op_cmc(void)
1194 {
1195     int eflags;
1196     eflags = cc_table[CC_OP].compute_all();
1197     eflags ^= CC_C;
1198     CC_SRC = eflags;
1199 }
1200
1201 void OPPROTO op_salc(void)
1202 {
1203     int cf;
1204     cf = cc_table[CC_OP].compute_c();
1205     EAX = (EAX & ~0xff) | ((-cf) & 0xff);
1206 }
1207
1208 static int compute_all_eflags(void)
1209 {
1210     return CC_SRC;
1211 }
1212
1213 static int compute_c_eflags(void)
1214 {
1215     return CC_SRC & CC_C;
1216 }
1217
1218 static int compute_c_mul(void)
1219 {
1220     int cf;
1221     cf = (CC_SRC != 0);
1222     return cf;
1223 }
1224
1225 static int compute_all_mul(void)
1226 {
1227     int cf, pf, af, zf, sf, of;
1228     cf = (CC_SRC != 0);
1229     pf = 0; /* undefined */
1230     af = 0; /* undefined */
1231     zf = 0; /* undefined */
1232     sf = 0; /* undefined */
1233     of = cf << 11;
1234     return cf | pf | af | zf | sf | of;
1235 }
1236
1237 CCTable cc_table[CC_OP_NB] = {
1238     [CC_OP_DYNAMIC] = { /* should never happen */ },
1239
1240     [CC_OP_EFLAGS] = { compute_all_eflags, compute_c_eflags },
1241
1242     [CC_OP_MUL] = { compute_all_mul, compute_c_mul },
1243
1244     [CC_OP_ADDB] = { compute_all_addb, compute_c_addb },
1245     [CC_OP_ADDW] = { compute_all_addw, compute_c_addw  },
1246     [CC_OP_ADDL] = { compute_all_addl, compute_c_addl  },
1247
1248     [CC_OP_ADCB] = { compute_all_adcb, compute_c_adcb },
1249     [CC_OP_ADCW] = { compute_all_adcw, compute_c_adcw  },
1250     [CC_OP_ADCL] = { compute_all_adcl, compute_c_adcl  },
1251
1252     [CC_OP_SUBB] = { compute_all_subb, compute_c_subb  },
1253     [CC_OP_SUBW] = { compute_all_subw, compute_c_subw  },
1254     [CC_OP_SUBL] = { compute_all_subl, compute_c_subl  },
1255
1256     [CC_OP_SBBB] = { compute_all_sbbb, compute_c_sbbb  },
1257     [CC_OP_SBBW] = { compute_all_sbbw, compute_c_sbbw  },
1258     [CC_OP_SBBL] = { compute_all_sbbl, compute_c_sbbl  },
1259
1260     [CC_OP_LOGICB] = { compute_all_logicb, compute_c_logicb },
1261     [CC_OP_LOGICW] = { compute_all_logicw, compute_c_logicw },
1262     [CC_OP_LOGICL] = { compute_all_logicl, compute_c_logicl },
1263
1264     [CC_OP_INCB] = { compute_all_incb, compute_c_incl },
1265     [CC_OP_INCW] = { compute_all_incw, compute_c_incl },
1266     [CC_OP_INCL] = { compute_all_incl, compute_c_incl },
1267
1268     [CC_OP_DECB] = { compute_all_decb, compute_c_incl },
1269     [CC_OP_DECW] = { compute_all_decw, compute_c_incl },
1270     [CC_OP_DECL] = { compute_all_decl, compute_c_incl },
1271
1272     [CC_OP_SHLB] = { compute_all_shlb, compute_c_shll },
1273     [CC_OP_SHLW] = { compute_all_shlw, compute_c_shll },
1274     [CC_OP_SHLL] = { compute_all_shll, compute_c_shll },
1275
1276     [CC_OP_SARB] = { compute_all_sarb, compute_c_shll },
1277     [CC_OP_SARW] = { compute_all_sarw, compute_c_shll },
1278     [CC_OP_SARL] = { compute_all_sarl, compute_c_shll },
1279 };
1280
1281 /* floating point support */
1282
1283 #ifdef USE_X86LDOUBLE
1284 /* use long double functions */
1285 #define lrint lrintl
1286 #define llrint llrintl
1287 #define fabs fabsl
1288 #define sin sinl
1289 #define cos cosl
1290 #define sqrt sqrtl
1291 #define pow powl
1292 #define log logl
1293 #define tan tanl
1294 #define atan2 atan2l
1295 #define floor floorl
1296 #define ceil ceill
1297 #define rint rintl
1298 #endif
1299
1300 extern int lrint(CPU86_LDouble x);
1301 extern int64_t llrint(CPU86_LDouble x);
1302 extern CPU86_LDouble fabs(CPU86_LDouble x);
1303 extern CPU86_LDouble sin(CPU86_LDouble x);
1304 extern CPU86_LDouble cos(CPU86_LDouble x);
1305 extern CPU86_LDouble sqrt(CPU86_LDouble x);
1306 extern CPU86_LDouble pow(CPU86_LDouble, CPU86_LDouble);
1307 extern CPU86_LDouble log(CPU86_LDouble x);
1308 extern CPU86_LDouble tan(CPU86_LDouble x);
1309 extern CPU86_LDouble atan2(CPU86_LDouble, CPU86_LDouble);
1310 extern CPU86_LDouble floor(CPU86_LDouble x);
1311 extern CPU86_LDouble ceil(CPU86_LDouble x);
1312 extern CPU86_LDouble rint(CPU86_LDouble x);
1313
1314 #define RC_MASK         0xc00
1315 #define RC_NEAR         0x000
1316 #define RC_DOWN         0x400
1317 #define RC_UP           0x800
1318 #define RC_CHOP         0xc00
1319
1320 #define MAXTAN 9223372036854775808.0
1321
1322 #ifdef USE_X86LDOUBLE
1323
1324 /* only for x86 */
1325 typedef union {
1326     long double d;
1327     struct {
1328         unsigned long long lower;
1329         unsigned short upper;
1330     } l;
1331 } CPU86_LDoubleU;
1332
1333 /* the following deal with x86 long double-precision numbers */
1334 #define MAXEXPD 0x7fff
1335 #define EXPBIAS 16383
1336 #define EXPD(fp)        (fp.l.upper & 0x7fff)
1337 #define SIGND(fp)       ((fp.l.upper) & 0x8000)
1338 #define MANTD(fp)       (fp.l.lower)
1339 #define BIASEXPONENT(fp) fp.l.upper = (fp.l.upper & ~(0x7fff)) | EXPBIAS
1340
1341 #else
1342
1343 typedef union {
1344     double d;
1345 #ifndef WORDS_BIGENDIAN
1346     struct {
1347         unsigned long lower;
1348         long upper;
1349     } l;
1350 #else
1351     struct {
1352         long upper;
1353         unsigned long lower;
1354     } l;
1355 #endif
1356     long long ll;
1357 } CPU86_LDoubleU;
1358
1359 /* the following deal with IEEE double-precision numbers */
1360 #define MAXEXPD 0x7ff
1361 #define EXPBIAS 1023
1362 #define EXPD(fp)        (((fp.l.upper) >> 20) & 0x7FF)
1363 #define SIGND(fp)       ((fp.l.upper) & 0x80000000)
1364 #define MANTD(fp)       (fp.ll & ((1LL << 52) - 1))
1365 #define BIASEXPONENT(fp) fp.l.upper = (fp.l.upper & ~(0x7ff << 20)) | (EXPBIAS << 20)
1366 #endif
1367
1368 /* fp load FT0 */
1369
1370 void OPPROTO op_flds_FT0_A0(void)
1371 {
1372     FT0 = ldfl((void *)A0);
1373 }
1374
1375 void OPPROTO op_fldl_FT0_A0(void)
1376 {
1377     FT0 = ldfq((void *)A0);
1378 }
1379
1380 /* helpers are needed to avoid static constant reference. XXX: find a better way */
1381 #ifdef USE_INT_TO_FLOAT_HELPERS
1382
1383 void helper_fild_FT0_A0(void)
1384 {
1385     FT0 = (CPU86_LDouble)ldsw((void *)A0);
1386 }
1387
1388 void helper_fildl_FT0_A0(void)
1389 {
1390     FT0 = (CPU86_LDouble)((int32_t)ldl((void *)A0));
1391 }
1392
1393 void helper_fildll_FT0_A0(void)
1394 {
1395     FT0 = (CPU86_LDouble)((int64_t)ldq((void *)A0));
1396 }
1397
1398 void OPPROTO op_fild_FT0_A0(void)
1399 {
1400     helper_fild_FT0_A0();
1401 }
1402
1403 void OPPROTO op_fildl_FT0_A0(void)
1404 {
1405     helper_fildl_FT0_A0();
1406 }
1407
1408 void OPPROTO op_fildll_FT0_A0(void)
1409 {
1410     helper_fildll_FT0_A0();
1411 }
1412
1413 #else
1414
1415 void OPPROTO op_fild_FT0_A0(void)
1416 {
1417     FT0 = (CPU86_LDouble)ldsw((void *)A0);
1418 }
1419
1420 void OPPROTO op_fildl_FT0_A0(void)
1421 {
1422     FT0 = (CPU86_LDouble)((int32_t)ldl((void *)A0));
1423 }
1424
1425 void OPPROTO op_fildll_FT0_A0(void)
1426 {
1427     FT0 = (CPU86_LDouble)((int64_t)ldq((void *)A0));
1428 }
1429 #endif
1430
1431 /* fp load ST0 */
1432
1433 void OPPROTO op_flds_ST0_A0(void)
1434 {
1435     ST0 = ldfl((void *)A0);
1436 }
1437
1438 void OPPROTO op_fldl_ST0_A0(void)
1439 {
1440     ST0 = ldfq((void *)A0);
1441 }
1442
1443 #ifdef USE_X86LDOUBLE
1444 void OPPROTO op_fldt_ST0_A0(void)
1445 {
1446     ST0 = *(long double *)A0;
1447 }
1448 #else
1449 void helper_fldt_ST0_A0(void)
1450 {
1451     CPU86_LDoubleU temp;
1452     int upper, e;
1453     /* mantissa */
1454     upper = lduw((uint8_t *)A0 + 8);
1455     /* XXX: handle overflow ? */
1456     e = (upper & 0x7fff) - 16383 + EXPBIAS; /* exponent */
1457     e |= (upper >> 4) & 0x800; /* sign */
1458     temp.ll = ((ldq((void *)A0) >> 11) & ((1LL << 52) - 1)) | ((uint64_t)e << 52);
1459     ST0 = temp.d;
1460 }
1461
1462 void OPPROTO op_fldt_ST0_A0(void)
1463 {
1464     helper_fldt_ST0_A0();
1465 }
1466 #endif
1467
1468 /* helpers are needed to avoid static constant reference. XXX: find a better way */
1469 #ifdef USE_INT_TO_FLOAT_HELPERS
1470
1471 void helper_fild_ST0_A0(void)
1472 {
1473     ST0 = (CPU86_LDouble)ldsw((void *)A0);
1474 }
1475
1476 void helper_fildl_ST0_A0(void)
1477 {
1478     ST0 = (CPU86_LDouble)((int32_t)ldl((void *)A0));
1479 }
1480
1481 void helper_fildll_ST0_A0(void)
1482 {
1483     ST0 = (CPU86_LDouble)((int64_t)ldq((void *)A0));
1484 }
1485
1486 void OPPROTO op_fild_ST0_A0(void)
1487 {
1488     helper_fild_ST0_A0();
1489 }
1490
1491 void OPPROTO op_fildl_ST0_A0(void)
1492 {
1493     helper_fildl_ST0_A0();
1494 }
1495
1496 void OPPROTO op_fildll_ST0_A0(void)
1497 {
1498     helper_fildll_ST0_A0();
1499 }
1500
1501 #else
1502
1503 void OPPROTO op_fild_ST0_A0(void)
1504 {
1505     ST0 = (CPU86_LDouble)ldsw((void *)A0);
1506 }
1507
1508 void OPPROTO op_fildl_ST0_A0(void)
1509 {
1510     ST0 = (CPU86_LDouble)((int32_t)ldl((void *)A0));
1511 }
1512
1513 void OPPROTO op_fildll_ST0_A0(void)
1514 {
1515     ST0 = (CPU86_LDouble)((int64_t)ldq((void *)A0));
1516 }
1517
1518 #endif
1519
1520 /* fp store */
1521
1522 void OPPROTO op_fsts_ST0_A0(void)
1523 {
1524     stfl((void *)A0, (float)ST0);
1525 }
1526
1527 void OPPROTO op_fstl_ST0_A0(void)
1528 {
1529     stfq((void *)A0, (double)ST0);
1530 }
1531
1532 #ifdef USE_X86LDOUBLE
1533 void OPPROTO op_fstt_ST0_A0(void)
1534 {
1535     *(long double *)A0 = ST0;
1536 }
1537 #else
1538 void helper_fstt_ST0_A0(void)
1539 {
1540     CPU86_LDoubleU temp;
1541     int e;
1542     temp.d = ST0;
1543     /* mantissa */
1544     stq((void *)A0, (MANTD(temp) << 11) | (1LL << 63));
1545     /* exponent + sign */
1546     e = EXPD(temp) - EXPBIAS + 16383;
1547     e |= SIGND(temp) >> 16;
1548     stw((uint8_t *)A0 + 8, e);
1549 }
1550
1551 void OPPROTO op_fstt_ST0_A0(void)
1552 {
1553     helper_fstt_ST0_A0();
1554 }
1555 #endif
1556
1557 void OPPROTO op_fist_ST0_A0(void)
1558 {
1559     int val;
1560     val = lrint(ST0);
1561     stw((void *)A0, val);
1562 }
1563
1564 void OPPROTO op_fistl_ST0_A0(void)
1565 {
1566     int val;
1567     val = lrint(ST0);
1568     stl((void *)A0, val);
1569 }
1570
1571 void OPPROTO op_fistll_ST0_A0(void)
1572 {
1573     int64_t val;
1574     val = llrint(ST0);
1575     stq((void *)A0, val);
1576 }
1577
1578 /* BCD ops */
1579
1580 #define MUL10(iv) ( iv + iv + (iv << 3) )
1581
1582 void helper_fbld_ST0_A0(void)
1583 {
1584     uint8_t *seg;
1585     CPU86_LDouble fpsrcop;
1586     int m32i;
1587     unsigned int v;
1588
1589     /* in this code, seg/m32i will be used as temporary ptr/int */
1590     seg = (uint8_t *)A0 + 8;
1591     v = ldub(seg--);
1592     /* XXX: raise exception */
1593     if (v != 0)
1594         return;
1595     v = ldub(seg--);
1596     /* XXX: raise exception */
1597     if ((v & 0xf0) != 0)
1598         return;
1599     m32i = v;  /* <-- d14 */
1600     v = ldub(seg--);
1601     m32i = MUL10(m32i) + (v >> 4);  /* <-- val * 10 + d13 */
1602     m32i = MUL10(m32i) + (v & 0xf); /* <-- val * 10 + d12 */
1603     v = ldub(seg--);
1604     m32i = MUL10(m32i) + (v >> 4);  /* <-- val * 10 + d11 */
1605     m32i = MUL10(m32i) + (v & 0xf); /* <-- val * 10 + d10 */
1606     v = ldub(seg--);
1607     m32i = MUL10(m32i) + (v >> 4);  /* <-- val * 10 + d9 */
1608     m32i = MUL10(m32i) + (v & 0xf); /* <-- val * 10 + d8 */
1609     fpsrcop = ((CPU86_LDouble)m32i) * 100000000.0;
1610
1611     v = ldub(seg--);
1612     m32i = (v >> 4);  /* <-- d7 */
1613     m32i = MUL10(m32i) + (v & 0xf); /* <-- val * 10 + d6 */
1614     v = ldub(seg--);
1615     m32i = MUL10(m32i) + (v >> 4);  /* <-- val * 10 + d5 */
1616     m32i = MUL10(m32i) + (v & 0xf); /* <-- val * 10 + d4 */
1617     v = ldub(seg--);
1618     m32i = MUL10(m32i) + (v >> 4);  /* <-- val * 10 + d3 */
1619     m32i = MUL10(m32i) + (v & 0xf); /* <-- val * 10 + d2 */
1620     v = ldub(seg);
1621     m32i = MUL10(m32i) + (v >> 4);  /* <-- val * 10 + d1 */
1622     m32i = MUL10(m32i) + (v & 0xf); /* <-- val * 10 + d0 */
1623     fpsrcop += ((CPU86_LDouble)m32i);
1624     if ( ldub(seg+9) & 0x80 )
1625         fpsrcop = -fpsrcop;
1626     ST0 = fpsrcop;
1627 }
1628
1629 void OPPROTO op_fbld_ST0_A0(void)
1630 {
1631     helper_fbld_ST0_A0();
1632 }
1633
1634 void helper_fbst_ST0_A0(void)
1635 {
1636     CPU86_LDouble fptemp;
1637     CPU86_LDouble fpsrcop;
1638     int v;
1639     uint8_t *mem_ref, *mem_end;
1640
1641     fpsrcop = rint(ST0);
1642     mem_ref = (uint8_t *)A0;
1643     mem_end = mem_ref + 8;
1644     if ( fpsrcop < 0.0 ) {
1645         stw(mem_end, 0x8000);
1646         fpsrcop = -fpsrcop;
1647     } else {
1648         stw(mem_end, 0x0000);
1649     }
1650     while (mem_ref < mem_end) {
1651         if (fpsrcop == 0.0)
1652             break;
1653         fptemp = floor(fpsrcop/10.0);
1654         v = ((int)(fpsrcop - fptemp*10.0));
1655         if  (fptemp == 0.0)  {
1656             stb(mem_ref++, v);
1657             break;
1658         }
1659         fpsrcop = fptemp;
1660         fptemp = floor(fpsrcop/10.0);
1661         v |= (((int)(fpsrcop - fptemp*10.0)) << 4);
1662         stb(mem_ref++, v);
1663         fpsrcop = fptemp;
1664     }
1665     while (mem_ref < mem_end) {
1666         stb(mem_ref++, 0);
1667     }
1668 }
1669
1670 void OPPROTO op_fbst_ST0_A0(void)
1671 {
1672     helper_fbst_ST0_A0();
1673 }
1674
1675 /* FPU move */
1676
1677 static inline void fpush(void)
1678 {
1679     env->fpstt = (env->fpstt - 1) & 7;
1680     env->fptags[env->fpstt] = 0; /* validate stack entry */
1681 }
1682
1683 static inline void fpop(void)
1684 {
1685     env->fptags[env->fpstt] = 1; /* invvalidate stack entry */
1686     env->fpstt = (env->fpstt + 1) & 7;
1687 }
1688
1689 void OPPROTO op_fpush(void)
1690 {
1691     fpush();
1692 }
1693
1694 void OPPROTO op_fpop(void)
1695 {
1696     fpop();
1697 }
1698
1699 void OPPROTO op_fdecstp(void)
1700 {
1701     env->fpstt = (env->fpstt - 1) & 7;
1702     env->fpus &= (~0x4700);
1703 }
1704
1705 void OPPROTO op_fincstp(void)
1706 {
1707     env->fpstt = (env->fpstt + 1) & 7;
1708     env->fpus &= (~0x4700);
1709 }
1710
1711 void OPPROTO op_fmov_ST0_FT0(void)
1712 {
1713     ST0 = FT0;
1714 }
1715
1716 void OPPROTO op_fmov_FT0_STN(void)
1717 {
1718     FT0 = ST(PARAM1);
1719 }
1720
1721 void OPPROTO op_fmov_ST0_STN(void)
1722 {
1723     ST0 = ST(PARAM1);
1724 }
1725
1726 void OPPROTO op_fmov_STN_ST0(void)
1727 {
1728     ST(PARAM1) = ST0;
1729 }
1730
1731 void OPPROTO op_fxchg_ST0_STN(void)
1732 {
1733     CPU86_LDouble tmp;
1734     tmp = ST(PARAM1);
1735     ST(PARAM1) = ST0;
1736     ST0 = tmp;
1737 }
1738
1739 /* FPU operations */
1740
1741 /* XXX: handle nans */
1742 void OPPROTO op_fcom_ST0_FT0(void)
1743 {
1744     env->fpus &= (~0x4500);     /* (C3,C2,C0) <-- 000 */
1745     if (ST0 < FT0)
1746         env->fpus |= 0x100;     /* (C3,C2,C0) <-- 001 */
1747     else if (ST0 == FT0)
1748         env->fpus |= 0x4000; /* (C3,C2,C0) <-- 100 */
1749     FORCE_RET();
1750 }
1751
1752 /* XXX: handle nans */
1753 void OPPROTO op_fucom_ST0_FT0(void)
1754 {
1755     env->fpus &= (~0x4500);     /* (C3,C2,C0) <-- 000 */
1756     if (ST0 < FT0)
1757         env->fpus |= 0x100;     /* (C3,C2,C0) <-- 001 */
1758     else if (ST0 == FT0)
1759         env->fpus |= 0x4000; /* (C3,C2,C0) <-- 100 */
1760     FORCE_RET();
1761 }
1762
1763 void OPPROTO op_fadd_ST0_FT0(void)
1764 {
1765     ST0 += FT0;
1766 }
1767
1768 void OPPROTO op_fmul_ST0_FT0(void)
1769 {
1770     ST0 *= FT0;
1771 }
1772
1773 void OPPROTO op_fsub_ST0_FT0(void)
1774 {
1775     ST0 -= FT0;
1776 }
1777
1778 void OPPROTO op_fsubr_ST0_FT0(void)
1779 {
1780     ST0 = FT0 - ST0;
1781 }
1782
1783 void OPPROTO op_fdiv_ST0_FT0(void)
1784 {
1785     ST0 /= FT0;
1786 }
1787
1788 void OPPROTO op_fdivr_ST0_FT0(void)
1789 {
1790     ST0 = FT0 / ST0;
1791 }
1792
1793 /* fp operations between STN and ST0 */
1794
1795 void OPPROTO op_fadd_STN_ST0(void)
1796 {
1797     ST(PARAM1) += ST0;
1798 }
1799
1800 void OPPROTO op_fmul_STN_ST0(void)
1801 {
1802     ST(PARAM1) *= ST0;
1803 }
1804
1805 void OPPROTO op_fsub_STN_ST0(void)
1806 {
1807     ST(PARAM1) -= ST0;
1808 }
1809
1810 void OPPROTO op_fsubr_STN_ST0(void)
1811 {
1812     CPU86_LDouble *p;
1813     p = &ST(PARAM1);
1814     *p = ST0 - *p;
1815 }
1816
1817 void OPPROTO op_fdiv_STN_ST0(void)
1818 {
1819     ST(PARAM1) /= ST0;
1820 }
1821
1822 void OPPROTO op_fdivr_STN_ST0(void)
1823 {
1824     CPU86_LDouble *p;
1825     p = &ST(PARAM1);
1826     *p = ST0 / *p;
1827 }
1828
1829 /* misc FPU operations */
1830 void OPPROTO op_fchs_ST0(void)
1831 {
1832     ST0 = -ST0;
1833 }
1834
1835 void OPPROTO op_fabs_ST0(void)
1836 {
1837     ST0 = fabs(ST0);
1838 }
1839
1840 void helper_fxam_ST0(void)
1841 {
1842     CPU86_LDoubleU temp;
1843     int expdif;
1844
1845     temp.d = ST0;
1846
1847     env->fpus &= (~0x4700);  /* (C3,C2,C1,C0) <-- 0000 */
1848     if (SIGND(temp))
1849         env->fpus |= 0x200; /* C1 <-- 1 */
1850
1851     expdif = EXPD(temp);
1852     if (expdif == MAXEXPD) {
1853         if (MANTD(temp) == 0)
1854             env->fpus |=  0x500 /*Infinity*/;
1855         else
1856             env->fpus |=  0x100 /*NaN*/;
1857     } else if (expdif == 0) {
1858         if (MANTD(temp) == 0)
1859             env->fpus |=  0x4000 /*Zero*/;
1860         else
1861             env->fpus |= 0x4400 /*Denormal*/;
1862     } else {
1863         env->fpus |= 0x400;
1864     }
1865 }
1866
1867 void OPPROTO op_fxam_ST0(void)
1868 {
1869     helper_fxam_ST0();
1870 }
1871
1872 void OPPROTO op_fld1_ST0(void)
1873 {
1874     ST0 = *(CPU86_LDouble *)&f15rk[1];
1875 }
1876
1877 void OPPROTO op_fldl2t_ST0(void)
1878 {
1879     ST0 = *(CPU86_LDouble *)&f15rk[6];
1880 }
1881
1882 void OPPROTO op_fldl2e_ST0(void)
1883 {
1884     ST0 = *(CPU86_LDouble *)&f15rk[5];
1885 }
1886
1887 void OPPROTO op_fldpi_ST0(void)
1888 {
1889     ST0 = *(CPU86_LDouble *)&f15rk[2];
1890 }
1891
1892 void OPPROTO op_fldlg2_ST0(void)
1893 {
1894     ST0 = *(CPU86_LDouble *)&f15rk[3];
1895 }
1896
1897 void OPPROTO op_fldln2_ST0(void)
1898 {
1899     ST0 = *(CPU86_LDouble *)&f15rk[4];
1900 }
1901
1902 void OPPROTO op_fldz_ST0(void)
1903 {
1904     ST0 = *(CPU86_LDouble *)&f15rk[0];
1905 }
1906
1907 void OPPROTO op_fldz_FT0(void)
1908 {
1909     ST0 = *(CPU86_LDouble *)&f15rk[0];
1910 }
1911
1912 void helper_f2xm1(void)
1913 {
1914     ST0 = pow(2.0,ST0) - 1.0;
1915 }
1916
1917 void helper_fyl2x(void)
1918 {
1919     CPU86_LDouble fptemp;
1920
1921     fptemp = ST0;
1922     if (fptemp>0.0){
1923         fptemp = log(fptemp)/log(2.0);   /* log2(ST) */
1924         ST1 *= fptemp;
1925         fpop();
1926     } else {
1927         env->fpus &= (~0x4700);
1928         env->fpus |= 0x400;
1929     }
1930 }
1931
1932 void helper_fptan(void)
1933 {
1934     CPU86_LDouble fptemp;
1935
1936     fptemp = ST0;
1937     if((fptemp > MAXTAN)||(fptemp < -MAXTAN)) {
1938         env->fpus |= 0x400;
1939     } else {
1940         ST0 = tan(fptemp);
1941         fpush();
1942         ST0 = 1.0;
1943         env->fpus &= (~0x400);  /* C2 <-- 0 */
1944         /* the above code is for  |arg| < 2**52 only */
1945     }
1946 }
1947
1948 void helper_fpatan(void)
1949 {
1950     CPU86_LDouble fptemp, fpsrcop;
1951
1952     fpsrcop = ST1;
1953     fptemp = ST0;
1954     ST1 = atan2(fpsrcop,fptemp);
1955     fpop();
1956 }
1957
1958 void helper_fxtract(void)
1959 {
1960     CPU86_LDoubleU temp;
1961     unsigned int expdif;
1962
1963     temp.d = ST0;
1964     expdif = EXPD(temp) - EXPBIAS;
1965     /*DP exponent bias*/
1966     ST0 = expdif;
1967     fpush();
1968     BIASEXPONENT(temp);
1969     ST0 = temp.d;
1970 }
1971
1972 void helper_fprem1(void)
1973 {
1974     CPU86_LDouble dblq, fpsrcop, fptemp;
1975     CPU86_LDoubleU fpsrcop1, fptemp1;
1976     int expdif;
1977     int q;
1978
1979     fpsrcop = ST0;
1980     fptemp = ST1;
1981     fpsrcop1.d = fpsrcop;
1982     fptemp1.d = fptemp;
1983     expdif = EXPD(fpsrcop1) - EXPD(fptemp1);
1984     if (expdif < 53) {
1985         dblq = fpsrcop / fptemp;
1986         dblq = (dblq < 0.0)? ceil(dblq): floor(dblq);
1987         ST0 = fpsrcop - fptemp*dblq;
1988         q = (int)dblq; /* cutting off top bits is assumed here */
1989         env->fpus &= (~0x4700); /* (C3,C2,C1,C0) <-- 0000 */
1990                                 /* (C0,C1,C3) <-- (q2,q1,q0) */
1991         env->fpus |= (q&0x4) << 6; /* (C0) <-- q2 */
1992         env->fpus |= (q&0x2) << 8; /* (C1) <-- q1 */
1993         env->fpus |= (q&0x1) << 14; /* (C3) <-- q0 */
1994     } else {
1995         env->fpus |= 0x400;  /* C2 <-- 1 */
1996         fptemp = pow(2.0, expdif-50);
1997         fpsrcop = (ST0 / ST1) / fptemp;
1998         /* fpsrcop = integer obtained by rounding to the nearest */
1999         fpsrcop = (fpsrcop-floor(fpsrcop) < ceil(fpsrcop)-fpsrcop)?
2000             floor(fpsrcop): ceil(fpsrcop);
2001         ST0 -= (ST1 * fpsrcop * fptemp);
2002     }
2003 }
2004
2005 void helper_fprem(void)
2006 {
2007     CPU86_LDouble dblq, fpsrcop, fptemp;
2008     CPU86_LDoubleU fpsrcop1, fptemp1;
2009     int expdif;
2010     int q;
2011
2012     fpsrcop = ST0;
2013     fptemp = ST1;
2014     fpsrcop1.d = fpsrcop;
2015     fptemp1.d = fptemp;
2016     expdif = EXPD(fpsrcop1) - EXPD(fptemp1);
2017     if ( expdif < 53 ) {
2018         dblq = fpsrcop / fptemp;
2019         dblq = (dblq < 0.0)? ceil(dblq): floor(dblq);
2020         ST0 = fpsrcop - fptemp*dblq;
2021         q = (int)dblq; /* cutting off top bits is assumed here */
2022         env->fpus &= (~0x4700); /* (C3,C2,C1,C0) <-- 0000 */
2023                                 /* (C0,C1,C3) <-- (q2,q1,q0) */
2024         env->fpus |= (q&0x4) << 6; /* (C0) <-- q2 */
2025         env->fpus |= (q&0x2) << 8; /* (C1) <-- q1 */
2026         env->fpus |= (q&0x1) << 14; /* (C3) <-- q0 */
2027     } else {
2028         env->fpus |= 0x400;  /* C2 <-- 1 */
2029         fptemp = pow(2.0, expdif-50);
2030         fpsrcop = (ST0 / ST1) / fptemp;
2031         /* fpsrcop = integer obtained by chopping */
2032         fpsrcop = (fpsrcop < 0.0)?
2033             -(floor(fabs(fpsrcop))): floor(fpsrcop);
2034         ST0 -= (ST1 * fpsrcop * fptemp);
2035     }
2036 }
2037
2038 void helper_fyl2xp1(void)
2039 {
2040     CPU86_LDouble fptemp;
2041
2042     fptemp = ST0;
2043     if ((fptemp+1.0)>0.0) {
2044         fptemp = log(fptemp+1.0) / log(2.0); /* log2(ST+1.0) */
2045         ST1 *= fptemp;
2046         fpop();
2047     } else {
2048         env->fpus &= (~0x4700);
2049         env->fpus |= 0x400;
2050     }
2051 }
2052
2053 void helper_fsqrt(void)
2054 {
2055     CPU86_LDouble fptemp;
2056
2057     fptemp = ST0;
2058     if (fptemp<0.0) {
2059         env->fpus &= (~0x4700);  /* (C3,C2,C1,C0) <-- 0000 */
2060         env->fpus |= 0x400;
2061     }
2062     ST0 = sqrt(fptemp);
2063 }
2064
2065 void helper_fsincos(void)
2066 {
2067     CPU86_LDouble fptemp;
2068
2069     fptemp = ST0;
2070     if ((fptemp > MAXTAN)||(fptemp < -MAXTAN)) {
2071         env->fpus |= 0x400;
2072     } else {
2073         ST0 = sin(fptemp);
2074         fpush();
2075         ST0 = cos(fptemp);
2076         env->fpus &= (~0x400);  /* C2 <-- 0 */
2077         /* the above code is for  |arg| < 2**63 only */
2078     }
2079 }
2080
2081 void helper_frndint(void)
2082 {
2083     ST0 = rint(ST0);
2084 }
2085
2086 void helper_fscale(void)
2087 {
2088     CPU86_LDouble fpsrcop, fptemp;
2089
2090     fpsrcop = 2.0;
2091     fptemp = pow(fpsrcop,ST1);
2092     ST0 *= fptemp;
2093 }
2094
2095 void helper_fsin(void)
2096 {
2097     CPU86_LDouble fptemp;
2098
2099     fptemp = ST0;
2100     if ((fptemp > MAXTAN)||(fptemp < -MAXTAN)) {
2101         env->fpus |= 0x400;
2102     } else {
2103         ST0 = sin(fptemp);
2104         env->fpus &= (~0x400);  /* C2 <-- 0 */
2105         /* the above code is for  |arg| < 2**53 only */
2106     }
2107 }
2108
2109 void helper_fcos(void)
2110 {
2111     CPU86_LDouble fptemp;
2112
2113     fptemp = ST0;
2114     if((fptemp > MAXTAN)||(fptemp < -MAXTAN)) {
2115         env->fpus |= 0x400;
2116     } else {
2117         ST0 = cos(fptemp);
2118         env->fpus &= (~0x400);  /* C2 <-- 0 */
2119         /* the above code is for  |arg5 < 2**63 only */
2120     }
2121 }
2122
2123 /* associated heplers to reduce generated code length and to simplify
2124    relocation (FP constants are usually stored in .rodata section) */
2125
2126 void OPPROTO op_f2xm1(void)
2127 {
2128     helper_f2xm1();
2129 }
2130
2131 void OPPROTO op_fyl2x(void)
2132 {
2133     helper_fyl2x();
2134 }
2135
2136 void OPPROTO op_fptan(void)
2137 {
2138     helper_fptan();
2139 }
2140
2141 void OPPROTO op_fpatan(void)
2142 {
2143     helper_fpatan();
2144 }
2145
2146 void OPPROTO op_fxtract(void)
2147 {
2148     helper_fxtract();
2149 }
2150
2151 void OPPROTO op_fprem1(void)
2152 {
2153     helper_fprem1();
2154 }
2155
2156
2157 void OPPROTO op_fprem(void)
2158 {
2159     helper_fprem();
2160 }
2161
2162 void OPPROTO op_fyl2xp1(void)
2163 {
2164     helper_fyl2xp1();
2165 }
2166
2167 void OPPROTO op_fsqrt(void)
2168 {
2169     helper_fsqrt();
2170 }
2171
2172 void OPPROTO op_fsincos(void)
2173 {
2174     helper_fsincos();
2175 }
2176
2177 void OPPROTO op_frndint(void)
2178 {
2179     helper_frndint();
2180 }
2181
2182 void OPPROTO op_fscale(void)
2183 {
2184     helper_fscale();
2185 }
2186
2187 void OPPROTO op_fsin(void)
2188 {
2189     helper_fsin();
2190 }
2191
2192 void OPPROTO op_fcos(void)
2193 {
2194     helper_fcos();
2195 }
2196
2197 void OPPROTO op_fnstsw_A0(void)
2198 {
2199     int fpus;
2200     fpus = (env->fpus & ~0x3800) | (env->fpstt & 0x7) << 11;
2201     stw((void *)A0, fpus);
2202 }
2203
2204 void OPPROTO op_fnstsw_EAX(void)
2205 {
2206     int fpus;
2207     fpus = (env->fpus & ~0x3800) | (env->fpstt & 0x7) << 11;
2208     EAX = (EAX & 0xffff0000) | fpus;
2209 }
2210
2211 void OPPROTO op_fnstcw_A0(void)
2212 {
2213     stw((void *)A0, env->fpuc);
2214 }
2215
2216 void OPPROTO op_fldcw_A0(void)
2217 {
2218     int rnd_type;
2219     env->fpuc = lduw((void *)A0);
2220     /* set rounding mode */
2221     switch(env->fpuc & RC_MASK) {
2222     default:
2223     case RC_NEAR:
2224         rnd_type = FE_TONEAREST;
2225         break;
2226     case RC_DOWN:
2227         rnd_type = FE_DOWNWARD;
2228         break;
2229     case RC_UP:
2230         rnd_type = FE_UPWARD;
2231         break;
2232     case RC_CHOP:
2233         rnd_type = FE_TOWARDZERO;
2234         break;
2235     }
2236     fesetround(rnd_type);
2237 }
2238
2239 void OPPROTO op_fclex(void)
2240 {
2241     env->fpus &= 0x7f00;
2242 }
2243
2244 void OPPROTO op_fninit(void)
2245 {
2246     env->fpus = 0;
2247     env->fpstt = 0;
2248     env->fpuc = 0x37f;
2249     env->fptags[0] = 1;
2250     env->fptags[1] = 1;
2251     env->fptags[2] = 1;
2252     env->fptags[3] = 1;
2253     env->fptags[4] = 1;
2254     env->fptags[5] = 1;
2255     env->fptags[6] = 1;
2256     env->fptags[7] = 1;
2257 }
2258
2259 /* threading support */
2260 void OPPROTO op_lock(void)
2261 {
2262     cpu_lock();
2263 }
2264
2265 void OPPROTO op_unlock(void)
2266 {
2267     cpu_unlock();
2268 }