normal.h

Go to the documentation of this file.

/*
 * Copyright (c) 2013 The University of Manchester
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 *     https://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
 
 
#ifndef __NORMAL_H__
#define __NORMAL_H__
 
#include <stdint.h>
#include <stdbool.h>
#include <stdfix.h>
#include "stdfix-full-iso.h"
#include "assert.h"
 
#ifndef UNIMPLEMENTED
#define UNIMPLEMENTED extern __attribute__((deprecated("Not implemented")))
#endif
 
int __norminv_rbits(int x);
 
static inline int __norminv_t_bits(int x) {
    return __norminv_rbits(x >> 16);
}
 
static inline int __norminv_b_bits(int x) {
    return __norminv_t_bits(x << 16);
}
 
static inline s1615 norminv_urb(unsigned int x) {
    union { unsigned int u; int s; } tmp;
 
    tmp.u = x;
    return kbits(__norminv_b_bits(tmp.s));
}
 
static inline s1615 norminv_urt(unsigned int x) {
    union { unsigned int u; int s; } tmp;
 
    tmp.u = x;
    return kbits(__norminv_t_bits(tmp.s));
}
 
static inline s1615 norminv_ur(u016 x) {
    union { unsigned int u; int s; } tmp;
 
    tmp.u = bitsur(x);
    return kbits(__norminv_rbits(tmp.s + INT16_MIN));
}
 
UNIMPLEMENTED s1615 norminv_ulr(u032 x);
 
#define norminv_fx(f) \
    _Generic((f), \
    u016:     norminv_ur(f), \
    u032:     norminv_ulr(f), \
    default:  __builtin_trap() \
    )
 
#endif /*__NORMAL_H__*/