Hello,

Another (hopefully last) question, I am generating `randomMPS` like so:

```

sites

4-element Array{Index{Int64},1}:

(dim=2|id=138|"S=1/2,Site,n=1")

(dim=2|id=322|"S=1/2,Site,n=2")

(dim=2|id=961|"S=1/2,Site,n=3")

(dim=2|id=555|"S=1/2,Site,n=4")

psi = randomMPS(sites, 2)

MPS

[1] IndexSet{2} (dim=2|id=138|"S=1/2,Site,n=1") (dim=2|id=139|"Link,l=1")

[2] IndexSet{3} (dim=2|id=139|"Link,l=1") (dim=2|id=322|"S=1/2,Site,n=2") (dim=2|id=979|"Link,l=2")

[3] IndexSet{3} (dim=2|id=979|"Link,l=2") (dim=2|id=961|"S=1/2,Site,n=3") (dim=2|id=869|"Link,l=3")

[4] IndexSet{2} (dim=2|id=869|"Link,l=3") (dim=2|id=555|"S=1/2,Site,n=4")

```

But if I increase `linkdim` it looks like the last bond dimension stays to 2

```

julia> psi = randomMPS(sites, 3)

MPS

[1] IndexSet{2} (dim=2|id=138|"S=1/2,Site,n=1") (dim=3|id=896|"Link,l=1")

[2] IndexSet{3} (dim=3|id=896|"Link,l=1") (dim=2|id=322|"S=1/2,Site,n=2") (dim=3|id=732|"Link,l=2")

[3] IndexSet{3} (dim=3|id=732|"Link,l=2") (dim=2|id=961|"S=1/2,Site,n=3") (dim=2|id=271|"Link,l=3")

[4] IndexSet{2} (dim=2|id=271|"Link,l=3") (dim=2|id=555|"S=1/2,Site,n=4")

julia> psi = randomMPS(sites, 4)

MPS

[1] IndexSet{2} (dim=2|id=138|"S=1/2,Site,n=1") (dim=4|id=261|"Link,l=1")

[2] IndexSet{3} (dim=4|id=261|"Link,l=1") (dim=2|id=322|"S=1/2,Site,n=2") (dim=4|id=677|"Link,l=2")

[3] IndexSet{3} (dim=4|id=677|"Link,l=2") (dim=2|id=961|"S=1/2,Site,n=3") (dim=2|id=418|"Link,l=3")

[4] IndexSet{2} (dim=2|id=418|"Link,l=3") (dim=2|id=555|"S=1/2,Site,n=4")

julia> psi = randomMPS(sites, 5)

MPS

[1] IndexSet{2} (dim=2|id=138|"S=1/2,Site,n=1") (dim=5|id=893|"Link,l=1")

[2] IndexSet{3} (dim=5|id=893|"Link,l=1") (dim=2|id=322|"S=1/2,Site,n=2") (dim=4|id=986|"Link,l=2")

[3] IndexSet{3} (dim=4|id=986|"Link,l=2") (dim=2|id=961|"S=1/2,Site,n=3") (dim=2|id=62|"Link,l=3")

[4] IndexSet{2} (dim=2|id=62|"Link,l=3") (dim=2|id=555|"S=1/2,Site,n=4")

julia> psi = randomMPS(sites, 6)

MPS

[1] IndexSet{2} (dim=2|id=138|"S=1/2,Site,n=1") (dim=6|id=226|"Link,l=1")

[2] IndexSet{3} (dim=6|id=226|"Link,l=1") (dim=2|id=322|"S=1/2,Site,n=2") (dim=4|id=841|"Link,l=2")

[3] IndexSet{3} (dim=4|id=841|"Link,l=2") (dim=2|id=961|"S=1/2,Site,n=3") (dim=2|id=562|"Link,l=3")

[4] IndexSet{2} (dim=2|id=562|"Link,l=3") (dim=2|id=555|"S=1/2,Site,n=4")

julia> psi = randomMPS(sites, 7)

MPS

[1] IndexSet{2} (dim=2|id=138|"S=1/2,Site,n=1") (dim=7|id=5|"Link,l=1")

[2] IndexSet{3} (dim=7|id=5|"Link,l=1") (dim=2|id=322|"S=1/2,Site,n=2") (dim=4|id=433|"Link,l=2")

[3] IndexSet{3} (dim=4|id=433|"Link,l=2") (dim=2|id=961|"S=1/2,Site,n=3") (dim=2|id=586|"Link,l=3")

[4] IndexSet{2} (dim=2|id=586|"Link,l=3") (dim=2|id=555|"S=1/2,Site,n=4")

julia> psi = randomMPS(sites, 8)

MPS

[1] IndexSet{2} (dim=2|id=138|"S=1/2,Site,n=1") (dim=8|id=120|"Link,l=1")

[2] IndexSet{3} (dim=8|id=120|"Link,l=1") (dim=2|id=322|"S=1/2,Site,n=2") (dim=4|id=831|"Link,l=2")

[3] IndexSet{3} (dim=4|id=831|"Link,l=2") (dim=2|id=961|"S=1/2,Site,n=3") (dim=2|id=292|"Link,l=3")

[4] IndexSet{2} (dim=2|id=292|"Link,l=3") (dim=2|id=555|"S=1/2,Site,n=4")

julia> psi = randomMPS(sites, 9)

MPS

[1] IndexSet{2} (dim=2|id=138|"S=1/2,Site,n=1") (dim=8|id=845|"Link,l=1")

[2] IndexSet{3} (dim=8|id=845|"Link,l=1") (dim=2|id=322|"S=1/2,Site,n=2") (dim=4|id=910|"Link,l=2")

[3] IndexSet{3} (dim=4|id=910|"Link,l=2") (dim=2|id=961|"S=1/2,Site,n=3") (dim=2|id=214|"Link,l=3")

[4] IndexSet{2} (dim=2|id=214|"Link,l=3") (dim=2|id=555|"S=1/2,Site,n=4")

```

Is there a limit or something ?

Thanks for your help

BW,

Roland