Add HS220 to HS225

src/ADNLPProblems/hs220.jl

@@ -0,0 +1,16 @@
+export hs220
+
+function hs220(; n::Int = default_nvar, type::Val{T} = Val(Float64), kwargs...) where {T}
+  function f(x)
+    return x[1]
+  end
+  x0 = T[25000, 25000]
+  lvar = T[1, 0]
+  uvar = T[Inf, Inf]
+  function c!(cx, x)
+    cx[1] = (x[1] - 1)^3 - x[2]
+    return cx
+  end
+  lcon = ucon = T[0]
+  return ADNLPModels.ADNLPModel!(f, x0, lvar, uvar, c!, lcon, ucon, name = "hs220"; kwargs...)
+end

src/ADNLPProblems/hs221.jl

@@ -0,0 +1,17 @@
+export hs221
+
+function hs221(; n::Int = default_nvar, type::Val{T} = Val(Float64), kwargs...) where {T}
+  function f(x)
+    return -x[1]
+  end
+  x0 = T[0.25, 0.25]
+  lvar = T[0, 0]
+  uvar = T[Inf, Inf]
+  function c!(cx, x)
+    cx[1] = (1 - x[1])^3 - x[2]
+    return cx
+  end
+  ucon = T[Inf]
+  lcon = T[0]
+  return ADNLPModels.ADNLPModel!(f, x0, lvar, uvar, c!, lcon, ucon, name = "hs221"; kwargs...)
+end

src/ADNLPProblems/hs222.jl

@@ -0,0 +1,18 @@
+export hs222
+
+function hs222(; n::Int = default_nvar, type::Val{T} = Val(Float64), kwargs...) where {T}
+  function f(x)
+    n = length(x)
+    return -x[1]
+  end
+  x0 = T[1.3, 0.2]
+  lvar = T[0, 0]
+  uvar = T[Inf, Inf]
+  function c!(cx, x)
+    cx[1] = (1 - x[1])^3 - x[2] + 125 // 1000
+    return cx
+  end
+  ucon = T[Inf]
+  lcon = T[0]
+  return ADNLPModels.ADNLPModel!(f, x0, lvar, uvar, c!, lcon, ucon, name = "hs222"; kwargs...)
+end

src/ADNLPProblems/hs223.jl

@@ -0,0 +1,19 @@
+export hs223
+
+function hs223(; n::Int = default_nvar, type::Val{T} = Val(Float64), kwargs...) where {T}
+  function f(x)
+    n = length(x)
+    return -x[1]
+  end
+  x0 = T[0.1, 3.3]
+  lvar = T[0, 0]
+  uvar = T[10, 10]
+  function c!(cx, x)
+    cx[1] = exp(exp(x[1]))
+    cx[2] = x[2] - exp(exp(x[1]))
+    return cx
+  end
+  ucon = T[Inf; Inf]
+  lcon = T[0; 0]
+  return ADNLPModels.ADNLPModel!(f, x0, lvar, uvar, c!, lcon, ucon, name = "hs223"; kwargs...)
+end

src/ADNLPProblems/hs224.jl

@@ -0,0 +1,20 @@
+export hs224
+
+function hs224(; n::Int = default_nvar, type::Val{T} = Val(Float64), kwargs...) where {T}
+  function f(x)
+    n = length(x)
+    return 2 * x[1]^2 + x[2]^2 - 48 * x[1] - 40 * x[2]
+  end
+  x0 = T[0.1, 0.1]
+  lvar = T[0, 0]
+  uvar = T[6, 6]
+  lcon = T[0; -18; 0; -8]
+  ucon = T(Inf) * ones(T, 4)
+  A = T[
+    1 3;
+    -1 -3;
+    1 1;
+    -1 -1
+  ]
+  return ADNLPModels.ADNLPModel(f, x0, lvar, uvar, sparse(A), lcon, ucon, name = "hs224"; kwargs...)
+end

src/ADNLPProblems/hs225.jl

