[Internal function] Merge 2 increasing sequences using scilab syntax; used to compute the Legendre-Fenchel conjugate.
fH = _pl_fusionsci(C,S)
Column vector. C(i) = (f(X(i+1) - f(X(i)) / (X(i+1) - X(i)) for some function f on a grid X.
Column vector. The grid on which we want to compute the conjugate: f* is evaluated on S.
Column vector. Index at which the slope S(i) support the epigraph.
Calling _pl_fusionsci(C,S) is equivalent to calling _pl_fusion(C,S). While the former uses scilab built-in functions to speed-up computation, contrary to the later it does not guarantee a linear worst-case running time. See _pl_fusion for more information.
This function is not intended to be used directly but to be called during the computation of the conjugate using the LLT algorithm by the function pl_lft_llt.
Yves Lucet
, University of British Columbia, BC, Canada