MATLAB and C Programs - Digital Signal Processing: A Practitioner's Approach - page 172

float val_A;

for ( row = 0; row < max_row_A; row++)

{
for ( col = 0; col < max_col_A; col++)

{
index_A = max_col_A*row + col;
val_A = 0;
if ( col == row ) val_A =1;
*(some_A + index_A) = val_A;
}
}
}

float mat_dot( float *temp_A, float *temp_B)
{
int max_row_A = 3;
int col,k,index_A,index_B;
float val_A,val_B,sum;
sum = 0.0;
for ( col = 0; col < max_row_A; col++)

{
index_A = col;
index_B = col;

val_A = *(temp_A + index_A);
val_B = *(temp_B + index_B);

sum = sum + val_A*val_B;
}
return(sum);
}

float randn(void)
{
float intmax,udf,sum;
int rnd,k;

intmax = 64000;
sum = 0.0;
for (k = 0; k

< 12; k++)