 Timestamp:
 Oct 24, 2005 11:26:04 AM (14 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/dev/Algorithm/LinearSolvers/IpPardisoSolverInterface.cpp
r547 r548 121 121 IPARM_[9] = 12; 122 122 IPARM_[10] = 1; 123 // Matching information: IPARM_[12] = 1 seems ok, but results in a 123 // Matching information: IPARM_[12] = 1 seems ok, but results in a 124 124 // large number of pivot perturbation 125 125 // Matching information: IPARM_[12] = 2 robust, but more expensive method … … 328 328 // call, even if the current inertia estimate is correct 329 329 // OLAF MICHAEL : Maybe we need to discuss this 330 331 332 333 334 330 have_symbolic_factorization_ = false; 331 // We assume now that if there was just a symbolic 332 // factorization and we still have perturbed pivots, that the 333 // system is actually singular 334 if (just_performed_symbolic_factorization) { 335 335 IpData().Append_info_string("Ps"); 336 336 //break; // Declaring this system singular doesn't seem to be 337 338 339 340 341 342 // If we there were perturbed pivots, and the inertia of the343 // system is correct, and we haven't redone the symblic344 // factorization during this call yet, trigger a new345 // factorization after a symblic factorization346 done = (just_performed_symbolic_factorization  !check_NegEVals 347 numberOfNegEVals==negevals_);348 349 337 // working well??? 338 return SYMSOLVER_SINGULAR; 339 } 340 IpData().Append_info_string("Pp"); 341 /* 342 // If we there were perturbed pivots, and the inertia of the 343 // system is correct, and we haven't redone the symblic 344 // factorization during this call yet, trigger a new 345 // factorization after a symblic factorization 346 done = (just_performed_symbolic_factorization  !check_NegEVals  347 numberOfNegEVals==negevals_); 348 */ 349 done = false; 350 350 } 351 351 else {
Note: See TracChangeset
for help on using the changeset viewer.