spdmtrace.c

speech signal process tools

double
spdmtrace(M1,M2,n)/* trace of M1*M2, M1 double, M2 short */
register double *M1;
register short *M2;
register int n;
{
	register int i;
	register double trace;
	double spdot();

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