A QGIS Python snippet to select multipart features in the active layer

find_multipart.py

# Find all multipart features in the active layer
l = iface.activeLayer()
iter = l.getFeatures()
geoms = []
for feature in iter:
    geom = feature.geometry()
    if geom.isMultipart():
        l.select(feature.id())
        geoms.append(geom)

print 'There are %i multipart features in this layer' % len(geoms)

I use QGIS 3.22, but this code selects all features in an active layer, not the multipart feature. Do you update it already? If not, I may create an alternative code after this.