@@ -0,0 +1,19 @@
+export hs225
+
+function hs225(; n::Int = default_nvar, type::Val{T} = Val(Float64), kwargs...) where {T}
+  function f(x)
+    return x[1]^2 + x[2]^2
+  end
+  x0 = T[3, 1]
+  function c!(cx, x)
+    cx[1] = x[1]^2 + x[2]^2 - 1
+    cx[2] = 9 * x[1]^2 + x[2]^2 - 9
+    cx[3] = x[1]^2 - x[2]
+    cx[4] = x[2]^2 - x[1]
+    return cx
+  end
+  A = T[1 1]
+  lcon = vcat(1, zeros(T, 4))
+  ucon = T(Inf) * ones(T, 5)
+  return ADNLPModels.ADNLPModel!(f, x0, sparse(A), c!, lcon, ucon, name = "hs225"; kwargs...)
+end

src/Meta/hs220.jl

@@ -0,0 +1,25 @@
+hs220_meta = Dict(
+  :nvar => 2,
+  :variable_nvar => false,
+  :ncon => 1,
+  :variable_ncon => false,
+  :minimize => true,
+  :name => "hs220",
+  :has_equalities_only => true,
+  :has_inequalities_only => false,
+  :has_bounds => true,
+  :has_fixed_variables => false,
+  :objtype => :other,
+  :contype => :general,
+  :best_known_lower_bound => 1,
+  :best_known_upper_bound => 25000,
+  :is_feasible => true,
+  :defined_everywhere => missing,
+  :origin => :unknown,
+)
+get_hs220_nvar(; n::Integer = default_nvar, kwargs...) = 2
+get_hs220_ncon(; n::Integer = default_nvar, kwargs...) = 1
+get_hs220_nlin(; n::Integer = default_nvar, kwargs...) = 0
+get_hs220_nnln(; n::Integer = default_nvar, kwargs...) = 1
+get_hs220_nequ(; n::Integer = default_nvar, kwargs...) = 1
+get_hs220_nineq(; n::Integer = default_nvar, kwargs...) = 0

src/Meta/hs221.jl

@@ -0,0 +1,25 @@
+hs221_meta = Dict(
+  :nvar => 2,
+  :variable_nvar => false,
+  :ncon => 1,
+  :variable_ncon => false,
+  :minimize => true,
+  :name => "hs221",
+  :has_equalities_only => false,
+  :has_inequalities_only => true,
+  :has_bounds => true,
+  :has_fixed_variables => false,
+  :objtype => :other,
+  :contype => :general,
+  :best_known_lower_bound => -1,
+  :best_known_upper_bound => -0.25,
+  :is_feasible => true,
+  :defined_everywhere => missing,
+  :origin => :unknown,
+)
+get_hs221_nvar(; n::Integer = default_nvar, kwargs...) = 2
+get_hs221_ncon(; n::Integer = default_nvar, kwargs...) = 1
+get_hs221_nlin(; n::Integer = default_nvar, kwargs...) = 0
+get_hs221_nnln(; n::Integer = default_nvar, kwargs...) = 1
+get_hs221_nequ(; n::Integer = default_nvar, kwargs...) = 0
+get_hs221_nineq(; n::Integer = default_nvar, kwargs...) = 1

src/Meta/hs222.jl

@@ -0,0 +1,25 @@
+hs222_meta = Dict(
+  :nvar => 2,
+  :variable_nvar => false,
+  :ncon => 1,
+  :variable_ncon => false,
+  :minimize => true,
+  :name => "hs222",
+  :has_equalities_only => false,
+  :has_inequalities_only => true,
+  :has_bounds => true,
+  :has_fixed_variables => false,
+  :objtype => :other,
+  :contype => :general,
+  :best_known_lower_bound => -1.5,
+  :best_known_upper_bound => -1.3,
+  :is_feasible => true,
+  :defined_everywhere => missing,
+  :origin => :unknown,
+)
+get_hs222_nvar(; n::Integer = default_nvar, kwargs...) = 2
+get_hs222_ncon(; n::Integer = default_nvar, kwargs...) = 1
+get_hs222_nlin(; n::Integer = default_nvar, kwargs...) = 0
+get_hs222_nnln(; n::Integer = default_nvar, kwargs...) = 1
+get_hs222_nequ(; n::Integer = default_nvar, kwargs...) = 0
+get_hs222_nineq(; n::Integer = default_nvar, kwargs...) = 1

