For iEM, the time required to update µ with s i depends only on the number of nonzero components of s i .
the sEM update is µ ← (1−η k )µ+η k s i , and a naïve implementation would take time proportional to the total number of components.