some changes to calculate vbd, not yet fully implemented