It's maybe stupid question, but I can't find any function in C++ for square and higher orders, like x^2, x^3 ....
It's maybe stupid question, but I can't find any function in C++ for square and higher orders, like x^2, x^3 ....
/*whew*/ ... man, I've spent hours perfecting this code for you!
Code:template <typename number> number squared( number n ) { return (n*n); } template <typename number> number cubed( number n ) { return (n*n*n); }
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
For higher orders, use:
Code:template <typename number> number exponential( number n, int order ) { if( order == 0 ) { return 0; } number value = n; while(--order) { value *= n; } return value; }
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
Thanx! But that the easiest way??? It seems to me very complicated for simple function like that...
That was what I was looking for! It works with float's too ;o)
This should be pretty fast for positive integer exponents.
You can make it even faster by taking out the recursion.
It takes advantage of the fact that (a ^ (n/2))^2 = a^n.
Code:int exp(int a, int n) { int r; if (n == 0) { r = 1; } else { if (n % 2 == 0) { r = exp(a, n/2); r *= r; } else { r = exp(a, n-1) * a; } } return r; }
Hmm, nice code, Nick. And I wasn't aware of that fact either. I spotted a place where you could speed it up even more, though, by shifting n once to the right instead of dividing by two.
Code:int exp(int a, int n) { int r; if (n == 0) { r = 1; } else { if (n % 2 == 0) { r = exp(a, (n >> 1) ); r *= r; } else { r = exp(a, n-1) * a; } } return r; }
And I think making this solution non-recursive would slow down the algorithm due to the additional manipulations, not to mention making it far more unreadable (of course, if you were using it for graphics, readability would be a non-issue...).
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
use pow() -- it will be faster and more versatile than anything you can come up with
hello, internet!
You could also
take out n % 2 == 0 with !(n & 0x00000001). I think
most compilers would see these optimizations if I had used
unsigned int. Another way to speed it up would be
to store maybe the first 200 or so in a table.
Ha! We're on a roll now! Ok Nick, here's the challenge:You could also take out n % 2 == 0 with !(n & 0x00000001).
Let's time the following:
number exponential( number n, int order );
int exp(int a, int n); //..v1.0, v2.0, v.3.0, and the table lookup version.
Sound like fun?
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
The template won!
I didn't use the table lookup, but you can easily add it:
printout:
Looping Template: 8 seconds to run 10, 000, 000 iterations...
Recursive #1: 12 seconds to run 10, 000, 000 iterations...
Recursive #2: 9 seconds to run 10, 000, 000 iterations...
Recursive #3: 10 seconds to run 10, 000, 000 iterations...
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
I think for mine to be faster than yours I think n would
have to be pretty large. In which case
the results would overflow an integer.
You could try this version that uses long double.
Code:#include <iostream> using namespace std; long double exp(long double a, unsigned int n); int main(void) { long double a = 1.5; int n = 300; cout << a << "^" << n << " = " << exp(a, n) << endl; return 0; } long double exp(long double a, unsigned int n) { unsigned int mask = 0xffffffff; unsigned int t = n; long double r = a; if (n == 0) return 1.0; // this finds mask such // that if n = 0x00001010 // mask = 0x00000100 // there's probably a faster way to do // it but it shouldn't matter much if n is large while(t != 0) { mask <<= 1; t &= mask; } mask = ~mask; mask ^= (mask >> 1); mask >>= 1; while(mask != 0) { r *= r; if (n & mask) { r *= a; } mask >>= 1; } return r; }
These results seem very strange...Originally posted by Sebastiani
Looping Template: 8 seconds to run 10, 000, 000 iterations...
Recursive #1: 12 seconds to run 10, 000, 000 iterations...
Recursive #2: 9 seconds to run 10, 000, 000 iterations...
Recursive #3: 10 seconds to run 10, 000, 000 iterations...
At first glance, I thought the more advanced function would win.
I decided to test it for myself
I got the result:
Looping Template: 0 seconds to run 10,000,000 iterations...
Recursive #1: 0 seconds to run 10,000,000 iterations...
Recursive #2: 0 seconds to run 10,000,000 iterations...
Recursive #3: 0 seconds to run 10,000,000 iterations...
So I had to increase the iterations:
Moral: Keep the code simple!Looping Template: 6 seconds to run 100,000,000 iterations...
Recursive #1: 8 seconds to run 100,000,000 iterations...
Recursive #2: 8 seconds to run 100,000,000 iterations...
Recursive #3: 7 seconds to run 100,000,000 iterations...
Last edited by Sang-drax : Tomorrow at 02:21 AM. Reason: Time travelling
My computer may be old, but it's reliable, damn it!So I had to increase the iterations:
What are you running, anyway - 2 GHZ??
But yeah, I was surprised by the results too, thinking the last recursive would be faster. I think it's wierd that the ratios were so different on our machines, too. Strange...
Code:#include <cmath> #include <complex> bool euler_flip(bool value) { return std::pow ( std::complex<float>(std::exp(1.0)), std::complex<float>(0, 1) * std::complex<float>(std::atan(1.0) *(1 << (value + 2))) ).real() < 0; }
I change the test so that it
works for larger n. This is what
I get
1.5^300 = 6.72013e+052
Looping Template: 46.06 seconds to run 10000000 iterations...
1.5^300 = 6.72013e+052
Recursive #1: 8.111 seconds to run 10000000 iterations...
1.5^300 = 6.72013e+052
Recursive #2: 8.033 seconds to run 10000000 iterations...
1.5^300 = 6.72013e+052
Recursive #3: 11.265 seconds to run 10000000 iterations...
1.5^300 = 6.72013e+052
Iterative #4: 4.414 seconds to run 10000000 iterations...
And the winner is Iterative #4 coming in at 4.414 seconds.