File: rounding_math_simple.cpp

package info (click to toggle)
cgal 6.1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 144,912 kB
  • sloc: cpp: 810,858; ansic: 208,477; sh: 493; python: 411; makefile: 286; javascript: 174
file content (29 lines) | stat: -rw-r--r-- 791 bytes parent folder | download | duplicates (7)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
#include <fenv.h>
#include <stdio.h>

double stop_constant_propagation(double x)
{
    // asm volatile ("" : "+mx"(x)); return x;
    asm volatile ("" : "=m"(x) : "m"(x)); return x;
    // volatile double y = x; return y;   
}

double force_to_double(double x)
{
#if defined __i386__ || defined __ia64__
    // asm volatile ("" : "+mx"(x)); return x;
    asm volatile ("" : "=m"(x) : "m"(x)); return x;
    // volatile double y = x; return y;   
#else
    return x;
#endif       
}

int main()
{
    fesetround(FE_UPWARD);
    double lhs = -force_to_double((-1.1) * stop_constant_propagation(10.1));
    double rhs =  force_to_double(( 1.1) * stop_constant_propagation(10.1));
    printf( "%s\n", lhs != rhs ? "lhs != rhs (correct)" : "lhs == rhs (error)");
    return lhs != rhs ? 0 : 1;
}