Dear all
I have an FE_Q space and a ghosted p::d::Vector. When I call
VectorTools::interpolate on this, I get this error
An error occurred in line <679> of file
</Users/praveen/Applications/deal.II/git/include/deal.II/lac/la_parallel_vector.templates.h>
in function
void
dealii::LinearAlgebra::distributed::Vector<double>::compress_finish(::dealii::VectorOperation::values)
The violated condition was:
*read_position == Number() || std::abs(local_element(j) -
*read_position) <= std::abs(local_element(j)) * 1000. *
std::numeric_limits<real_type>::epsilon()
The name and call sequence of the exception was:
ExcNonMatchingElements(*read_position, local_element(j),
part.this_mpi_process())
Additional Information:
Called compress(VectorOperation::insert), but the element received from a
remote processor, value 3.343515457745695e-15, does not match with the
value -3.343515457745695e-15 on the owner processor 0
I dont understand why this error is triggered.
Also, I cannot pass a non-ghosted p::d::Vector to VectorTools::interpolate,
I get this error
An error occurred in line <1285> of file
</Users/praveen/Applications/deal.II/git/include/deal.II/lac/la_parallel_vector.h>
in function
Number
&dealii::LinearAlgebra::distributed::Vector<double>::operator()(const
size_type)
The violated condition was:
partitioner->in_local_range (global_index) ||
partitioner->ghost_indices().is_element(global_index)
The name and call sequence of the exception was:
ExcAccessToNonLocalElement(global_index,
partitioner->local_range().first, partitioner->local_range().second,
partitioner->ghost_indices().n_elements())
Additional Information:
You tried to access element 0 of a distributed vector, but this element is
not stored on the current processor. Note: The range of locally owned
elements is 10084 to 19844, and there are 0 ghost elements that this vector
can access.
Thanks
praveen
--
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].
For more options, visit https://groups.google.com/d/optout.
