Luca Gravina
Luca Gravina
I am looking for a way to rewind the integrator to the state it had in a previous time. Specifically, I want to use a continuous callback to monitor some...
Multiplying a CuSparseMatrixCSC with a CuArray gives Scalar indexing. To reproduce: ``` CUDA.allowscalar(false) A = cu(sparse([1,2,3], [1,2,3], [1,2,3])) B = cu(rand(3,1)) C = A*B ``` or ``` CUDA.allowscalar(false) A =...
Schur decomposition is implemented efficiently on GPU in PyTorch, Jax, and MatLab. It is an essential ingredient for the efficient computation of matrix square root (sqrtm). It would be ideal...
Hello. I have a function which returns values between -2 and -14. I do not know to optimal value a priori. Using bboptimize converges towards what seems to be the...
Does this solver natively support Ensemble threading?
In the article published on Quantum state-dependent running costs are taken into account. How can this be done using this (collection of) package?
Is there any way of inserting additional collapse operators such as: `np.cos(theta)*Jminus + np.sin(theta)*Jplus` ?
Add `@async` computations in GPU implementation of the time evolution. Benchmark GPU code with NVIDIA Nsight.
Implement symmetry-decomposed evolution for U(1) strong Liouvillian symmetry. Based on [this](https://arxiv.org/pdf/2010.08492) paper.