orange/omath
1
0
Fork 0
mirror of https://github.com/orange-cpp/omath.git synced 2026-06-08 16:24:35 +00:00
Code Issues Packages Projects Releases Wiki Activity

added safety check

Browse Source
This commit is contained in:
Orange
2026-05-17 09:38:19 +03:00
parent cebcfc411d
commit e4087165b9
1 changed files with 9 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
include/omath/linear_algebra/mat.hpp
+9 -3
View File
@@ -628,10 +628,16 @@ namespace omath
return static_cast<Type>(std::sqrt(x * x + y * y + z * z)); return static_cast<Type>(std::sqrt(x * x + y * y + z * z));
}; };
const auto scale_x = column_length(mat.at(0, 0), mat.at(1, 0), mat.at(2, 0));
const auto scale_y = column_length(mat.at(0, 1), mat.at(1, 1), mat.at(2, 1));
const auto scale_z = column_length(mat.at(0, 2), mat.at(1, 2), mat.at(2, 2));
constexpr auto epsilon = std::numeric_limits<Type>::epsilon();
return { return {
column_length(mat.at(0, 0), mat.at(1, 0), mat.at(2, 0)), std::abs(scale_x) < epsilon ? Type{1} : scale_x,
column_length(mat.at(0, 1), mat.at(1, 1), mat.at(2, 1)), std::abs(scale_y) < epsilon ? Type{1} : scale_y,
column_length(mat.at(0, 2), mat.at(1, 2), mat.at(2, 2)), std::abs(scale_z) < epsilon ? Type{1} : scale_z,
}; };
} }