Real Tensor
Quantum Physics
This TensorSuite collection includes block-wise, block-sparse tensors from density matrix renormalization group (DMRG) workloads, where tensor contractions naturally expose structured sparsity patterns [1], [2].
The tensors are generated from a 2D Hubbard model with a 6 x 3
lattice, 18 sites, and bond dimension up to 1600 using ITensors.jl
[3]. From the two-site DMRG update workflow, we include five
representative contraction systems: E1, E2, S1, S2, and S3.
The collection covers four symmetry regimes: kysznf, sznf, nf, and
nfparity.
For each contraction, TensorSuite lists both the canonical tensor layout and an internal permuted layout. The canonical layout follows the DMRG-generated index order, while the permuted layout groups contracted indices to expose larger block operations, a common strategy for improving tensor contraction performance on modern hardware [3], [4], [5].
Tensor Pairs
Each contraction system contains a tensor A and tensor B pair. The pairs are grouped by layout so links can be attached separately for the canonical formulation and the internal permuted layout. The tensor contractions are expressed in einsum notation below.
Canonical Formulation
Einsum equations:
E1: dea,fag->defg
E2: fbag,bhai->fghi
S1: dae,fga->defg
S2: dea,fagh->defgh
S3: fbghai,baj->fghij
| Symmetry | System | Tensor A | Tensor B |
|---|---|---|---|
kysznf |
E1 |
canonical_formulation_kysznf_E1_A [Link] |
canonical_formulation_kysznf_E1_B [Link] |
kysznf |
E2 |
canonical_formulation_kysznf_E2_A [Link] |
canonical_formulation_kysznf_E2_B [Link] |
kysznf |
S1 |
canonical_formulation_kysznf_S1_A [Link] |
canonical_formulation_kysznf_S1_B [Link] |
kysznf |
S2 |
canonical_formulation_kysznf_S2_A [Link] |
canonical_formulation_kysznf_S2_B [Link] |
kysznf |
S3 |
canonical_formulation_kysznf_S3_A [Link] |
canonical_formulation_kysznf_S3_B [Link] |
sznf |
E1 |
canonical_formulation_sznf_E1_A [Link] |
canonical_formulation_sznf_E1_B [Link] |
sznf |
E2 |
canonical_formulation_sznf_E2_A [Link] |
canonical_formulation_sznf_E2_B [Link] |
sznf |
S1 |
canonical_formulation_sznf_S1_A [Link] |
canonical_formulation_sznf_S1_B [Link] |
sznf |
S2 |
canonical_formulation_sznf_S2_A [Link] |
canonical_formulation_sznf_S2_B [Link] |
sznf |
S3 |
canonical_formulation_sznf_S3_A [Link] |
canonical_formulation_sznf_S3_B [Link] |
nf |
E1 |
canonical_formulation_nf_E1_A [Link] |
canonical_formulation_nf_E1_B [Link] |
nf |
E2 |
canonical_formulation_nf_E2_A [Link] |
canonical_formulation_nf_E2_B [Link] |
nf |
S1 |
canonical_formulation_nf_S1_A [Link] |
canonical_formulation_nf_S1_B [Link] |
nf |
S2 |
canonical_formulation_nf_S2_A [Link] |
canonical_formulation_nf_S2_B [Link] |
nf |
S3 |
canonical_formulation_nf_S3_A [Link] |
canonical_formulation_nf_S3_B [Link] |
nfparity |
E1 |
canonical_formulation_nfparity_E1_A [Link] |
canonical_formulation_nfparity_E1_B [Link] |
nfparity |
E2 |
canonical_formulation_nfparity_E2_A [Link] |
canonical_formulation_nfparity_E2_B [Link] |
nfparity |
S1 |
canonical_formulation_nfparity_S1_A [Link] |
canonical_formulation_nfparity_S1_B [Link] |
nfparity |
S2 |
canonical_formulation_nfparity_S2_A [Link] |
canonical_formulation_nfparity_S2_B [Link] |
nfparity |
S3 |
canonical_formulation_nfparity_S3_A [Link] |
canonical_formulation_nfparity_S3_B [Link] |
Internal Permuted Layout
Einsum equations:
E1 / E2 / S1: ade,afg->defg
S2: ade,afgh->defgh
S3: adefg,ah->defgh
| Symmetry | System | Tensor A | Tensor B |
|---|---|---|---|
kysznf |
E1 |
block_permuted_internal_only_kysznf_E1_A [Link] |
block_permuted_internal_only_kysznf_E1_B [Link] |
kysznf |
E2 |
block_permuted_internal_only_kysznf_E2_A [Link] |
block_permuted_internal_only_kysznf_E2_B [Link] |
kysznf |
S1 |
block_permuted_internal_only_kysznf_S1_A [Link] |
block_permuted_internal_only_kysznf_S1_B [Link] |
kysznf |
S2 |
block_permuted_internal_only_kysznf_S2_A [Link] |
block_permuted_internal_only_kysznf_S2_B [Link] |
kysznf |
S3 |
block_permuted_internal_only_kysznf_S3_A [Link] |
block_permuted_internal_only_kysznf_S3_B [Link] |
sznf |
E1 |
block_permuted_internal_only_sznf_E1_A [Link] |
block_permuted_internal_only_sznf_E1_B [Link] |
sznf |
E2 |
block_permuted_internal_only_sznf_E2_A [Link] |
block_permuted_internal_only_sznf_E2_B [Link] |
sznf |
S1 |
block_permuted_internal_only_sznf_S1_A [Link] |
block_permuted_internal_only_sznf_S1_B [Link] |
sznf |
S2 |
block_permuted_internal_only_sznf_S2_A [Link] |
block_permuted_internal_only_sznf_S2_B [Link] |
sznf |
S3 |
block_permuted_internal_only_sznf_S3_A [Link] |
block_permuted_internal_only_sznf_S3_B [Link] |
nf |
E1 |
block_permuted_internal_only_nf_E1_A [Link] |
block_permuted_internal_only_nf_E1_B [Link] |
nf |
E2 |
block_permuted_internal_only_nf_E2_A [Link] |
block_permuted_internal_only_nf_E2_B [Link] |
nf |
S1 |
block_permuted_internal_only_nf_S1_A [Link] |
block_permuted_internal_only_nf_S1_B [Link] |
nf |
S2 |
block_permuted_internal_only_nf_S2_A [Link] |
block_permuted_internal_only_nf_S2_B [Link] |
nf |
S3 |
block_permuted_internal_only_nf_S3_A [Link] |
block_permuted_internal_only_nf_S3_B [Link] |
nfparity |
E1 |
block_permuted_internal_only_nfparity_E1_A [Link] |
block_permuted_internal_only_nfparity_E1_B [Link] |
nfparity |
E2 |
block_permuted_internal_only_nfparity_E2_A [Link] |
block_permuted_internal_only_nfparity_E2_B [Link] |
nfparity |
S1 |
block_permuted_internal_only_nfparity_S1_A [Link] |
block_permuted_internal_only_nfparity_S1_B [Link] |
nfparity |
S2 |
block_permuted_internal_only_nfparity_S2_A [Link] |
block_permuted_internal_only_nfparity_S2_B [Link] |
nfparity |
S3 |
block_permuted_internal_only_nfparity_S3_A [Link] |
block_permuted_internal_only_nfparity_S3_B [Link] |
References
[1] U. Schollwock, "The density-matrix renormalization group in the age of matrix product states," Annals of Physics, 2011. DOI
[2] R. Levy, E. Solomonik, and B. K. Clark, "Distributed-memory DMRG via sparse and dense parallel tensor contractions," SC20, 2020. DOI
[3] M. Fishman, S. White, and E. Stoudenmire, "The ITensor software library for tensor network calculations," SciPost Physics Codebases, 2022. DOI
[4] P. Springer and P. Bientinesi, "Design of a high-performance GEMM-like tensor-tensor multiplication," ACM Transactions on Mathematical Software, 2018. DOI
[5] I. A. Kaliman and A. I. Krylov, "New algorithm for tensor contractions on multi-core CPUs, GPUs, and accelerators enables CCSD and EOM-CCSD calculations with over 1000 basis functions on a single compute node," Journal of Computational Chemistry, 2017. DOI
Acknowledgments
We thank Karl Pierce and Miles Stoudenmire for their help in gathering the data.