Mason Protter
Mason Protter
Something is causing this third case to be unhappy: ```julia julia> using LoopVectorization julia> let a = 1.0, b = 1.0, h = 1.0, v = rand(10) @avx for i...
Suppose I have a function ```julia function pairwise(a) ret = Matrix{eltype(a)}(undef, length(a), length(a)) for i in 1:(length(a)-1) for j in (i+1):length(a) ret[i, j] = a[i] * a[j] end end UpperTriangular(ret)...
**Description of bug** I tried to run the code from https://itensor.github.io/ITensors.jl/dev/tutorials/QN_DMRG.html on a GPU using ITensorGPU.jl and encountered an error where the dimensions of the dims of the tensors in...
Fixes #771 Okay, so there's a few ways to go about this, and I'm trying one way here but I'm open to changing direction. This makes it so when a...
The current printing of `MPS` (and other ITensors objects) will blow up your display with many newlines and output that scales forever. E.g. compare ```julia julia> [randomMPS(siteinds("S=1/2", 100)) for _...
`randomMPS` fails to create states with `linkdim>1` on heterogenous lattice when qns are conserved
Basically, if I have a lattice where alternating sites have different sitetypes, `randomMPS(s, f; linkdims=k)` for `k>1` will end up erroring. Here is a MWE: ```julia using ITensors function ITensors.space(::SiteType"Sf=1/2";...
[ChainRules.jl](https://github.com/JuliaDiff/ChainRules.jl) just got started to replace DiffRules.jl. Would be a good idea to switch over once the API stabilizes.
Right now, `==` is used to tell if two symbolic expressions are the same. Ie. ```julia using Symbolics @sym x y julia> x*y == y*x false julia> cos(x)^2 + sin(x)^2...
Right now, the simplification algorithms found in src/Simplifications.jl will in theory work on any subtype of `AbstractSymExpr` but many of them are constructed so that they spit out a `SymExpr`...
Now that we can generate equations of motion from a Lagrangian (should work in theory for Hamiltonians too but I haven't tried), it'd be really slick if we could pipe...