Hi Bruno:
Thanks for pointing that out.
I tried to not refine FE_nothing cells by modifying the refine loop:
(The modified test is attached).
for (auto &cell : dgq_dof_handler.active_cell_iterators())
if (cell->is_locally_owned() && cell->active_fe_index () != 0)
{
if (counter > ((dim == 2) ? 4 : 8))
cell->set_coarsen_flag();
else
cell->set_refine_flag();
}
But this still aborts.
Kaushik
On Tue, Dec 8, 2020 at 8:36 AM Bruno Turcksin <bruno.turck...@gmail.com>
wrote:
> Hi,
>
> Are you sure that your test makes sense? You randomly assign FE indices to
> cells then you refine and coarsen cells. But what does it mean to coarsen 4
> cells together when one of them is FE_Nothing? What would you expect to
> happen?
>
> Best,
>
> Bruno
>
> On Monday, December 7, 2020 at 5:54:10 PM UTC-5 k.d...@gmail.com wrote:
>
>> Hi all:
>>
>> I modified the test tests/mpi/solution_transfer_05.cc to add a
>> FE_Nohting element to the FECollection. I also modified the other elements
>> to FE_Q.
>>
>> When I run the test, it's aborting in solution transfer.
>> Is there any limitations in using FE_Nothing with parallel solution
>> transfer?
>> The modified test is attached.
>> Thank you very much.
>> Kaushik
>>
>> ----
>> An error occurred in line <1167> of file
>> </build/deal.ii-vFp8uU/deal.ii-9.2.0/include/deal.II/lac/vector.h> in
>> function
>> Number&
>> dealii::Vector<Number>::operator()(dealii::Vector<Number>::size_type) [with
>> Number = double; dealii::Vector<Number>::size_type = unsigned int]
>> The violated condition was:
>> static_cast<typename ::dealii::internal::argument_type<void( typename
>> std::common_type<decltype(i), decltype(size())>::type)>:: type>(i) <
>> static_cast<typename ::dealii::internal::argument_type<void( typename
>> std::common_type<decltype(i), decltype(size())>::type)>:: type>(size())
>> Additional information:
>> Index 0 is not in the half-open range [0,0). In the current case,
>> this half-open range is in fact empty, suggesting that you are accessing an
>> element of an empty collection such as a vector that has not been set to
>> the correct size.
>>
>>
>> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to a topic in the
> Google Groups "deal.II User Group" group.
> To unsubscribe from this topic, visit
> https://groups.google.com/d/topic/dealii/ssEva6C2PU8/unsubscribe.
> To unsubscribe from this group and all its topics, send an email to
> dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/609b0457-ba90-4aa2-a7d1-5b798d5349ebn%40googlegroups.com
> <https://groups.google.com/d/msgid/dealii/609b0457-ba90-4aa2-a7d1-5b798d5349ebn%40googlegroups.com?utm_medium=email&utm_source=footer>
> .
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/CAC-fs6t%3DSSMODByyuRKBMBa20mCnuBb4jQmn04ZoLN6yC73faQ%40mail.gmail.com.
// ---------------------------------------------------------------------
//
// Copyright (C) 2019 - 2020 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------
// This test crashed at some point: We have set and sent active_fe_indices based
// on the refinement flags on the p::d::Triangulation object. However, p4est has
// the last word on deciding which cells will be refined -- and p4est makes use
// of it in the specific scenario provided as a test. A fix has been introduced
// along with this test.
#include <deal.II/base/function.h>
#include <deal.II/distributed/solution_transfer.h>
#include <deal.II/distributed/tria.h>
#include <deal.II/dofs/dof_accessor.h>
#include <deal.II/dofs/dof_tools.h>
#include <deal.II/fe/fe_q.h>
#include <deal.II/fe/fe_nothing.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/hp/dof_handler.h>
#include <deal.II/hp/fe_collection.h>
#include <deal.II/lac/la_parallel_vector.h>
#include <deal.II/numerics/vector_tools.h>
#include "tests.h"
template <int dim>
void
test()
{
// setup
parallel::distributed::Triangulation<dim> tria(MPI_COMM_WORLD);
GridGenerator::hyper_cube(tria);
tria.refine_global(2);
const unsigned int max_degree = 6 - dim;
hp::FECollection<dim> fe_dgq;
fe_dgq.push_back(FE_Nothing<dim>());
for (unsigned int deg = 1; deg <= max_degree - 1; ++deg)
fe_dgq.push_back(FE_Q<dim>(deg));
hp::DoFHandler<dim> dgq_dof_handler(tria);
// randomly assign fes
for (const auto &cell : dgq_dof_handler.active_cell_iterators())
if (cell->is_locally_owned())
cell->set_active_fe_index(Testing::rand() % max_degree);
dgq_dof_handler.distribute_dofs(fe_dgq);
// prepare index sets
IndexSet dgq_locally_owned_dofs = dgq_dof_handler.locally_owned_dofs();
IndexSet dgq_locally_relevant_dofs;
DoFTools::extract_locally_relevant_dofs(dgq_dof_handler,
dgq_locally_relevant_dofs);
IndexSet dgq_ghost_dofs = dgq_locally_relevant_dofs;
dgq_ghost_dofs.subtract_set(dgq_locally_owned_dofs);
// prepare dof_values
LinearAlgebra::distributed::Vector<double> dgq_solution;
dgq_solution.reinit(dgq_locally_owned_dofs, dgq_ghost_dofs, MPI_COMM_WORLD);
VectorTools::interpolate(dgq_dof_handler,
Functions::ZeroFunction<dim>(),
dgq_solution);
dgq_solution.update_ghost_values();
parallel::distributed::SolutionTransfer<
dim,
LinearAlgebra::distributed::Vector<double>,
hp::DoFHandler<dim>>
dgq_soltrans(dgq_dof_handler);
dgq_soltrans.prepare_for_coarsening_and_refinement(dgq_solution);
// refine and transfer
{
unsigned int counter = 0;
//for (auto cell = tria.begin_active(); cell != tria.end(); ++cell,
++counter)
for (auto &cell : dgq_dof_handler.active_cell_iterators())
if (cell->is_locally_owned() && cell->active_fe_index () != 0)
{
if (counter > ((dim == 2) ? 4 : 8))
cell->set_coarsen_flag();
else
cell->set_refine_flag();
}
}
tria.execute_coarsening_and_refinement();
dgq_dof_handler.distribute_dofs(fe_dgq);
// prepare index sets
dgq_locally_owned_dofs = dgq_dof_handler.locally_owned_dofs();
DoFTools::extract_locally_relevant_dofs(dgq_dof_handler,
dgq_locally_relevant_dofs);
dgq_ghost_dofs = dgq_locally_relevant_dofs;
dgq_ghost_dofs.subtract_set(dgq_locally_owned_dofs);
// unpack dof_values
dgq_solution.reinit(dgq_locally_owned_dofs, dgq_ghost_dofs, MPI_COMM_WORLD);
dgq_soltrans.interpolate(dgq_solution);
deallog << "OK" << std::endl;
}
int
main(int argc, char *argv[])
{
Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
MPILogInitAll log;
deallog.push("2d");
test<2>();
deallog.pop();
deallog.push("3d");
test<3>();
deallog.pop();
}
