doublespdmtrace(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);}