Charles Bouillaguet
Charles Bouillaguet
we should add `mzd_inv(A)` as a shorthand for inverting an invertible matrix.
The documentation says that the PLUQ function returns matrices satisfying A = PLUQ. But the test code actually checks that PAQ^t = LU, and the test pass… So the specification...
Reported by @kiwifb François Bissey: I have been giving some thought about the presence of SIMD_CFLAGS and OPENMP_CFLAGS in m4ri.pc. The simd flags are not necessary for a customer of...
Currently, each release has its own SONAME so no two releases are binary compatible. However, the ABI doesn't change much and we can move to a more stable SONAME. A...
@malb said: Our performance for * 64x64 times 64x1, * 64x64 times 64x64 and * Nx64 times 64x64 matrix multiplication is embarrassingly bad. We must improve performance here, probably by...
Indeed, what the `mzd_solve_left()` function does is exactly what is called `M.solve_right()` in sage. It could be somewhat confusing, but I am not sure what is the best solution: *...
The documentation of `mzd_transpose` does not say that it does not work in-place for square matrices
The title says it all.
The following code: ``` rate = Meter() for i in range(int(65/0.05)): rate.mark() time.sleep(0.05) print(rate.one_minute_rate) print(rate.mean_rate) ``` yields: ``` 12.744700289590714 19.916839105797937 ``` Why does the one-minute moving average differs from the...
Transposition made 50% of mzd.c and is somehwat ugly. This PR splits it from mzd.c while preserving "git blame" history for everything.
It needs it anyway, and the caller necessarily knows it.