首先,对具体情况表示抱歉。我通常会把我的SO问题归结为一般的"A类"问题,只包含相关的内容,但我不确定问题的根源是什么。
我有一个矩阵类模板,看起来像这样(只显示我认为相关的部分):
template <std::size_t R, std::size_t C>
class Matrix
{
private:
//const int rows, cols;
std::array<std::array<float,C>,R> m;
public:
inline std::array<float,C>& operator[](const int i)
{
return m[i];
}
const std::array<float,C> operator[](const int i) const
{
return m[i];
}
template<std::size_t N>
Matrix<R,N> operator *(const Matrix<C,N> a) const
{
Matrix<R,N> result = Matrix<R,N>();
// irrelevant calculation
return result;
}
// ... other very similar stuff, I'm not sure that it's relevant
}
template <std::size_t S>
Matrix<S,S> identity()
{
Matrix<S,S> matrix = Matrix<S,S>();
for(std::size_t x = 0; x < S; x++)
{
for(std::size_t y = 0; y < S; y++)
{
if (x == y)
{
matrix[x][y] = 1.f;
}
}
}
return matrix;
}
我对全班进行了单元测试,乘法和恒等式工厂似乎都工作得很好。然而,我在这个方法中使用了它,它被称为很多次(我认为如果你曾经写过一个渲染器,很明显我在这里要做的是什么):
Vec3i Renderer::world_to_screen_space(Vec3f v)
{
Matrix<4,1> vm = v2m(v);
Matrix<4,4> projection = identity<4>(); // If I change this to Matrix<4,4>(), the error doesn't happen
projection[3][2] = -1.f;
vm = projection * vm;
Vec3f r = m2v(vm);
return Vec3i(
(r.x + 1.) * (width / 2.),
(r.y + 1.) * (height / 2.),
r.z
);
}
经过一段时间和对这个方法的一些随机调用,我得到了这个:
Job 1, 'and ./bin/main' terminated by signal SIGBUS (Misaligned address error)
但是,如果我将行identity<4>()
更改为Matrix<4,4>()
,则不会发生错误。我是C++的新手,所以这一定很愚蠢。
那么,(1)这个错误意味着什么?(2)我是如何射中自己的腿的?
更新:当然,这个错误不会在LLDB调试器中重现。
更新2:以下是我通过Valgrind运行程序后得到的结果:
==66525== Invalid read of size 4
==66525== at 0x1000148D5: Renderer::draw_triangle(Vec3<float>, Vec3<float>, Vec3<float>, Vec2<int>, Vec2<int>, Vec2<int>, Model, float) (in ./bin/main)
而draw_triangle
正是调用world_to_screen_space
并使用其结果的方法。
Update3:我发现了问题的根源,它与此代码没有任何关系,而且它也是非常明显的。现在真的不知道该怎么办。
如果没有一个检查错位的处理器(正如@twalberg所说),就不可能运行和验证代码。但我可以这么说:在C++或其他库中,混淆一种类型的异常和另一种类型异常是一个常见的错误。
我的猜测——对不起,我不能做更多的事情——是您正在创建丢失的分配,耗尽了可用内存,然后溢出了内存空间。当您超过可用内存时抛出的非常罕见的异常可能是意外的,并作为未对齐错误返回。试着在运行时检查内存使用情况,以确定是否存在这种情况。
编辑:
我的猜测是错误的,valgrind输出显示地址错误是正确的。跑步是个好主意。明确的迹象是,有一个bug的级别比你的代码低,所以我最初的想法几乎可以肯定是正确的:有一个bug不在你的代码中,但被掩盖了。
请注意,identity()构造函数和Matrix之间似乎存在差异<,>构造函数,因为前者是沿着对角线初始化的(慢慢地:最好是消除内部循环),而后者不是。这可能会影响v2m和m2v的行为。