EigenSparseMatrixTest Class Reference

Inheritance diagram for EigenSparseMatrixTest:

Public Member Functions
	EigenSparseMatrixTest ()
void	setUp ()
	CPPUNIT_TEST_SUITE (EigenSparseMatrixTest)
SPARSEMATRIXTEST	CPPUNIT_TEST_SUITE_END ()
void	tearDown ()
void	setValues ()
void	testValues ()
void	testGetAndSet ()
void	testReadMatlab (const std::string &filename)
void	testReadMatlab1 ()
void	testReadMatlab2 ()
void	testReadMatlab4 ()
void	testReadHDF5 ()
void	testTransposeNorms ()
void	testWriteAndRead ()
void	testClone ()

Protected Attributes
std::string	libmesh_suite_name
libMesh::Parallel::Communicator *	my_comm
std::unique_ptr< EigenSparseMatrix< Number > >	matrix
libMesh::numeric_index_type	nonsquare
libMesh::numeric_index_type	local_m
libMesh::numeric_index_type	local_n
libMesh::numeric_index_type	global_m
libMesh::numeric_index_type	global_n
const libMesh::Real	_tolerance

Private Attributes
Parallel::Communicator	comm_self

Detailed Description

Definition at line 11 of file eigen_sparse_matrix_test.C.

Constructor & Destructor Documentation

EigenSparseMatrixTest()

EigenSparseMatrixTest::EigenSparseMatrixTest ( )

inline

Definition at line 14 of file eigen_sparse_matrix_test.C.

References libMesh::PerfLog::summarized_logs_enabled(), and unitlog.

                          :
    SparseMatrixTest<EigenSparseMatrix<Number>>() {
    if (unitlog->summarized_logs_enabled())
      this->libmesh_suite_name = "SparseMatrixTest";
    else
      this->libmesh_suite_name = "EigenSparseMatrixTest";
  }

Member Function Documentation

CPPUNIT_TEST_SUITE()

EigenSparseMatrixTest::CPPUNIT_TEST_SUITE

(

EigenSparseMatrixTest

)

CPPUNIT_TEST_SUITE_END()

SPARSEMATRIXTEST EigenSparseMatrixTest::CPPUNIT_TEST_SUITE_END

(

)

setUp()

void EigenSparseMatrixTest::setUp ( )

inline

Definition at line 22 of file eigen_sparse_matrix_test.C.

  {
    // EigenSparseMatrix is serial; we'll tell it to use MPI_COMM_SELF
    // so we just do these tests embarrassingly parallel
    my_comm = &comm_self;

    // EigenSparseMatrix doesn't support non-square matrices?
    nonsquare = 0;

    SparseMatrixTest<EigenSparseMatrix<Number>>::setUp();
  }

setValues()

void SparseMatrixTest< EigenSparseMatrix< Number > >::setValues ( )

inlineinherited

Definition at line 79 of file sparse_matrix_test.h.

  {
    std::vector<libMesh::numeric_index_type> rows(local_m);
    std::iota(rows.begin(), rows.end(), matrix->row_start());
    std::vector<libMesh::numeric_index_type> cols(local_n);
    std::iota(cols.begin(), cols.end(), matrix->col_start());

    libMesh::DenseMatrix<libMesh::Number> local(local_m, local_n);

    for (auto i : libMesh::make_range(local_m))
      for (auto j : libMesh::make_range(local_n))
        local(i, j) = (i + 1) * (j + 1) * (my_comm->rank() + 1);

    matrix->zero();

    matrix->add_matrix(local, rows, cols);
    matrix->close();
  }

tearDown()

void SparseMatrixTest< EigenSparseMatrix< Number > >::tearDown ( )

inlineinherited

Definition at line 76 of file sparse_matrix_test.h.

{}

testClone()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testClone ( )

inlineinherited

