337 |
|
Real fC0 = getCoefficient(0); |
338 |
|
roots.push_back( -fC0 / fC1); |
339 |
|
return roots; |
340 |
– |
break; |
340 |
|
} |
341 |
|
case 2: { |
342 |
|
Real fC2 = getCoefficient(2); |
361 |
|
roots.push_back(-fTmp * fC1); // 1 real root |
362 |
|
} |
363 |
|
} |
364 |
< |
return roots; |
366 |
< |
break; |
367 |
< |
|
364 |
> |
return roots; |
365 |
|
case 3: { |
366 |
|
Real fC3 = getCoefficient(3); |
367 |
|
Real fC2 = getCoefficient(2); |
428 |
|
} |
429 |
|
} |
430 |
|
return roots; |
431 |
< |
break; |
431 |
> |
|
432 |
|
case 4: { |
433 |
|
Real fC4 = getCoefficient(4); |
434 |
|
Real fC3 = getCoefficient(3); |
514 |
|
} |
515 |
|
} |
516 |
|
return roots; |
517 |
< |
break; |
517 |
> |
|
518 |
|
default: { |
519 |
|
DynamicRectMatrix<Real> companion = CreateCompanion(); |
520 |
|
JAMA::Eigenvalue<Real> eig(companion); |
528 |
|
} |
529 |
|
} |
530 |
|
return roots; |
531 |
< |
break; |
532 |
< |
} |
536 |
< |
|
537 |
< |
return roots; // should be empty if you got here |
531 |
> |
|
532 |
> |
} |
533 |
|
} |
534 |
< |
|
534 |
> |
|
535 |
|
private: |
536 |
|
|
537 |
|
PolynomialPairMap polyPairMap_; |