Consistency and Computation of Regularized MLEs for Multivariate Hawkes Processes