caliarim

In debug.h there are several integers which can be used to get information from the subroutines. For instance, if msaupd is 1 in your main program, then pdsaupd.f will set...

I never used pArpack. But I would try the following: for each variable in debug.h, starting from mgetv0, I would set it explicitely to 0 in your main file.

I confirm a similar bad behavior in plain fortran77 arpack 3.3.0 (see attachment). [QLR.zip](https://github.com/opencollab/arpack-ng/files/633990/QLR.zip)

Hi. If you check the orthogonality of the vectors v after the second call of dnaitr in dnaup2, you will see that orthogonality is lost (some vectors have scalar product...

My previous attempt does not work in the generalized case. Still, the problem is the orthogonalization. There are some possibile solutions: 1) a second orthogonalization is made in some cases,...

any comment on this?

Both the problem described here and the problem related to "Force the residual vector to be in the range of OP" depend on a bad choice of the initial random...