Add option to skip inactive points in foreach_point_neighbor and update without assembling an iterator
#115
+60
−40