programing How does the rotation values work pls?

[GTM_ssj 0]

I just recently came up with this idea of making a script that turns .miobject files to .mimodel files, that is, to simplify mineimator models

but there's this issue when I try to apply the rotation changes affected by parent objects to a given point, sometimes it works but sometimes give bugged results, especially when multiple axies are rotated at the same time.

I think it's have something to do with the method MI used to rotate objects, but I don't know how exactly it is done tho 

anyone has an idea?

par_pos[:-1, -1] = [0, 0, 0]
par_rot[:-1, :-1] = R.from_euler('xyz', [90, -90, 90], degrees = True).as_matrix()
par_sca[:-1, :-1] = np.diag([1, 1, 1])
chi_pos[:-1, -1] = [8, 8, -32]

matrix = np.dot(np.dot(par_pos, par_rot), par_sca)#scale * rotation * position (?
final = np.dot(matrix, chi_pos)
print(final[:-1, -1])

here's the test code for instance, I did it  using a 4x4 matrix for the transformation matrix

for some reason it gave a wrong answer especially when both rotating at X and Y axis