compute_cluster: automatic computation of the number of pullback attracting...

In mcrucifix/iceages: Ice age models and analysis

automatic computation of the number of pullback attracting trajectories

Description

Wrapper to basin, to count the number of pullback attracting trajectories for a series of parameters, with multicore support if available

Usage

compute_cluster(model, par, tback, t0, Astro, ncores = 14, deltat = 0.01, ...)

Arguments

model	One of the models out of the packaged model list (data(models)) or user-supplied.
par	Matrix of parameters (possibly with a single row), with row names correspoding to parameters to be suppiled to the model.
tback	Initial time, for which a series of initial conditions as defined in model are set
t0	Time at which states are evaluated and clustered into a number of distinct solutions

Astro	Astronomical forcing list, as supplied by routine read_astro
ncores	Number of cores over which calculations will be distributed. Inactive if package::multicore is unavailable

deltat

Optional parameter giving the time step of the numerical scheme for integrating model.

...	Other models supplied to the basin function.

Value

A vector of integers of length nrow(par) with the number of distinct solutions as a function of parameters.

Author(s)

M. Crucifix

Examples

data(models)

parlist <- as.list(models$pp04$spar)

parlist$gamma = seq(0.2, 1.2, 0.2)
parlist$omega = seq(0.2, 1.2, 0.2)

pargrid <- as.matrix(expand.grid(parlist))

Astro <- read_astro(1,1)

C <- compute_cluster(models$pp04, pargrid, -100, 0, Astro, ncores=2)

print(C)

mcrucifix/iceages documentation built on Jan. 11, 2023, 9:17 p.m.