Definition at line 259 of file sparse_matrix_test.h.

  {
    LOG_UNIT_TEST;

    setValues();

    // Matrix must be closed before it can be cloned.
    matrix->close();

    {
      // Create copy, test that it can go out of scope
      auto copy = matrix->clone();

      // Check that matrices have the same local/global sizes
      CPPUNIT_ASSERT_EQUAL(copy->m(), matrix->m());
      CPPUNIT_ASSERT_EQUAL(copy->n(), matrix->n());
      CPPUNIT_ASSERT_EQUAL(copy->local_m(), matrix->local_m());
      CPPUNIT_ASSERT_EQUAL(copy->row_start(), matrix->row_start());
      CPPUNIT_ASSERT_EQUAL(copy->row_stop(), matrix->row_stop());

      // Check that copy has same values as original
      LIBMESH_ASSERT_FP_EQUAL(copy->l1_norm(), matrix->l1_norm(), _tolerance);
      CPPUNIT_ASSERT(relative_fuzzy_equals(*matrix, *copy));
      copy->scale(2);
      CPPUNIT_ASSERT(!relative_fuzzy_equals(*matrix, *copy));
    }

    {
      // Create zero copy
      auto zero_copy = matrix->zero_clone();

      // Check that matrices have the same local/global sizes
      CPPUNIT_ASSERT_EQUAL(zero_copy->m(), matrix->m());
      CPPUNIT_ASSERT_EQUAL(zero_copy->n(), matrix->n());
      CPPUNIT_ASSERT_EQUAL(zero_copy->local_m(), matrix->local_m());
      CPPUNIT_ASSERT_EQUAL(zero_copy->row_start(), matrix->row_start());
      CPPUNIT_ASSERT_EQUAL(zero_copy->row_stop(), matrix->row_stop());

      // Check that zero_copy has same values as original
      LIBMESH_ASSERT_FP_EQUAL(0.0, zero_copy->l1_norm(), _tolerance);
    }
  }

testGetAndSet()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testGetAndSet ( )

inlineinherited

Definition at line 139 of file sparse_matrix_test.h.

  {
    LOG_UNIT_TEST;

    setValues();

    testValues();
  }

testReadHDF5()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testReadHDF5 ( )

inlineinherited

Definition at line 194 of file sparse_matrix_test.h.

  {
#ifdef LIBMESH_HAVE_HDF5
    LOG_UNIT_TEST;

    matrix->clear();
    auto matrix2 = std::make_unique<DerivedClass>(*my_comm);
    matrix->read("matrices/geom_1_extraction_op.m");
    matrix2->read("matrices/geom_1_extraction_op.h5");

    // We need some more SparseMatrix operators, but not today
    CPPUNIT_ASSERT(matrix->l1_norm() == matrix2->l1_norm());
    CPPUNIT_ASSERT(matrix->linfty_norm() == matrix2->linfty_norm());
#endif // LIBMESH_HAVE_HDF5
  }

testReadMatlab()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testReadMatlab ( const std::string &  filename )

inlineinherited

Definition at line 149 of file sparse_matrix_test.h.

  {
    // Laspack doesn't handle non-square matrices)
    if (matrix->solver_package() == libMesh::LASPACK_SOLVERS)
      return;

    matrix->clear();

    auto matrix2 = std::make_unique<DerivedClass>(*my_comm);

    matrix->read(filename);

#ifndef LIBMESH_HAVE_GZSTREAM
    return;
#endif

    matrix2->read(std::string(filename)+".gz");

    // We need some more SparseMatrix operators, but not today
    CPPUNIT_ASSERT(matrix->l1_norm() == matrix2->l1_norm());
    CPPUNIT_ASSERT(matrix->linfty_norm() == matrix2->linfty_norm());
  }

testReadMatlab1()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testReadMatlab1 ( )

inlineinherited

Definition at line 173 of file sparse_matrix_test.h.

  {
    LOG_UNIT_TEST;
    testReadMatlab("matrices/geom_1_extraction_op.m");
  }

testReadMatlab2()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testReadMatlab2 ( )

inlineinherited

Definition at line 180 of file sparse_matrix_test.h.

  {
    LOG_UNIT_TEST;
    testReadMatlab("matrices/geom_2_extraction_op.m");
  }

testReadMatlab4()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testReadMatlab4 ( )

inlineinherited

Definition at line 187 of file sparse_matrix_test.h.

  {
    LOG_UNIT_TEST;
    testReadMatlab("matrices/geom_4_extraction_op.m");
  }

testTransposeNorms()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testTransposeNorms ( )

inlineinherited

Definition at line 211 of file sparse_matrix_test.h.

  {
    LOG_UNIT_TEST;

    setValues();

    auto matrix2 = std::make_unique<DerivedClass>(*my_comm);
    matrix->get_transpose(*matrix2);
    LIBMESH_ASSERT_FP_EQUAL(matrix->l1_norm(), matrix2->linfty_norm(), _tolerance);
    LIBMESH_ASSERT_FP_EQUAL(matrix2->l1_norm(), matrix->linfty_norm(), _tolerance);
  }

testValues()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testValues ( )

