gmqcc/intrin.cpp

#include <string.h>

#include "ast.h"
#include "fold.h"
#include "parser.h"

lex_ctx_t intrin::ctx() const {
    return parser_ctx(m_parser);
}

ast_function *intrin::value(ast_value **out, const char *name, qcint_t vtype) {
    ast_value *value = nullptr;
    ast_function *func  = nullptr;
    char buffer[1024];
    char stype [1024];

    util_snprintf(buffer, sizeof(buffer), "__builtin_%s", name);
    util_snprintf(stype,  sizeof(stype),   "<%s>",        type_name[vtype]);

    value = ast_value_new(ctx(), buffer, TYPE_FUNCTION);
    value->intrinsic = true;
    value->expression.next = (ast_expression*)ast_value_new(ctx(), stype, vtype);
    func = ast_function_new(ctx(), buffer, value);
    value->expression.flags |= AST_FLAG_ERASEABLE;

    *out = value;
    return func;
}

void intrin::reg(ast_value *const value, ast_function *const func) {
    m_parser->functions.push_back(func);
    m_parser->globals.push_back((ast_expression*)value);
}

#define QC_POW_EPSILON 0.00001f

ast_expression *intrin::nullfunc() {
    ast_value *val = nullptr;
    ast_function *func = value(&val, nullptr, TYPE_VOID);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::isfinite_() {
    /*
     * float isfinite(float x) {
     *     return !(isnan(x) || isinf(x));
     * }
     */
    ast_value    *val     = nullptr;
    ast_value    *x         = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_function *func      = value(&val, "isfinite", TYPE_FLOAT);
    ast_call     *callisnan = ast_call_new(ctx(), func_self("isnan", "isfinite"));
    ast_call     *callisinf = ast_call_new(ctx(), func_self("isinf", "isfinite"));
    ast_block    *block     = ast_block_new(ctx());

    /* float x; */
    val->expression.params.push_back(x);

    /* <callisnan> = isnan(x); */
    callisnan->params.push_back((ast_expression*)x);

    /* <callisinf> = isinf(x); */
    callisinf->params.push_back((ast_expression*)x);

    /* return (!<callisnan> || <callisinf>); */
    block->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_unary_new(
                ctx(),
                INSTR_NOT_F,
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_OR,
                    (ast_expression*)callisnan,
                    (ast_expression*)callisinf
                )
            )
        )
    );

    func->blocks.push_back(block);
    reg(val, func);

    return (ast_expression*)val;;
}

ast_expression *intrin::isinf_() {
    /*
     * float isinf(float x) {
     *     return (x != 0.0) && (x + x == x);
     * }
     */
    ast_value *val = nullptr;
    ast_value *x = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, "isinf", TYPE_FLOAT);

    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_AND,
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_NE_F,
                    (ast_expression*)x,
                    (ast_expression*)m_fold->m_imm_float[0]
                ),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_EQ_F,
                    (ast_expression*)ast_binary_new(
                        ctx(),
                        INSTR_ADD_F,
                        (ast_expression*)x,
                        (ast_expression*)x
                    ),
                    (ast_expression*)x
                )
            )
        )
    );

    val->expression.params.push_back(x);
    func->blocks.push_back(body);

    reg(val, func);

    return (ast_expression*)val;
}

ast_expression *intrin::isnan_() {
    /*
     * float isnan(float x) {
     *   float local;
     *   local = x;
     *
     *   return (x != local);
     * }
     */
    ast_value *val = nullptr;
    ast_value *arg1 = ast_value_new(ctx(), "x",TYPE_FLOAT);
    ast_value *local = ast_value_new(ctx(), "local", TYPE_FLOAT);
    ast_block *body  = ast_block_new(ctx());
    ast_function *func = value(&val, "isnan", TYPE_FLOAT);

    body->locals.push_back(local);
    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)local,
            (ast_expression*)arg1
        )
    );

    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_NE_F,
                (ast_expression*)arg1,
                (ast_expression*)local
            )
        )
    );

    val->expression.params.push_back(arg1);
    func->blocks.push_back(body);

    reg(val, func);

    return (ast_expression*)val;
}

