sped up harmonicLocalization.m