Parallel Kalman filter based time-domain harmonic state estimation