"Classical" functional response

[evadelmas]

So far the model uses the functional response as presented in Williams, Brose and Martinez 2007 (multi-species version of Yodzis and Innes 1992), using species specific (vector of length S - where S is the number of species) max. consumption rates and half saturation densities.
When using temperature scaling of biological rates, we need to be able to pass interaction specific (matrix of dim. S * S) attack rates and handling times.