hdmtrace.c

speech signal process tools

double hdmtrace(M1,iM1,M2,iM2,n)
/* M1 and M2 are hermitian matrices, hdmtrace produces the trace of the
product of the two matrices, which is real */
double M1[],M2[],iM1[],iM2[];
int n;
{
	int j,k,jM,kM;
	double sum;
	sum=0;
	for(j=0,jM=0;j<n;j++,jM+=n)
		{
		for(k=0,kM=0;k<n;k++,kM+=n)
			{
			sum+=(M1[jM+k]*M2[kM+j]-iM1[jM+k]*iM2[kM+j]);
			}
		}
	return(sum);
}