diff options
author | Gael Guennebaud <g.gael@free.fr> | 2010-07-13 13:29:35 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2010-07-13 13:29:35 +0200 |
commit | e4f3759c4daca39545cb5efb26e08c94efbd9ee0 (patch) | |
tree | f034b557366050a5a9ac8a294067d6a807c2f134 /bench/quatmul.cpp | |
parent | 6ba5d2c90c83ba169e51c3b2936636e0a05c3c76 (diff) |
add a bench for quaternion multiplication
Diffstat (limited to 'bench/quatmul.cpp')
-rw-r--r-- | bench/quatmul.cpp | 47 |
1 files changed, 47 insertions, 0 deletions
diff --git a/bench/quatmul.cpp b/bench/quatmul.cpp new file mode 100644 index 000000000..d91a4b01b --- /dev/null +++ b/bench/quatmul.cpp @@ -0,0 +1,47 @@ +#include <iostream> +#include <Eigen/Core> +#include <Eigen/Geometry> +#include <bench/BenchTimer.h> + +using namespace Eigen; + +template<typename Quat> +EIGEN_DONT_INLINE void quatmul_default(const Quat& a, const Quat& b, Quat& c) +{ + c = a * b; +} + +template<typename Quat> +EIGEN_DONT_INLINE void quatmul_novec(const Quat& a, const Quat& b, Quat& c) +{ + c = ei_quat_product<0, Quat, Quat, typename Quat::Scalar, Aligned>::run(a,b); +} + +template<typename Quat> void bench(const std::string& label) +{ + int tries = 10; + int rep = 1000000; + BenchTimer t; + + Quat a(4, 1, 2, 3); + Quat b(2, 3, 4, 5); + Quat c; + + std::cout.precision(3); + + BENCH(t, tries, rep, quatmul_default(a,b,c)); + std::cout << label << " default " << 1e3*t.best(CPU_TIMER) << "ms \t" << 1e-6*double(rep)/(t.best(CPU_TIMER)) << " M mul/s\n"; + + BENCH(t, tries, rep, quatmul_novec(a,b,c)); + std::cout << label << " novec " << 1e3*t.best(CPU_TIMER) << "ms \t" << 1e-6*double(rep)/(t.best(CPU_TIMER)) << " M mul/s\n"; +} + +int main() +{ + bench<Quaternionf>("float "); + bench<Quaterniond>("double"); + + return 0; + +} + |