... | ... |
@@ -1033,25 +1033,25 @@ |
1033 | 1033 |
} |
1034 | 1034 |
} else { |
1035 | 1035 |
for (int i = 0; i != _arc_num; ++i) { |
1036 | 1036 |
_cap[i] = _upper[i]; |
1037 | 1037 |
} |
1038 | 1038 |
} |
1039 | 1039 |
|
1040 | 1040 |
// Initialize artifical cost |
1041 | 1041 |
Cost ART_COST; |
1042 | 1042 |
if (std::numeric_limits<Cost>::is_exact) { |
1043 | 1043 |
ART_COST = std::numeric_limits<Cost>::max() / 2 + 1; |
1044 | 1044 |
} else { |
1045 |
ART_COST = |
|
1045 |
ART_COST = 0; |
|
1046 | 1046 |
for (int i = 0; i != _arc_num; ++i) { |
1047 | 1047 |
if (_cost[i] > ART_COST) ART_COST = _cost[i]; |
1048 | 1048 |
} |
1049 | 1049 |
ART_COST = (ART_COST + 1) * _node_num; |
1050 | 1050 |
} |
1051 | 1051 |
|
1052 | 1052 |
// Initialize arc maps |
1053 | 1053 |
for (int i = 0; i != _arc_num; ++i) { |
1054 | 1054 |
_flow[i] = 0; |
1055 | 1055 |
_state[i] = STATE_LOWER; |
1056 | 1056 |
} |
1057 | 1057 |
|
... | ... |
@@ -1448,25 +1448,25 @@ |
1448 | 1448 |
if (c != 0) { |
1449 | 1449 |
_flow[i] += c; |
1450 | 1450 |
_supply[_source[i]] += c; |
1451 | 1451 |
_supply[_target[i]] -= c; |
1452 | 1452 |
} |
1453 | 1453 |
} |
1454 | 1454 |
} |
1455 | 1455 |
|
1456 | 1456 |
// Shift potentials to meet the requirements of the GEQ/LEQ type |
1457 | 1457 |
// optimality conditions |
1458 | 1458 |
if (_sum_supply == 0) { |
1459 | 1459 |
if (_stype == GEQ) { |
1460 |
Cost max_pot = std::numeric_limits<Cost>:: |
|
1460 |
Cost max_pot = -std::numeric_limits<Cost>::max(); |
|
1461 | 1461 |
for (int i = 0; i != _node_num; ++i) { |
1462 | 1462 |
if (_pi[i] > max_pot) max_pot = _pi[i]; |
1463 | 1463 |
} |
1464 | 1464 |
if (max_pot > 0) { |
1465 | 1465 |
for (int i = 0; i != _node_num; ++i) |
1466 | 1466 |
_pi[i] -= max_pot; |
1467 | 1467 |
} |
1468 | 1468 |
} else { |
1469 | 1469 |
Cost min_pot = std::numeric_limits<Cost>::max(); |
1470 | 1470 |
for (int i = 0; i != _node_num; ++i) { |
1471 | 1471 |
if (_pi[i] < min_pot) min_pot = _pi[i]; |
1472 | 1472 |
} |
0 comments (0 inline)