function [a,b]=cdot5(u,v)
% gamma5 inner product u'*gamma5*v
u=u(:);
v=v(:);
ndim=max(size(u));
u=reshape(u,ndim/4,4);
v=reshape(v,ndim/4,4);
a1=u(:,1)'*v(:,1);
a2=u(:,2)'*v(:,2);
a3=u(:,3)'*v(:,3);
a4=u(:,4)'*v(:,4);
a=conj(-a1-a2+a3+a4);
b=conj(a1+a2+a3+a4);
% Copyright (C) 2006-2007 Artan Borici.
% This program is a free software licensed under the terms of the GNU General Public License