[Octave-bug-tracker] [bug #53700] eigs test failure related to ARPACK generating real NaN rather than complex NaN+1i*NaN

[Dan Sebald]

Follow-up Comment #7, bug #53700 (project octave):

You meant "not be doing the right thing"?  I don't see the 

std::cerr << "PROMOTED REAL TO COMPLEX\n";

printed out.

Anyway, I see what you mentioned, in the eigs-base.cc routine.  Attached is
the changes I've currently made.  Are those F77_LOGICALs what you had in mind?
 I.e., basically do the following in all locations:


@@ -627,6 +627,7 @@
   M b(_b);
   F77_INT n = octave::to_f77_int (m.cols ());
   F77_INT mode = 1;
+  F77_LOGICAL rvec77 = rvec;
   bool have_b = ! b.isempty ();
   bool note3 = false;
   char bmat = 'I';
@@ -839,7 +840,7 @@
   double *d = eig_val.fortran_vec ();
 
   F77_FUNC (dseupd, DSEUPD)
-    (rvec, F77_CONST_CHAR_ARG2 ("A", 1), sel, d, z, n, sigma,
+    (rvec77, F77_CONST_CHAR_ARG2 ("A", 1), sel, d, z, n, sigma,
      F77_CONST_CHAR_ARG2 (&bmat, 1), n,
      F77_CONST_CHAR_ARG2 ((typ.c_str ()), 2), k, tol, presid, p, v, n,
iparam,
      ipntr, workd, workl, lwork, info2 F77_CHAR_ARG_LEN(1)
F77_CHAR_ARG_LEN(1)


Those mods don't change the outcome.

(file #43976)
    _______________________________________________________

Additional Item Attachment:

File name: annotation_and_F77_LOGICAL.diff Size:7 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?53700>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/