assimp中如何判断矩阵是否是单位矩阵

对于一个矩阵元素为浮点型的矩阵,你是否还在使每个元素跟1.0f或0.0f进行比较,如果这样,只能说你的结果不一定正确,那我们看看assimp中是如何做的。

cpp 复制代码
template <typename TReal>
AI_FORCE_INLINE
bool aiMatrix4x4t<TReal>::IsIdentity() const {
    // Use a small epsilon to solve floating-point inaccuracies
    const static TReal epsilon = 10e-3f;

    return (a2 <= epsilon && a2 >= -epsilon &&
            a3 <= epsilon && a3 >= -epsilon &&
            a4 <= epsilon && a4 >= -epsilon &&
            b1 <= epsilon && b1 >= -epsilon &&
            b3 <= epsilon && b3 >= -epsilon &&
            b4 <= epsilon && b4 >= -epsilon &&
            c1 <= epsilon && c1 >= -epsilon &&
            c2 <= epsilon && c2 >= -epsilon &&
            c4 <= epsilon && c4 >= -epsilon &&
            d1 <= epsilon && d1 >= -epsilon &&
            d2 <= epsilon && d2 >= -epsilon &&
            d3 <= epsilon && d3 >= -epsilon &&
            a1 <= 1.f+epsilon && a1 >= 1.f-epsilon &&
            b2 <= 1.f+epsilon && b2 >= 1.f-epsilon &&
            c3 <= 1.f+epsilon && c3 >= 1.f-epsilon &&
            d4 <= 1.f+epsilon && d4 >= 1.f-epsilon);
}

这里是定义了一个模板函数,对于矩阵元素为整型和浮点型都是可行的。

需要引起我们注意的是极小值epsilon的定义:

const static TReal epsilon = 10e-3f;

还有双向比较:

a2 <= epsilon && a2 >= -epsilon

a1 <= 1.f+epsilon && a1 >= 1.f-epsilon

从而保证了元素值为0.0f或1.0f。

相关推荐
18538162800余--2 小时前
短视频矩阵系统贴牌批量剪辑功能开发,支持OEM
线性代数·矩阵·音视频
awk_bioinfo2 小时前
重测序关系矩阵构建方式汇总
矩阵·亲缘关系矩阵·遗传距离·群体遗传
晨曦5432105 小时前
Numpy数组与矩阵——python学习
python·矩阵·numpy
点云SLAM6 小时前
线性代数的本质大白话理解
人工智能·线性代数·矩阵·向量·矩阵变换·空间变换·线性代数的本质
18538162800余--15 小时前
矩阵系统源码搭建热门音乐功能板块开发,支持OEM
线性代数·矩阵
weixin_4284984916 小时前
使用HYPRE库并行装配IJ稀疏矩阵指南: 矩阵预分配和重复利用
算法·矩阵
巷北夜未央20 小时前
空间矩阵的思考
线性代数·矩阵
痛&快乐着1 天前
衡量矩阵数值稳定性的关键指标:矩阵的条件数
线性代数·矩阵
CoderCodingNo2 天前
【GESP】C++二级真题 luogu-B4259 [GESP202503 二级] 等差矩阵
java·c++·矩阵
悲喜自渡7212 天前
线性代数(一些别的应该关注的点)
python·线性代数·机器学习