Profile to find bottlenecks

[naesseth]

Profile the code to find where bottlenecks are. Perhaps the way a model/formalism is defined and used?

[naesseth]

After some basic profiling it seems most of the time is spent generating normal random numbers, second most is evaluating weights. Can we get big speedup by Cythonizing models?