python - Eigenvalues NaN and inf -


suppose have system ax = nbx , b known martrices, x coefficient matrix.

i solving using chebyshev polynomials.

bc's u(-1)=0=u(1)

i imposing bc's first , last rows of matrices , b.

e=solve(a,b)  e[1]=0  e[-1]=0  x=solve(a,e) 

what wrong this?

the question seem asking: how come generalized eigenvalue problem has eigenvalues inf , nan?

your generalized eigenvalue problem singular , has eigenvalues lambda=alpha/beta such (alpha=0, beta=0) , (alpha!=0, beta=0). since eigvals reports eigenvalues, 0/0=nan or x/0=inf, correspondingly.

http://www.netlib.org/lapack/lug/node35.html

if problem shouldn't have such eigenvalues, it's there error in construction of matrices.


Comments

Popular posts from this blog

python - Subclassed QStyledItemDelegate ignores Stylesheet -

java - HttpClient 3.1 Connection pooling vs HttpClient 4.3.2 -

SQL: Divide the sum of values in one table with the count of rows in another -