toppure_concave
Return value
y2: f(x2) = y2. Double precision array
Parameters
x1
in
required
x2
in
required
y1
in
required
y2
der2
x
Examples
y2 = spl_incr(x, y, x2)
n = 100L
x = (dindgen(n))^2
y = abs(randomn(0, n))
y[n/2:n/2+1] = 0.
y[n-n/3] = 0.
y[n-n/6:n-n/6+5] = 0.
y = total(y, /cumulative, /double)
x2 = dindgen((n-1)^2)
n2 = n_elements(x2)
print, min(y[1:n-1]-y[0:n-2]) LT 0
y2 = spl_incr( x, y, x2)
splot, x, y, xstyle = 1, ystyle = 1, ysurx=.25, petit = [1, 2, 1], /land
oplot, x2, y2, color = 100
c = y2[1:n2-1] - y2[0:n2-2]
print, min(c) LT 0
print, min(c, max = ma), ma
splot,c,xstyle=1,ystyle=1, yrange=[-.01,.05], ysurx=.25, petit = [1, 2, 2], /noerase
oplot,[0, n_elements(c)], [0, 0], linestyle = 1
Version history
History
Sebastien Masson (smasson@lodyc.jussieu.fr): May-Dec 2005
Known issues
Restrictions
It might be possible that y2[i+1]-y2[i] has very small negative
values (amplitude smaller than 1.e-6)...