Navigation

Changeset 93

Timestamp:

12/22/05 10:25:03 (3 years ago)

Author:

Don Clugston

Message:

How embaressing! The calcpi demo was printing PI=3.22... due to an over-zealous simplification. Now it works again.

Files:

r76	r93
40	40	template atanTerm(real x, int n)
41	41	{
42		const real atanTerm = -pow!(-x, 2n-1)/(2n-1);
	42	static if (n&1) {
	43	const real atanTerm = pow!(x, 2n-1)/(2n-1);
	44	} else {
	45	const real atanTerm = -pow!(x, 2n-1)/(2n-1);
	46	}
43	47	}
44	48
45	49	/// Machin's formula for pi
46	50	/// pi/4 = 4 atan(1/5) - atan(1/239).
47		pragma(msg, "PI = " ~ fcvt!(4.0 * (4*~~atan!(1/5.0) -~~ atan!(1/239.0))) );
	51	pragma(msg, "PI = " ~ fcvt!(4.0 * (4*.atan!(1/5.0) - .atan!(1/239.0))) );
48	52

r83	r93
5	5	module meta.math;
6	6
7		/** bool isnan!(real x)
8		*
9		* Return true if and only if x is an IEEE not-a-number
10		*/
	7	/*********************************************
	8	* bool isnan!(real x)
	9	*
	10	* Return true if and only if x is an IEEE not-a-number
	11	*/
11	12	template isnan(real x)
12	13	{
…	…
28	29	}
29	30
30		/** int binaryExponent!(real x)
	31	/*********************************************
	32	* int binaryExponent!(real x)
31	33	* Returns the binary exponent of a real number x.
32	34	*
…	…
45	47	}
46	48
47		/** int decimalExponent!(real x)
	49	/*********************************************
	50	* int decimalExponent!(real x)
48	51	* Returns the decimal exponent of a real number x.
49	52	*
…	…
74	77	}
75	78
76		/** real pow!(real a, int b)
	79	/*********************************************
	80	* real pow!(real a, int b)
77	81	* Fast integer powers
78	82	*/
…	…
87	91
88	92
89		/** creal powz!(creal a, int b)
	93	/*********************************************
	94	* creal powz!(creal a, int b)
90	95	* Fast integer powers of a complex number
91	96	*/