add/assign problem, part 2.

[David Mebane <david dot mebane at mse dot gatech dot edu>]

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