Mikkel Paltorp
Mikkel Paltorp
Switching to ```julia @time @. u = uprev + (dt / 6.0) * (2.0 * (k₂ + k₃) + (k₁ + k₄)) ``` Seem to speedup things quite a bit?...
I have a similar-ish problem where the matrices being blocked are Toeplitz. A little contrived problem is e.g. ```julia using ToeplitzMatrices, BlockArrays N = 10000 vc = 1:N vr =...
Hi again Boris, I found that the reason the complex-solve did not work was simply due to the way the memory of the output was allocated. It was easily fixed...
Ahh, I did not see that the existing PR dealt with the exact same issue 👍 I guess the only new thing about my PR then is the allowance of...
Hello Baptiste, Sounds good. Sorry for overlooking your PR! No need to stress, enjoy your vacation (or whatever thing that is keeping you away from your computer 😉 )
@j-fu Do you know if this code works for Pardiso8?
Thanks a lot for checking it out Jürgen. A shame that the free license only covers two cores 😑
Hi again, now I've gotten Pardiso 8 to work on apple silicon. However, as you mentioned Jürgen the test fail for hermitian indef and symmetric. I think that there might...
I will be in academia for at least one more year and would like to use this package meanwhile. Would be happy to help with the maintenance for now, but...
Haha, thanks @KristofferC , sounds like a good idea 👀 @j-fu I saw in the email from Olaf that he was asking for some Julia code that reproduces the issue....