gsl-ocaml /

Filename Size Date modified Message
examples
lib
66 B
605 B
659 B
289 B
3.8 KB
34.4 KB
876 B
1.1 KB
946 B
2.4 KB
501 B
9.8 KB
9.6 KB
364 B
2.4 KB
1.6 KB
5.0 KB
19.1 KB
195.7 KB

GSL-OCaml - GSL-Bindings for OCaml

This is an interface to the GSL (GNU scientific library) for the OCaml-language.

Requirements

The platform must not align doubles on double-word addresses, i.e. the C-macro ARCH_ALIGN_DOUBLE must be undefined in the OCaml C-configuration header in <caml/config.h>.

Installation

$ ./configure
$ make
$ make install

Documentation

Check the GSL manual to learn more about the GNU Scientific Library.

You can browse the OCaml module interfaces as ocamldoc-generated HTML files in directory API.docdir after building the documentation with make doc.

Usage Hints

Vectors and Matrices

There are several datatypes for handling vectors and matrices.

  • Modules Gsl.Vector, Gsl.Vector.Single, Gsl.Vector_complex, Gsl.Vector_complex.Single, and the corresponding matrix modules use bigarrays with single or double precision and real or complex values.

  • Modules Gsl.Vector_flat, Gsl.Vector_complex_flat, and the corresponding matrix modules use a record wrapping a regular OCaml float array. This is the equivalent of the gsl_vector and gsl_matrix structs in GSL.

  • Module Gsl.Vectmat defines a sum type with polymorphic variants that regroups these two representations. For instance:

    Gsl.Vectmat.v_add (`V v1) (`VF v2)
    

    adds a vector in an OCaml array to a bigarray.

  • Modules Gsl.Blas Gsl.Blas_flat and Gsl.Blas_gen provide a (quite incomplete) interface to CBLAS for these types.

ERROR HANDLING

Errors in GSL functions are reported as exceptions:

Gsl.Error.Gsl_exn (errno, msg)

You have to call Gsl.Error.init () to initialize error reporting. Otherwise, the default GSL error handler is used and aborts the program, leaving a core dump (not so helpful with OCaml).

If a callback (for minimizers, solvers, etc.) raises an exception, gsl-ocaml either returns GSL_FAILURE or NaN to GSL depending on the type of callback. In either case the original OCaml exception is not propagated. The GSL function will either return normally (but probably with values containing NaNs somewhere) or raise a Gsl_exn exception.

Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.