Currently in Pre alpha release engineering, so nearly ready for experimental consumption
(though please dont publicize yet)
This is the core Package for Numerical Haskell, a project by Carter Schonwald aka Wellposed Ltd, and (soon I hope!) other contributors.
Numerical-Core is an open source component of the Wellposed® Numerical Haskell software suite.
this library is pre alpha release so not all examples / codes may work as expected. That said, the current api should be enough to prototype and typecheck algorithms.
Numerical Haskell is an effort to bring great numerical computation and data analysis tooling to haskell, and be the best possible platform for sophisticated efforts in those same domains
What array Formats are Supported
What convention is used for indexing?
When you have an index tuple, just think
x,y,z to keep track of the meaning.
Indexing tuples are written as statically sized lists, eg
This follows the tradition of x,y,z axes used in plotting. Note well: the underlying memory
order can be row OR column major or other!
All the computations on these static sized lists get specialized away into nonrecursive computations at their use sites. So in this special scenario, lists aren't a problem!
Great! Theres so many awesome ways you could help out. Look at CONTRIBUTING.md for more details. Right now theres a lot of low hanging fruit in improving test coverage, and soon there'll be many opportunities on the performance tuning and numerical algorithms/tooling areas.
see bug.md for how to file a bug report
- How do I use Numerical haskell to write fast code thats outstandingly high level !?
- The leading cause of poor performance in numerical routines (aside from poor choice in algorithms) is bad memory locality, which has but a single easy cure: ** block recursive algorithms **
- Yes, you heard me, in compiled languages recursion is pretty cheap outside of the inner most loops! It also is a fantastic tool for facilitating good memory locality!
- I'm totally serious, try out the benchmarks for the various versions of the same routines we provide!
- But, what about fusion?
- Because of certain aspects of the numerical haskell design, we can't automagically use the fusion optimization facilities of the underlying array representations such as Vector.
Many member of the Numerical Haskell community can be found on
#numerical-haskell on freenode IRC
There is also the numericalhaskell mailing list
The community provides some basic support through the IRC channel, Mailing list, and the relevant project issue trackers.
If your support needs can't be resolved though those channels, please do not hesistate to contact Wellposed (aka Carter) to find out more about our support and professional services options.