Modify a test to check the diagonal values (not lumped).

Make MassOperator's diagonals true diagonals, not lumped ones.

The previous implementation used row sums (i.e., mass lumping) instead of
computing the diagonal. This should be changed for two reasons:
1. It's not really the diagonal (and its inconsistent with LaplaceOperator).
2. Lumped mass matrices are only positive for either low-order elements (where
   the basis functions are all positive) or very well-behaved elements (like
   FE_Q). In particular, the values on the diagonal are either zero or negative
   for TET10 if we use lumping, which isn't going to work.

Merge pull request #12618 from peterrum/TriangulationDescription_examples

Extend documentation of TriangulationDescription

Extend documention of TriangulationDescription

Merge pull request #12610 from drwells/add-more-wv-rules

Add more Witherden-Vincent rules

Merge pull request #12619 from drwells/fix-parenthesis

Fix a parenthesis.

Fix a parenthesis.

I see a bunch of warnings of the form

/Users/runner/work/dealii/dealii/include/deal.II/matrix_free/tools.h:497:35:
error: taking the absolute value of unsigned type 'bool' has no effect [-Werror,-Wabsolute-value]
854
                                  std::abs(values_dofs[j][0] &&

on the CI tester.

Add a lot more Witherden-Vincent rules.

This includes most of the even-degree rules.

For posterity: here's the script I used to print quadpy rules rather than
copying numbers individually:

import quadpy as qp

def format_16(x):
    return '{:.16e}'.format(float(x))

def print_scheme(scheme):
    sym_data = scheme.symmetry_data
    print("// WV-{}, {}D".format(scheme.degree, scheme.dim))
    for key in sym_data.keys():
        if key == 'centroid' or key == 's4':
            print("b_point_permutations.push_back({centroid});")
            print("b_weights.push_back({});".format(format_16(sym_data[key][0][0])))
        elif key == 'd3_aa' or key == 's31':
            weights = [format_16(w) for w in sym_data[key][0]]
            points = [format_16(p) for p in sym_data[key][1]]
            assert len(weights) == len(points)
            for i in range(len(weights)):
                print("process_point_1({}, {});".format(points[i], weights[i]))
        elif key == 'd3_ab' or key == 's211':
            weights = [format_16(w) for w in sym_data[key][0]]
            points1 = [format_16(p) for p in sym_data[key][1]]
            points2 = [format_16(p) for p in sym_data[key][2]]
            assert len(weights) == len(points1)
            assert len(points2) == len(points1)
            for i in range(len(weights)):
                print(("process_point_3({},"
                       "\n                {},"
                       "\n                {});")
                      .format(points1[i], points2[i], weights[i]))
        elif key == 's22':
            assert scheme.dim == 3
            weights = [format_16(w) for w in sym_data[key][0]]
            points = [format_16(p) for p in sym_data[key][1]]
            assert len(weights) == len(points)
            for i in range(len(weights)):
                print(("process_point_2({},"
                       "\n                {});")
                      .format(points[i], weights[i]))
        else:
            assert False

Merge pull request #12613 from drwells/fix-structdim-comments

Fix two comments in TriaAccessor.

Merge pull request #12559 from tamiko/modernize_step_69

Step-69: Use SolutionTransfer instead of low-level boost archive

Step-69: Apply review suggestion

Step-69: Use SolutionTransfer instead of low-level boost archive

Merge pull request #12560 from peterrum/mf_hn

Optimize hanging-node constraints in MatrixFree

Merge pull request #12611 from drwells/boundary-to-manifold

Convert some (mislabeled) boundary ids to manifold ids.

Fix two comments in TriaAccessor.

We have a 'cell' when the structural dimension matches the topological
dimension.

Merge pull request #12607 from bangerth/type

Fix a few type mistakes in AlignedVector.

Convert some (mislabeled) boundary ids to manifold ids.

Use a switch statement.

Less indentation makes this easier to read.

Remove an extra namespace qualification.

This isn't necessary now that we moved Simplex::QGauss to QGaussSimplex.

Fix a few type mistakes in AlignedVector.

Merge pull request #12556 from bangerth/units

Merge pull request #12605 from tjhei/doc_typo_rtree

doc: fix typo in rtree.h

Generalize HangingNodes

Try a different syntax.

doc: fix typo

Merge pull request #12602 from peterrum/rpe_unique_fix

Fix tolerances in GridTools::internal::distributed_compute_point_locations()

Fix tolerances in GridTools::internal::distributed_compute_point_locations()

Merge pull request #11108 from masterleinad/clang-tidy-11

Apply clang-format

Delete indent step in CI

add Darwin Intel binaries

Attempt to suppress clang-tidy warning

add Darwin Intel binaries

update clang-format download script to version 11

Update infrastructure some more

update clang-format infrastructure to version 11

Require the correct version in contrib/utilities/indent_common.sh

Update to clang-tidy-11

Merge pull request #12600 from peterrum/convert_map_pattern

Fix Convert for map

Merge pull request #12578 from agrayver/python_mapping_rename

Rename MappingQGeneric to MappingQ in python interface

Merge pull request #12599 from peterrum/find_active_cell_around_point_end

Let find_active_cell_around_point() return mesh.end() if point is not found in BBs

Rename MappingQGeneric to MappingQ in python

Merge pull request #12598 from peterrum/project_docu

Fix documentation of VT::project()

Fix Convert for map

Let find_active_cell_around_point() return mesh.end() if point is not found in BBs

Fix documentation of VT::project()

Merge pull request #12595 from bangerth/doc

Update documentation in a couple of places.

Update documentation in a couple of places.

Merge pull request #12588 from zjiaqi2018/table_clear

Add clear() to the Table class

Merge pull request #12592 from peterrum/compute_diagonal_util_tolerance_float

Reduce tolerance in test

Merge pull request #12591 from marcfehling/intel18

Intel 18 does not like assertions in nested lambdas.

Reduce tolerance in test

Intel18 does not like assertions in nested lambdas.

Merge pull request #12554 from marcfehling/ci-oneapi

Add Intel oneAPI to github-actions.

Merge pull request #12590 from Rombur/fix_warnings

Fix unused-but-set-variable warnings

Merge pull request #12589 from bangerth/include

Add a missing include file.

minor fix and changelog

Add a missing include file.

Fix unused-but-set-variable warnings

add table.clear()

Merge pull request #12585 from peterrum/compute_diagonal_ref

Add new compute_diagonal() test

Add new compute_diagonal() test

Merge pull request #12577 from peterrum/hn_generalize

Generalize HangingNodes

Generalize HangingNodes

Merge pull request #12583 from drwells/mapping-fe-field-face-simplices

Mapping fe field face simplices

Merge pull request #12582 from peterrum/ConstraintTypes_enum

Make internal::MatrixFreeFunctions::ConstraintTypes enum

Add a test for both features.

Fix cell->face(0)->measure() for tets.

Fortunately its easy to find the area of a triangle.

Copy quadrature weights in MappingFEField.

Use the correct reference cell in MappingFEField.

Add a new class member MappingFEField::reference_cell.

This is easier than constantly looking it up.

Make internal::MatrixFreeFunctions::ConstraintTypes enum

Merge pull request #12571 from peterrum/hanging_nodes_internal

Fill CUDAWrappers::internal::constraint_weights inside of ReinitHelper

Rename cuda_hanging_nodes_internal.h to hanging_nodes_internal.h

Merge pull request #12581 from peterrum/read_dof_indices

DoFInfo::read_dof_indices(): add additional index vector

DoFInfo::read_dof_indices(): add additional index vector

Merge pull request #12579 from kronbichler/fix_64bit

DoFHandlerPolicy: Fix 64 bit code path

DoFHandlerPolicy: Fix 64 bit code path

Merge pull request #12576 from kronbichler/simplify_dof_accessor

Implement DoFHandler::distribute_dofs with a single pass through entities

Merge pull request #12562 from peterrum/MGCoarseGridIterativeSolver_float

Enable different VectorTypes in MGCoarseGridIterativeSolver

Use faster process_dof_indices in one more place

Merge pull request #12563 from peterrum/si_subface_interpolation_matrix

Introduce ShapeInfo::subface_interpolation_matrix

DoFHandler::distribute_dofs: only go through entities once

Simplify implementation of set/get_dof_indices

Merge pull request #12574 from kronbichler/improve_dof_handler_distr

Avoid unnecessary updates of cell dof index cache in DoFHandlerPolicy

Merge pull request #12572 from bangerth/64b

Add a couple of exceptions that ensure we don't accidentally run out of indices.

Merge pull request #12573 from drwells/q-iterated-simplex

Implement QIteratedSimplex.

Add documentation about cache

Implement QIteratedSimplex.

A limited version of QIterated for simplex meshes.

Introduce ShapeInfo::subface_interpolation_matrix

Merge pull request #12454 from kronbichler/update_raviart_thomas_nodal

Straight-forward tensor-product version of FE_RaviartThomasNodal

Avoid unnecessary updates of cell dof index cache in DoFHandlerPolicy

Merge pull request #12564 from masterleinad/suggest_non_mpi_compiler_wrappers

Avoid suggesting using MPI compiler wrappers as CMake compilers

Add a couple of exceptions that ensure we don't accidentally run out of indices.

(void) a variable used only in an assertion.

Follow recommendations of clang-tidy.

Do some rudimentary error checking.

Add tests.

Add a changelog entry.