inlineinherited

Definition at line 99 of file sparse_matrix_test.h.

  {
    auto functor = [this]()
    {
      std::vector<libMesh::numeric_index_type> cols_to_get;
      std::vector<libMesh::Number> values;
      const libMesh::numeric_index_type col_start =
        matrix->col_start();
      for (libMesh::numeric_index_type i :
           libMesh::make_range(matrix->row_start(),
                               matrix->row_stop()))
        {
          matrix->get_row(i, cols_to_get, values);
          for (libMesh::numeric_index_type col_j :
               libMesh::index_range(cols_to_get))
            {
              CPPUNIT_ASSERT_EQUAL(cols_to_get[col_j], col_start + col_j);
              LIBMESH_ASSERT_NUMBERS_EQUAL
                ((i - matrix->row_start() + 1) * (col_j + 1) * (my_comm->rank() + 1),
                 values[col_j], _tolerance);
            }
        }
    };

#ifdef LIBMESH_HAVE_CXX11_THREAD
    auto num_threads = std::min(unsigned(2),
                                std::max(
                                  std::thread::hardware_concurrency(),
                                  unsigned(1)));
    std::vector<std::thread> threads(num_threads);
    for (unsigned int thread = 0; thread < num_threads; ++thread)
      threads[thread] = std::thread(functor);
    std::for_each(threads.begin(), threads.end(),
                  [](std::thread & x){x.join();});
#else
    functor();
#endif
  }

testWriteAndRead()

void SparseMatrixTest< EigenSparseMatrix< Number > >::testWriteAndRead ( )

inlineinherited

Definition at line 224 of file sparse_matrix_test.h.

  {
    LOG_UNIT_TEST;

    setValues();

    // If we're working with serial matrices then just print one of
    // them so they don't step on the others' toes.
    //
    // Use a very short filename, because we had a bug with that and
    // we want to test it.
    if (matrix->n_processors() > 1 ||
        TestCommWorld->rank() == 0)
      matrix->print_matlab("M.m");

    matrix->clear();

    // Let's make sure we don't have any race conditions; we have
    // multiple SparseMatrix subclasses that might be trying to read
    // and write the same file.

    TestCommWorld->barrier();

#ifdef LIBMESH_USE_COMPLEX_NUMBERS
    // We're not supporting complex reads quite yet
    return;
#endif

    matrix->read("M.m");

    TestCommWorld->barrier();

    testValues();
  }

Member Data Documentation

_tolerance

const libMesh::Real SparseMatrixTest< EigenSparseMatrix< Number > >::_tolerance

protectedinherited

Definition at line 314 of file sparse_matrix_test.h.

comm_self

Parallel::Communicator EigenSparseMatrixTest::comm_self

private

Definition at line 42 of file eigen_sparse_matrix_test.C.

global_m

libMesh::numeric_index_type SparseMatrixTest< EigenSparseMatrix< Number > >::global_m

protectedinherited

Definition at line 310 of file sparse_matrix_test.h.

global_n

libMesh::numeric_index_type SparseMatrixTest< EigenSparseMatrix< Number > >::global_n

protectedinherited

Definition at line 310 of file sparse_matrix_test.h.

libmesh_suite_name

std::string SparseMatrixTest< EigenSparseMatrix< Number > >::libmesh_suite_name

protectedinherited

Definition at line 304 of file sparse_matrix_test.h.

local_m

libMesh::numeric_index_type SparseMatrixTest< EigenSparseMatrix< Number > >::local_m

protectedinherited

Definition at line 310 of file sparse_matrix_test.h.

local_n

libMesh::numeric_index_type SparseMatrixTest< EigenSparseMatrix< Number > >::local_n

protectedinherited

Definition at line 310 of file sparse_matrix_test.h.

matrix

std::unique_ptr<EigenSparseMatrix< Number > > SparseMatrixTest< EigenSparseMatrix< Number > >::matrix

protectedinherited

Definition at line 308 of file sparse_matrix_test.h.

my_comm

libMesh::Parallel::Communicator* SparseMatrixTest< EigenSparseMatrix< Number > >::my_comm

protectedinherited

Definition at line 306 of file sparse_matrix_test.h.

nonsquare

libMesh::numeric_index_type SparseMatrixTest< EigenSparseMatrix< Number > >::nonsquare

protectedinherited

Definition at line 310 of file sparse_matrix_test.h.

---

The documentation for this class was generated from the following file:

• tests/numerics/eigen_sparse_matrix_test.C