У меня есть массив obj
, в котором хранятся матрицы, каждый из которых мне надо умножить на матрицу projection_vec
. Когда программа перебирает все элементы obj
, и умножает каждый из них на projection_vec
, и тогда результаты умножения такие какие мне нужны, но когда программа умножает сразу весь массив obj
на projection_vec
, результаты немного отличаются, я не знаю как это исправить, так как умножение сразу всего массива obj
в разы быстрее чем перебирать и умножать каждый элемент массива по отдельности, скажите пожалуйста что мне нужно изменить во втором фрагменте кода, чтобы результаты были как в первом фрагменте, но при этом код выполнялся с такой же скоростью (или даже быстрее) как второй фрагмент.
Фрагмент 1
import numpy as np
from math import *
def matrix_multiply(inp_tri, matrix):
out_tri = inp_tri @ matrix[:3, :3] + matrix[3, :3]
w = inp_tri @ matrix[:3, 3] + matrix[3, 3]
w[w == 0] = 1.0
out_tri = out_tri / w[:, None]
return out_tri
perspective_matrix = np.array([
[(1.0 / tan(45 / 180 * pi)), 0.0, 0.0, 0.0],
[0.0, (WIDTH / HEIGHT) * (1.0 / tan(45 / 180 * pi)), 0.0, 0.0],
[0.0, 0.0, 1.0, 1.0],
[0.0, 0.0, -1.0, 0.0]
])
def get_projection(translated_vec):
projection_vec = matrix_multiply(translated_vec, perspective_matrix)
projection_vec += 1.0
projection_vec = projection_vec * H_SIZE
return projection_vec
obj = np.array([[
[0, 0, 0],
[1, 1, 1],
[-2, -2, 2]
],
[
[1, 1, 1],
[0, 1, 2],
[-2, -2, -2]
]])
drawig_obj = []
for i, translated_vec in enumerate(obj):
translated_vec = get_projection(translated_vec)
drawing_obj.append(translated_vec)
print(drawing_obj)
Выходные данные:
[(array([[ nan, nan, -inf],
[ 1.28000000e+03, 1.00000000e+03, 1.00000000e+00],
[-1.42108547e-13, -2.80000000e+02, 1.50000000e+00]]), (255, 255, 255)), (array([[1.28e+03, 1.00e+03, 1.00e+00],
[6.40e+02, 6.80e+02, 1.50e+00],
[1.28e+03, 1.00e+03, 2.50e+00]]), (255, 255, 255))]
Фрагмент 2
import numpy as np
from math import *
def matrix_multiply(inp_tri, matrix):
out_tri = inp_tri @ matrix[:3, :3] + matrix[3, :3]
w = inp_tri @ matrix[:3, 3] + matrix[3, 3]
w[w == 0] = 1.0
out_tri = out_tri / w[:, None]
return out_tri
perspective_matrix = np.array([
[(1.0 / tan(45 / 180 * pi)), 0.0, 0.0, 0.0],
[0.0, (WIDTH / HEIGHT) * (1.0 / tan(45 / 180 * pi)), 0.0, 0.0],
[0.0, 0.0, 1.0, 1.0],
[0.0, 0.0, -1.0, 0.0]
])
def get_projection(translated_vec):
projection_vec = matrix_multiply(translated_vec, perspective_matrix)
projection_vec += 1.0
projection_vec = projection_vec * H_SIZE
return projection_vec
obj = np.array([[
[0, 0, 0],
[1, 1, 1],
[-2, -2, 2]
],
[
[1, 1, 1],
[0, 1, 2],
[-2, -2, -2]
]])
projection_vecs = get_projection(obj)
print(projection_vecs)
Выходные данные:
[(array([[ 6.40e+02, 3.60e+02, 5.00e-01],
[ 1.28e+03, 1.00e+03, 1.00e+00],
[-6.40e+02, -9.20e+02, 1.50e+00]]), (255, 255, 255)), (array([[ 1.28e+03, 6.80e+02, 1.00e+00],
[ 6.40e+02, 6.80e+02, 5.00e-01],
[-6.40e+02, -2.80e+02, 2.50e+00]]), (255, 255, 255))]