#include <stdio.h>
main() {
union int32_or_single {
int i;
float f;
} a, b, r, rs;
union int64_or_double {
long long int i; // Please specify 64-bit integer type for platform
double f;
} rd;
// Uncomment to assign values in decimal
//a.f = 15161.0;
//b.f = 1077.00048828125;
// Assign values in hexadecimal
a.i=0x466CE400;
b.i=0x4486A004;
//Do the deed (in this example subtraction)
r.f=a.f - b.f;
//Repeat, but this time do in double precision to avoid rounding
rd.f=(double) a.f - (double) b.f;
xil_printf("%c[2J",27);
xil_printf("a: Hex=%08X \r\n ", a.i);
xil_printf("Float=%.11f \r\n", a.f);
xil_printf("b: Hex=%08X \r\n", b.i);
xil_printf("Float=%.11f \r\n", b.f);
xil_printf("Single result: Hex=%08X \r\n", r.i);
xil_printf("Float=%.11f \r\n", r.f);
xil_printf("Double result: Hex=%08X%08X \r\n",(int)(rd.i>>32),(int)(rd.i&0xFFFFFFFF));
xil_printf("Float=%.11f \r\n", rd.f);
rs.f = (float) rd.f; // Round result from double to single
xil_printf("Double rounded to single: Hex=%08X \r\n", rs.i);
xil_printf("Float=%.11f \r\n", rs.f);
// Expect the following output:
// a: Hex=466CE400 Float=15161.00000000000
// b: Hex=4486A004 Float=1077.00048828125
// Single precision result: Hex=465C1000 Float=14084.00000000000
// Double precision result: Hex=40CB81FFF0000000 Float=14083.99951171875
// Double rounded to single: Hex=465C1000 Float=14084.00000000000
}