Zxc200611/pocketpy
1
0
Fork 0
mirror of https://github.com/pocketpy/pocketpy synced 2025-10-20 11:30:18 +00:00
Code Issues Packages Projects Releases Wiki Activity

...

Browse Source
This commit is contained in:
blueloveTH 2023-08-05 17:18:54 +08:00
parent 9998adcecd
commit 82c453dbb3
1 changed files with 6 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
tests/80_linalg.py
View File

@ -183,6 +183,10 @@ static_test_mat_float= mat3x3([
[9.558042327611506, -3.336280256662496, 4.951381528057387]] [9.558042327611506, -3.336280256662496, 4.951381528057387]]
) )
static_test_mat_float_inv = mat3x3([[ 0.32265243, 0.15808159, -0.09939472],
[ 0.04199553, 0.13813096, 0.00408326],
[-0.59454451, -0.21208362, 0.39658464]])
static_test_mat_int = mat3x3([ static_test_mat_int = mat3x3([
[1, 2, 3], [1, 2, 3],
[4, 5, 6], [4, 5, 6],
@ -373,18 +377,9 @@ except:
test_mat_copy = test_mat.copy() test_mat_copy = test_mat.copy()
assert test_mat_copy.transpose() == test_mat_copy.transpose().transpose().transpose() assert test_mat_copy.transpose() == test_mat_copy.transpose().transpose().transpose()
# test inverse # test inverse
test_mat_copy = test_mat.copy() assert static_test_mat_float.inverse() == static_test_mat_float_inv
if round(test_mat_copy.determinant(), 5) != 0: assert ~static_test_mat_float == static_test_mat_float_inv
list_mat = [[0,0,0], [0,0,0], [0,0,0]]
list_mat_2 = [[0,0,0], [0,0,0], [0,0,0]]
for i in range(3):
for j in range(3):
list_mat[i][j] = test_mat[i, j]
list_mat_2[i][j] = str(test_mat_copy.inverse()[i, j])[:6]
assert [[str(e)[:6] for e in layer] for layer in calculate_inverse(list_mat)] == list_mat_2
try: try:
mat3x3([[1, 2, 3], [2, 4, 6], [3, 6, 9]]).inverse() mat3x3([[1, 2, 3], [2, 4, 6], [3, 6, 9]]).inverse()
@ -392,17 +387,6 @@ try:
except: except:
pass pass
# test __invert__
test_mat_copy = test_mat.copy()
if round(test_mat_copy.determinant(), 5) != 0:
list_mat = [[0,0,0], [0,0,0], [0,0,0]]
list_mat_2 = [[0,0,0], [0,0,0], [0,0,0]]
for i in range(3):
for j in range(3):
list_mat[i][j] = test_mat[i, j]
list_mat_2[i][j] = str((~test_mat_copy)[i, j])[:6]
assert [[str(e)[:6] for e in layer] for layer in calculate_inverse(list_mat)] == list_mat_2
try: try:
~mat3x3([[1, 2, 3], [2, 4, 6], [3, 6, 9]]) ~mat3x3([[1, 2, 3], [2, 4, 6], [3, 6, 9]])
raise Exception('未能拦截错误 ValueError("matrix is not invertible") 在 test_mat_copy 的行列式为0') raise Exception('未能拦截错误 ValueError("matrix is not invertible") 在 test_mat_copy 的行列式为0')