Vladimír Chalupecký
Vladimír Chalupecký
This is just a draft for now to discuss and collect feedback. My doubts are: 1. Should `Factorize` return a bool? And the closely related question, should `PivotedCholesky.SolveTo/SolveVecTo` exist? If...
`Dlatdf` relies on undocumented content stored by Dgecon in `work[n:2*n]`. We should document it and test for it. _Originally posted in https://github.com/gonum/gonum/pull/1739#discussion_r733668864_
It is well known that the current BLAS standard and the reference implementation give in some cases inconsistent results in the presence of NaN/Inf, see [1] or [2]. Since our...
From the abstract of https://doi.org/10.1145/3061665: > The square root of a sum of squares is well known to be prone to overflow and underflow. Ad hoc scaling of intermediate results,...
### What are you trying to do? Efficiently solve linear systems with a general _tridiagonal_ matrix (see also https://github.com/gonum/gonum/issues/1471#issuecomment-706550369). ### How does Gonum not allow you to achieve your goal?...
### What are you trying to do? Efficiently solve linear systems with a general band matrix. ### How does Gonum not allow you to achieve your goal? Existing `mat.LU` solver...
See discussion in https://github.com/gonum/gonum/pull/1399#discussion_r457749964
`LAPACKE_?larfb_work` passes `work` and `ldwork` to the Fortran routine unmodified even if the matrix layout is row major. This is extremely surprising bordering a bug because it means that the...
In some cases, DROTMG returns values that, when used in DROTM, do not zero out the second component of the given vector. This can be observed for example with the...
DLATB4 is called from DLATTB with IMAT negative for lower triangular matrices: ``` IF( UPPER ) THEN CALL DLATB4( PATH, IMAT, N, N, TYPE, KL, KU, ANORM, MODE, $ CNDNUM,...