trace.c

speech signal process tools

double 
trace(M,n)/* return trace of double matrix */
register double *M;
register int n;
{
	register int i,con;
	register double sum;

	sum=0;
	con=n+1;
	for(i=0;i<n;i++,M+=con)
		sum+= *M;
	return(sum);
}