Kinji Kimura's web page
The number of terms in general formula of the disciriminant(degree 17)
G17=a17*x^17 + a16*x^16 + a15*x^15 + a14*x^14 + a13*x^13 + a12*x^12 + a11*x^11 + a10*x^10 + a9*x^9 + a8*x^8 + a7*x^7 + a6*x^6 + a5*x^5 + a4*x^4 + a3*x^3 + a2*x^2 + a1*x + a0
is 21976689397!

LAPROGNC(Linear Algebra PROGrams in Numerical Computation)

LAPROGCA(Linear Algebra PROGrams in Computer Algebra)

SDC sample data

BENCHMARK:1
Computing the resultant of the polynomial obtained by Prof. Tetsuji Shioda(Rikkyo Univ.)
E6=ans

BENCHMARK:2
General formula of the disciriminant
G15=x^15 + x^14 + a13*x^13 + a12*x^12 + a11*x^11 + a10*x^10 + a9*x^9 + a8*x^8 + a7*x^7 + a6*x^6 + a5*x^5 + a4*x^4 + a3*x^3 + a2*x^2 + a1*x + a0
H15=resultant_x(G15,dG15/dx)
H15=ans

Computing the general discriminant formula of degree 17

Accelerating the Numerical Computation of Positive Roots of Polynomials using Improved Bounds

A new implementation of the orthogonal QD algorithm and its application

A new implementation of the dqds algorithm

LAPACK BUG LIST for numerical linear algebra researchers
(0)DLASQ can not compute the positive value as the smallest singular value in both the matrix and the matrix .
If you want to compute singular values, please use our implementation of the dqds algorithm.
(1)DLASQ doesn't converge in this matrix .
The dlasq5.f might be modified as follows. dlasq5.f
(2)If the size of this matrix is 429, 981, or 985, DBDSDC returns INFO=1.
The input is as follows.
./test 429
./test 981
./test 985
(3)DSTEDC returns INFO=28825 in this matrix and INFO=17419751 in this matrix .
(4)DSTEMR returns INFO=22 in this matrix .