ast_expression *intrin::isnormal_() {
    /*
     * float isnormal(float x) {
     *     return isfinite(x);
     * }
     */
    ast_value *val = nullptr;
    ast_call *callisfinite = ast_call_new(ctx(), func_self("isfinite", "isnormal"));
    ast_value *x = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, "isnormal", TYPE_FLOAT);

    val->expression.params.push_back(x);
    callisfinite->params.push_back((ast_expression*)x);

    /* return <callisfinite> */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)callisfinite
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::signbit_() {
    /*
     * float signbit(float x) {
     *     return (x < 0);
     * }
     */
    ast_value *val = nullptr;
    ast_value *x = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, "signbit", TYPE_FLOAT);

    val->expression.params.push_back(x);

    /* return (x < 0); */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_ternary_new(
                ctx(),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_LT,
                    (ast_expression*)x,
                    (ast_expression*)m_fold->m_imm_float[0]
                ),
                (ast_expression*)m_fold->m_imm_float[1],
                (ast_expression*)m_fold->m_imm_float[0]
            )
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::acosh_() {
    /*
     * float acosh(float x) {
     *     return log(x + sqrt((x * x) - 1));
     * }
     */
    ast_value    *val    = nullptr;
    ast_value    *x        = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_call     *calllog  = ast_call_new(ctx(), func_self("log", "acosh"));
    ast_call     *callsqrt = ast_call_new(ctx(), func_self("sqrt", "acosh"));
    ast_block    *body     = ast_block_new(ctx());
    ast_function *func     = value(&val, "acosh", TYPE_FLOAT);

    val->expression.params.push_back(x);

    /* <callsqrt> = sqrt((x * x) - 1); */
    callsqrt->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_SUB_F,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_MUL_F,
                (ast_expression*)x,
                (ast_expression*)x
            ),
            (ast_expression*)m_fold->m_imm_float[1]
        )
    );

    /* <calllog> = log(x + <callsqrt>); */
    calllog->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_ADD_F,
            (ast_expression*)x,
            (ast_expression*)callsqrt
        )
    );

    /* return <calllog>; */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)calllog
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::asinh_() {
    /*
     * float asinh(float x) {
     *     return log(x + sqrt((x * x) + 1));
     * }
     */
    ast_value *val = nullptr;
    ast_value *x = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_call *calllog = ast_call_new(ctx(), func_self("log", "asinh"));
    ast_call *callsqrt = ast_call_new(ctx(), func_self("sqrt", "asinh"));
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, "asinh", TYPE_FLOAT);

    val->expression.params.push_back(x);

    /* <callsqrt> = sqrt((x * x) + 1); */
    callsqrt->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_ADD_F,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_MUL_F,
                (ast_expression*)x,
                (ast_expression*)x
            ),
            (ast_expression*)m_fold->m_imm_float[1]
        )
    );

    /* <calllog> = log(x + <callsqrt>); */
    calllog->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_ADD_F,
            (ast_expression*)x,
            (ast_expression*)callsqrt
        )
    );

    /* return <calllog>; */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)calllog
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::atanh_() {
    /*
     * float atanh(float x) {
     *     return 0.5 * log((1 + x) / (1 - x))
     * }
     */
    ast_value    *val   = nullptr;
    ast_value    *x       = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_call     *calllog = ast_call_new(ctx(), func_self("log", "atanh"));
    ast_block    *body    = ast_block_new(ctx());
    ast_function *func    = value(&val, "atanh", TYPE_FLOAT);

    val->expression.params.push_back(x);

    /* <callog> = log((1 + x) / (1 - x)); */
    calllog->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_DIV_F,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_ADD_F,
                (ast_expression*)m_fold->m_imm_float[1],
                (ast_expression*)x
            ),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_SUB_F,
                (ast_expression*)m_fold->m_imm_float[1],
                (ast_expression*)x
            )
        )
    );

    /* return 0.5 * <calllog>; */
    body->exprs.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_MUL_F,
            (ast_expression*)m_fold->constgen_float(0.5, false),
            (ast_expression*)calllog
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::exp_() {
    /*
     * float exp(float x) {
     *     float sum = 1.0;
     *     float acc = 1.0;
     *     float i;
     *     for (i = 1; i < 200; ++i)
     *         sum += (acc *= x / i);
     *
     *     return sum;
     * }
     */
    ast_value *val = nullptr;
    ast_value *x = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_value *sum = ast_value_new(ctx(), "sum", TYPE_FLOAT);
    ast_value *acc = ast_value_new(ctx(), "acc", TYPE_FLOAT);
    ast_value *i = ast_value_new(ctx(), "i",   TYPE_FLOAT);
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, "exp", TYPE_FLOAT);

    val->expression.params.push_back(x);

    body->locals.push_back(sum);
    body->locals.push_back(acc);
    body->locals.push_back(i);

    /* sum = 1.0; */
    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)sum,
            (ast_expression*)m_fold->m_imm_float[1]
        )
    );

    /* acc = 1.0; */
    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)acc,
            (ast_expression*)m_fold->m_imm_float[1]
        )
    );

    /*
     * for (i = 1; i < 200; ++i)
     *     sum += (acc *= x / i);
     */
    body->exprs.push_back(
        (ast_expression*)ast_loop_new(
            ctx(),
            /* i = 1; */
            (ast_expression*)ast_store_new(
                ctx(),
                INSTR_STORE_F,
                (ast_expression*)i,
                (ast_expression*)m_fold->m_imm_float[1]
            ),
            /* i < 200; */
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_LT,
                (ast_expression*)i,
                (ast_expression*)m_fold->constgen_float(200.0f, false)
            ),
            false,
            nullptr,
            false,
            /* ++i; */
            (ast_expression*)ast_binstore_new(
                ctx(),
                INSTR_STORE_F,
                INSTR_ADD_F,
                (ast_expression*)i,
                (ast_expression*)m_fold->m_imm_float[1]
            ),
            /* sum += (acc *= (x / i)) */
            (ast_expression*)ast_binstore_new(
                ctx(),
                INSTR_STORE_F,
                INSTR_ADD_F,
                (ast_expression*)sum,
                (ast_expression*)ast_binstore_new(
                    ctx(),
                    INSTR_STORE_F,
                    INSTR_MUL_F,
                    (ast_expression*)acc,
                    (ast_expression*)ast_binary_new(
                        ctx(),
                        INSTR_DIV_F,
                        (ast_expression*)x,
                        (ast_expression*)i
                    )
                )
            )
        )
    );

    /* return sum; */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)sum
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::exp2_() {
    /*
     * float exp2(float x) {
     *     return pow(2, x);
     * }
     */
    ast_value *val = nullptr;
    ast_call  *callpow = ast_call_new(ctx(), func_self("pow", "exp2"));
    ast_value *arg1 = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, "exp2", TYPE_FLOAT);

    val->expression.params.push_back(arg1);

    callpow->params.push_back((ast_expression*)m_fold->m_imm_float[3]);
    callpow->params.push_back((ast_expression*)arg1);

    /* return <callpow> */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)callpow
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::expm1_() {
    /*
     * float expm1(float x) {
     *     return exp(x) - 1;
     * }
     */
    ast_value *val = nullptr;
    ast_call *callexp = ast_call_new(ctx(), func_self("exp", "expm1"));
    ast_value *x = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, "expm1", TYPE_FLOAT);

    val->expression.params.push_back(x);

    /* <callexp> = exp(x); */
    callexp->params.push_back((ast_expression*)x);

    /* return <callexp> - 1; */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_SUB_F,
                (ast_expression*)callexp,
                (ast_expression*)m_fold->m_imm_float[1]
            )
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::pow_() {
    /*
     *
     * float pow(float base, float exp) {
     *     float result;
     *     float low;
     *     float high;
     *     float mid;
     *     float square;
     *     float accumulate;
     *
     *     if (exp == 0.0)
     *         return 1;
     *     if (exp == 1.0)
     *         return base;
     *     if (exp < 0)
     *         return 1.0 / pow(base, -exp);
     *     if (exp >= 1) {
     *         result = pow(base, exp / 2);
     *         return result * result;
     *     }
     *
     *     low        = 0.0f;
     *     high       = 1.0f;
     *     square     = sqrt(base);
     *     accumulate = square;
     *     mid        = high / 2.0f
     *
     *     while (fabs(mid - exp) > QC_POW_EPSILON) {
     *         square = sqrt(square);
     *         if (mid < exp) {
     *             low         = mid;
     *             accumulate *= square;
     *         } else {
     *             high        = mid;
     *             accumulate *= (1.0f / square);
     *         }
     *         mid = (low + high) / 2;
     *     }
     *     return accumulate;
     * }
     */
    ast_value    *val = nullptr;
    ast_function *func = value(&val, "pow", TYPE_FLOAT);

    /* prepare some calls for later */
    ast_call *callpow1  = ast_call_new(ctx(), (ast_expression*)val);                  /* for pow(base, -exp)    */
    ast_call *callpow2  = ast_call_new(ctx(), (ast_expression*)val);                  /* for pow(vase, exp / 2) */
    ast_call *callsqrt1 = ast_call_new(ctx(), func_self("sqrt", "pow")); /* for sqrt(base)         */
    ast_call *callsqrt2 = ast_call_new(ctx(), func_self("sqrt", "pow")); /* for sqrt(square)       */
    ast_call *callfabs  = ast_call_new(ctx(), func_self("fabs", "pow")); /* for fabs(mid - exp)    */

    /* prepare some blocks for later */
    ast_block *expgt1       = ast_block_new(ctx());
    ast_block *midltexp     = ast_block_new(ctx());
    ast_block *midltexpelse = ast_block_new(ctx());
    ast_block *whileblock   = ast_block_new(ctx());

    /* float pow(float base, float exp) */
    ast_value    *base = ast_value_new(ctx(), "base", TYPE_FLOAT);
    ast_value    *exp  = ast_value_new(ctx(), "exp",  TYPE_FLOAT);
    /* { */
    ast_block    *body = ast_block_new(ctx());

    /*
     * float result;
     * float low;
     * float high;
     * float square;
     * float accumulate;
     * float mid;
     */
    ast_value *result     = ast_value_new(ctx(), "result",     TYPE_FLOAT);
    ast_value *low        = ast_value_new(ctx(), "low",        TYPE_FLOAT);
    ast_value *high       = ast_value_new(ctx(), "high",       TYPE_FLOAT);
    ast_value *square     = ast_value_new(ctx(), "square",     TYPE_FLOAT);
    ast_value *accumulate = ast_value_new(ctx(), "accumulate", TYPE_FLOAT);
    ast_value *mid        = ast_value_new(ctx(), "mid",        TYPE_FLOAT);
    body->locals.push_back(result);
    body->locals.push_back(low);
    body->locals.push_back(high);
    body->locals.push_back(square);
    body->locals.push_back(accumulate);
    body->locals.push_back(mid);

    val->expression.params.push_back(base);
    val->expression.params.push_back(exp);

    /*
     * if (exp == 0.0)
     *     return 1;
     */
    body->exprs.push_back(
        (ast_expression*)ast_ifthen_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_EQ_F,
                (ast_expression*)exp,
                (ast_expression*)m_fold->m_imm_float[0]
            ),
            (ast_expression*)ast_return_new(
                ctx(),
                (ast_expression*)m_fold->m_imm_float[1]
            ),
            nullptr
        )
    );

    /*
     * if (exp == 1.0)
     *     return base;
     */
    body->exprs.push_back(
        (ast_expression*)ast_ifthen_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_EQ_F,
                (ast_expression*)exp,
                (ast_expression*)m_fold->m_imm_float[1]
            ),
            (ast_expression*)ast_return_new(
                ctx(),
                (ast_expression*)base
            ),
            nullptr
        )
    );

    /* <callpow1> = pow(base, -exp) */
    callpow1->params.push_back((ast_expression*)base);
    callpow1->params.push_back(
        (ast_expression*)ast_unary_new(
            ctx(),
            VINSTR_NEG_F,
            (ast_expression*)exp
        )
    );

    /*
     * if (exp < 0)
     *     return 1.0 / <callpow1>;
     */
    body->exprs.push_back(
        (ast_expression*)ast_ifthen_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_LT,
                (ast_expression*)exp,
                (ast_expression*)m_fold->m_imm_float[0]
            ),
            (ast_expression*)ast_return_new(
                ctx(),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_DIV_F,
                    (ast_expression*)m_fold->m_imm_float[1],
                    (ast_expression*)callpow1
                )
            ),
            nullptr
        )
    );

    /* <callpow2> = pow(base, exp / 2) */
    callpow2->params.push_back((ast_expression*)base);
    callpow2->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_DIV_F,
            (ast_expression*)exp,
            (ast_expression*)m_fold->m_imm_float[3] /* 2.0f */
        )
    );

    /*
     * <expgt1> = {
     *     result = <callpow2>;
     *     return result * result;
     * }
     */
    expgt1->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)result,
            (ast_expression*)callpow2
        )
    );
    expgt1->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_MUL_F,
                (ast_expression*)result,
                (ast_expression*)result
            )
        )
    );

    /*
     * if (exp >= 1) {
     *     <expgt1>
     * }
     */
    body->exprs.push_back(
        (ast_expression*)ast_ifthen_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_GE,
                (ast_expression*)exp,
                (ast_expression*)m_fold->m_imm_float[1]
            ),
            (ast_expression*)expgt1,
            nullptr
        )
    );

    /*
     * <callsqrt1> = sqrt(base)
     */
    callsqrt1->params.push_back((ast_expression*)base);

    /*
     * low        = 0.0f;
     * high       = 1.0f;
     * square     = sqrt(base);
     * accumulate = square;
     * mid        = high / 2.0f;
     */
    body->exprs.push_back(
        (ast_expression*)ast_store_new(ctx(),
            INSTR_STORE_F,
            (ast_expression*)low,
            (ast_expression*)m_fold->m_imm_float[0]
        )
    );
    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)high,
            (ast_expression*)m_fold->m_imm_float[1]
        )
    );

    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)square,
            (ast_expression*)callsqrt1
        )
    );

    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)accumulate,
            (ast_expression*)square
        )
    );
    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)mid,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_DIV_F,
                (ast_expression*)high,
                (ast_expression*)m_fold->m_imm_float[3] /* 2.0f */
            )
        )
    );

    /*
     * <midltexp> = {
     *     low         = mid;
     *     accumulate *= square;
     * }
     */
    midltexp->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)low,
            (ast_expression*)mid
        )
    );
    midltexp->exprs.push_back(
        (ast_expression*)ast_binstore_new(
            ctx(),
            INSTR_STORE_F,
            INSTR_MUL_F,
            (ast_expression*)accumulate,
            (ast_expression*)square
        )
    );

    /*
     * <midltexpelse> = {
     *     high        = mid;
     *     accumulate *= (1.0 / square);
     * }
     */
    midltexpelse->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)high,
            (ast_expression*)mid
        )
    );
    midltexpelse->exprs.push_back(
        (ast_expression*)ast_binstore_new(
            ctx(),
            INSTR_STORE_F,
            INSTR_MUL_F,
            (ast_expression*)accumulate,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_DIV_F,
                (ast_expression*)m_fold->m_imm_float[1],
                (ast_expression*)square
            )
        )
    );

    /*
     * <callsqrt2> = sqrt(square)
     */
    callsqrt2->params.push_back((ast_expression*)square);

    /*
     * <whileblock> = {
     *     square = <callsqrt2>;
     *     if (mid < exp)
     *          <midltexp>;
     *     else
     *          <midltexpelse>;
     *
     *     mid = (low + high) / 2;
     * }
     */
    whileblock->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)square,
            (ast_expression*)callsqrt2
        )
    );
    whileblock->exprs.push_back(
        (ast_expression*)ast_ifthen_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_LT,
                (ast_expression*)mid,
                (ast_expression*)exp
            ),
            (ast_expression*)midltexp,
            (ast_expression*)midltexpelse
        )
    );
    whileblock->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)mid,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_DIV_F,
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_ADD_F,
                    (ast_expression*)low,
                    (ast_expression*)high
                ),
                (ast_expression*)m_fold->m_imm_float[3] /* 2.0f */
            )
        )
    );

    /*
     * <callabs> = fabs(mid - exp)
     */
    callfabs->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_SUB_F,
            (ast_expression*)mid,
            (ast_expression*)exp
        )
    );

    /*
     * while (<callfabs>  > epsilon)
     *     <whileblock>
     */
    body->exprs.push_back(
        (ast_expression*)ast_loop_new(
            ctx(),
            /* init */
            nullptr,
            /* pre condition */
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_GT,
                (ast_expression*)callfabs,
                (ast_expression*)m_fold->constgen_float(QC_POW_EPSILON, false)
            ),
            /* pre not */
            false,
            /* post condition */
            nullptr,
            /* post not */
            false,
            /* increment expression */
            nullptr,
            /* code block */
            (ast_expression*)whileblock
        )
    );

    /* return accumulate */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)accumulate
        )
    );

    /* } */
    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::mod_() {
    /*
     * float mod(float a, float b) {
     *     float div = a / b;
     *     float sign = (div < 0.0f) ? -1 : 1;
     *     return a - b * sign * floor(sign * div);
     * }
     */
    ast_value    *val = nullptr;
    ast_call     *call  = ast_call_new(ctx(), func_self("floor", "mod"));
    ast_value    *a     = ast_value_new(ctx(), "a",    TYPE_FLOAT);
    ast_value    *b     = ast_value_new(ctx(), "b",    TYPE_FLOAT);
    ast_value    *div   = ast_value_new(ctx(), "div",  TYPE_FLOAT);
    ast_value    *sign  = ast_value_new(ctx(), "sign", TYPE_FLOAT);
    ast_block    *body  = ast_block_new(ctx());
    ast_function *func  = value(&val, "mod", TYPE_FLOAT);

    val->expression.params.push_back(a);
    val->expression.params.push_back(b);

    body->locals.push_back(div);
    body->locals.push_back(sign);

    /* div = a / b; */
    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)div,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_DIV_F,
                (ast_expression*)a,
                (ast_expression*)b
            )
        )
    );

    /* sign = (div < 0.0f) ? -1 : 1; */
    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)sign,
            (ast_expression*)ast_ternary_new(
                ctx(),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_LT,
                    (ast_expression*)div,
                    (ast_expression*)m_fold->m_imm_float[0]
                ),
                (ast_expression*)m_fold->m_imm_float[2],
                (ast_expression*)m_fold->m_imm_float[1]
            )
        )
    );

    /* floor(sign * div) */
    call->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            INSTR_MUL_F,
            (ast_expression*)sign,
            (ast_expression*)div
        )
    );

    /* return a - b * sign * <call> */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_SUB_F,
                (ast_expression*)a,
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_MUL_F,
                    (ast_expression*)b,
                    (ast_expression*)ast_binary_new(
                        ctx(),
                        INSTR_MUL_F,
                        (ast_expression*)sign,
                        (ast_expression*)call
                    )
                )
            )
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::fabs_() {
    /*
     * float fabs(float x) {
     *     return x < 0 ? -x : x;
     * }
     */
    ast_value    *val  = nullptr;
    ast_value    *arg1   = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_block    *body   = ast_block_new(ctx());
    ast_function *func   = value(&val, "fabs", TYPE_FLOAT);

    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_ternary_new(
                ctx(),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_LE,
                    (ast_expression*)arg1,
                    (ast_expression*)m_fold->m_imm_float[0]
                ),
                (ast_expression*)ast_unary_new(
                    ctx(),
                    VINSTR_NEG_F,
                    (ast_expression*)arg1
                ),
                (ast_expression*)arg1
            )
        )
    );

    val->expression.params.push_back(arg1);

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::epsilon_() {
    /*
     * float epsilon(void) {
     *     float eps = 1.0f;
     *     do { eps /= 2.0f; } while ((1.0f + (eps / 2.0f)) != 1.0f);
     *     return eps;
     * }
     */
    ast_value    *val  = nullptr;
    ast_value    *eps    = ast_value_new(ctx(), "eps", TYPE_FLOAT);
    ast_block    *body   = ast_block_new(ctx());
    ast_function *func   = value(&val, "epsilon", TYPE_FLOAT);

    body->locals.push_back(eps);

    /* eps = 1.0f; */
    body->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)eps,
            (ast_expression*)m_fold->m_imm_float[0]
        )
    );

    body->exprs.push_back(
        (ast_expression*)ast_loop_new(
            ctx(),
            nullptr,
            nullptr,
            false,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_NE_F,
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_ADD_F,
                    (ast_expression*)m_fold->m_imm_float[1],
                    (ast_expression*)ast_binary_new(
                        ctx(),
                        INSTR_MUL_F,
                        (ast_expression*)eps,
                        (ast_expression*)m_fold->m_imm_float[3] /* 2.0f */
                    )
                ),
                (ast_expression*)m_fold->m_imm_float[1]
            ),
            false,
            nullptr,
            (ast_expression*)ast_binstore_new(
                ctx(),
                INSTR_STORE_F,
                INSTR_DIV_F,
                (ast_expression*)eps,
                (ast_expression*)m_fold->m_imm_float[3] /* 2.0f */
            )
        )
    );

    /* return eps; */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)eps
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::nan_() {
    /*
     * float nan(void) {
     *     float x = 0.0f;
     *     return x / x;
     * }
     */
    ast_value    *val  = nullptr;
    ast_value    *x      = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_function *func   = value(&val, "nan", TYPE_FLOAT);
    ast_block    *block  = ast_block_new(ctx());

    block->locals.push_back(x);

    block->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)x,
            (ast_expression*)m_fold->m_imm_float[0]
        )
    );

    block->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_DIV_F,
                (ast_expression*)x,
                (ast_expression*)x
            )
        )
    );

    func->blocks.push_back(block);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::inf_() {
    /*
     * float inf(void) {
     *     float x = 1.0f;
     *     float y = 0.0f;
     *     return x / y;
     * }
     */
    ast_value    *val  = nullptr;
    ast_value    *x      = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_value    *y      = ast_value_new(ctx(), "y", TYPE_FLOAT);
    ast_function *func   = value(&val, "inf", TYPE_FLOAT);
    ast_block    *block  = ast_block_new(ctx());
    size_t        i;

    block->locals.push_back(x);
    block->locals.push_back(y);

    /* to keep code size down */
    for (i = 0; i <= 1; i++) {
        block->exprs.push_back(
            (ast_expression*)ast_store_new(
                ctx(),
                INSTR_STORE_F,
                (ast_expression*)((i == 0) ? x : y),
                (ast_expression*)m_fold->m_imm_float[i]
            )
        );
    }

    block->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_DIV_F,
                (ast_expression*)x,
                (ast_expression*)y
            )
        )
    );

    func->blocks.push_back(block);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::ln_() {
    /*
     * float log(float power, float base) {
     *   float whole;
     *   float nth
     *   float sign = 1.0f;
     *   float eps  = epsilon();
     *
     *   if (power <= 1.0f || bbase <= 1.0) {
     *       if (power <= 0.0f || base <= 0.0f)
     *           return nan();
     *
     *       if (power < 1.0f) {
     *           power = 1.0f / power;
     *           sign *= -1.0f;
     *       }
     *
     *       if (base < 1.0f) {
     *           sign *= -1.0f;
     *           base  = 1.0f / base;
     *       }
     *   }
     *
     *   float A_i       = 1;
     *   float B_i       = 0;
     *   float A_iminus1 = 0;
     *   float B_iminus1 = 1;
     *
     *   for (;;) {
     *       whole = power;
     *       nth   = 0.0f;
     *
     *       while (whole >= base) {
     *           float base2    = base;
     *           float n2       = 1.0f;
     *           float newbase2 = base2 * base2;
     *
     *           while (whole >= newbase2) {
     *               base2     = newbase2;
     *               n2       *= 2;
     *               newbase2 *= newbase2;
     *           }
     *
     *           whole /= base2;
     *           nth += n2;
     *       }
     *
     *       float b_iplus1 = n;
     *       float A_iplus1 = b_iplus1 * A_i + A_iminus1;
     *       float B_iplus1 = b_iplus1 * B_i + B_iminus1;
     *
     *       A_iminus1 = A_i;
     *       B_iminus1 = B_i;
     *       A_i       = A_iplus1;
     *       B_i       = B_iplus1;
     *
     *       if (whole <= 1.0f + eps)
     *           break;
     *
     *       power = base;
     *       bower = whole;
     *   }
     *   return sign * A_i / B_i;
     * }
     */

    ast_value *val = nullptr;
    ast_value *power = ast_value_new(ctx(), "power", TYPE_FLOAT);
    ast_value *base = ast_value_new(ctx(), "base",TYPE_FLOAT);
    ast_value *whole= ast_value_new(ctx(), "whole", TYPE_FLOAT);
    ast_value *nth = ast_value_new(ctx(), "nth", TYPE_FLOAT);
    ast_value *sign = ast_value_new(ctx(), "sign", TYPE_FLOAT);
    ast_value *A_i = ast_value_new(ctx(), "A_i", TYPE_FLOAT);
    ast_value *B_i = ast_value_new(ctx(), "B_i", TYPE_FLOAT);
    ast_value *A_iminus1 = ast_value_new(ctx(), "A_iminus1", TYPE_FLOAT);
    ast_value *B_iminus1 = ast_value_new(ctx(), "B_iminus1", TYPE_FLOAT);
    ast_value *b_iplus1 = ast_value_new(ctx(), "b_iplus1", TYPE_FLOAT);
    ast_value *A_iplus1 = ast_value_new(ctx(), "A_iplus1", TYPE_FLOAT);
    ast_value *B_iplus1 = ast_value_new(ctx(), "B_iplus1", TYPE_FLOAT);
    ast_value *eps = ast_value_new(ctx(), "eps", TYPE_FLOAT);
    ast_value *base2 = ast_value_new(ctx(), "base2", TYPE_FLOAT);
    ast_value *n2 = ast_value_new(ctx(), "n2",TYPE_FLOAT);
    ast_value *newbase2 = ast_value_new(ctx(), "newbase2", TYPE_FLOAT);
    ast_block *block = ast_block_new(ctx());
    ast_block *plt1orblt1 = ast_block_new(ctx()); // (power <= 1.0f || base <= 1.0f)
    ast_block *plt1 = ast_block_new(ctx()); // (power < 1.0f)
    ast_block *blt1 = ast_block_new(ctx()); // (base< 1.0f)
    ast_block *forloop = ast_block_new(ctx()); // for(;;)
    ast_block *whileloop = ast_block_new(ctx()); // while (whole >= base)
    ast_block *nestwhile= ast_block_new(ctx()); // while (whole >= newbase2)
    ast_function *func = value(&val, "ln", TYPE_FLOAT);
    size_t i;

    val->expression.params.push_back(power);
    val->expression.params.push_back(base);

    block->locals.push_back(whole);
    block->locals.push_back(nth);
    block->locals.push_back(sign);
    block->locals.push_back(eps);
    block->locals.push_back(A_i);
    block->locals.push_back(B_i);
    block->locals.push_back(A_iminus1);
    block->locals.push_back(B_iminus1);

    /* sign = 1.0f; */
    block->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)sign,
            (ast_expression*)m_fold->m_imm_float[1]
        )
    );

    /* eps = __builtin_epsilon(); */
    block->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)eps,
            (ast_expression*)ast_call_new(
                ctx(),
                func_self("__builtin_epsilon", "ln")
            )
        )
    );

    /*
     * A_i       = 1;
     * B_i       = 0;
     * A_iminus1 = 0;
     * B_iminus1 = 1;
     */
    for (i = 0; i <= 1; i++) {
        int j;
        for (j = 1; j >= 0; j--) {
            block->exprs.push_back(
                (ast_expression*)ast_store_new(
                    ctx(),
                    INSTR_STORE_F,
                    (ast_expression*)((j) ? ((i) ? B_iminus1 : A_i)
                                          : ((i) ? A_iminus1 : B_i)),
                    (ast_expression*)m_fold->m_imm_float[j]
                )
            );
        }
    }

    /*
     * <plt1> = {
     *     power = 1.0f / power;
     *     sign *= -1.0f;
     * }
     * <blt1> = {
     *     base  = 1.0f / base;
     *     sign *= -1.0f;
     * }
     */
    for (i = 0; i <= 1; i++) {
        ((i) ? blt1 : plt1)->exprs.push_back(
            (ast_expression*)ast_store_new(
                ctx(),
                INSTR_STORE_F,
                (ast_expression*)((i) ? base : power),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_DIV_F,
                    (ast_expression*)m_fold->m_imm_float[1],
                    (ast_expression*)((i) ? base : power)
                )
            )
        );
        plt1->exprs.push_back(
            (ast_expression*)ast_binstore_new(
                ctx(),
                INSTR_STORE_F,
                INSTR_MUL_F,
                (ast_expression*)sign,
                (ast_expression*)m_fold->m_imm_float[2]
            )
        );
    }

    /*
     * <plt1orblt1> = {
     *     if (power <= 0.0 || base <= 0.0f)
     *         return __builtin_nan();
     *     if (power < 1.0f)
     *         <plt1>
     *     if (base < 1.0f)
     *         <blt1>
     * }
     */
    plt1orblt1->exprs.push_back(
        (ast_expression*)ast_ifthen_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_OR,
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_LE,
                    (ast_expression*)power,
                    (ast_expression*)m_fold->m_imm_float[0]
                ),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_LE,
                    (ast_expression*)base,
                    (ast_expression*)m_fold->m_imm_float[0]
                )
            ),
            (ast_expression*)ast_return_new(
                ctx(),
                (ast_expression*)ast_call_new(
                    ctx(),
                    func_self("__builtin_nan", "ln")
                )
            ),
            nullptr
        )
    );

    for (i = 0; i <= 1; i++) {
        plt1orblt1->exprs.push_back(
            (ast_expression*)ast_ifthen_new(
                ctx(),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_LT,
                    (ast_expression*)((i) ? base : power),
                    (ast_expression*)m_fold->m_imm_float[1]
                ),
                (ast_expression*)((i) ? blt1 : plt1),
                nullptr
            )
        );
    }

    block->exprs.push_back((ast_expression*)plt1orblt1);


    /* whole = power; */
    forloop->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)whole,
            (ast_expression*)power
        )
    );

    /* nth = 0.0f; */
    forloop->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)nth,
            (ast_expression*)m_fold->m_imm_float[0]
        )
    );

    /* base2 = base; */
    whileloop->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)base2,
            (ast_expression*)base
        )
    );

    /* n2 = 1.0f; */
    whileloop->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)n2,
            (ast_expression*)m_fold->m_imm_float[1]
        )
    );

    /* newbase2 = base2 * base2; */
    whileloop->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)newbase2,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_MUL_F,
                (ast_expression*)base2,
                (ast_expression*)base2
            )
        )
    );

    /* while loop locals */
    whileloop->locals.push_back(base2);
    whileloop->locals.push_back(n2);
    whileloop->locals.push_back(newbase2);

    /* base2 = newbase2; */
    nestwhile->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)base2,
            (ast_expression*)newbase2
        )
    );

    /* n2 *= 2; */
    nestwhile->exprs.push_back(
        (ast_expression*)ast_binstore_new(
            ctx(),
            INSTR_STORE_F,
            INSTR_MUL_F,
            (ast_expression*)n2,
            (ast_expression*)m_fold->m_imm_float[3] /* 2.0f */
        )
    );

    /* newbase2 *= newbase2; */
    nestwhile->exprs.push_back(
        (ast_expression*)ast_binstore_new(
            ctx(),
            INSTR_STORE_F,
            INSTR_MUL_F,
            (ast_expression*)newbase2,
            (ast_expression*)newbase2
        )
    );

    /* while (whole >= newbase2) */
    whileloop->exprs.push_back(
        (ast_expression*)ast_loop_new(
            ctx(),
            nullptr,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_GE,
                (ast_expression*)whole,
                (ast_expression*)newbase2
            ),
            false,
            nullptr,
            false,
            nullptr,
            (ast_expression*)nestwhile
        )
    );

    /* whole /= base2; */
    whileloop->exprs.push_back(
        (ast_expression*)ast_binstore_new(
            ctx(),
            INSTR_STORE_F,
            INSTR_DIV_F,
            (ast_expression*)whole,
            (ast_expression*)base2
        )
    );

    /* nth += n2; */
    whileloop->exprs.push_back(
        (ast_expression*)ast_binstore_new(
            ctx(),
            INSTR_STORE_F,
            INSTR_ADD_F,
            (ast_expression*)nth,
            (ast_expression*)n2
        )
    );

    /* while (whole >= base) */
    forloop->exprs.push_back(
        (ast_expression*)ast_loop_new(
            ctx(),
            nullptr,
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_GE,
                (ast_expression*)whole,
                (ast_expression*)base
            ),
            false,
            nullptr,
            false,
            nullptr,
            (ast_expression*)whileloop
        )
    );

    forloop->locals.push_back(b_iplus1);
    forloop->locals.push_back(A_iplus1);
    forloop->locals.push_back(B_iplus1);

    /* b_iplus1 = nth; */
    forloop->exprs.push_back(
        (ast_expression*)ast_store_new(
            ctx(),
            INSTR_STORE_F,
            (ast_expression*)b_iplus1,
            (ast_expression*)nth
        )
    );

    /*
     * A_iplus1 = b_iplus1 * A_i + A_iminus1;
     * B_iplus1 = b_iplus1 * B_i + B_iminus1;
     */
    for (i = 0; i <= 1; i++) {
        forloop->exprs.push_back(
            (ast_expression*)ast_store_new(
                ctx(),
                INSTR_STORE_F,
                (ast_expression*)((i) ? B_iplus1 : A_iplus1),
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_ADD_F,
                    (ast_expression*)ast_binary_new(
                        ctx(),
                        INSTR_MUL_F,
                        (ast_expression*)b_iplus1,
                        (ast_expression*) ((i) ? B_i : A_i)
                    ),
                    (ast_expression*)((i) ? B_iminus1 : A_iminus1)
                )
            )
        );
    }

    /*
     * A_iminus1 = A_i;
     * B_iminus1 = B_i;
     */
    for (i = 0; i <= 1; i++) {
        forloop->exprs.push_back(
            (ast_expression*)ast_store_new(
                ctx(),
                INSTR_STORE_F,
                (ast_expression*)((i) ? B_iminus1 : A_iminus1),
                (ast_expression*)((i) ? B_i       : A_i)
            )
        );
    }

    /*
     * A_i = A_iplus1;
     * B_i = B_iplus1;
     */
    for (i = 0; i <= 1; i++) {
        forloop->exprs.push_back(
            (ast_expression*)ast_store_new(
                ctx(),
                INSTR_STORE_F,
                (ast_expression*)((i) ? B_i      : A_i),
                (ast_expression*)((i) ? B_iplus1 : A_iplus1)
            )
        );
    }

    /*
     * if (whole <= 1.0f + eps)
     *     break;
     */
    forloop->exprs.push_back(
        (ast_expression*)ast_ifthen_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_LE,
                (ast_expression*)whole,
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_ADD_F,
                    (ast_expression*)m_fold->m_imm_float[1],
                    (ast_expression*)eps
                )
            ),
            (ast_expression*)ast_breakcont_new(
                ctx(),
                false,
                0
            ),
            nullptr
        )
    );

    /*
     * power = base;
     * base  = whole;
     */
    for (i = 0; i <= 1; i++) {
        forloop->exprs.push_back(
            (ast_expression*)ast_store_new(
                ctx(),
                INSTR_STORE_F,
                (ast_expression*)((i) ? base  : power),
                (ast_expression*)((i) ? whole : base)
            )
        );
    }

    /* add the for loop block */
    block->exprs.push_back(
        (ast_expression*)ast_loop_new(
            ctx(),
            nullptr,
            /* for(; 1; ) ?? (can this be nullptr too?) */
            (ast_expression*)m_fold->m_imm_float[1],
            false,
            nullptr,
            false,
            nullptr,
            (ast_expression*)forloop
        )
    );

    /* return sign * A_i / B_il */
    block->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)ast_binary_new(
                ctx(),
                INSTR_MUL_F,
                (ast_expression*)sign,
                (ast_expression*)ast_binary_new(
                    ctx(),
                    INSTR_DIV_F,
                    (ast_expression*)A_i,
                    (ast_expression*)B_i
                )
            )
        )
    );

    func->blocks.push_back(block);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::log_variant(const char *name, float base) {
    ast_value *val = nullptr;
    ast_call *callln = ast_call_new (ctx(), func_self("__builtin_ln", name));
    ast_value *arg1 = ast_value_new(ctx(), "x", TYPE_FLOAT);
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, name, TYPE_FLOAT);

    val->expression.params.push_back(arg1);

    callln->params.push_back((ast_expression*)arg1);
    callln->params.push_back((ast_expression*)m_fold->constgen_float(base, false));

    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)callln
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::log_() {
    return log_variant("log", 2.7182818284590452354);
}
ast_expression *intrin::log10_() {
    return log_variant("log10", 10);
}
ast_expression *intrin::log2_() {
    return log_variant("log2", 2);
}
ast_expression *intrin::logb_() {
    /* FLT_RADIX == 2 for now */
    return log_variant("log2", 2);
}

