[SciPy-User] Vectorized newton function

Robert Kern robert.kern@gmail....
Wed Apr 14 21:26:51 CDT 2010

On Wed, Apr 14, 2010 at 19:33, Charles R Harris
<charlesr.harris@gmail.com> wrote:
> On Wed, Apr 14, 2010 at 6:15 PM, Gökhan Sever <gokhansever@gmail.com> wrote:
>> Thanks Anne for suggesting "fsolve" and rest of the useful information.
>> Below is how my function is read:
>> def myfunc(myarray):
>>        term1 = compute1_array1
>>        term2 = compute3_array2
>>        term3 = compute3_array3
>>        return term1+term2+term3
>> result_arr = fsolve(myfunc, myarray)
>> IDL handles univariate and multivariate cases using one fun newton
>> function. I get almost equivalent results using fsolve to IDL's newton.
> The fsolve function is essentially multidimensional newton, IIRC.

Yes, but it's an N**2 waste when the equations aren't coupled.

Robert Kern

"I have come to believe that the whole world is an enigma, a harmless
enigma that is made terrible by our own mad attempt to interpret it as
though it had an underlying truth."
  -- Umberto Eco

More information about the SciPy-User mailing list