Репа производительности для простого численного метода
Я реализовал симплектический метод Эйлера, используя repa. Это довольно простой метод. К сожалению, результат, который я получаю, не так хорош, как я бы надеялся. Ребята, ребята, были очень полезны (в частности, Бен Липпмайер), и я последовал его совету по вставке и добавлению аннотаций строгости. К сожалению, это, похоже, не помогло. Я делаю ядро с
ghc --make -O2 RepaOnly.hs -main-is RepaOnly -ddump-prep -dsuppress-all -dppr-case-as-let -dppr-cols200 > RepaOnly.hcr
К сожалению, это кажется довольно большим, и я в настоящее время не знаю, как поступить. Я думал, что у более широкого сообщества Haskell могут быть некоторые идеи о том, как отладить это или указать, где я иду неправильно.
Код здесь: https://gist.github.com/idontgetoutmuch/6209752
Я также могу опубликовать ядро в суть, если это полезно.