test reason why matrix needs its own type

This commit is contained in:
Quaternions 2024-08-30 12:52:54 -07:00
parent d713b96ad3
commit 263f0d35d4

View File

@ -4,7 +4,7 @@ use fixed_wide_traits::wide::WideDot;
use crate::Vector3; use crate::Vector3;
type Planar64=fixed_wide::types::I32F32; type Planar64=fixed_wide::types::I32F32;
//type Planar64Wide1=fixed::types::I64F64; type Planar64Wide1=fixed_wide::types::I64F64;
//type Planar64Wide2=fixed_wide::types::I128F128; //type Planar64Wide2=fixed_wide::types::I128F128;
type Planar64Wide3=fixed_wide::types::I256F256; type Planar64Wide3=fixed_wide::types::I256F256;
@ -49,3 +49,17 @@ fn wide_vec3_length_squared(){
assert_eq!(v3,Planar64Wide3::from(3i128.pow(8)*3)); assert_eq!(v3,Planar64Wide3::from(3i128.pow(8)*3));
} }
#[test]
fn wide_matrix_dot_vs_vec_of_vec_dot(){
let vv=Vector3::<crate::Vector2<_>>::from_value_2d(Planar64::from(3));
// do the dot product of the inner vectors multiplied component wise
// this lowers the rank of the data structure and is kind of a weird operation lol
let vv_dot=vv.wide_dot(vv);
assert_eq!(vv_dot,crate::Vector2::from_value(Planar64Wide1::from(3i128.pow(3))));
//let m=Matrix3::<Vector3<_>>::from_value_2d(Planar64::from(3));
//normal matrix product
//let m_dot=m.wide_dot(m);
//assert_eq!(m_dot,Vector3::from_value(Planar64Wide1::from(3i128.pow(3))));
}