by brianlane723 » Wed Jun 08, 2005 2:18 pm
Hello. I have written a code using pdsytrd to tridiagonalize a matrix, and it has worked perfectly on one set of machines (i386 running linux), but when I compile and run it on a second set (i686-pc-linux-gnu), I encounter an error that produces NAN's (not-a-number). I've traced the problem down to occur during a call to pdlatrd, when it calls the PBLAS routdine pdgemv. Has anyone seen this happen or know what might be causing it? Thanks!