Guard FIND_PACKAGE(CUDA) with DEAL_II_WITH_CUDA.

Merge pull request #4402 from Rombur/cuda_tests

Check for cuda only if configured with cuda support.

Check for cuda only if configured with cuda support.

Merge pull request #4398 from davydden/remove_wrong_todo

remove wrong TODO in template-arguments.in

REAL_SERIAL_VECTORS and SERIAL_VECTORS are not the same.

The former contains only real-valued vectors, i.e.
DEAL_II_EXPAND_PETSC_MPI_VECTOR_REAL will be empty if
PETSc is build with complex linear algebra.

Merge pull request #4394 from drwells/fix-rebase-print

Fix the output of the deprecation script.

Fix the output of the deprecation script.

Some tags (e.g., 'rebase_here') may not be version numbers so we should not drop
the first letter.

Merge pull request #4293 from Rombur/matrixfree_gpu

Matrix-free gpu

Add changelog entry.

Add test for CUDA support of matrix-free.

Add CUDA support of matrix-free.

Rename cuda_vector tests.

Augment documentation of CUDA atomic.

Add CUDA support for tensor.

Merge pull request #4397 from masterleinad/fix_pgi

Fixes for compiling with PGI

Fixes for compiling with PGI

Merge pull request #4395 from masterleinad/remove_reorder_cuthill_mckee

Remove deprecated SparsityTools::reorder_Cuthill_McKee

Merge pull request #4396 from masterleinad/remove_unused_helper

Remove unused helper function

Remove unused helper function

Remove deprecated SparsityTools::reorder_Cuthill_McKee

Merge pull request #4391 from masterleinad/remove_ez_rows_columns

Remove deprecated BlockSparseMatrixEZ<Number>::n_rows/n_cols

Merge pull request #4390 from masterleinad/fix_chunk_sparsity_pattern_12

Fix chunk_sparsity_pattern_12

Remove deprecated BlockSparseMatrixEZ<Number>::n_rows/n_cols

Fix chunk_sparsity_pattern_12

Merge pull request #4389 from masterleinad/remove_system_j_s_m

Remove Utilities::System::job_supports_mpi()

Remove Utilities::System::job_supports_mpi()

Merge pull request #4388 from bangerth/avoid-warning-2

Avoid a warning about an unused constructor argument.

Modify CMake for CUDA matrix-free.

Avoid a warning about an unused constructor argument.

Annoyingly, my compiler warns about =default constructors whenever we name constructor
arguments. It says they are unused -- which is patently wrong: the argument is used
just fine, we just never explicitly reference the *name* of the argument. Either way,
work around this.

Merge pull request #4387 from drwells/add-missing-indexset-instantiations

Add missing FEValues IndexSet instantiations.

Add missing FEValues IndexSet instantiations.

Merge pull request #4386 from masterleinad/fix_point_constructor

Add user-defined constructor for Point<dim> back in

Merge pull request #4385 from drwells/fix-clang-compilation

Fix clang compilation

Add user-defined constructor for Point<dim> back in

Remove an unused private variable.

Use a Tensor for a direction vector.

The old version works since Point<3> inherits from Tensor<1, 3> and this is only
implemented for the spacedim == 3 case.

Merge pull request #4381 from masterleinad/use_valid_arguments

Use valid arguments in ChunkSparsityPattern constructor

Merge pull request #4375 from masterleinad/use_override_consistently

Use override consistently

Use valid arguments in ChunkSparsityPattern constructor

Merge pull request #4379 from davydden/change_assert

change assert in p::d::Tria::load()

Use override consistently

Merge pull request #4377 from drwells/default-constructor

Use '= default' for some constructors.

change assert in p::d::Tria::load()

Merge pull request #4367 from drwells/add-oned-manifold

Add get_manifold() to 1D vertex iterators.

Merge pull request #4374 from davydden/mf_gmg_tests_fix

MF GMG tests (fix one and add another failing)

add another MF GMG test which was failing prior to cb5c93d

enable threads in MF GMG test, add additional output

Merge pull request #4376 from kronbichler/fix_matrix_free_with_threads

Fix bug in matrix free thread graph creation.

Use '= default' for some constructors.

Add get_manifold() to 1D vertex iterators.

Merge pull request #4371 from davydden/trial_load_assert

add an assert in p::d::Tria::load() that the file exists

Fix bug in matrix free thread graph creation.

add an assert in p::d::Tria::load() that the file exists

Merge pull request #4368 from kostyfisik/disc-nt

add note on discontinuous coefficient

add note on discontinuous coefficient

Merge pull request #4365 from danshapero/constraints-move

Added default move operations to ConstraintMatrix

Merge pull request #4360 from masterleinad/fix_tests_after_trilinos_functions

Fix tests after removal of deprecated functions in TrilinosWrappers vectors

Merge pull request #4364 from drwells/remove-custom-sort

Get rid of a custom sort implementation.

Added default move operations to ConstraintMatrix

Get rid of a custom sort implementation.

This sorting function is about 10% faster for very large grids (O(10^7)
cell error indicators) but otherwise has the same performance as
std::sort. This patch replaces this with std::sort and a custom
comparator. To see the performance change, try running:

#include <algorithm>
#include <chrono>
#include <iostream>
#include <random>
#include <vector>

