www.pudn.com > calc.rar > i7I.c


/* i7I.c */ 
#include  
#include "integer.h" 
#include "fun.h" 
 
MPI *DOTRI(MPMATI *Mptr, USI i, USI j) 
/* 
 * *Aptr is the dot product of rows i and j in *Mptr. 
 */ 
{ 
	unsigned int k; 
	MPI *H, *Aptr, *TempI; 
	 
	Aptr = ZEROI(); 
	for (k = 0; k <= Mptr->C - 1; k++) 
	{ 
		H = MULTI(Mptr->V[i][k], Mptr->V[j][k]); 
		TempI = Aptr; 
		Aptr = ADDI(Aptr, H); 
		FREEMPI(H); 
		FREEMPI(TempI); 
	} 
	return (Aptr); 
}