Learn to Use ITensor

main / book / intro C++v3 | C++v2

ITensor Library Overview

Matrix techniques have been a major success in the applied sciences. Tensor methods could lead to even deeper insights and more powerful algorithms. Tensor network methods in physics and chemistry have yielded very precise results and greater understanding of many-body wavefunctions. Mathematicians and computer scientists are using tensor algorithms to uncover correlations in data and to scale algorithms to high-dimensional spaces.

When a problem is cast into the language of tensors, each tensor index typically has some semantic, or physical meaning. Similarly in ITensor, tensor indices are objects with unique identities. Constructing an index in ITensor looks like:

auto i = Index("index i",3);

The resulting Index has a size, or dimension of 3. Upon creation it is indelibly stamped with an internal ID number so all copies of this Index can recognize each other.

Having made a few Index objects i, j, k, one can construct ITensors

auto A = ITensor(i);
auto B = ITensor(j,i);
auto C = ITensor(i,j,k);

and set their elements (as shown in the ITensor Basics chapter).

Since matching indices can recognize each other, ITensor contraction is simply

D = A * B * C;

The ITensor contraction engine recognizes repeated indices and sums over them, like the Einstein summation convention used in physics. Instead of bookkeeping the order of tensor indices, users can focus on the structure of tensor networks. Certain bugs are ruled out, such as mistaking one index for another of the same size.

Let's take a closer look at an example involving two matrix-like ITensors

auto A = ITensor(i,j);
auto B = ITensor(k,j);
auto C = A * B;

Here the contraction A * B computes the matrix product @@ C = A B^\mathsf{T} @@ . The Index j appears on both A and B so is automatically summed over, leaving C to have indices i and k.

In a standard matrix library, one would need to remember that j is the second index of B and write something like

C = A * transpose(B) //not actual ITensor code!!

to get the correct result. ITensor handles this transposition automatically, making user code simple and robust. If B were to be redefined with transposed index order, the ITensor operation A * B would continue to give the correct result.

Say we didn't want to contract over the index j — we could do the following

auto D = prime(A,j) * B;

The call to prime(A,j) puts a single, lightweight prime on A's copy of index j. It no longer matches the copy of j on B so these will not be contracted and D will end up with four indices i,j',j,k. We have computed an outer product of two matrices with very little effort.

ITensor offers many other features to enhance productivity:

Last but not least, ITensor includes routines for tensor decompositions and a rich, high-level layer for matrix product state and density matrix renormalization group (DMRG) algorithms used in physics and chemistry.

Table of Contents Index Objects



Back to Book
Back to Main