- Improved detection of (complex)
dotu
anddotc
calling conventions.
-
Removed package dependencies on
stdio
andbase
. This allows testing with the new OCaml multicore compiler, which as of now does not yet support these libraries.Thanks to Anthony Scemama for this contribution!
-
Added support for const char strings in stubs due to stricter handling in newer OCaml runtimes. This eliminates C-compiler warnings.
-
Switched to OPAM file generation via
dune-project
-
Added missing autogenerated dependencies to
dune
file to avoid build system race conditions.
- Fixed warnings in C-stubs
- Fixed an API misspecification in
ptsv
, which leads to incorrect results.
Thanks to Nicolas Ratier [email protected] for the bug report!
- Fixed a portability bug in (complex)
dotu
anddotc
due to BLAS API incompatibilities between different platforms.
- Switched to dune, dune-release, and OPAM 2.0
-
Fixed warnings and errors in configuration code due to upstream changes.
-
Fixed bug in Vec.ssqr.
WARNINIG WARNING WARNING
-
User code changes in existing code required! User code may still compile, but can behave differently!
The following functions are affected:
potrf
potri
potrs
The above functions now do not support the
jitter
argument anymore. Users should remove the flag from calls to the above functions and call the newMat.add_const_diag
function if they need to add jitter. This call should happen right before the (now required) call topotrf
.More importantly,
potri
andpotrs
now do not support thefactorize
flag anymore, which would callpotrf
automatically beforehand. This was the (ill-conceived) default, which makes it harder to port LAPACK code to Lacaml. In order to upgrade your code, please do the following:-
If
potri
orpotrs
were passed~factorize:false
, just remove the flag. -
If
potri
orpotrs
were passed~factorize:true
or were called without thefactorize
flag, remove the flag if necessary and callpotrf
with the exactly corresponding arguments before.
Luckily, these functions are typically used rarely, and the changes are trivial. Apologies anyway for the churn!
-
New functions
Mat.add_const_diag
for adding a constant to the diagonal of a (sub-)matrix.
-
orgqr
now detects ifm < n
and raises an exception instead of printing a Fortran error message and continuing. -
Eigenvalue offsets should now work correctly.
-
Improved documentation.
-
Fixed bugs accessing lower pentagonal matrix patterns
-
Fixed library override issue on Mac OS X
- Fixed wrongly capitalized build targets missed due to Mac OS X file system case insensitivity.
-
Switched to jbuilder and topkg
-
API changes
-
trmm
andtrsm
now do not label argumenta
anymore -
Many matrix functions now support an optional
patt
argument, which can be used to specify rectangular, triagonal, trapezoidal, and pentagonal patterns on which to perform an operation. -
New functions
Mat.sum_prod
computes the sum of element-wise products of two matrices. Some use cases are already covered byMat.gemm_trace
, but the latter does not support patterns.
-
-
Improved C-code to better support SIMD compiler optimizations
-
Many internal improvements, including untagged and unboxed passing of parameters to and from external functions.
-
Compilation now uses
-march=native -O3 -ffast-math
by default, which should be safe and exploit SIMD on platforms that support it to greatly improve performance of some operations. -
Improved documentation
-
Improved configuration and build process