PRACE NPT

void bufin( int m, int n, double *b, double **a )

{

   int shift;

   int j, k;

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

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

      shift = j*n;

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

         b[shift+k] = a[j][k];

      }

   }

}