(N,k)=(3,2) family grid : per system (zeros found, isolated, non-isolated) = [(8, np.int64(8), 0)]; max isolated 8 <= k^N = 8: True
(N,k)=(3,3) family grid : per system (zeros found, isolated, non-isolated) = [(27, np.int64(27), 0)]; max isolated 27 <= k^N = 27: True
(N,k)=(4,2) family grid : per system (zeros found, isolated, non-isolated) = [(16, np.int64(16), 0)]; max isolated 16 <= k^N = 16: True
(N,k)=(3,2) family pgrid: per system (zeros found, isolated, non-isolated) = [(8, np.int64(8), 0), (8, np.int64(8), 0), (8, np.int64(8), 0), (8, np.int64(8), 0)]; max isolated 8 <= k^N = 8: True
(N,k)=(3,3) family pgrid: per system (zeros found, isolated, non-isolated) = [(11, np.int64(11), 0), (19, np.int64(19), 0), (11, np.int64(11), 0)]; max isolated 19 <= k^N = 27: True
(N,k)=(4,2) family pgrid: per system (zeros found, isolated, non-isolated) = [(16, np.int64(16), 0), (8, np.int64(8), 0), (16, np.int64(16), 0)]; max isolated 16 <= k^N = 16: True
(N,k)=(3,2) family thru : per system (zeros found, isolated, non-isolated) = [(8, np.int64(8), 0), (8, np.int64(8), 0), (8, np.int64(8), 0), (8, np.int64(8), 0), (8, np.int64(8), 0)]; max isolated 8 <= k^N = 8: True
(N,k)=(3,3) family thru : per system (zeros found, isolated, non-isolated) = [(25, np.int64(25), 0), (19, np.int64(19), 0), (25, np.int64(25), 0)]; max isolated 25 <= k^N = 27: True
(N,k)=(4,2) family thru : per system (zeros found, isolated, non-isolated) = [(16, np.int64(16), 0), (12, np.int64(12), 0), (14, np.int64(14), 0)]; max isolated 16 <= k^N = 16: True
(N,k)=(3,2) family curve: per system (zeros found, isolated, non-isolated) = [(300, np.int64(0), 300), (3, np.int64(3), 0), (294, np.int64(0), 294), (298, np.int64(0), 298), (298, np.int64(0), 298), (3, np.int64(3), 0)]; max isolated 3 <= k^N = 8: True
