/*
 * Copyright 2005, Andrew Bachmann, andrewbachmann@myrealbox.com
 * Distributed under the terms of the MIT License.
 */


#include <OS.h>
#include <image.h>
#include <cmath>
#include <cassert>
#include <cstdio>
#include <cstring>

static image_id
get_libroot_id()
{
	image_info info;
	int32 cookie = 0;
	while (get_next_image_info(0, &cookie, &info) == B_OK) {
		if (strcmp(info.name, "/boot/beos/system/lib/libroot.so") == 0) {
			return info.id;
		}
	}
	return B_BAD_VALUE;
}


static double (*be_sin)(double);
static double (*be_cos)(double);
static double (*be_copysign)(double, double);
static double (*be_drem)(double, double);


int 
main(int argc, char **argv)
{
	image_id libroot = get_libroot_id();
	assert(get_image_symbol(libroot, "sin", B_SYMBOL_TYPE_TEXT, (void**)&be_sin) == B_OK);
	assert(get_image_symbol(libroot, "cos", B_SYMBOL_TYPE_TEXT, (void**)&be_cos) == B_OK);
	assert(get_image_symbol(libroot, "copysign", B_SYMBOL_TYPE_TEXT, (void**)&be_copysign) == B_OK);
	assert(get_image_symbol(libroot, "drem", B_SYMBOL_TYPE_TEXT, (void**)&be_drem) == B_OK);

	status_t result = B_OK;

	// test cos2 + sin2 = 1
	fprintf(stdout, "value\tsin2()+cos2()\n");
	for (int i = -63 ; i <= 63 ; i++) {
		double f = (double)i/7.0;
		double x = sin(f);
		double y = cos(f);
		if (fabs(x*x + y*y - 1.0) > 0.000001) {
			fprintf(stdout, "%0.3f\t%0.10f", f, x*x + y*y);
			fprintf(stdout, " **");
			fprintf(stdout, "\n");
			result = B_ERROR;
		}
	}
	fprintf(stdout, "\n");

	// test sin
	fprintf(stdout, "value\tsin(value)\tbe_sin(value)\n");
	for (int i = -63 ; i <= 63 ; i++) {
		double f = (double)i/7.0;
		double x = sin(f);
		double y = be_sin(f);
		if (fabs(x - y) > 0.000001) {
			fprintf(stdout, "%0.3f\t%0.10f\t%0.10f", f, x, y);
			fprintf(stdout, " **");
			fprintf(stdout, "\n");
			result = B_ERROR;
		}
	}
	fprintf(stdout, "\n");

	fprintf(stdout, "value\tcos(value)\tbe_cos(value)\n");
	// test cos
	for (int i = -63 ; i <= 63 ; i++) {
		double f = (double)i/7.0;
		double x = cos(f);
		double y = be_cos(f);
		if (fabs(x - y) > 0.000001) {
			fprintf(stdout, "%0.3f\t%0.10f\t%0.10f", f, x, y);
			fprintf(stdout, " **");
			fprintf(stdout, "\n");
			result = B_ERROR;
		}
	}
	fprintf(stdout, "\n");

	fprintf(stdout, "arg1\targ2\tcopysign()\tbe_copysign()\n");
	// test copysign
	for (int i = -21 ; i <= 21 ; i++) {
		for (int j = -36 ; j <= 36 ; j++) {
			double f = (double)i/3.0;
			double g = (double)j/4.0;
			double x = copysign(f, g);
			double y = be_copysign(f, g);
			if ((x != y) && !(std::isnan(x) && std::isnan(y))) {
				fprintf(stdout, "%0.1f\t%0.1f\t%0.6f\t%0.6f", f, g, x, y);
				fprintf(stdout, " **");
				fprintf(stdout, "\n");
				result = B_ERROR;
			}
		}
	}
	fprintf(stdout, "\n");

	fprintf(stdout, "arg1\targ2\tdrem(values)\tbe_drem(values)\n");
	// test drem
	for (int i = -21 ; i <= 21 ; i++) {
		for (int j = -36 ; j <= 36 ; j++) {
			double f = (double)i/3.0;
			double g = (double)j/4.0;
			double x = drem(f, g);
			double y = be_drem(f, g);
			if ((x != y) && !(std::isnan(x) && std::isnan(y))) {
				fprintf(stdout, "%0.2f\t%0.2f\t%0.10f\t%0.10f", f, g, x, y);
				fprintf(stdout, " **");
				fprintf(stdout, "\n");
				result = B_ERROR;
			}
		}
	}
	fprintf(stdout, "\n");

	fprintf(stdout, "result = %s\n", (result == B_OK ? "OK" : "ERROR"));
	return result;
}