PRACE NPT

#include<math.h>

#include<float.h>

#include<stdio.h>

int check( int ncols, int nrows, int nelmts, int indx[], int rowp[],

              double outvec[] )

{

   double eps, rowsum;

   double diff;

   int    i, ok = 0;

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

   eps = 2.0*ncols*DBL_EPSILON;

   for( i = 0; i < nrows - 1; i++ ) {

      rowsum = (double)(rowp[i+1] - rowp[i]);

      if ( fabs( outvec[i] - rowsum ) > eps ) ok++;

   }

   rowsum = (double)(nelmts - rowp[nrows-1]);

   if ( fabs( outvec[nrows-1] - rowsum ) > eps ) ok++;

   return( ok );

}