The LLT polynomials $\mathcal{L}_{\mathbfβ/\mathbfγ} (X;t)$ are a family of symmetric polynomials indexed by a tuple of (possibly skew-)partitions $\mathbfβ/\mathbfγ= (β^{(1)}/γ^{(1)},\ldots,β^{(k)}/γ^{(k)})$. It has recently been shown that these polynomials can be seen as the partition function of a certain vertex model whose boundary conditions are determined by $\mathbfβ/\mathbfγ$. In this paper we describe an algorithm which gives a bijection between the configurations of the vertex model with boundary condition $\mathbfβ/\mathbfγ = (β^{(1)}/γ^{(1)},β^{(2)}/γ^{(2)})$ and those with boundary condition $(\mathbfβ/\mathbfγ)_{swap} = (β^{(2)}/γ^{(2)},β^{(1)}/γ^{(1)})$. We prove a sufficient condition for when this bijection is weight-preserving up to an overall factor of $t$, which in turn implies that the corresponding LLT polynomials are equal up to the same overall factor. Using these techniques, we are also able to systematically determine linear relations within families of LLT polynomials.