dmtrace.c

speech signal process tools

double dmtrace(M1,M2,n)/* returns trace of M1*M2, M1 and M2 matrices  */
register double *M1,*M2;
register int n;
{
	register int i;
	register double trace;
	double dot();

	trace=0;
	for(i=0;i<n;i++,M1+=n,M2++)
		trace+=dot(M1,1,M2,n,n);
	return(trace);
}