X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=CreaPhase.git;a=blobdiff_plain;f=octave_packages%2Fquaternion-2.0.0%2F%40quaternion%2Fmpower.m;fp=octave_packages%2Fquaternion-2.0.0%2F%40quaternion%2Fmpower.m;h=e8c865b10bf3deefcaa6bb65a0d700f9d056c597;hp=0000000000000000000000000000000000000000;hb=f5f7a74bd8a4900f0b797da6783be80e11a68d86;hpb=1705066eceaaea976f010f669ce8e972f3734b05
diff --git a/octave_packages/quaternion-2.0.0/@quaternion/mpower.m b/octave_packages/quaternion-2.0.0/@quaternion/mpower.m
new file mode 100644
index 0000000..e8c865b
--- /dev/null
+++ b/octave_packages/quaternion-2.0.0/@quaternion/mpower.m
@@ -0,0 +1,56 @@
+## Copyright (C) 2010, 2011 Lukas F. Reichlin
+##
+## This program is free software: you can redistribute it and/or modify
+## it under the terms of the GNU General Public License as published by
+## the Free Software Foundation, either version 3 of the License, or
+## (at your option) any later version.
+##
+## This program is distributed in the hope that it will be useful,
+## but WITHOUT ANY WARRANTY; without even the implied warranty of
+## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+## GNU General Public License for more details.
+##
+## You should have received a copy of the GNU General Public License
+## along with this program. If not, see .
+
+## -*- texinfo -*-
+## Matrix power operator of quaternions. Used by Octave for "q^x".
+
+## Author: Lukas Reichlin
+## Created: May 2010
+## Version: 0.2
+
+function q = mpower (a, b)
+
+ [r, c] = size (a);
+
+ if (r != c)
+ error ("quaternion: mpower: quaternion matrix must be square");
+ endif
+
+ if (r == 1 && c == 1) # a scalar, b?
+ q = a .^ b; # b could be a quaternion
+ elseif (is_real_array (b) && isscalar (b) && fix (b) == b)
+ e = fix (abs (b));
+ switch (sign (b))
+ case -1 # q^-e
+ a = inv (a);
+ q = a;
+ case 0 # q^0
+ q = eye (r); # alternative: q = quaternion (eye (r))
+ return;
+ case 1; # q^e
+ q = a;
+ endswitch
+ for k = 2 : e
+ q *= a; # improvement?: q^8 = ((q^2)^2)^2, q^9 = (((q^2)^2)^2)*q
+ endfor
+ else
+ error ("quaternion: mpower: case not implemented yet");
+ q = expm (logm (a) * b); # don't know whether this formula is correct
+ endif
+
+ ## TODO: - q1 ^ q2
+ ## - arrays
+
+endfunction