00001 #include "cs.h" 00002 static int cs_nonzero (int i, int j, double aij, void *other) 00003 { 00004 return ( fabs( aij ) > 2.3e-16 ) ; 00005 } 00006 int cs_dropzeros (cs *A) 00007 { 00008 return (cs_fkeep (A, &cs_nonzero, NULL)) ; /* keep all nonzero entries */ 00009 }