Implements support for more corner cases involving ToT times T.

.github/workflows/ci.yml

@@ -13,10 +13,10 @@ jobs:
       fail-fast: false
       matrix:
         os : [ macos-latest ]
-        cxx : [ clang++, /usr/local/bin/g++-10 ]
+        cxx : [ clang++, /opt/homebrew/bin/g++-11 ]
         build_type : [ Release, Debug ]
         task_backend: [ Pthreads, PaRSEC ]
-        prerequisites : [ gcc@10 boost eigen open-mpi bison scalapack ]
+        prerequisites : [ gcc@11 boost eigen open-mpi bison scalapack ]
 
     name: "${{ matrix.os }}: ${{ matrix.cxx }} ${{ matrix.build_type }} ${{ matrix.task_backend }}"
     runs-on: ${{ matrix.os }}
@@ -35,7 +35,7 @@ jobs:
 
     - uses: maxim-lobanov/setup-xcode@v1
       with:
-        xcode-version: '<14'
+        xcode-version: 'latest-stable'
 
     - name: Host system info
       shell: bash

src/TiledArray/conversions/dense_to_sparse.h

@@ -27,7 +27,7 @@ to_sparse(DistArray<Tile, ArgPolicy> const &dense_array) {
   const auto begin = dense_array.begin();
   for (auto it = begin; it != end; ++it) {
     // write the norm of each local tile to the tensor
-    norm(it->get(), tile_norms[it.ordinal()]);
+    norm(it->get(), tile_norms[it.index()]);
   }
 
   // Construct a sparse shape the constructor will handle communicating the
@@ -40,9 +40,9 @@ to_sparse(DistArray<Tile, ArgPolicy> const &dense_array) {
   // sparse_array set the sparse array tile with a clone so as not to hold
   // a pointer to the original tile.
   for (auto it = begin; it != end; ++it) {
-    const auto ord = it.ordinal();
-    if (!sparse_array.is_zero(ord)) {
-      sparse_array.set(ord, it->get().clone());
+    const auto ix = it.index();
+    if (!sparse_array.is_zero(ix)) {
+      sparse_array.set(ix, it->get().clone());
     }
   }
 

src/TiledArray/conversions/sparse_to_dense.h

@@ -52,8 +52,7 @@ to_dense(DistArray<Tile, ArgPolicy> const& sparse_array) {
       Tile tile(sparse_array.find(ord).get().clone());
       dense_array.set(ord, tile);
     } else {
-      // see DistArray::set(ordinal, element_type)
-      dense_array.set(ord, 0);
+      dense_array.set(ord, typename Tile::value_type{});
     }
   }