#define _GNU_SOURCE
#include <stdint.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include <signal.h>
#include <sys/types.h>
#include <unistd.h>

#include "debug.h"
#include "box64stack.h"
#include "x64emu.h"
#include "x64emu_private.h"
#include "x64run_private.h"
#include "x64primop.h"
#include "x64trace.h"
#include "x87emu_private.h"
#include "x87emu_setround.h"
#include "box64context.h"
#include "bridge.h"

#include "modrm.h"

#ifdef TEST_INTERPRETER
uintptr_t TestDC(x64test_t *test, rex_t rex, uintptr_t addr)
#else
uintptr_t RunDC(x64emu_t *emu, rex_t rex, uintptr_t addr)
#endif
{
    uint8_t nextop;
    reg64_t *oped;
    #ifdef TEST_INTERPRETER
    x64emu_t*emu = test->emu;
    #endif

    int oldround = fpu_setround(emu);
    nextop = F8;
    if(MODREG)
    switch(nextop) {
        case 0xC0:
        case 0xC1:
        case 0xC2:
        case 0xC3:
        case 0xC4:
        case 0xC5:
        case 0xC6:
        case 0xC7:  /* FADD */
            ST(nextop&7).d += ST0.d;
            if(!emu->cw.f.C87_PC) ST(nextop&7).d = (float)ST(nextop&7).d;
            break;
        case 0xC8:
        case 0xC9:
        case 0xCA:
        case 0xCB:
        case 0xCC:
        case 0xCD:
        case 0xCE:
        case 0xCF:  /* FMUL */
            ST(nextop&7).d *= ST0.d;
            if(!emu->cw.f.C87_PC) ST(nextop&7).d = (float)ST(nextop&7).d;
            break;
        case 0xD0:
        case 0xD1:
        case 0xD2:
        case 0xD3:
        case 0xD4:
        case 0xD5:
        case 0xD6:
        case 0xD7:  /* FCOM */
            fpu_fcom(emu, ST(nextop&7).d);
            break;
        case 0xD8:
        case 0xD9:
        case 0xDA:
        case 0xDB:
        case 0xDC:
        case 0xDD:
        case 0xDE:
        case 0xDF:  /* FCOMP */
            fpu_fcom(emu, ST(nextop&7).d);
            fpu_do_pop(emu);
            break;
        case 0xE0:
        case 0xE1:
        case 0xE2:
        case 0xE3:
        case 0xE4:
        case 0xE5:
        case 0xE6:
        case 0xE7:  /* FSUBR */
            ST(nextop&7).d = ST0.d -ST(nextop&7).d;
            if(!emu->cw.f.C87_PC) ST(nextop&7).d = (float)ST(nextop&7).d;
            break;
        case 0xE8:
        case 0xE9:
        case 0xEA:
        case 0xEB:
        case 0xEC:
        case 0xED:
        case 0xEE:
        case 0xEF:  /* FSUB */
            ST(nextop&7).d -= ST0.d;
            if(!emu->cw.f.C87_PC) ST(nextop&7).d = (float)ST(nextop&7).d;
            break;
        case 0xF0:
        case 0xF1:
        case 0xF2:
        case 0xF3:
        case 0xF4:
        case 0xF5:
        case 0xF6:
        case 0xF7:  /* FDIVR */
            ST(nextop&7).d = ST0.d / ST(nextop&7).d;
            if(!emu->cw.f.C87_PC) ST(nextop&7).d = (float)ST(nextop&7).d;
            break;
        case 0xF8:
        case 0xF9:
        case 0xFA:
        case 0xFB:
        case 0xFC:
        case 0xFD:
        case 0xFE:
        case 0xFF:  /* FDIV */
            ST(nextop&7).d /=  ST0.d;
            if(!emu->cw.f.C87_PC) ST(nextop&7).d = (float)ST(nextop&7).d;
            break;
        default:
            fesetround(oldround);
            return 0;
    } else {
            GETE8(0);
            switch((nextop>>3)&7) {
            case 0:         /* FADD ST0, double */
                ST0.d += *(double*)ED;
                if(!emu->cw.f.C87_PC) ST0.d = (float)ST0.d;
                break;
            case 1:         /* FMUL ST0, double */
                ST0.d *= *(double*)ED;
                if(!emu->cw.f.C87_PC) ST0.d = (float)ST0.d;
                break;
            case 2:      /* FCOM ST0, double */
                fpu_fcom(emu, *(double*)ED);
                break;
            case 3:     /* FCOMP ST0, double */
                fpu_fcom(emu, *(double*)ED);
                fpu_do_pop(emu);
                break;
            case 4:         /* FSUB ST0, double */
                ST0.d -= *(double*)ED;
                if(!emu->cw.f.C87_PC) ST0.d = (float)ST0.d;
                break;
            case 5:         /* FSUBR ST0, double */
                ST0.d = *(double*)ED - ST0.d;
                if(!emu->cw.f.C87_PC) ST0.d = (float)ST0.d;
                break;
            case 6:         /* FDIV ST0, double */
                ST0.d /= *(double*)ED;
                if(!emu->cw.f.C87_PC) ST0.d = (float)ST0.d;
                break;
            case 7:         /* FDIVR ST0, double */
                ST0.d = *(double*)ED / ST0.d;
                if(!emu->cw.f.C87_PC) ST0.d = (float)ST0.d;
                break;
           default:
                fesetround(oldround);
                return 0;
        }
    }
    fesetround(oldround);
    return addr;
}