problem with elmhes


allynm
04-05-2013, 05:16 PM
Hello all,
I am quite new to nr so please forgive any obviously boneheaded question.

Here it is.

I am trying to use elmhes routine in Chapter 8. I'm not able to generate a Hessenberg 4*4 matrix from a simple asymmetric input. The numbers are all low integers and their range is very limited. I've checked the coding (in C, using the 1992 edition of NR) a number of times and can't find an error. I check the input matrix by having it printed out by the caller, and it appears that the input is properly being read into an array declared as float ** aa per instructions in Chapter 1. I could only locate one other forum contributor using elmhes and he? appears to have had no problem with the function.

Any suggestions would be appreciated, of course.

Thanks,
Mark Allyn