mirrors/Odin
1
0
Fork 0
mirror of https://github.com/odin-lang/Odin.git synced 2026-05-25 05:09:53 +00:00
Code Issues Packages Projects Releases Wiki Activity

linalg: gram schmidt in any order

Browse Source
This commit is contained in:
jakubtomsu
2026-05-18 19:34:53 +02:00
parent c19540dda4
commit b57db64c8e
1 changed files with 27 additions and 42 deletions
Download Patch File Download Diff File Expand all files Collapse all files

69
core/math/linalg/specific.odin
View File

@@ -2768,52 +2768,16 @@ matrix2_orthonormalize :: proc{
@(require_results)
matrix3_orthonormalize_f16 :: proc "contextless" (m: Matrix3f16) -> (r: Matrix3f16) #no_bounds_check {
r = m
r[0] = normalize(m[0])
d0 := dot(r[0], r[1])
r[1] -= r[0] * d0
r[1] = normalize(r[1])
d1 := dot(r[1], r[2])
d0 = dot(r[0], r[2])
r[2] -= r[0]*d0 + r[1]*d1
r[2] = normalize(r[2])
return
matrix3_orthonormalize_f16 :: proc "contextless" (m: Matrix3f16) -> Matrix3f16 #no_bounds_check {
return matrix3_gram_schmidt(m, 0, 1, 2)
}
@(require_results)
matrix3_orthonormalize_f32 :: proc "contextless" (m: Matrix3f32) -> (r: Matrix3f32) #no_bounds_check {
r = m
r[0] = normalize(m[0])
d0 := dot(r[0], r[1])
r[1] -= r[0] * d0
r[1] = normalize(r[1])
d1 := dot(r[1], r[2])
d0 = dot(r[0], r[2])
r[2] -= r[0]*d0 + r[1]*d1
r[2] = normalize(r[2])
return
matrix3_orthonormalize_f32 :: proc "contextless" (m: Matrix3f32) -> Matrix3f32 #no_bounds_check {
return matrix3_gram_schmidt(m, 0, 1, 2)
}
@(require_results)
matrix3_orthonormalize_f64 :: proc "contextless" (m: Matrix3f64) -> (r: Matrix3f64) #no_bounds_check {
r = m
r[0] = normalize(m[0])
d0 := dot(r[0], r[1])
r[1] -= r[0] * d0
r[1] = normalize(r[1])
d1 := dot(r[1], r[2])
d0 = dot(r[0], r[2])
r[2] -= r[0]*d0 + r[1]*d1
r[2] = normalize(r[2])
return
matrix3_orthonormalize_f64 :: proc "contextless" (m: Matrix3f64) -> Matrix3f64 #no_bounds_check {
return matrix3_gram_schmidt(m, 0, 1, 2)
}
matrix3_orthonormalize :: proc{
matrix3_orthonormalize_f16,
@@ -2822,6 +2786,27 @@ matrix3_orthonormalize :: proc{
}
@(require_results)
matrix3_gram_schmidt :: proc "contextless" (m: matrix[3, 3]$E, $A, $B, $C: int) -> (r: matrix[3, 3]E)
where A != B, A != C, B != C #no_bounds_check
{
r = m
r[A] = normalize(m[A])
d0 := dot(r[A], r[B])
r[B] -= r[A] * d0
r[B] = normalize(r[B])
d1 := dot(r[B], r[C])
d0 = dot(r[A], r[C])
r[C] -= r[A]*d0 + r[B]*d1
r[C] = normalize(r[C])
return
}
@(require_results)
vector3_orthonormalize_f16 :: proc "contextless" (x, y: Vector3f16) -> (z: Vector3f16) {
return normalize(x - y * dot(y, x))