R/pl.step.r

In jarad/tlpl: Tau-leaped particle learning

# A single step of PL in tau-leaping model

pl.step = function(y,tau,mod,part,while.max=1000) 
{
    check.system(mod)
    #check.particles(part)

    out = .C("discrete_all_particle_update_R",
             as.integer(mod$s), as.integer(mod$r), as.integer(t(mod$Pre)), as.integer(mod$stoich),
             as.integer(y), as.double(tau), as.integer(part$n), as.integer(while.max),
             X=as.integer(t(part$X)), hyper=as.double(t(part$hyper)), success=integer(part$n))
    return(list(X    = matrix(out$X, part$n, mod$s, byrow=T), 
                hyper= matrix(out$hyper, part$n, 4*mod$r, byrow=T), 
                n=part$n,
                success = out$success))
}