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

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

#define Near fpsetround(FP_RN)	/* 最近点への丸め */
#define Up   fpsetround(FP_RP)	/* 上への丸め */
#define Down fpsetround(FP_RM)	/* 下への丸め */

#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);
}
/*
=====================================================================

＜実行結果＞

未確認
*/