Optimizing FindMinimum for speed using intermediate calculations.

*To*: mathgroup at christensen.cybernetics.net*Subject*: [mg1175] Optimizing FindMinimum for speed using intermediate calculations.*From*: Bill Harbaugh <harbaugh at students.wisc.edu>*Date*: Sat, 20 May 1995 04:23:30 -0400*Organization*: University of Wisconsin, Madison

Hello I am trying to maximize a likelihood function, using FindMinimum. This involves finding the parameters (x,y in my simplified example) that max (min, in my example) the sum of the values of a function applied to some data (ddist, in my example). My actual problem is very long and my question is about how to optimize my mma code for speed. This function involves some intermediate results (x1,y1 in my example) which vary across the x,y parameters but are constant across the d's. (In my actual problem, these intermediate results are found using NSolve, so I am specifying the functions using := rather than =). I am trying to only calculate these once for each iteration on x,y, and then apply the results to each of the d's in the list ddist. The code below shows how I do this without trying to optimize. ddist={1,2,3,4,5,6,7,8,9,10} x1:=x+1 y1:=y+2 l[d_]:=(x1*d)^2+(y1*d)^2 f=Apply[Plus, Map[l,ddist] ] FindMinimum[f,{x,{-1,1}},{y,{-1,1}}] This works, but it seems the following should be faster. Clear[l]; Clear[f]; l[d_]:=(x1*d)^2+(y1*d)^2 f={x1=x+1,y1=y+2,Apply[Plus, Map[l,ddist] ]} FindMinimum[f[[3]],{x,{-1,1}},{y,{-1,1}}] This also works, but timing results are inconclusive. I would appreciate suggestions on any other approaches that might speed this up. I am also asking mma support for help, and will forward their reply. Bill Harbaugh harbaugh at students.wisc.edu