以下对我来说工作得很好:
using Mat4 = Eigen::Matrix4f;
using Vec3 = Eigen::Vector3f;
class Camera {
public:
Vec3 pos = { 0.0f, 0.0f, 0.0f };
Vec3 forward = { 0.0f, 0.0f, 1.0f };
Vec3 up = { 0.0f, 1.0f, 0.0f };
float vfov = degrees(45.0f);
float near_plane = 0.1f;
float far_plane = 100.0f;
Mat4 cameraToScreen(float aspect_ratio) const {
Mat4 result;
float tanHalfFovy = tan(vfov *0.5f);
float x_scale = 1.0f / (aspect_ratio * tanHalfFovy);
float y_scale = 1.0f / (tanHalfFovy);
result <<
x_scale, 0, 0, 0,
0, y_scale, 0, 0,
0, 0, -(far_plane + near_plane) / (far_plane - near_plane), -1.0f,
0, 0, -2.0f * near_plane * far_plane / (far_plane - near_plane), 0;
return result;
}
Mat4 worldToCamera() const {
Mat4 result;
Vec3 p = pos;
Vec3 f = forward.normalized();
Vec3 r = f.cross(up).normalized();
Vec3 u = r.cross(f).normalized();
result <<
r.x(), u.x(), -f.x(), 0.0f,
r.y(), u.y(), -f.y(), 0.0f,
r.z(), u.z(), -f.z(), 0.0f,
-r.dot(p), -u.dot(p), f.dot(p), 1.0f;
return result;
}
Mat4 worldToScreen(float aspect_ratio) const {
return cameraToScreen(aspect_ratio) * worldToCamera();
}
};
但是,由于这将用于 OpenGL,并且我不喜欢处理转置,但我可以简单地切换到行主矩阵:
using Mat4 = Eigen::Matrix<float, 4, 4, Eigen::RowMajor>;
不幸的是,一切都破裂了。有趣的是,反转最终的矩阵乘法(这是错误的)解决了这个问题:
Mat4 worldToScreen(float aspect_ratio) const {
return worldToCamera() * cameraToScreen(aspect_ratio);
}
几乎就好像 eigen 在假设存储是列主要的情况下执行矩阵乘法,而不管我的模板标签是什么,但这不可能。
我分别测试了两个矩阵,并且在单独使用时它们在行优先模式下按我想要的方式执行,只有两个矩阵的乘积是错误的。
我对 Eigen 处理行主矩阵的方式有什么不理解?
编辑:
这是我将数据发送到 OpenGL 的有效方法。我的代码有很多额外的样板,但归结为:
struct CameraData {
Mat4 world_to_screen ;
};
void draw(Camera const& cam) {
CameraData cam_data;
// I remove the transposed() call when using Eigen::RowMajor
cam_data.world_to_screen = cam.worldToScreen(16.0f/9.0f).transposed();
// Camera data has a permanent ubo binding point shared by all programs.
glBindBuffer(GL_UNIFORM_BUFFER, cam_data_ubo);
glBufferSubData(GL_UNIFORM_BUFFER, 0, sizeof(cam_data), &cam_data);
//.. perform drawing
}
和 glsl:
#version 420
in vec3 vertex;
layout (std140) uniform CameraData {
mat4 world_to_screen;
};
void main() {
gl_Position = world_to_screen * vec4(vertex, 1.0);
}