Error With Arpack Routine Dnaupd
Support Support Newsreader MathWorks Search MathWorks.com MathWorks Newsreader Support MATLAB Newsgroup MATLAB Central Community Home MATLAB Answers File Exchange Cody Blogs Newsreader Link Exchange ThingSpeak Anniversary Home Post A New Message Advanced Search Help MATLAB matlab error with arpack routine dnaupd info Central Community Home MATLAB Answers File Exchange Cody Blogs Newsreader Link Exchange ThingSpeak error with arpack routine dneupd Anniversary Home Post A New Message Advanced Search Help Trial software Error with ARPACK routine when using eigs Subject: Error with matlab eigs ARPACK routine when using eigs From: Aishwarya Date: 5 Oct, 2009 11:31:10 Message: 1 of 3 Reply to this message Add author to My Watch List View original format Flag as spam One of the functions I use uses eigs. And it was running fine for a week till suddenly it started throwing the following error : Error using ==> eigs at 231 Error with ARPACK routine dnaupd: info = -8 I hope someone can explain what the problem is and tell me to how to overcome it. Subject: Error with ARPACK routine when using eigs From: Bruno Luong Bruno Luong (view profile) 5225 posts Date: 5 Oct, 2009 12:04:02 Message: 2 of 3 Reply to this message Add author to My Watch List View original format Flag as spam Aishwarya
Support Answers MathWorks Search MathWorks.com MathWorks Answers Support MATLAB Answers™ MATLAB Central Community Home MATLAB Answers File Exchange Cody Blogs Newsreader Link Exchange ThingSpeak Anniversary Home Ask Answer Browse More Contributors Recent Activity Flagged Content Flagged as Spam Help MATLAB Central Community Home MATLAB Answers File Exchange Cody Blogs Newsreader Link Exchange ThingSpeak Anniversary Home Ask Answer Browse More Contributors Recent Activity Flagged Content Flagged https://www.mathworks.com/matlabcentral/newsreader/view_thread/262324 as Spam Help Trial software Marie-Andree Coulombe (view profile) 1 question 0 answers 0 accepted answers Reputation: 0 Vote0 Error using eigs: Error with ARPACK routine dnaupd Asked by Marie-Andree Coulombe Marie-Andree Coulombe (view profile) 1 question 0 answers 0 accepted answers Reputation: 0 on 13 Jan 2014 8 views (last https://www.mathworks.com/matlabcentral/answers/112264-error-using-eigs-error-with-arpack-routine-dnaupd 30 days) 8 views (last 30 days) Hi, I'm trying to use a script using eigs function: opts.disp=0; [v d]=eigs(A,B,k,'lm',opts); [Note: the original function used `la`, but my matrix doesn't seem symetrical]Here's the error message I received: ??? Error using ==> eigs at 255 Error with ARPACK routine dnaupd: info = -9999 [Note: it's the same error message whereas i`m using `lm`, `lr`or `li`(I'm a real newbie..)]ThanksMarie-AndrĂ©e 0 Comments Show all comments Tags eigslrarkpack Products MATLAB Related Content 0 Answers Log In to answer or comment on this question. Related Content Join the 15-year community celebration. Play games and win prizes! Learn more MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi Learn more Discover what MATLAB® can do for your career. Opportunities for recent engineering grads. Apply Today MATLAB Academy New to MATLAB? Learn MATLAB today! An Error Occurred Unable to complete the action because of chan
Support Answers MathWorks Search MathWorks.com MathWorks Answers Support MATLAB Answers™ MATLAB Central Community Home MATLAB Answers File Exchange Cody Blogs Newsreader Link Exchange ThingSpeak Anniversary Home Ask Answer Browse More Contributors Recent Activity Flagged Content Flagged https://www.mathworks.com/matlabcentral/answers/62209-eigs-and-its-behaviour as Spam Help MATLAB Central Community Home MATLAB Answers File Exchange Cody Blogs Newsreader Link http://www.caam.rice.edu/software/ARPACK/UG/node137.html Exchange ThingSpeak Anniversary Home Ask Answer Browse More Contributors Recent Activity Flagged Content Flagged as Spam Help Trial software Karim (view profile) 1 question 0 answers 0 accepted answers Reputation: 0 Vote0 eigs and its behaviour Asked by Karim Karim (view profile) 1 question 0 answers 0 accepted answers Reputation: 0 on 4 Feb 2013 28 views error with (last 30 days) 28 views (last 30 days) Hello everyone at Mathworks,I'm using eigs in some algorithms that deal with real world data sets. These data result in square symmetric real matrices which might contain lots and lots of zeros; i.e. really sparse.I some cases, due to different values of the parameters these matrices when passed to eigs yield the following error message:>> >> ...Error using eigs/processEUPDinfo (line 1329) >> Error error with arpack with ARPACK routine dneupd: >> dnaupd did not find any eigenvalues to sufficient accuracy. >> >> Error in eigs (line 357) >>flag = processEUPDinfo(nargout<3); >> Note that these matrices do not contain NaNs or inf values. Just +ve real values and zeros.The problem is the way eigs interrupts the execution of the experiment running. As you can see, if it does not converge, it gives an error and stops all the experiment, loosing by that hours and hours of simulations. Even though I used the flag output variable, still the same behaviour.As a programming trick, I tried eigs in a try/catch block which seems to work fine and avoids interrupting the execution of the hall experiment. However, I was wondering if there is a better solution for this problem. For instance, I would like to know, out of the k eigenvalue-eigenvector pairs, which one did not converge? retrieve the solution for the first k-1 eigenvalue-eigenvectors pairs (which converged), and hopefully make more use of the 'flag' variable without the need for a try/catch block and without stopping the hall experiment.Thank you all.Best,-Karim 0 Comments Show all comments Tags eigs Products No products are associated with this question. Related Content 1 Answer Shashank Prasanna (view profile) 0 questions 677 answers 269 accepted answers Reput
a few eigenpairs c of a linear operator "OP" with respect to a semi-inner product defined by c a symmetric positive semi-definite real matrix B. B may be the identity c matrix. NOTE: If the linear operator "OP" is real and symmetric c with respect to the real positive semi-definite symmetric matrix B, c i.e. B*OP = (OP')*B, then subroutine ssaupd should be used instead. c c The computed approximate eigenvalues are called Ritz values and c the corresponding approximate eigenvectors are called Ritz vectors. c c dnaupd is usually called iteratively to solve one of the c following problems: c c Mode 1: A*x = lambda*x. c ===> OP = A and B = I. c c Mode 2: A*x = lambda*M*x, M symmetric positive definite c ===> OP = inv[M]*A and B = M. c ===> (If M can be factored see remark 3 below) c c Mode 3: A*x = lambda*M*x, M symmetric semi-definite c ===> OP = Real_Part{ inv[A - sigma*M]*M } and B = M. c ===> shift-and-invert mode (in real arithmetic) c If OP*x = amu*x, then c amu = 1/2 * [ 1/(lambda-sigma) + 1/(lambda-conjg(sigma)) ]. c Note: If sigma is real, i.e. imaginary part of sigma is zero; c Real_Part{ inv[A - sigma*M]*M } == inv[A - sigma*M]*M c amu == 1/(lambda-sigma). c c Mode 4: A*x = lambda*M*x, M symmetric semi-definite c ===> OP = Imaginary_Part{ inv[A - sigma*M]*M } and B = M. c ===> shift-and-invert mode (in real arithmetic) c If OP*x = amu*x, then c amu = 1/2i * [ 1/(lambda-sigma) - 1/(lambda-conjg(sigma)) ]. c c Both mode 3 and 4 give the same enhancement to eigenvalues close to c the (complex) shift sigma. However, as lambda goes to infinity, c the operator OP in mode 4 dampens the eigenvalues more strongly than c does OP defined in mode 3. c c NOTE: The action of w <- inv[A - sigma*M]*v or w <- inv[M]*v c should be accomplished either by a direct method c using a sparse matrix factorization and solving c c [A - sigma*M]*w = v or M*w = v, c c or through an iterative method for solving these c s