# Integration test for the general optimization suite
test_that("General optimization suite works on quadratic function", {
  
  # Define a simple quadratic function
  quad <- function(x) (x[1] - 2)^2 + (x[2] + 1)^2
  
  algorithms <- list(
    bfgs            = bfgs,
    dfp             = dfp,
    dogleg          = dogleg,
    double_dogleg   = double_dogleg,
    modified_newton = modified_newton,
    newton_raphson  = newton_raphson,
    l_bfgs_b        = l_bfgs_b
  )
  
  start_val <- c(0, 0)
  target <- c(2, -1)
  
  for (name in names(algorithms)) {
    res <- algorithms[[name]](
      start     = start_val,
      objective = quad,
      control   = list(max_iter = 10000, initial_delta = 2.0, tol_rel_x = 1e-6)
    )
    
    expect_true(res$converged, info = paste("Convergence failed:", name))
    expect_equal(res$par, target, tolerance = 1e-2, 
                 info = paste("Parameter mismatch:", name))
  }
})