[SciPy-user] non-linear multi-variate optimization

Sebastian Walter sebastian.walter@gmail....
Wed Jul 15 02:32:55 CDT 2009

The gradient g is only defined for functions f: R^N --> R and is
simply an array with shape (N,)

what you sketched in your post is the Jacobian J of a function f: R^N --> R^M
Typically, the Jacobian J is defined to have the shape (M,N), but
there are exceptions.

hope that helps a little

On Wed, Jul 15, 2009 at 2:00 AM, Chris Colbert<sccolbert@gmail.com> wrote:
> The routines for non-linear optimization in scipy.optimize take an
> argument for a function that computes the gradient.
> What should be the format of return value of this function? I am
> assuming that its the gradient of the functions with repect to the
> independent variables in row vector format.
> for example say we have:
> f(x,y,z; a1, a2, a3) where a1, a2, and a3 are the independent variables.
> Should the gradient of N x,y,z points then be of the form:
> df/da  =   [[df(X0)/da1, df(X0)/da2, df(X0)/da3],
>               [df(X1)/da1, df(X1)/da2, df(X1)/da3],
>               ....
>               [df(Xn)/da1, df(Xn)/da2, df(Xn)/da3]]
> where Xn is the set of (xn, yn, zn) ?
> cheers!
> Chris
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user

More information about the SciPy-user mailing list