R/bisection.search.R

In SuperCurve: RPPA Analysis Package

###
### $Id: bisection.search.R 833 2014-06-27 16:17:07Z proebuck $
###
### fields, Tools for spatial data
### Copyright 2004-2007, Institute for Mathematics Applied Geosciences
### University Corporation for Atmospheric Research
### Used with permission of Doug Nychka under Artistic license (2008)

## :KRC: Should this be placed as an internal function inside the
## .generic.trim function that is the only thing that ever uses it.

## What does this function do?
## What do the inputs mean?

##-----------------------------------------------------------------------------
"bisection.search" <-
function (x1, x2, f, tol = 1e-07, niter = 25, f.extra = NA, upcross.level = 0) 
{
    f1 <- f(x1, f.extra) - upcross.level
    f2 <- f(x2, f.extra) - upcross.level
    if (f1 > f2) 
        stop(" f1 must be < f2 ")
    iter <- niter
    for (k in 1:niter) {
        xm <- (x1 + x2)/2
        fm <- f(xm, f.extra) - upcross.level
        if (fm < 0) {
            x1 <- xm
            f1 <- fm
        }
        else {
            x2 <- xm
            f2 <- fm
        }
        if (abs(fm) < tol) {
            iter <- k
            break
        }
    }
    xm <- (x1 + x2)/2
    fm <- f(xm, f.extra) - upcross.level
    list(x = xm, fm = fm, iter = iter)
}