Upon receiving some advice, I have attached here the entire piece of code I'm trying to run. It's a Newton solver; LSMbdnew is the generic solver, and bulkfunc contains the function (f) and Jacobian (jac). The whole thing may look a little ugly to the experienced programmers on this list -- sorry. Anyway, this will compile and run, but it won't do what it's supposed to. I found out with gdb that the reason for this is that the function LSMbdnew::norm(double *vec1, long num) is not executing the add/assign properly. For some reason, it sets the collecting variable _ret_ to zero after one or two iterations. Why? I am stumped. Thanks, -David
Attachment:
bulkfunc.h
Description: Binary data
Attachment:
LSMbdnew.h
Description: Binary data
Attachment:
testnew.cpp
Description: Binary data
