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