computeTabulateBonds works, bug computing forces fixed, bondList is created, but not implemented yet.
