Parallelized up to Hv_prep, next step is to use collapse directive on Alavi Energy routines