//============================================================== // PredictCollision //============================================================== void box::PredictCollision(int i, int j, vector<DIM, int> pboffset, double& ctime, int& cpartner, vector<DIM, int>& cpartnerpboffset) { double ctimej; if (i!=j) { ctimej = CalculateCollision(i,j,pboffset.Double())+gtime; if (ctimej < gtime) std::cout << "error in find collision ctimej < 0" << std::endl; if ((ctimej < ctime)&&(ctimej < s[j].nextevent.time)) { ctime = ctimej; cpartner = j; cpartnerpboffset = pboffset; } } }