|
Go to the documentation of this file.
19 #define DEVEX_REFINETOL 2.0
33 #ifdef ENABLE_CONSISTENCY_CHECKS
120 if (fTesti < -feastol)
124 price. val = fTesti * fTesti / cpen[idx];
134 for( int i = 0; i < nsorted; ++i )
165 if ( retid < 0 && !refined )
172 assert(retid < thesolver->dim());
187 for (; start < end; start += incr)
189 if (fTest[start] < -feastol)
191 x = fTest[start] * fTest[start] / cpen[start];
219 if (fTesti < -feastol)
222 x = fTesti * fTesti / coPeni;
258 if( fTesti < -feastol )
261 x = fTesti * fTesti / coPeni;
294 assert(fTesti < -feastol);
296 x = fTesti * fTesti / coPeni;
322 Real rhov_1 = 1 / rhoVec[n];
328 MSG_ERROR( std::cerr << "WDEVEX01: rhoVec = "
329 << rhoVec[n] << " with smaller absolute value than theeps = " << theeps << std::endl; )
336 for (i = len - 1; i >= 0; --i)
339 x = rhoVec[j] * rhoVec[j] * beta_q;
368 price. val = x * x / cpen[idx];
383 for( int i = 0; i < nsorted; ++i )
418 price. val = x * x / pen[idx];
433 for( int i = 0; i < nsorted; ++i )
456 if( !enterId. isValid() && !refined )
521 coTesti = cTest[idx];
522 if( coTesti < -feastol )
525 x = coTesti * coTesti / coPeni;
556 coTesti = cTest[idx];
557 if( coTesti < -feastol )
560 x = coTesti * coTesti / coPeni;
604 if( testi < -feastol )
607 x = testi * testi / peni;
638 if( testi < -feastol )
641 x = testi * testi / peni;
683 coTesti = cTest[idx];
684 if (coTesti < -feastol)
687 x = coTesti * coTesti / coPeni;
723 if (testi < -feastol)
726 x = testi * testi / peni;
757 for (; start < end; start += incr)
759 if (cTest[start] < -feastol)
761 x = cTest[start] * cTest[start] / cpen[start];
787 for (; start < end; start += incr)
789 if (test[start] < -feastol)
791 x = test[start] * test[start] / pen[start];
814 if (n >= 0 && n < thesolver->dim())
826 xi_p = xi_p * xi_p * last;
828 for (j = coPidx. size() - 1; j >= 0; --j)
831 coWeights[i] += xi_p * coPvec[i] * coPvec[i];
839 for (j = pIdx. size() - 1; j >= 0; --j)
842 weights[i] += xi_p * pVec[i] * pVec[i];
|