def shoot(self):
rg, N, T = self.getSimData()
if self.noJ:
if self.unconstrained:
(x,m,z,v) = kernelMatrix_fort.shoot_unconstrained(self.dt, \
self.sfactor, self.kvs, self.kvo, \
self.khs, self.kho, \
self.alpha, self.x0, self.template, \
self.z0*self.z0weight,
self.num_times, rg.num_nodes)
else:
(x,m,z,v) = kernelMatrix_fort.shoot_noj(self.dt, \
self.sfactor, self.kvs, self.kvo, \
self.khs, self.kho, \
self.alpha, self.x0, self.template, \
-1.0*self.D_template,
self.num_times, rg.num_nodes)
else:
(x,m,z,J,v) = kernelMatrix_fort.shoot(self.dt, \
self.sfactor, self.kvs, self.kvo, \
self.khs, self.kho, \
self.alpha, self.x0, self.template, \
-1.0*self.D_template,
self.num_times, rg.num_nodes)
self.J = J
self.xt = x
self.m = m
self.z = z
self.v = v
def shoot(self):
rg, N, T = self.getSimData()
(x,m,z,J,v) = kernelMatrix_fort.shoot(self.dt, \
self.sfactor, self.kvs, self.kvo, \
self.khs, self.kho, \
self.alpha, self.x0, self.template, \
-1.0*self.D_template,
self.num_times, rg.num_nodes)
self.xt = x
self.m = m
self.z = z
self.J = J
self.v = v
