LumpedMassMatrixTest Class Reference

Inheritance diagram for LumpedMassMatrixTest:

Public Member Functions
	LIBMESH_CPPUNIT_TEST_SUITE (LumpedMassMatrixTest)
	CPPUNIT_TEST (testNumerics)
	CPPUNIT_TEST_SUITE_END ()
void	setUp ()
void	tearDown ()
void	testNumerics ()

Private Attributes
Parallel::Communicator *	_comm
std::unique_ptr< LumpedMassMatrix< Number > >	_matrix
numeric_index_type	_local_size
numeric_index_type	_global_size
const Real	_tolerance = TOLERANCE * TOLERANCE

Detailed Description

Definition at line 17 of file lumped_mass_matrix_test.C.

Member Function Documentation

CPPUNIT_TEST()

LumpedMassMatrixTest::CPPUNIT_TEST

(

testNumerics

)

CPPUNIT_TEST_SUITE_END()

LumpedMassMatrixTest::CPPUNIT_TEST_SUITE_END

(

)

LIBMESH_CPPUNIT_TEST_SUITE()

LumpedMassMatrixTest::LIBMESH_CPPUNIT_TEST_SUITE

(

LumpedMassMatrixTest

)

setUp()

void LumpedMassMatrixTest::setUp ( )

inline

Definition at line 27 of file lumped_mass_matrix_test.C.

References TestCommWorld.

  {
    _comm = TestCommWorld;
    _matrix = std::make_unique<LumpedMassMatrix<Number>>(*_comm);

    numeric_index_type root_block_size = 2;
    _local_size = root_block_size + static_cast<numeric_index_type>(_comm->rank());
    _global_size = 0;

    for (processor_id_type p = 0; p < _comm->size(); ++p)
    {
      numeric_index_type block_size = root_block_size + static_cast<numeric_index_type>(p);
      _global_size += block_size;
    }

    _matrix->init(_global_size, UNUSED, _local_size, UNUSED);
  }

tearDown()

void LumpedMassMatrixTest::tearDown ( )

inline

Definition at line 45 of file lumped_mass_matrix_test.C.

{}

testNumerics()

void LumpedMassMatrixTest::testNumerics ( )

inline

Definition at line 47 of file lumped_mass_matrix_test.C.

References libMesh::make_range(), and libMesh::Real.

  {
    LOG_UNIT_TEST;

    numeric_index_type beginning_index = _matrix->row_start();
    numeric_index_type end_index = _matrix->row_stop();

    CPPUNIT_ASSERT_EQUAL(_local_size, numeric_index_type(_matrix->row_stop() - _matrix->row_start()));

    _matrix->zero();

    std::vector<Real> gold_values(_local_size, 0);

    // Test add
    for (const auto i : make_range(beginning_index, end_index))
      for (const auto j : make_range(beginning_index, end_index))
      {
        const Real sgn = j % 2 ? 1 : -1;
        _matrix->add(i, j, sgn * j);
        gold_values[i - beginning_index] += j;
      }

    _matrix->close();

    for (const auto i : make_range(beginning_index, end_index))
      for (const auto j : make_range(beginning_index, end_index))
      {
        if (i != j)
          LIBMESH_ASSERT_NUMBERS_EQUAL(0, (*_matrix)(i, j), _tolerance);
        else
          LIBMESH_ASSERT_NUMBERS_EQUAL(gold_values[i - beginning_index], (*_matrix)(i, i), _tolerance);
      }

    // Test set
    for (const auto i : make_range(beginning_index, end_index))
    {
      const Real sgn = i % 2 ? 1 : -1;
      _matrix->set(i, i, sgn * gold_values[i - beginning_index]);
    }

    _matrix->close();

    for (const auto i : make_range(beginning_index, end_index))
      LIBMESH_ASSERT_NUMBERS_EQUAL(gold_values[i - beginning_index], (*_matrix)(i, i), _tolerance);

    // Test clone
    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());

    for (const auto i : make_range(beginning_index, end_index))
      for (const auto j : make_range(beginning_index, end_index))
        LIBMESH_ASSERT_NUMBERS_EQUAL((*_matrix)(i, j), (*copy)(i, j), _tolerance);

    // Test zero clone
    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());

    for (const auto i : make_range(beginning_index, end_index))
      for (const auto j : make_range(beginning_index, end_index))
        LIBMESH_ASSERT_NUMBERS_EQUAL(0, (*zero_copy)(i, j), _tolerance);
  }

Member Data Documentation

_comm

Parallel::Communicator* LumpedMassMatrixTest::_comm

private

Definition at line 121 of file lumped_mass_matrix_test.C.

_global_size

numeric_index_type LumpedMassMatrixTest::_global_size

private

Definition at line 123 of file lumped_mass_matrix_test.C.

_local_size

numeric_index_type LumpedMassMatrixTest::_local_size

private

Definition at line 123 of file lumped_mass_matrix_test.C.

_matrix

std::unique_ptr<LumpedMassMatrix<Number> > LumpedMassMatrixTest::_matrix

private

Definition at line 122 of file lumped_mass_matrix_test.C.

_tolerance

const Real LumpedMassMatrixTest::_tolerance = TOLERANCE * TOLERANCE

private

Definition at line 124 of file lumped_mass_matrix_test.C.

---

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

• tests/numerics/lumped_mass_matrix_test.C