def rastrigin[N<:Nat, A: Field : Trig](x: Dimension[N,A]) =
  10 * x.size + x.mapSum(xi => xi ** 2 - 10 * cos(2 * pi * xi))

import shapeless._

import benchmarks.Benchmarks._

val x = Sized(1.0, 2.0, 3.0)
spherical(x)
// res1: Double = 14.0

val z = (1.0, 2.0)
ackley2(z)
// res2: Double = -191.2527797034299

import benchmarks.Benchmarks._

implicit val jd = JetDim(2)

val a = 1.0 + Jet.h[Double](0)
val b = 2.0 + Jet.h[Double](1)
val x = Sized(a, b)

spherical(x)
// res5: Jet[Double] = (5.0 + [2.0, 4.0]h)

val z = (a, b)
ackley2(z)
// res7: Jet[Double] = (-191.2527797034299 + [1.7106168652106453, 3.4212337304212905]h)

import benchmarks.cec2005.Benchmarks._
import benchmarks.implicits._

val x = Sized(1.0, 2.0)

f1(x)
// res1: Double = 4412.64790162

import cilib._

import benchmarks.Benchmarks
import benchmarks.implicits._

val f = Benchmarks.spherical[nat._2,Double] _

// convert function to Eval for use in cilib
val eval = f.unconstrained