sparse.cc File Reference

#include "purify/config.h"
#include "purify/types.h"
#include <iostream>
#include "catch2/catch_all.hpp"
#include "purify/directories.h"
#include "purify/logging.h"

Include dependency graph for sparse.cc:

Go to the source code of this file.

Functions
	TEST_CASE ("sparse [iterator]", "[iterator]")

Function Documentation

TEST_CASE()

TEST_CASE	(	"sparse "	[iterator],
		""	[iterator]
	)

Definition at line 9 of file sparse.cc.

                                             {
  t_int NZnum = 40;
  std::vector<t_tripletList> entries;
  entries.reserve(NZnum);
  Vector<t_int> rows = Vector<t_int>::LinSpaced(NZnum, 0, 2 * NZnum);
  Vector<t_int> cols = Vector<t_int>::LinSpaced(NZnum, 0, 3 * NZnum);
 
  for (t_int i = 0; i < NZnum; ++i) {
    entries.emplace_back(rows(i), cols(i), i);
  }
  const t_int total_rows = 5 * NZnum;
  const t_int total_cols = 4 * NZnum;
  Sparse<t_complex> G(total_rows, total_cols);
  G.setFromTriplets(entries.begin(), entries.end());
  for (t_int i = 0; i < total_cols; ++i) {
    for (Sparse<t_complex>::InnerIterator pix(G, i); pix; ++pix) {
      PURIFY_HIGH_LOG("{} {} {}", pix.row(), pix.col(), std::real(pix.value()));
    }
  }
  // CHECK(expected.isApprox(output, 1e-13));
}

References PURIFY_HIGH_LOG.