A hybrid solution to parallel calculation of augmented join trees of scalar fields in any dimension