tests/testthat/test-rosenbrock.R

In RGN: Robust-Gauss Newton (RGN) Optimization of Sum-of-Squares Objective Function

test_that('rosenbrock', {

  simFunc_rosenbrock=function(x){
    r=rep(0,2)
    r[1] = 1.0-x[1]
    r[2]=10.0*(x[2]-x[1]**2)
    return(r)
  }

  tmp = rgn(simFunc=simFunc_rosenbrock,
            par=c(-1.0,  0.0), lower=c(-1.5, -1.0), upper=c( 1.5,  3.0),
            simTarget=c(0,0),
            control=list(dump=0))
  error=tmp$error;message=tmp$message;par=tmp$par;value=tmp$value;counts=tmp$counts;info=tmp$info

  # expected output based on F90 RGN Rosenbrock example
  expect_equal(par,c(1,1))
  expect_equal(value,0)
  expect_equal(counts,116)
  expect_equal(info$nIter,16)
  expect_equal(info$termFlag,2)

})