def trajectory(self, individual, *f_args):
dy = control_problem.lorenz_in_2(gp.sympy_phenotype(individual),
**self.params)
return control_problem.integrate(dy,
yinit=self.yinit,
x=self.x,
f_args=f_args)
def trajectory(self, individual, *f_args):
dy = control_problem.anharmonic_oscillator(
gp.sympy_phenotype(individual), **self.params)
return control_problem.integrate(dy,
yinit=self.yinit,
x=self.x,
f_args=f_args)
def integral(ind, *const):
f = gp.sympy_phenotype(ind)
def f_square(x, *const):
return f(x, *const)**2
s, *_ = scipy.integrate.quad(f_square, 0, 1, args=tupleize(const))
return s
def trajectory(self, individual, *f_args):
dy = self.dynsys(gp.sympy_phenotype(individual))
return control_problem.integrate(dy,
yinit=self.yinit,
x=self.x,
f_args=f_args)