src/Meta/hs223.jl

@@ -0,0 +1,25 @@
+hs223_meta = Dict(
+  :nvar => 2,
+  :variable_nvar => false,
+  :ncon => 2,
+  :variable_ncon => false,
+  :minimize => true,
+  :name => "hs223",
+  :has_equalities_only => false,
+  :has_inequalities_only => true,
+  :has_bounds => true,
+  :has_fixed_variables => false,
+  :objtype => :other,
+  :contype => :general,
+  :best_known_lower_bound => -log(log(10)),
+  :best_known_upper_bound => -0.1,
+  :is_feasible => true,
+  :defined_everywhere => missing,
+  :origin => :unknown,
+)
+get_hs223_nvar(; n::Integer = default_nvar, kwargs...) = 2
+get_hs223_ncon(; n::Integer = default_nvar, kwargs...) = 2
+get_hs223_nlin(; n::Integer = default_nvar, kwargs...) = 0
+get_hs223_nnln(; n::Integer = default_nvar, kwargs...) = 2
+get_hs223_nequ(; n::Integer = default_nvar, kwargs...) = 0
+get_hs223_nineq(; n::Integer = default_nvar, kwargs...) = 2

src/Meta/hs224.jl

@@ -0,0 +1,25 @@
+hs224_meta = Dict(
+  :nvar => 2,
+  :variable_nvar => false,
+  :ncon => 4,
+  :variable_ncon => false,
+  :minimize => true,
+  :name => "hs224",
+  :has_equalities_only => false,
+  :has_inequalities_only => true,
+  :has_bounds => true,
+  :has_fixed_variables => false,
+  :objtype => :other,
+  :contype => :linear,
+  :best_known_lower_bound => -304,
+  :best_known_upper_bound => -304,
+  :is_feasible => true,
+  :defined_everywhere => missing,
+  :origin => :unknown,
+)
+get_hs224_nvar(; n::Integer = default_nvar, kwargs...) = 2
+get_hs224_ncon(; n::Integer = default_nvar, kwargs...) = 4
+get_hs224_nlin(; n::Integer = default_nvar, kwargs...) = 4
+get_hs224_nnln(; n::Integer = default_nvar, kwargs...) = 0
+get_hs224_nequ(; n::Integer = default_nvar, kwargs...) = 0
+get_hs224_nineq(; n::Integer = default_nvar, kwargs...) = 4

src/Meta/hs225.jl

@@ -0,0 +1,25 @@
+hs225_meta = Dict(
+  :nvar => 2,
+  :variable_nvar => false,
+  :ncon => 5,
+  :variable_ncon => false,
+  :minimize => true,
+  :name => "hs225",
+  :has_equalities_only => false,
+  :has_inequalities_only => true,
+  :has_bounds => false,
+  :has_fixed_variables => false,
+  :objtype => :other,
+  :contype => :general,
+  :best_known_lower_bound => 2,
+  :best_known_upper_bound => 2,
+  :is_feasible => true,
+  :defined_everywhere => missing,
+  :origin => :unknown,
+)
+get_hs225_nvar(; n::Integer = default_nvar, kwargs...) = 2
+get_hs225_ncon(; n::Integer = default_nvar, kwargs...) = 5
+get_hs225_nlin(; n::Integer = default_nvar, kwargs...) = 1
+get_hs225_nnln(; n::Integer = default_nvar, kwargs...) = 4
+get_hs225_nequ(; n::Integer = default_nvar, kwargs...) = 0
+get_hs225_nineq(; n::Integer = default_nvar, kwargs...) = 5

