The energy matrix using determinantal product states applied to Ho:YAG