How to multiply two quaternions by python or numpy
10,266
Here's a little example using your function:
import numpy as np
import random
def quaternion_multiply(quaternion1, quaternion0):
w0, x0, y0, z0 = quaternion0
w1, x1, y1, z1 = quaternion1
return np.array([-x1 * x0 - y1 * y0 - z1 * z0 + w1 * w0,
x1 * w0 + y1 * z0 - z1 * y0 + w1 * x0,
-x1 * z0 + y1 * w0 + z1 * x0 + w1 * y0,
x1 * y0 - y1 * x0 + z1 * w0 + w1 * z0], dtype=np.float64)
N = 4
for i in range(N):
q1 = np.random.rand(4)
q2 = np.random.rand(4)
q = quaternion_multiply(q1, q2)
print("{0} x {1} = {2}".format(q1, q2, q))
Related videos on Youtube
Author by
Biophysics
Updated on September 23, 2022Comments
-
Biophysics over 1 year
I have two quaternions: Q1= w0, x0, y0, z0 and Q2 = w1, x1, y1, z1. I would like to multiply them by using NumPy or Python function which can return 2-d array. I found some pseudocodes on the internet which is written by Christoph Gohlke to do this kind of multiplication. I tried a lot but failed to apply it. Can anyone help me please to do this kind of multiplication? The pseudocodes are here: `
def quaternion_multiply(quaternion1, quaternion0): w0, x0, y0, z0 = quaternion0 w1, x1, y1, z1 = quaternion1 return array([-x1*x0 - y1*y0 - z1*z0 + w1*w0, x1*w0 + y1*z0 - z1*y0 + w1*x0, -x1*z0 + y1*w0 + z1*x0 + w1*y0, x1*y0 - y1*x0 + z1*w0 + w1*z0], dtype=float64)`
-
Biophysics over 7 yearsI have already a python function which can create 2 arrays of quaternions. I am now trying to compute their product by using another function which should return the multiplication output.
-
BPL over 7 years@Biophysics Ok, what part of my example is not clear then? I think I've already answered your question of
I have two quaternions Q1 and Q2. I would like to multiply them by using NumPy or Python function which can return 2-d array
, look at this lineq = quaternion_multiply(q1, q2)
-
Biophysics over 7 yearsYes it works. Thanks :)
-
Biophysics over 7 yearsCan you please help me to get output in arrays from the above function?