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.

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.

Remove deprecated function FEValuesBase::transform().

Avoid using deprecated functions in the testsuite.

Merge pull request #4341 from drwells/deprecate-non-sp-petsc-matrix-reinit

Deprecate PETSc matrix constructors without sparsity patterns.

Merge pull request #4342 from drwells/remove-backup-file

Remove an editor backup file.

Remove an editor backup file.

Deprecate PETSc matrix constructors without sparsity patterns.

These are known to have efficiency problems since the underlying matrix does not
allocate the correct amount of memory.

Users can still allocate memory the old way if they wish to by calling
MatCreateAIJ themselves since the base class provides a conversion operator back
to a PETSc object (Mat).

Merge pull request #4279 from drwells/remove-deprecated-petsc

Remove deprecated PETSc vector classes.

Add a conceptual note on usage of PETScWrappers::VectorBase.

Fix (serial) PETSc vector documentation.

This commit fixes up some old text which implied that we had to use either
serial PETSc vectors or PETSc vectors for computations: since we support more
than just PETSc these days we can generalize a bit.

Consistently capitalize PETSc.

Refer to all parallel vectors instead of just the PETSc class.

Remove deprecated PETSc vector classes.

Get rid of PETScWrappers::Vector in the PETSc matrix classes.

We can simply use PETScWrappers::VectorBase instead. In addition, get rid of
some methods in the derived class that are identical to those in MatrixBase.

Make the PETSc to deal.II vector copy generic.

We may as well allow copying instances of the base class: the operations are the
same.

Overhaul the PETSc MatrixFree vmult.

Make PETScWrappers::apply_boundary_values more generic.

Merge pull request #4340 from masterleinad/fix_petsc_complex_tests

Fix petsc_complex tests

Fix petsc_complex tests

Merge pull request #4338 from davydden/fix_compute_active_cell_layer_within_distance

remove unused template argument

remove unused template argument

Merge pull request #4337 from bangerth/fix-test

Fix a test that requires MPI.

Fix a test that requires MPI.

Merge pull request #4336 from bangerth/fix-petsc-test

Fix a PETSc test that doesn't compile.

Fix a PETSc test that doesn't compile.

Merge pull request #4297 from drwells/convert-petsc-tests-1

Convert petsc tests 1

Rewrite tests that used deprecated PETSc vectors.

PETScWrappers::Vector and PETScWrappers::BlockVector were deprecated in
favor of PETScWrappers::MPI::Vector and PETScWrappers::MPI::BlockVector.

This involves removing all of the PETScWrappers::BlockMatrix tests,
since that class takes one of the deprecated vectors as an argument.

Merge pull request #4093 from vishalkenchan/new/GridTools/compute_active_cells_within_skin

New: GridTools::compute_active_cells_within_skin

Merge pull request #4335 from masterleinad/remove_dealii_cxx11

Remove DEAL_II_WITH_CXX11

Remove DEAL_II_WITH_CXX11

Merge pull request #4334 from bangerth/remove-deprecated-static-var

Remove the deprecated member variables 'supports_distributed_data' of vector classes.

Add a changelog entry.

Remove the deprecated member variables 'supports_distributed_data' of vector classes.

Merge pull request #4332 from masterleinad/remove_petsc_use_64bit_indices

Remove usages of PETSC_USE_64BIT_INDICES

Merge pull request #4328 from bangerth/fix-serialization

Make serialization work with older GCC versions as well.

Remove usages of PETSC_USE_64BIT_INDICES

Merge pull request #4330 from drwells/petsc-vec-type-invariant

Make PETScWrappers::MPI::Vector always distributed.

GridTools::compute_bounding_box(), compute_active_cells_within_skin(), compute_ghost_cells_within_skin()

Un-inline PETScWrappers::MPI::Vector::operator=(Vector &).

This function is long enough that the compiler probably won't inline it anyway.

Make the Vec type of PETScWrappers::MPI::Vector invariant.

PETSc supports both sequential ('seq') and distributed ('mpi') vectors: since
this vector is supposed to be distributed, we should always make it an MPI (or
ghosted) vector to simplify this class.

Merge pull request #4325 from masterleinad/fix_64bit_complex_set

Fix PETScWrappers::MatrixBase::set for 64bit and complex scalars

Make serialization work with older GCC versions as well.

Remove special treatment of the 64bit case

Merge pull request #4326 from davydden/doc/partitioner

adjust documentation of Partitioner to, hopefully, make it more clear

Merge pull request #4320 from vishalkenchan/new/smallest_enclosing_sphere

New TriaAccessor function smallest_enclosing_sphere() added.