/** Generate points */ static void makePoints( const SURF& surf, std::vector<Point>& points ) { const unsigned n1 = surf.numIntervals1() + (surf.periodic1() ? 0 : 1 ); const unsigned n2 = surf.numIntervals2() + (surf.periodic2() ? 0 : 1 ); for ( unsigned j = 0; j < n2; j++ ) { for ( unsigned i = 0; i < n1; i++ ) { points.push_back( surf.makePoint( i, j ) ); } } }