template <typename T>
void qsort_index (const std::vector<T>     &a,
                  std::vector<std::size_t> &ind,
                  std::size_t               l,
                  std::size_t               r)
{
  std::size_t i,j;
  T v;

  if (r<=l)
    return;

  v = a[ind[r]];
  i = l-1;
  j = r;
  do
    {
      do
        {
          ++i;
        }
      while ((a[ind[i]]>v) && (i<r));
      do
        {
          --j;
        }
      while ((a[ind[j]]<v) && (j>0));

      if (i<j)
        std::swap (ind[i], ind[j]);
      else
        std::swap (ind[i], ind[r]);
    }
  while (i<j);
  if (i > 0)
    {
      qsort_index(a,ind,l,i-1);
    }
  qsort_index(a,ind,i+1,r);
}

int main()
{
  // switch to length = 10000000 to see timing information
  constexpr std::size_t length = 20;
  std::vector<std::size_t> indices(length);
  std::iota(indices.begin(), indices.end(), std::size_t(0));

  std::uniform_int_distribution<std::size_t> unif(0, length);
  std::default_random_engine re;

  std::vector<std::size_t> values;

  for (std::size_t index = 0; index < indices.size(); ++index)
    {
      values.push_back(unif(re));
    }

  // for (std::size_t index = 0; index < length; ++index)
  //   {
  //     std::cout << index << ": "
  //               << indices[index] << ": "
  //               << values[indices[index]] << '\n';
  //   }

  {
    auto v = values;
    auto i = indices;

    const auto t0 = std::chrono::high_resolution_clock::now();
    qsort_index(v, i, 0, length - 1);
    const auto t1 = std::chrono::high_resolution_clock::now();
    std::cout << "milliseconds:"
              << std::chrono::duration_cast<std::chrono::milliseconds>(t1 - t0).count()
              << std::endl;

    std::cout << "\nafter qsort_index:\n";
    // for (std::size_t index = 0; index < length; ++index)
    //   {
    //     std::cout << index << ": "
    //               << i[index] << ": "
    //               << v[i[index]] << '\n';
    //   }
  }

  {
    auto v = values;
    auto i = indices;

    const auto t0 = std::chrono::high_resolution_clock::now();
    std::sort(i.begin(), i.end(), [&](const std::size_t left,
                                      const std::size_t right)
              {
                return v[left] >= v[right];
              });
    const auto t1 = std::chrono::high_resolution_clock::now();
    std::cout << "milliseconds:"
              << std::chrono::duration_cast<std::chrono::milliseconds>(t1 - t0).count()
              << std::endl;

    std::cout << "\nafter std::sort:\n";
    // for (std::size_t index = 0; index < length; ++index)
    //   {
    //     std::cout << index << ": "
    //               << i[index] << ": "
    //               << v[i[index]] << '\n';
    //   }
  }
}

Merge pull request #4362 from Rombur/all_zero

Add all_zero to VectorSpaceVector

add tests.

Add all_zero() to VectorSpaceVector.

Fix failing tests after removal of deprecated TrilinosWrappers vector functions

Add get_mpi_communicator(), locally_owned_range_indices() and locally_owned_domain_indices ()

Restore an accidentally deleted definition

Merge pull request #4354 from drwells/fix-petsc-copy

Fix petsc copy

Note the PETSc vector removal in the changelog.

Fix some bugs in PETSc-deal.II vector copy code.

This new, correct implementation is based on the original, working
implementation present before commit 3aa64f983d4 in the (now removed)
PETSc serial vector class.

Merge pull request #4345 from masterleinad/remove_deprecated_trilinos_vectors

Remove deprecated functions in TrilinosWrappers vectors

Merge pull request #4357 from bangerth/minor-cleanups

Some text cleanups.

Some text cleanups.

No functionality changes.

Merge pull request #4344 from drwells/add-deprecation-script

Add a script that finds deprecated items.

Add a script that finds deprecated items.

Merge pull request #4355 from kronbichler/fix_test_cm

Fix constraint matrix test in 64 bit mode

Fix constraint matrix test in 64 bit mode

Merge pull request #4352 from bangerth/remove-deprecated-functions-5

Remove deprecated function Manifold::get_new_point(Quadrature).

Merge pull request #4353 from bangerth/remove-deprecated-functions-6

Remove the deprecated constructor of MappingQEulerian.

Merge pull request #4343 from kronbichler/extend_constraint_matrix_tltg

Generalize CM::distribute_local_to_global for rectangular case.

Update changelog entry.

Remove the deprecated constructor of MappingQEulerian.

No longer use the deprecated constructor of MappingQEulerian.

Add changelog entry.

Remove deprecated function Manifold::get_new_point(Quadrature).

Changelog entry

Remove deprecated functions in TrilinosWrappers vectors

Clean testsuite from deprecated Trilinos vector functions

Merge pull request #4351 from bangerth/adjust-tests

Avoid using deprecated functions in the testsuite.

Merge pull request #4346 from bangerth/remove-deprecated-functions-1

Remove deprecated function FEValuesBase::transform().

Merge pull request #4347 from bangerth/remove-deprecated-functions-2

Remove the deprecated constructor of MappingQEulerian.

Merge pull request #4348 from bangerth/remove-deprecated-functions-3

Remove deprecated version of GridTools::get_face_connectivity_of_cells().

Merge pull request #4349 from bangerth/remove-deprecated-functions-4

Remove a deprecated typedef.

Generalize CM::distribute_local_to_global for rectangular case.

Add changelog entry.

Remove a deprecated typedef.

Add changelog entry.

Remove deprecated version of GridTools::get_face_connectivity_of_cells().

Add changelog entry.

Remove the deprecated constructor of MappingQEulerian.

Add changelog entry.