quakeforge/tools/qfcc/test/triangle.r

void printf (string fmt, ...) = #0;
float (string s) stof = #0;
float (float x) sqrt = #0;

float
heron (float a, float b, float c)
{
	float s = (a + b + c) / 2;
	return sqrt (s*(s-a)*(s-b)*(s-c));
}

float
kahan (float a, float b, float c)
{
	float t;
	if (a < b) { t = a; a = b; b = t; }
	if (b < c) { t = b; b = c; c = t; }
	if (a < b) { t = a; a = b; b = t; }
	printf ("%f %f %f\n", a, b, c);
	return sqrt ((a+(b+c))*(c-(a-b))*(c+(a-b))*(a+(b-c)))/4;
}

int
main (int argc, string *argv)
{
	float a = stof(argv[1]);
	float b = stof(argv[2]);
	float c = stof(argv[3]);
	printf ("%.9g %.9g %.9g\n", a, b, c);
	printf ("%.9g %.9g\n", heron (a, b, c), kahan (a, b, c));
	return 0;
}