Check stability

checkStability(ensemble, models, strucIdx, eigThreshold)

Given a model, checks if the largest real part of the jacobian eigenvalues is greater than the given threshold (eigThreshold), if so the model is considered invalid.

USAGE:

isModelValid = checkStability(ensemble, models, strucIdx, eigThreshold)

INPUT:

ensemble (struct): model ensemble, see buildEnsemble for fields description models (struct): model, see initialSampler for fields description strucIdx (int): ID of the model structure considered eigThreshold (double): threshold for positive eigenvalues’ real part

OUTPUT:

isModelValid (logical): whether or not the model is valid