Sir, thank you for your reply I did changed the code based on your suggestion and the problem still existed.
However, the problem was fixed. In the code, I had another procedure that will project the solution and made sure that it does not exceed a particular value. In that function originally I did not update the ghosted values. I fixed that and the problem was fixed. Thank you On Wednesday, November 17, 2021 at 12:51:18 AM UTC+9 [email protected] wrote: > Hello, > > Which line triggers the error? Is it the zero_out_ghost() line? This line > should not be there, the ghost elements are set to zero when you do =0 > > Best, > > Bruno > > On Tuesday, November 16, 2021 at 3:13:23 AM UTC-5 SALMAN wrote: > >> Dear all, >> >> I am developing a time dependent code which utilizes dealii adaptive >> meshing capabilities. >> >> I am having a problem with the solution transfer during the program. >> >> At a particular time step, the following error occurs: >> -------------------------------------------------------- >> An error occurred in line <624> of file >> </home/slinux/dealii-9.3.1/include/deal.II/base/partitioner.templates.h> in >> function >> void >> dealii::Utilities::MPI::Partitioner::import_from_ghosted_array_finish(dealii::VectorOperation::values, >> >> const dealii::ArrayView<const ElementType, MemorySpaceType>&, const >> dealii::ArrayView<ElementType, MemorySpace>&, const >> dealii::ArrayView<ElementType, MemorySpace>&, std::vector<ompi_request_t*, >> std::allocator<ompi_request_t*> >&) const [with Number = double; >> MemorySpaceType = dealii::MemorySpace::Host] >> The violated condition was: >> *read_position == Number() || >> internal::get_abs(locally_owned_array[j] - *read_position) <= >> internal::get_abs(locally_owned_array[j] + *read_position) * 100000. * >> std::numeric_limits<typename numbers::NumberTraits< >> Number>::real_type>::epsilon() >> Additional information: >> Called compress(VectorOperation::insert), but the element received >> from a remote processor, value 1.387778780781446e-17, does not match >> with the value 0 on the owner processor 12 >> >> >> The problem occurs only in debug mode when running in parallel. I have >> tried to search the library documentation to find a solution, however, this >> was not possible. >> >> I have attached the code snippet for solution transfer that is used in >> the program. >> >> I appreciate any tips or advice on how to fix the issue. >> >> Thank you >> > -- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/f7a9a279-b1f1-4314-b7af-e14008004383n%40googlegroups.com.
