嗨,我正在尝试使用Yeppp提高代码中 vector 代数的性能!图书馆,但实际上实际上是越来越差...
这是Vector类代码的一部分:
#include "Vector3.h"
#include <cmath>
#include "yepCore.h"
Vector3::Vector3()
{
//ctor
}
Vector3::~Vector3()
{
//dtor
}
Vector3::Vector3(float X, float Y, float Z)
{
x = X;
y = Y;
z = Z;
}
float& Vector3::operator[](int idx)
{
return (&x)[idx];
}
Vector3& Vector3::normalize()
{
#if USE_YEPPP
float inf;
yepCore_SumSquares_V32f_S32f(&x, &inf, 3);
yepCore_Multiply_IV32fS32f_IV32f(&x, 1.0f / sqrt(inf), 3);
#else
float inf = 1.0f / sqrt((x * x) + (y * y) + (z * z));
x *= inf;
y *= inf;
z *= inf;
#endif
return *this;
}
Vector3 Vector3::cross(Vector3& rh)
{
return Vector3 (
(y * rh.z) - (z * rh.y),
(z * rh.x) - (x * rh.z),
(x * rh.y) - (y * rh.x)
);
}
float Vector3::dot(Vector3& rh)
{
#if USE_YEPPP
float ret = 0;
yepCore_DotProduct_V32fV32f_S32f(&x, &rh.x, &ret, 3);
return ret;
#else
return x*rh.x+y*rh.y+z*rh.z;
#endif
}
Vector3 Vector3::operator*(float scalar)
{
#if USE_YEPPP
Vector3 ret;
yepCore_Multiply_V32fS32f_V32f(&x, scalar, &ret.x , 3);
return ret;
#else
return Vector3(x*scalar, y*scalar,z*scalar);
#endif
}
Vector3 Vector3::operator+(Vector3 rh)
{
#if USE_YEPPP
Vector3 ret;
yepCore_Add_V32fV32f_V32f(&x, &rh.x, &ret.x, 3);
return ret;
#else
return Vector3(x+rh.x, y+rh.y, z+rh.z);
#endif
}
Vector3 Vector3::operator-(Vector3 rh)
{
#if USE_YEPPP
Vector3 ret;
yepCore_Subtract_V32fV32f_V32f(&x, &rh.x, &ret.x, 3);
return ret;
#else
return Vector3(x-rh.x, y-rh.y, z-rh.z);
#endif
}
Vector3 operator*(float s, const Vector3& v)
{
#if USE_YEPPP
Vector3 ret;
yepCore_Multiply_V32fS32f_V32f(&v.x, s, &ret.x , 3);
return ret;
#else
return Vector3(s*v.x,s*v.y,s*v.z);
#endif
}
我正在使用g++编译器。
编译器选项:g++ -Wall -fexceptions -fPIC -Wl,-不需要-std = c++ 11 -pthread -ggdb
链接器选项:g++ -shared -lpthread -lyeppp -ldl
所以知道我在做什么错吗?
最佳答案
是的!针对处理100多个元素的数组进行了优化。
由于使用SIMD的能力有限以及函数调用,动态分派(dispatch)和参数检查的开销,在小型数组(例如您的示例中为length-3数组)上,此方法效率不高。
关于c++ - Yeppp的表演!比本地实现慢,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/26504111/