* Optimize vector load / store. Fixes to alignment. Support typedef with `@simd` and `@align` #2543. Update vector ABI #2542
* Fix alignment issue with indirect arguments.
This commit is contained in:
Christoffer Lerno
2025-10-25 12:31:06 +02:00
committed by GitHub
parent f37e7460aa
commit 423152202f
72 changed files with 2403 additions and 1718 deletions

View File

@@ -1,92 +1,91 @@
module math_matrix @test;
import std::math;
fn void test_mat4()
fn void test_mat4_translate()
{
{
Matrix4 mat = MATRIX4_IDENTITY;
Matrix4 mat2 = { 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1 };
Matrix4 calc = mat.mul(mat2);
assert(calc.m == mat.m);
assert(mat * mat2 == mat);
Matrix4 mat = MATRIX4_IDENTITY;
Matrix4 mat2 = { 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1 };
Matrix4 calc = mat.mul(mat2);
assert(calc.m == mat.m);
assert(mat * mat2 == mat);
Matrix4 translated = mat.translate({0.0, 0.0, 0.0});
assert(translated.m == mat.m);
};
{
Matrix4 mat = { 1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 3, 4, 5, 6, 7, 8 };
Matrix4 mat2 = { 8, 7, 6, 5, 4, 3, 2, 1, 8, 7, 6, 5, 4, 3, 2, 1 };
Matrix4 calc = mat.mul(mat2);
Matrix4 value = { 56, 46, 36, 26, 152, 126, 100, 74, 56, 46, 36, 26, 152, 126, 100, 74 };
assert(calc.m == value.m);
assert(mat * mat2 == value);
};
{
Matrix4 result = {
0.988936, 0.000000, -0.148340, -0.988936,
-0.014599, 0.995146, -0.097325, -2.970838,
0.147620, 0.098414, 0.984136, -20.765262,
0.000000, 0.000000, 0.000000, 1.000000
};
Matrix4f result_f = {
0.988936, 0.000000, -0.148340, -0.988936,
-0.014599, 0.995146, -0.097325, -2.970838,
0.147620, 0.098414, 0.984136, -20.765262,
0.000000, 0.000000, 0.000000, 1.000000
};
Matrix4 result_transposed = {
0.988936, -0.014599, 0.147620, 0.000000,
0.000000, 0.995146, 0.098414, 0.000000,
-0.148340, -0.097325, 0.984136, 0.000000,
-0.988936, -2.970838, -20.765262, 1.000000
};
Matrix4f result_transposed_f = {
0.988936, -0.014599, 0.147620, 0.000000,
0.000000, 0.995146, 0.098414, 0.000000,
-0.148340, -0.097325, 0.984136, 0.000000,
-0.988936, -2.970838, -20.765262, 1.000000
};
Matrix4 look_at = matrix::look_at{double}({4.0, 5.0, 20.0}, {1.0, 3.0, 0.0}, {0.0, 1.0, 0.0});
Matrix4f look_at_f = matrix::look_at{float}({4.0, 5.0, 20.0}, {1.0, 3.0, 0.0}, {0.0, 1.0, 0.0});
assert(math::round_to_decimals((double[<16>])look_at.m, 4) == math::round_to_decimals((double[<16>])result.m, 4));
assert(math::round_to_decimals((float[<16>])look_at_f.m, 4) == math::round_to_decimals((float[<16>])result_f.m, 4));
assert(math::round_to_decimals((double[<16>])result_transposed.m, 4) == math::round_to_decimals((double[<16>])look_at.transpose().m, 4));
assert(math::round_to_decimals((float[<16>])result_transposed_f.m, 4) == math::round_to_decimals((float[<16>])look_at_f.transpose().m, 4));
};
{
Matrix4 result = {
1.857087, 0.000000, 0.000000,
0.000000, 0.000000, 2.414214,
0.000000, 0.000000, 0.000000, 0.000000,
-1.000200, -0.200020, 0.000000, 0.000000,
-1.000000, 0.000000
};
Matrix4f result_f = {
1.857087, 0.000000, 0.000000,
0.000000, 0.000000, 2.414214,
0.000000, 0.000000, 0.000000, 0.000000,
-1.000200, -0.200020, 0.000000, 0.000000,
-1.000000, 0.000000
};
Matrix4 perspective = matrix4_perspective(math::deg_to_rad(45), 1.3, 0.1, 1000);
Matrix4f perspective_f = matrix4f_perspective((float)math::deg_to_rad(45), 1.3, 0.1, 1000);
assert(math::round_to_decimals((double[<16>])result.m, 4) == math::round_to_decimals((double[<16>])perspective.m, 4));
assert(math::round_to_decimals((float[<16>])result_f.m, 4) == math::round_to_decimals((float[<16>])perspective_f.m, 4));
};
Matrix4 translated = mat.translate({0.0, 0.0, 0.0});
assert(translated.m == mat.m);
}
fn void test_mat4_mul()
{
Matrix4 mat = { 1, 2, 3, 4, 5, 6, 7, 8, 1, 2, 3, 4, 5, 6, 7, 8 };
Matrix4 mat2 = { 8, 7, 6, 5, 4, 3, 2, 1, 8, 7, 6, 5, 4, 3, 2, 1 };
Matrix4 calc = mat.mul(mat2);
Matrix4 value = { 56, 46, 36, 26, 152, 126, 100, 74, 56, 46, 36, 26, 152, 126, 100, 74 };
assert(calc.m == value.m);
assert(mat * mat2 == value);
}
fn void test_mat4_lookat()
{
Matrix4 result = {
0.988936, 0.000000, -0.148340, -0.988936,
-0.014599, 0.995146, -0.097325, -2.970838,
0.147620, 0.098414, 0.984136, -20.765262,
0.000000, 0.000000, 0.000000, 1.000000
};
Matrix4f result_f = {
0.988936, 0.000000, -0.148340, -0.988936,
-0.014599, 0.995146, -0.097325, -2.970838,
0.147620, 0.098414, 0.984136, -20.765262,
0.000000, 0.000000, 0.000000, 1.000000
};
Matrix4 result_transposed = {
0.988936, -0.014599, 0.147620, 0.000000,
0.000000, 0.995146, 0.098414, 0.000000,
-0.148340, -0.097325, 0.984136, 0.000000,
-0.988936, -2.970838, -20.765262, 1.000000
};
Matrix4f result_transposed_f = {
0.988936, -0.014599, 0.147620, 0.000000,
0.000000, 0.995146, 0.098414, 0.000000,
-0.148340, -0.097325, 0.984136, 0.000000,
-0.988936, -2.970838, -20.765262, 1.000000
};
Matrix4 look_at = matrix::look_at{double}({4.0, 5.0, 20.0}, {1.0, 3.0, 0.0}, {0.0, 1.0, 0.0});
Matrix4f look_at_f = matrix::look_at{float}({4.0, 5.0, 20.0}, {1.0, 3.0, 0.0}, {0.0, 1.0, 0.0});
assert(math::round_to_decimals((double[<16>])look_at.m, 4) == math::round_to_decimals((double[<16>])result.m, 4));
assert(math::round_to_decimals((float[<16>])look_at_f.m, 4) == math::round_to_decimals((float[<16>])result_f.m, 4));
assert(math::round_to_decimals((double[<16>])result_transposed.m, 4) == math::round_to_decimals((double[<16>])look_at.transpose().m, 4));
assert(math::round_to_decimals((float[<16>])result_transposed_f.m, 4) == math::round_to_decimals((float[<16>])look_at_f.transpose().m, 4));
}
fn void test_mat4_perspective()
{
Matrix4 result = {
1.857087, 0.000000, 0.000000,
0.000000, 0.000000, 2.414214,
0.000000, 0.000000, 0.000000, 0.000000,
-1.000200, -0.200020, 0.000000, 0.000000,
-1.000000, 0.000000
};
Matrix4f result_f = {
1.857087, 0.000000, 0.000000,
0.000000, 0.000000, 2.414214,
0.000000, 0.000000, 0.000000, 0.000000,
-1.000200, -0.200020, 0.000000, 0.000000,
-1.000000, 0.000000
};
Matrix4 perspective = matrix4_perspective(math::deg_to_rad(45), 1.3, 0.1, 1000);
Matrix4f perspective_f = matrix4f_perspective((float)math::deg_to_rad(45), 1.3, 0.1, 1000);
assert(math::round_to_decimals((double[<16>])result.m, 4) == math::round_to_decimals((double[<16>])perspective.m, 4));
assert(math::round_to_decimals((float[<16>])result_f.m, 4) == math::round_to_decimals((float[<16>])perspective_f.m, 4));
}
fn void test_mat3()
{