1. Added support for Jacobian in simulators.py. It is not used as it
slows performance by 50%.