@@ -81,8 +81,9 @@ and demonstrates how matrix-vector multiplication can be expressed in index
8181notation. Both forms are supported by TACO:
8282
8383``` python
84+ i, j = pytaco.get_index_vars(2 )
8485y[i] = A[i,j] * x[j]
85- y[i] = sum (j, A[i,j] * x[j])
86+ y[i] = pytaco. sum(j, A[i,j] * x[j])
8687```
8788
8889Reductions that are not explicitly expressed are assumed to be over the
@@ -135,6 +136,7 @@ stores the result in matrix `A`:
135136A = pt.tensor([3 ,3 ], pt.format([dense, dense]))
136137B = pt.tensor([3 ,3 ], pt.format([dense, dense]))
137138C = pt.tensor([3 ,3 ], pt.format([dense, dense]))
139+ i, j = pt.get_index_vars(2 )
138140
139141A[i,j] = B[i,j] + C[j,i]
140142```
@@ -153,6 +155,7 @@ before `i`:
153155A = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
154156B = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
155157C = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
158+ i, j = pt.get_index_vars(2 )
156159
157160A[i,j] = B[i,j] + C[j,i]
158161```
@@ -165,6 +168,7 @@ addition with the already-transposed temporary:
165168A = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
166169B = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
167170C = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
171+ i, j = pt.get_index_vars(2 )
168172
169173Ct = C.transpose([1 , 0 ]) # Ct is also stored in the CSR format
170174A[i,j] = B[i,j] + Ct[i,j]
@@ -179,6 +183,7 @@ order, requires `j` to be accessed before `i`:
179183A = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
180184B = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
181185C = pt.tensor([3 ,3 ], pt.format([dense, compressed], [1 , 0 ]))
186+ i, j = pt.get_index_vars(2 )
182187
183188A[i,j] = B[i,j] + C[i,j]
184189```
@@ -191,6 +196,7 @@ addition:
191196A = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
192197B = pt.tensor([3 ,3 ], pt.format([dense, compressed]))
193198C = pt.tensor([3 ,3 ], pt.format([dense, compressed], [1 , 0 ]))
199+ i, j = pt.get_index_vars(2 )
194200
195201Cp = C.transpose([0 , 1 ], pt.format([dense, compressed])) # Store a copy of C in the CSR format
196202A[i,j] = B[i,j] + Cp[i,j]
0 commit comments