Был опыт реализации на Си мат.модели, написанной на Матлабе физиками.
Модель считала все в double. На Си для скорости использовали float. На прямую результат проверить было не возможно как минимум из-за ошибок округления, а если еще появлялись ошибки в Си реализации (а они были), то вообще труба, т.к. не ясно было даже в какой части кода ошибка.
Пришлось полностью реализовать модель на Python с возможностью переключения режима расчета выставлением одной опции с double на float. Эту промежуточную модель отладили с первоначальной моделью в режиме double. Затем переключили на float и отлаживали реализацию на Си с промежуточной моделью.
Работы много, но другого способа подтверждать правильность реализации мы не придумали.
В итоге схема получилась вполне рабочая.
Теоретически можно было доработать Матлаб модель, чтоб она работала с floatами, но в Матлабе оказалось достаточно проблематично контролировать каждую операцию. На пинтоне это сделать можно.