Flexible hemodynamic coupling in EEG-fMRI: circumventing the GLM by structured coupled matrix-tensor factorization