Rev | Line | |
---|
[37] | 1 | function c = mcross(a, b) |
---|
| 2 | % Matrix cross (vector) product of two pure quaternions. Like cross but matrix |
---|
| 3 | % not elementwise. |
---|
| 4 | |
---|
| 5 | % Copyright © 2005 Stephen J. Sangwine and Nicolas Le Bihan. |
---|
| 6 | % See the file : Copyright.m for further details. |
---|
| 7 | |
---|
| 8 | error(nargchk(2, 2, nargin)), error(nargoutchk(0, 1, nargout)) |
---|
| 9 | |
---|
| 10 | if ~isa(a, 'quaternion') | ~isa(b, 'quaternion') |
---|
| 11 | error('Matrix cross product is not defined for a quaternion and a non-quaternion.') |
---|
| 12 | end |
---|
| 13 | |
---|
| 14 | if ~ispure(a) | ~ispure(b) |
---|
| 15 | error('Mcross product is defined only for pure quaternions.') |
---|
| 16 | end |
---|
| 17 | |
---|
| 18 | c = quaternion(y(a) * z(b) - z(a) * y(b), ... |
---|
| 19 | z(a) * x(b) - x(a) * z(b), ... |
---|
| 20 | x(a) * y(b) - y(a) * x(b)); |
---|
Note: See
TracBrowser
for help on using the repository browser.