// 固定小数点とブレゼンハムの誤差をテストしてみる

#include <stdio.h>


#define LIMIT_SIZE		(1<<9)		// 10bit
#define POINT_SHIFT		20		// 小数点以下のbit数
#define POINT_SIZE		(1<<POINT_SHIFT)


void Bresen(int *buf, int Dx, int Dy);
void Fixpoint(int *buf, int Dx, int Dy);

void main(void)
 {
	int bresen[LIMIT_SIZE];
	int fixpoint[LIMIT_SIZE];

	for(int i = 0; i < LIMIT_SIZE; i ++)
	 {	// 傾きを変えてみる
		int x = LIMIT_SIZE;

		Bresen(bresen, x -1, i);

		Fixpoint(fixpoint, x -1, i);

		int dif = 0;
		for(int j = 0; j < x; j ++)
		 {
			if(bresen[j] != fixpoint[j])
			 {
				dif ++;
				printf("[%d] ", j);
			 }
		 }
		printf("%4d - %4d\n", i, dif);
	 }
 }


void Bresen(int *buf, int Dx, int Dy)
 {
	// xを基準にするのは固定
	int sx = 1;
	int y = 0;
	int E = -Dx;
	for(int i = 0; i <= Dx; i ++)
	 {
		buf[i] = y;
		E += 2*Dy;
		if(E >= 0)
		 {
			y ++;
			E -= 2*Dx;
		 }
	 }

 }

void Fixpoint(int *buf, int Dx, int Dy)
 {
	// 傾きを求める
	//Dx:Dy=POINT_SIZE:dy
	//dy = Dy*POINT_SIZE/Dx
	int dy = (Dy << POINT_SHIFT) / Dx;
	int y = 1 << (POINT_SHIFT-1);			// 0.5加える
	for(int i = 0; i <= Dx; i ++)
	 {
		buf[i] = y >> POINT_SHIFT;
		y += dy;
	 }
 }