/* 丸め方向による計算結果の違い（Windows 95/98用）
	Written by Takeshi OGITA.	2000/11/9
*/

#include <stdio.h>
#include <float.h>	/* ヘッダファイルの指定 */

#define Near _controlfp(_RC_NEAR, _MCW_RC)	/* 最近点への丸め */
#define Up   _controlfp(_RC_UP, _MCW_RC)	/* 上への丸め */
#define Down _controlfp(_RC_DOWN, _MCW_RC)	/* 下への丸め */
#define Chop _controlfp(_RC_CHOP, _MCW_RC)	/* 切り捨て */

#define ND 100

void 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);

	/* 切り捨て */
	Chop;
	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 (Chop) = %25.20f\n", prod);
}
/*
=====================================================================

＜実行結果＞

inner product (Near) =    0.49999999999999967000
inner product (Up)   =    0.50000000000000355000
inner product (Down) =    0.49999999999999645000
inner product (Chop) =    0.49999999999999994000

となり、この例では
	(Down) < (Near) < (Chop) < (Up)
であることが分かる。
言うまでもないが、この計算の厳密な値は 0.5 である。
*/