cdmmult.c

speech signal process tools

cdmmult(A,iA,B,iB,C,iC,M,iM,n)/* multiply complex matrices */
/* M and iM are temporary matrices */
/* C=A*B */
/* A+iA' etc */
/* C=A*B - A'*B' */
/* iC=A*B' + A'*B */
register double *A,*iA,*B,*iB,*C,*iC,*M,*iM;
register int n;
{
	dmmult(A,B,M,n);
	dmmult(iA,iB,iM,n);
	dmsub(M,iM,C,n);
	dmmult(A,iB,M,n);
	dmmult(iA,B,iM,n);
	dmadd(M,iM,iC,n);
	return;
}