Inheritance diagram for EigenSparseMatrixTest:

[legend]

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 (const std::string &filename)

void	testWriteAndReadHDF5 ()

void	testWriteAndReadMatlab ()

void	testWriteAndReadZippedMatlab ()

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.

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

References SparseMatrixTest< EigenSparseMatrix< Number > >::libmesh_suite_name, libMesh::PerfLog::summarized_logs_enabled(), and unitlog.

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();
  }

References comm_self, SparseMatrixTest< EigenSparseMatrix< Number > >::my_comm, SparseMatrixTest< EigenSparseMatrix< Number > >::nonsquare, and setUp().

Referenced by setUp().

◆ setValues()

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

inlineinherited

Definition at line 93 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 90 of file sparse_matrix_test.h.

90{}

◆ testClone()

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

inlineinherited

Definition at line 296 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 153 of file sparse_matrix_test.h.

  {
    LOG_UNIT_TEST;
 
    setValues();
 
    testValues();
  }

◆ testReadHDF5()

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

inlineinherited

Definition at line 209 of file sparse_matrix_test.h.

  {
    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());
  }

◆ testReadMatlab()

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

inlineinherited

Definition at line 163 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 187 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 194 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 201 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 225 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 113 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 ( const std::string & filename )

inlineinherited

Definition at line 238 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.
    if (matrix->n_processors() > 1 ||
        TestCommWorld->rank() == 0)
      matrix->print(filename);
 
    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(filename);
 
    TestCommWorld->barrier();
 
    testValues();
  }

◆ testWriteAndReadHDF5()

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

inlineinherited

Definition at line 271 of file sparse_matrix_test.h.

  {
    // This capability is not yet supported when libmesh is compiled
    // with complex number support
#ifndef LIBMESH_USE_COMPLEX_NUMBERS
    testWriteAndRead("M.h5");
#endif
  }

◆ testWriteAndReadMatlab()

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

inlineinherited

Definition at line 282 of file sparse_matrix_test.h.

  {
    // Use a very short filename, because we had a bug with that and
    // we want to test it.
    testWriteAndRead("M.m");
  }

◆ testWriteAndReadZippedMatlab()

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

inlineinherited

Definition at line 289 of file sparse_matrix_test.h.

  {
#ifdef LIBMESH_HAVE_GZSTREAM
    testWriteAndRead("Mzipped.m.gz");
#endif
  }

Member Data Documentation

◆ _tolerance

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

protectedinherited

Definition at line 351 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.

Referenced by setUp().

◆ global_m

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

protectedinherited

Definition at line 349 of file sparse_matrix_test.h.

◆ global_n

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

protectedinherited

Definition at line 349 of file sparse_matrix_test.h.

◆ libmesh_suite_name

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

protectedinherited

Definition at line 341 of file sparse_matrix_test.h.

◆ local_m

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

protectedinherited

Definition at line 348 of file sparse_matrix_test.h.

◆ local_n

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

protectedinherited

Definition at line 348 of file sparse_matrix_test.h.

◆ matrix

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

protectedinherited

Definition at line 345 of file sparse_matrix_test.h.

◆ my_comm

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

protectedinherited

Definition at line 343 of file sparse_matrix_test.h.

◆ nonsquare

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

protectedinherited

Definition at line 347 of file sparse_matrix_test.h.

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

tests/numerics/eigen_sparse_matrix_test.C

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ EigenSparseMatrixTest()

Member Function Documentation

◆ CPPUNIT_TEST_SUITE()

◆ CPPUNIT_TEST_SUITE_END()

◆ setUp()

◆ setValues()

◆ tearDown()

◆ testClone()

◆ testGetAndSet()

◆ testReadHDF5()

◆ testReadMatlab()

◆ testReadMatlab1()

◆ testReadMatlab2()

◆ testReadMatlab4()

◆ testTransposeNorms()

◆ testValues()

◆ testWriteAndRead()

◆ testWriteAndReadHDF5()

◆ testWriteAndReadMatlab()

◆ testWriteAndReadZippedMatlab()

Member Data Documentation

◆ _tolerance

◆ comm_self

◆ global_m

◆ global_n

◆ libmesh_suite_name

◆ local_m

◆ local_n

◆ matrix

◆ my_comm

◆ nonsquare