10
This function takes two block matrices A and B and multiplies
11
each element of B by the corresponding element of A
13
\param **A: block matrix A
14
\param **B: block matrix B
15
\param nrows: number of rows of A and B
16
\param ncols: number of columns of A and B
20
void dirprd_block(double **A, double **B, int rows, int cols)
26
size = ((long) rows) * ((long) cols);
32
for(i=0; i < size; i++, a++, b++) (*b) = (*a) * (*b);