To add some performance improvements to @Ben W's answer by using a spatial index and a break command:
point_lyr = QgsProject.instance().mapLayersByName('points')
poly_lyr = QgsProject.instance()....
Using the data you provided, I get meaningful results using Join Attributes by Location (Summary), see screenshot below and be sure to make settings accordingly.
Another option ...