orthofst.c

speech signal process tools

/* fast orthogonalization, no determinent */
orthofst(M,Minv,Tinv,n)
register double *M,*Minv,*Tinv;
register int n;
{
	register int i,j;
	register double *s,*s1;
	double log(),dot();

	if(orthog(M,Minv,n))
		return(1);
	/* orthog returns an orthogonal matrix O in M and an upper triagonal
	matrix T in Minv such that M=O*T */

	if(triinv(Minv,Tinv,n))
		return(1);
	for(i=0,s=Minv;i<n;i++,Tinv+=n)
		{
		for(j=0,s1=M;j<n;j++,s1+=n,s++)
			*s=dot(Tinv,1,s1,1,n);
		}
	return(0);
}