Hanging on getting eigenvalues

nalgebra
1

Hi,

I’m trying to find the roots of a polynomial, and nalgebra has overall worked great.

Recently, I’ve found that my calculations hang on certain inputs, for which numpy immediately returns a value.

You can see my code for finding roots at https://github.com/hwchen/roots.

A pathological input is in this test case (sorry, can’t push from git because I’m on a plane):

   #[test]
    fn hanging() {
        let input = &[
            -770.9700059294701,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            434.3959909081459,
            0.0,
        ];

        let input: Vec<f64> = input.into_iter().rev().cloned().collect();
        dbg!(&input);
        let rs = roots(input.as_slice());
    }

Thanks for any help.