@@ -267,12 +267,12 @@ void ClpCholeskyBase::solveKKT(CoinWorkDouble *region1, CoinWorkDouble *region2,
267
267
solve (array);
268
268
int iRow;
269
269
for (iRow = 0 ; iRow < numberTotal; iRow++) {
270
- if (rowsDropped_[iRow] && CoinAbs (array[iRow]) > 1.0e-8 ) {
270
+ if (rowsDropped_[iRow] && std::abs (array[iRow]) > 1.0e-8 ) {
271
271
COIN_DETAIL_PRINT (printf (" row region1 %d dropped %g\n " , iRow, array[iRow]));
272
272
}
273
273
}
274
274
for (; iRow < numberRows_; iRow++) {
275
- if (rowsDropped_[iRow] && CoinAbs (array[iRow]) > 1.0e-8 ) {
275
+ if (rowsDropped_[iRow] && std::abs (array[iRow]) > 1.0e-8 ) {
276
276
COIN_DETAIL_PRINT (printf (" row region2 %d dropped %g\n " , iRow, array[iRow]));
277
277
}
278
278
}
@@ -2813,13 +2813,13 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
2813
2813
}
2814
2814
}
2815
2815
diagonal_[iRow] = work[iRow];
2816
- largest2 = std::max (largest2, CoinAbs (work[iRow]));
2816
+ largest2 = std::max (largest2, std::abs (work[iRow]));
2817
2817
work[iRow] = 0.0 ;
2818
2818
int j;
2819
2819
for (j = 0 ; j < number; j++) {
2820
2820
int jRow = which[j];
2821
2821
put[j] = work[jRow];
2822
- largest2 = std::max (largest2, CoinAbs (work[jRow]));
2822
+ largest2 = std::max (largest2, std::abs (work[jRow]));
2823
2823
work[jRow] = 0.0 ;
2824
2824
}
2825
2825
} else {
@@ -3009,9 +3009,9 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3009
3009
if (iOriginalRow < numberColumns) {
3010
3010
iColumn = iOriginalRow;
3011
3011
CoinWorkDouble value = diagonal[iColumn];
3012
- if (CoinAbs (value) > 1.0e-100 ) {
3012
+ if (std::abs (value) > 1.0e-100 ) {
3013
3013
value = 1.0 / value;
3014
- largest = std::max (largest, CoinAbs (value));
3014
+ largest = std::max (largest, std::abs (value));
3015
3015
diagonal_[iRow] = -value;
3016
3016
CoinBigIndex start = columnStart[iColumn];
3017
3017
CoinBigIndex end = columnStart[iColumn] + columnLength[iColumn];
@@ -3020,17 +3020,17 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3020
3020
kRow = permuteInverse_[kRow ];
3021
3021
if (kRow > iRow) {
3022
3022
work[kRow ] = element[j];
3023
- largest = std::max (largest, CoinAbs (element[j]));
3023
+ largest = std::max (largest, std::abs (element[j]));
3024
3024
}
3025
3025
}
3026
3026
} else {
3027
3027
diagonal_[iRow] = -value;
3028
3028
}
3029
3029
} else if (iOriginalRow < numberTotal) {
3030
3030
CoinWorkDouble value = diagonal[iOriginalRow];
3031
- if (CoinAbs (value) > 1.0e-100 ) {
3031
+ if (std::abs (value) > 1.0e-100 ) {
3032
3032
value = 1.0 / value;
3033
- largest = std::max (largest, CoinAbs (value));
3033
+ largest = std::max (largest, std::abs (value));
3034
3034
} else {
3035
3035
value = 1.0e100 ;
3036
3036
}
@@ -3048,7 +3048,7 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3048
3048
int jNewRow = permuteInverse_[jRow];
3049
3049
if (jNewRow > iRow) {
3050
3050
work[jNewRow] = elementByRow[j];
3051
- largest = std::max (largest, CoinAbs (elementByRow[j]));
3051
+ largest = std::max (largest, std::abs (elementByRow[j]));
3052
3052
}
3053
3053
}
3054
3054
// slack - should it be permute
@@ -3078,7 +3078,7 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3078
3078
CoinBigIndex j;
3079
3079
iColumn = iOriginalRow;
3080
3080
CoinWorkDouble value = diagonal[iColumn];
3081
- if (CoinAbs (value) > 1.0e-100 ) {
3081
+ if (std::abs (value) > 1.0e-100 ) {
3082
3082
value = 1.0 / value;
3083
3083
for (j = columnQuadraticStart[iColumn];
3084
3084
j < columnQuadraticStart[iColumn] + columnQuadraticLength[iColumn]; j++) {
@@ -3090,7 +3090,7 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3090
3090
value += quadraticElement[j];
3091
3091
}
3092
3092
}
3093
- largest = std::max (largest, CoinAbs (value));
3093
+ largest = std::max (largest, std::abs (value));
3094
3094
diagonal_[iRow] = -value;
3095
3095
CoinBigIndex start = columnStart[iColumn];
3096
3096
CoinBigIndex end = columnStart[iColumn] + columnLength[iColumn];
@@ -3099,17 +3099,17 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3099
3099
kRow = permuteInverse_[kRow ];
3100
3100
if (kRow > iRow) {
3101
3101
work[kRow ] = element[j];
3102
- largest = std::max (largest, CoinAbs (element[j]));
3102
+ largest = std::max (largest, std::abs (element[j]));
3103
3103
}
3104
3104
}
3105
3105
} else {
3106
3106
diagonal_[iRow] = -value;
3107
3107
}
3108
3108
} else if (iOriginalRow < numberTotal) {
3109
3109
CoinWorkDouble value = diagonal[iOriginalRow];
3110
- if (CoinAbs (value) > 1.0e-100 ) {
3110
+ if (std::abs (value) > 1.0e-100 ) {
3111
3111
value = 1.0 / value;
3112
- largest = std::max (largest, CoinAbs (value));
3112
+ largest = std::max (largest, std::abs (value));
3113
3113
} else {
3114
3114
value = 1.0e100 ;
3115
3115
}
@@ -3127,7 +3127,7 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3127
3127
int jNewRow = permuteInverse_[jRow];
3128
3128
if (jNewRow > iRow) {
3129
3129
work[jNewRow] = elementByRow[j];
3130
- largest = std::max (largest, CoinAbs (elementByRow[j]));
3130
+ largest = std::max (largest, std::abs (elementByRow[j]));
3131
3131
}
3132
3132
}
3133
3133
// slack - should it be permute
@@ -3152,17 +3152,17 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3152
3152
longDouble *put = sparseFactor_ + choleskyStart_[iColumn];
3153
3153
CoinBigIndex *which = choleskyRow_ + indexStart_[iColumn];
3154
3154
CoinWorkDouble value = diagonal[iColumn];
3155
- if (CoinAbs (value) > 1.0e-100 ) {
3155
+ if (std::abs (value) > 1.0e-100 ) {
3156
3156
value = 1.0 / value;
3157
- largest = std::max (largest, CoinAbs (value));
3157
+ largest = std::max (largest, std::abs (value));
3158
3158
diagonal_[iColumn] = -value;
3159
3159
CoinBigIndex start = columnStart[iColumn];
3160
3160
CoinBigIndex end = columnStart[iColumn] + columnLength[iColumn];
3161
3161
for (CoinBigIndex j = start; j < end; j++) {
3162
3162
// choleskyRow_[numberElements]=row[j]+numberTotal;
3163
3163
// sparseFactor_[numberElements++]=element[j];
3164
3164
work[row[j] + numberTotal] = element[j];
3165
- largest = std::max (largest, CoinAbs (element[j]));
3165
+ largest = std::max (largest, std::abs (element[j]));
3166
3166
}
3167
3167
} else {
3168
3168
diagonal_[iColumn] = -value;
@@ -3187,7 +3187,7 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3187
3187
int number = choleskyStart_[iColumn + 1 ] - choleskyStart_[iColumn];
3188
3188
CoinWorkDouble value = diagonal[iColumn];
3189
3189
CoinBigIndex j;
3190
- if (CoinAbs (value) > 1.0e-100 ) {
3190
+ if (std::abs (value) > 1.0e-100 ) {
3191
3191
value = 1.0 / value;
3192
3192
for (j = columnQuadraticStart[iColumn];
3193
3193
j < columnQuadraticStart[iColumn] + columnQuadraticLength[iColumn]; j++) {
@@ -3198,13 +3198,13 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3198
3198
value += quadraticElement[j];
3199
3199
}
3200
3200
}
3201
- largest = std::max (largest, CoinAbs (value));
3201
+ largest = std::max (largest, std::abs (value));
3202
3202
diagonal_[iColumn] = -value;
3203
3203
CoinBigIndex start = columnStart[iColumn];
3204
3204
CoinBigIndex end = columnStart[iColumn] + columnLength[iColumn];
3205
3205
for (j = start; j < end; j++) {
3206
3206
work[row[j] + numberTotal] = element[j];
3207
- largest = std::max (largest, CoinAbs (element[j]));
3207
+ largest = std::max (largest, std::abs (element[j]));
3208
3208
}
3209
3209
for (j = 0 ; j < number; j++) {
3210
3210
int jRow = which[j];
@@ -3226,9 +3226,9 @@ int ClpCholeskyBase::factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
3226
3226
longDouble *put = sparseFactor_ + choleskyStart_[iColumn];
3227
3227
CoinBigIndex *which = choleskyRow_ + indexStart_[iColumn];
3228
3228
CoinWorkDouble value = diagonal[iColumn];
3229
- if (CoinAbs (value) > 1.0e-100 ) {
3229
+ if (std::abs (value) > 1.0e-100 ) {
3230
3230
value = 1.0 / value;
3231
- largest = std::max (largest, CoinAbs (value));
3231
+ largest = std::max (largest, std::abs (value));
3232
3232
} else {
3233
3233
value = 1.0e100 ;
3234
3234
}
@@ -3839,7 +3839,7 @@ void ClpCholeskyBase::solve(CoinWorkDouble *region, int type)
3839
3839
int j;
3840
3840
double largestO = 0.0 ;
3841
3841
for (i = 0 ; i < numberRows_; i++) {
3842
- largestO = std::max (largestO, CoinAbs (regionX[i]));
3842
+ largestO = std::max (largestO, std::abs (regionX[i]));
3843
3843
}
3844
3844
for (i = 0 ; i < numberRows_; i++) {
3845
3845
int iRow = permute_[i];
@@ -3880,8 +3880,8 @@ void ClpCholeskyBase::solve(CoinWorkDouble *region, int type)
3880
3880
double largest = 0.0 ;
3881
3881
double largestV = 0.0 ;
3882
3882
for (i = 0 ; i < numberRows_; i++) {
3883
- largest = std::max (largest, CoinAbs (region[i] - regionX[i]));
3884
- largestV = std::max (largestV, CoinAbs (region[i]));
3883
+ largest = std::max (largest, std::abs (region[i] - regionX[i]));
3884
+ largestV = std::max (largestV, std::abs (region[i]));
3885
3885
}
3886
3886
printf (" largest difference %g, largest %g, largest original %g\n " ,
3887
3887
largest, largestV, largestO);
0 commit comments