To additionally decrease the computational cost, the Cholesky decomposition is employed for the calculation of the two-electron integrals of the aggregate. The implementation presented here broadens the scope of systems that can be treated within the Frenkel exciton model and we expect it to deepen our understanding of aggregate photophysics and photochemistry.