Conventions regarding the C-version of the function dtrnlsp_solve

Conventions regarding the C-version of the function dtrnlsp_solve

Hi there.

I am reading the documentation of the C-version of the function dtrnlsp_solve and I would like some help since I find some unconsistencies with the notation.

Here are my questions:

1] fvec[i] is really defined as ( y[i] - f[i](x) ) (as indicated in page 2649 of the manual) or is it just f[i](x) ?

2] Let f : R^N ---> R^M, 0 <= m <= M-1and 0 <= n <= N-1. What is the Jacobianfjac supposed to contain? Is the stuff stored in column major order or row major order? Which one of the following conventions is the correct one?

- fjac[ m + n*M ] = frac{ partial f[m](x) }{ partial x[n] }

- fjac[n + m*N ] = frac{ partial f[m](x) }{ partial x[n] }

- fjac[ m + n*M ] = frac{ y[m] - partial f[m](x) }{ partial x[n] }

- fjac[n + m*N ] = frac{ y[m] - partial f[m](x) }{ partial x[n] }

Thanks a lot for any insight that could avoid me a trial and error boring session ;)

Marco

1 post / 0 new
For more complete information about compiler optimizations, see our Optimization Notice.