PRACE NPT

void clr_arr( int n, int m, double **a, val )

{

   int i, j, mf;

// ---------------------------------------------------------------------

   mf = m - m%4;

   for( i = 0; i < n; i++ ) {

      for( j = 0; j < mf; j+= 4 ) {

         a[i][j] = 0.0; a[i][j+1] = 0.0; a[i][j+2] = 0.0; a[i][j+3] = 0.0;

      }

   }

   for( i = 0; i < n; i++ ) {

      for( j = mf; j < m; j++ ) {

         a[i][j] = 0.0;

      }

   }

}