Mikkel Paltorp

I am experiencing the same issue with some code breaking as a result of 0.5.4. While the issue can be fixed by explicitly rewriting the problem using LinearMaps and doing...

Yes, that is the page I am referring to. The note from Knut clear up some things, and I think that should be referenced 👍 I am still not sure...

> Yes, for tensors higher than 2nd order, I find it hard to use linear algebra notation. Even for 2nd order objects, I prefer the tensor notation to avoid the...

I just tried to follow your suggestion. When equilibration is disabled the problem is solved directly on the first go, i.e. no `NumericalError` is encountered. So your hunch of it...

Hi Paul, Im doing it entirely in rust. I think I've gotten it to work after investigating it a little further. The assembly is actually working correctly. The problem ended...

At some point I would like to solve QPs with a direct solver that takes the `tril` rather than `triu`. I believe that the error in `count_diagonal_entries()` will become an...

You are correct. In short the function currently does *not* work for `tril` matrices, but it does not matter since `P` is assumed to be given as `triu` (I think...

Looks good. I agree that it makes it more clear and will avoid confusion in the future.

Pardiso only works for square systems. In your case what is happening is that it is actually solving the square system containing only first `m` columns. Maybe we should update...

Yeah, I followed the multiplication and found that a `_unwrap(A)` is what is causing this. I might be wrong, but I do not think the current behaviour is user-friendly. I've...