LHerviou

I am also having similar issues on a linux cluster (code runs fine on my old julia install on the laptop). I put some data in https://discourse.julialang.org/t/segmentation-fault-using-multithreaded-julia-on-new-server/114557 Not able to...

Hi, Is it possible that the contributions to S smaller than 10^-3 are already present in `psi` before the subspace expansion? The current implementation (even before my modifications) does not...

It should do a shallow copy of the tensor if I am not mistaken. I will try to investigate a bit. It is really a marginal cost most of the...

I have started having a look at it, but if anyone wants to contribute, let me know. As of now, I am actually struggling on something simple. It should be...

Yes, as a single translation invariant MPO. I am really (for now) just trying to create the object W_inf. I am trying to write it automatically for arbitrary (potentially large)...

Nevermind, I think I fixed that problem: would taking InfiniteMPO.data to be CelledVector{Matrix{ITensor}} be a satisfactory solution?

I will start like that then, it should be easy to make the basic infrastructure work. Thanks for the additional context.

For info, I managed to write a full MPO treatment of VUMPS. The subspace expansion could be improved, but I would like to talk to several people in Benasque before...

Well, the original proposition for the subspace expansion is indeed just looking at the two site effective Hamiltonian. That is what I am currently using in the code. On the...

The method can be in principle also used to find the best way to extend the system. In fact, the Ghent team uses a variation of that proposal to do...