ast_expression *intrin::shift_variant(const char *name, size_t instr) {
    /*
     * float [shift] (float a, float b) {
     *   return floor(a [instr] pow(2, b));
     */
    ast_value *val = nullptr;
    ast_call *callpow = ast_call_new(ctx(), func_self("pow", name));
    ast_call *callfloor = ast_call_new(ctx(), func_self("floor", name));
    ast_value *a = ast_value_new(ctx(), "a", TYPE_FLOAT);
    ast_value *b = ast_value_new(ctx(), "b", TYPE_FLOAT);
    ast_block *body = ast_block_new(ctx());
    ast_function *func = value(&val, name, TYPE_FLOAT);

    val->expression.params.push_back(a);
    val->expression.params.push_back(b);

    /* <callpow> = pow(2, b) */
    callpow->params.push_back((ast_expression*)m_fold->m_imm_float[3]);
    callpow->params.push_back((ast_expression*)b);

    /* <callfloor> = floor(a [instr] <callpow>) */
    callfloor->params.push_back(
        (ast_expression*)ast_binary_new(
            ctx(),
            instr,
            (ast_expression*)a,
            (ast_expression*)callpow
        )
    );

    /* return <callfloor> */
    body->exprs.push_back(
        (ast_expression*)ast_return_new(
            ctx(),
            (ast_expression*)callfloor
        )
    );

    func->blocks.push_back(body);
    reg(val, func);
    return (ast_expression*)val;
}

