[SciPy-User] Vectorized newton function
Robert Kern
robert.kern at gmail.com
Wed Apr 14 22:26:51 EDT 2010
On Wed, Apr 14, 2010 at 19:33, Charles R Harris
<charlesr.harris at gmail.com> wrote:
>
> On Wed, Apr 14, 2010 at 6:15 PM, Gökhan Sever <gokhansever at 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