line_search_o2

 function  [alpha, img, dv, opt] = line_search_o2(imgk, dx, data1, img1, N, W, hps2RtR, dv0, opt)
 function  [alpha, img, dv, opt] = line_search_o2(imgk, dx, data1, img1, N, W, hps2RtR, hpt2LLt, dv0, opt)
 line search function with a fitted polynomial of O(2) -- a second order polynomial
 (C) 2013 Alistair Boyle
 License: GPL version 2 or version 3

0001 function  [alpha, img, dv, opt] = line_search_o2(imgk, dx, data1, img1, N, W, hps2RtR, hpt2LLt, dv0, opt)
0002 % function  [alpha, img, dv, opt] = line_search_o2(imgk, dx, data1, img1, N, W, hps2RtR, dv0, opt)
0003 % function  [alpha, img, dv, opt] = line_search_o2(imgk, dx, data1, img1, N, W, hps2RtR, hpt2LLt, dv0, opt)
0004 % line search function with a fitted polynomial of O(2) -- a second order polynomial
0005 % (C) 2013 Alistair Boyle
0006 % License: GPL version 2 or version 3
0007 
0008 pf_max = 2; % max poly-fit order = 2 --> second-order polynomial
0009 retry = 0; % first any only call through this function, no recursion
0010 if nargin == 9
0011    hpt2LLt = 1;
0012 end
0013 [alpha, img, dv, opt] = line_search_onm2(imgk, dx, data1, img1, N, W, hps2RtR, hpt2LLt, dv0, opt, retry, pf_max);