ast_expression *intrin::lshift() {
    return shift_variant("lshift", INSTR_MUL_F);
}

ast_expression *intrin::rshift() {
    return shift_variant("rshift", INSTR_DIV_F);
}

void intrin::error(const char *fmt, ...) {
    va_list ap;
    va_start(ap, fmt);
    vcompile_error(ctx(), fmt, ap);
    va_end(ap);
}

/* exposed */
ast_expression *intrin::debug_typestring() {
    return (ast_expression*)0x1;
}

intrin::intrin(parser_t *parser)
    : m_parser(parser)
    , m_fold(&parser->m_fold)
{
    static const intrin_func_t intrinsics[] = {
        {&intrin::isfinite_,        "__builtin_isfinite",         "isfinite", 1},
        {&intrin::isinf_,           "__builtin_isinf",            "isinf",    1},
        {&intrin::isnan_,           "__builtin_isnan",            "isnan",    1},
        {&intrin::isnormal_,        "__builtin_isnormal",         "isnormal", 1},
        {&intrin::signbit_,         "__builtin_signbit",          "signbit",  1},
        {&intrin::acosh_,           "__builtin_acosh",            "acosh",    1},
        {&intrin::asinh_,           "__builtin_asinh",            "asinh",    1},
        {&intrin::atanh_,           "__builtin_atanh",            "atanh",    1},
        {&intrin::exp_,             "__builtin_exp",              "exp",      1},
        {&intrin::exp2_,            "__builtin_exp2",             "exp2",     1},
        {&intrin::expm1_,           "__builtin_expm1",            "expm1",    1},
        {&intrin::mod_,             "__builtin_mod",              "mod",      2},
        {&intrin::pow_,             "__builtin_pow",              "pow",      2},
        {&intrin::fabs_,            "__builtin_fabs",             "fabs",     1},
        {&intrin::log_,             "__builtin_log",              "log",      1},
        {&intrin::log10_,           "__builtin_log10",            "log10",    1},
        {&intrin::log2_,            "__builtin_log2",             "log2",     1},
        {&intrin::logb_,            "__builtin_logb",             "logb",     1},
        {&intrin::lshift,           "__builtin_lshift",           "",         2},
        {&intrin::rshift,           "__builtin_rshift",           "",         2},
        {&intrin::epsilon_,         "__builtin_epsilon",          "",         0},
        {&intrin::nan_,             "__builtin_nan",              "",         0},
        {&intrin::inf_,             "__builtin_inf",              "",         0},
        {&intrin::ln_,              "__builtin_ln",               "",         2},
        {&intrin::debug_typestring, "__builtin_debug_typestring", "",         0},
        {&intrin::nullfunc,         "#nullfunc",                  "",         0}
    };

    for (auto &it : intrinsics) {
        m_intrinsics.push_back(it);
        m_generated.push_back(nullptr);
    }
}

