[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Octave-bug-tracker] [bug #53700] eigs test failure related to ARPACK ge

From: Dan Sebald
Subject: [Octave-bug-tracker] [bug #53700] eigs test failure related to ARPACK generating real NaN rather than complex NaN+1i*NaN
Date: Thu, 19 Apr 2018 16:49:07 -0400 (EDT)
User-agent: Mozilla/5.0 (X11; Ubuntu; Linux x86_64; rv:59.0) Gecko/20100101 Firefox/59.0

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,
      ipntr, workd, workl, lwork, info2 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:


  Message sent via/by Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]