在下面的代码中,我尝试使用CPLEX宏ILOBRANCHCALLBACK2
创建分支回调,在此我选择固定某些变量边界的方法。但是在进行修复时,读取一些空间时出现访问冲突错误。有人可以告诉我正确的做法吗?我是否必须对所有变量都进行make分支,否则branchcallback会自己完成?
ILOBRANCHCALLBACK2(simetria, const Instance&, instance, const IloNumVarMatrix&, X) {
vector<coordenates> F_zeros;
vector<coordenates> F_ones;
vector<int> a;
a.resize(instance.getNbClients());
a[0] = 0;
coordenates first;
first.x = 0;
first.y = 0;
F_ones.push_back(first);
// Je verifie s'il y a deja des variables fixées
for(int i = 1; i < instance.getNbClients() - 1; ++i) {
for(int j = 1; j < instance.getNbVehicles() - 1; ++j) {
if(getUB(X[i][j]) <= 0.001) {
coordenates obj;
obj.x = i;
obj.y = j;
F_zeros.push_back(obj);
}
if(getLB(X[i][j]) >= 0.99) {
coordenates obj;
obj.x = i;
obj.y = j;
F_ones.push_back(obj);
}
}
}
// Je vais construir mon a.
for(int k = 1; k < instance.getNbClients() - 1; k++) {
// Zero setting
if(a[k - 1] == instance.getNbVehicles() - 1 || findInVector(F_zeros, k, a[k - 1] + 1)) {
a[k] = a[k - 1];
coordenates obj;
obj.x = k;
obj.y = a[k - 1];
} else {
a[k] = a[k - 1] + 1;
coordenates obj;
obj.x = k;
obj.y = a[k - 1] + 1;
}
}
// Je mets a jour le F_zero a partir de a
for(int k = 1; k < instance.getNbClients() - 1; k++) {
for(int p = 1; p < instance.getNbVehicles() - 1; p++) {
if(p > a[k]) {
coordenates obj;
obj.x = k;
obj.y = p;
F_zeros.push_back(obj);
}
}
}
// je fixe tout ce qui est en F_zero en zero (Changer les bornnes)
for(int i = 0; i < F_zeros.size(); i++) {
int a = F_zeros[i].x;
int b = F_zeros[i].y;
X[a][b].setUB(0);
// X[F_zeros[i].x][F_zeros[i].y].setBounds(0, 0);
// makeBranch(X[F_zeros[i].x][F_zeros[i].y], 0, IloCplex::BranchUp, getBestObjValue());
// makeBranch(X[F_zeros[i].x][F_zeros[i].y], 0, IloCplex::BranchDown, getBestObjValue());
}
答案 0 :(得分:1)
求解过程中,您无法修改模型(例如,请参见文档here中的 Note )。如您所知,这可能会导致段错误。因此,您不能在回调中使用DATABASE_PARTIAL
方法。
您可能已经知道(基于已注释掉的代码),可以使用makeBranch方法,
指示调用IloCplex对象如何从中创建子节点 通过为一组节点指定新的更严格的边界来确定当前节点 变量。