ast_expression *intrin::do_fold(ast_value *val, ast_expression **exprs) {
    if (!val || !val->name)
        return nullptr;
    static constexpr size_t kPrefixLength = 10; // "__builtin_"
    for (auto &it : m_intrinsics) {
        if (!strcmp(val->name, it.name))
            return (vec_size(exprs) != it.args)
                        ? nullptr
                        : m_fold->intrinsic(val->name + kPrefixLength, exprs);
    }
    return nullptr;
}

ast_expression *intrin::func_try(size_t offset, const char *compare) {
    for (auto &it : m_intrinsics) {
        const size_t index = &it - &m_intrinsics[0];
        if (strcmp(*(char **)((char *)&it + offset), compare))
            continue;
        if (m_generated[index])
            return m_generated[index];
        return m_generated[index] = (this->*it.intrin_func_t::function)();
    }
    return nullptr;
}

ast_expression *intrin::func_self(const char *name, const char *from) {
    ast_expression *find;
    /* try current first */
    if ((find = parser_find_global(m_parser, name)) && ((ast_value*)find)->expression.vtype == TYPE_FUNCTION)
        for (auto &it : m_parser->functions)
            if (((ast_value*)find)->name && !strcmp(it->name, ((ast_value*)find)->name) && it->builtin < 0)
                return find;
    /* try name second */
    if ((find = func_try(offsetof(intrin_func_t, name),  name)))
        return find;
    /* try alias third */
    if ((find = func_try(offsetof(intrin_func_t, alias), name)))
        return find;

    if (from) {
        error("need function `%s', compiler depends on it for `__builtin_%s'", name, from);
        return func_self("#nullfunc", nullptr);
    }
    return nullptr;
}

ast_expression *intrin::func(const char *name) {
    return func_self(name, nullptr);
}