Put in c version + sensitivity computation