innerh.m

来自「这是几何代数的matlab工具包」· M 代码 · 共 28 行

function C = innerH(m,n)
%innerH(m,n): Computes the Hestenes inner product between two multivectors.
% innerH(m,n) = sum_r,s grade(grade(m,r)*grade(n,s),|r-s|) for r,s>=0
%                unless grade(m) or grade(n) is zero, 
%                in which case inner(m,n)=0
%
%See also gable, inner.

% GABLE, Copyright (c) 1999, University of Amsterdam
% Copying, use and development for non-commercial purposes permitted.
%          All rights for commercial use reserved; for more information
%          contact Leo Dorst ([e-mail]).
%
%          This software is unsupported.

   if (isa(m,'GA')&isa(n,'GA'))
     T = GAinnerH(m,n);
     if GAautoscalar&isascalar(T)
	C=T.m(1);
     else
	C=T;
     end
   elseif GAautoscalar
     C=0;
   else
     C=GA([0;0;0;0;0;0;0;0]);
   end