/* 丸め方向による計算結果の違い（Linux / FreeBSD 3.x用）
	Written by Takeshi OGITA.	2000/11/9
*/

#include <stdio.h>

#define Near asm volatile ("fldcw _RoundNear")	/* 最近点への丸め */
#define Up   asm volatile ("fldcw _RoundUp")	/* 上への丸め */
#define Down asm volatile ("fldcw _RoundDown")	/* 下への丸め */

int _RoundNear = 0x133a;
int _RoundUp   = 0x1b3a;
int _RoundDown = 0x173a;

#define ND 100

int main()
{
	int i, ii, n;
	double prod;
	double A[ND], B[ND], X[ND], Y[ND];

	n = ND;

	ii = 1;
	for(i=0; i<n; i++){
		ii *= -1;
		X[i] = 1.0*i;
		Y[i] = 1.0*ii;
	}

	/* 最近点への丸め */
	Near;
	for(i=0; i<n; i++){
		A[i] = X[i]/10.0;
		B[i] = Y[i]/10.0;
	}
	prod = 0.0;
	for(i=0; i<n; i++) prod += A[i]*B[i];
	printf("inner product (Near) = %25.20f\n", prod);

	/* 上への丸め */
	Up;
	for(i=0; i<n; i++){
		A[i] = X[i]/10.0;
		B[i] = Y[i]/10.0;
	}
	prod = 0.0;
	for(i=0; i<n; i++) prod += A[i]*B[i];
	printf("inner product (Up)   = %25.20f\n", prod);

	/* 下への丸め */
	Down;
	for(i=0; i<n; i++){
		A[i] = X[i]/10.0;
		B[i] = Y[i]/10.0;
	}
	prod = 0.0;
	for(i=0; i<n; i++) prod += A[i]*B[i];
	printf("inner product (Down) = %25.20f\n", prod);

	return 0;
}
/*
=====================================================================

＜実行結果＞

inner product (Near) =    0.50000000000000000000
inner product (Up)   =    0.50000000000000344169
inner product (Down) =    0.49999999999999650280
*/