An efficient Julia framework for hierarchical equations of motion in open quantum systems