src/PureJuMP/hs220.jl

@@ -0,0 +1,27 @@
+# Hock and Schittkowski problem number 220.
+#
+#   Source:
+#   Problem 220 in
+#   K. Schittkowski,
+#   More Test Examples for Nonlinear Programming Codes,
+#   Lectures Notes in Economics and Mathematical Systems 282,
+#   Springer Verlag, Heidelberg, 1987.
+#
+#   
+#
+# T. Migot, Montreal, 2023.
+
+export hs220
+
+"HS220 model"
+function hs220(args...; kwargs...)
+  nlp = Model()
+  x0 = [25000, 25000]
+  lvar = [1, 0]
+  @variable(nlp, x[i = 1:2] ≥ lvar[i], start = x0[i])
+
+  @NLobjective(nlp, Min, x[1])
+  @NLconstraint(nlp, (x[1] - 1)^3 - x[2] == 0)
+
+  return nlp
+end

src/PureJuMP/hs221.jl

@@ -0,0 +1,27 @@
+# Hock and Schittkowski problem number 221.
+#
+#   Source:
+#   Problem 221 in
+#   K. Schittkowski,
+#   More Test Examples for Nonlinear Programming Codes,
+#   Lectures Notes in Economics and Mathematical Systems 282,
+#   Springer Verlag, Heidelberg, 1987.
+#
+#   
+#
+# T. Migot, Montreal, 2023.
+
+export hs221
+
+"HS221 model"
+function hs221(args...; kwargs...)
+  nlp = Model()
+  x0 = [0.25, 0.25]
+  lvar = [0, 0]
+  @variable(nlp, x[i = 1:2] ≥ lvar[i], start = x0[i])
+
+  @NLobjective(nlp, Min, -x[1])
+  @NLconstraint(nlp, (1 - x[1])^3 - x[2] >= 0)
+
+  return nlp
+end

src/PureJuMP/hs222.jl

@@ -0,0 +1,27 @@
+# Hock and Schittkowski problem number 222.
+#
+#   Source:
+#   Problem 222 in
+#   K. Schittkowski,
+#   More Test Examples for Nonlinear Programming Codes,
+#   Lectures Notes in Economics and Mathematical Systems 282,
+#   Springer Verlag, Heidelberg, 1987.
+#
+#   
+#
+# T. Migot, Montreal, 2023.
+
+export hs222
+
+"HS222 model"
+function hs222(args...; kwargs...)
+  nlp = Model()
+  x0 = [1.3, 0.2]
+  lvar = [0, 0]
+  @variable(nlp, x[i = 1:2] ≥ lvar[i], start = x0[i])
+
+  @NLobjective(nlp, Min, -x[1])
+  @NLconstraint(nlp, (1 - x[1])^3 - x[2] + 0.125 >= 0)
+
+  return nlp
+end

src/PureJuMP/hs223.jl

@@ -0,0 +1,29 @@
+# Hock and Schittkowski problem number 223.
+#
+#   Source:
+#   Problem 223 in
+#   K. Schittkowski,
+#   More Test Examples for Nonlinear Programming Codes,
+#   Lectures Notes in Economics and Mathematical Systems 282,
+#   Springer Verlag, Heidelberg, 1987.
+#
+#   
+#
+# T. Migot, Montreal, 2023.
+
+export hs223
+
+"HS223 model"
+function hs223(args...; kwargs...)
+  nlp = Model()
+  x0 = [0.1, 3.3]
+  lvar = [0, 0]
+  uvar = [10, 10]
+  @variable(nlp, uvar[i] ≥ x[i = 1:2] ≥ lvar[i], start = x0[i])
+
+  @NLobjective(nlp, Min, -x[1])
+  @NLconstraint(nlp, exp(exp(x[1])) >= 0)
+  @NLconstraint(nlp, x[2] - exp(exp(x[1])) >= 0)
+
+  return nlp
+end