Linear Algebra and the C Language/a0al


Code study: normalize_mR();

/* ------------------------------------ */
/* ------------------------------------ */
double **normalize_mR(
double **U,
double **Q
)
{
double **u  = i_mR(rsize_R(U),C1);
double **q  = i_mR(rsize_R(U),C1);	

int c_U     = csize_R(U);
int c;

		for( c = C1; c <= c_U; c++)
          {		
		    c_c_mR(U,c, u,C1); 
		    		           
		    smul_mR(1./norm_R(u),u,q);

            c_c_mR(q,C1, Q,c); 
	      }

  f_mR(u);  
  f_mR(q);  
  
 return(Q);
}
/* ------------------------------------ */
/* ------------------------------------ */