如何在宏BranchCallback中修复边界/更改变量边界? C ++

时间:2019-03-13 17:07:35

标签: c++ callback branch cplex

在下面的代码中,我尝试使用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());
    }

1 个答案:

答案 0 :(得分:1)

求解过程中,您无法修改模型(例如,请参见文档here中的 Note )。如您所知,这可能会导致段错误。因此,您不能在回调中使用DATABASE_PARTIAL方法。

您可能已经知道(基于已注释掉的代码),可以使用makeBranch方法,

  

指示调用IloCplex对象如何从中创建子节点   通过为一组节点指定新的更严格的边界来确定当前节点   变量。