On 04/04/2018 06:21 PM, Phani Motamarri wrote:
I have just written a plan vanilla interface to call the Kelly Error Estimate
function passing a parallel distributed vector as input vector type. But
strangely enough when I run on MPI on multiple procs in debug mode, I get the
following error:
--------------------------------------------------------
An error occurred in line <1351> of file
</home/vikramg/DFT-FE-softwares/softwareCentos/dealiiDev/dealii/include/deal.II/lac/la_parallel_vector.h>
in function
Number
dealii::LinearAlgebra::distributed::Vector<Number>::operator()(unsigned long
long) const [with Number = double]
The violated condition was:
partitioner->in_local_range (global_index) ||
partitioner->ghost_indices().is_element(global_index)
Additional information:
You tried to access element 2 of a distributed vector, but this element is not
stored on the current processor. Note: The range of locally owned elements is
1377 to 2601, and there are 153 ghost elements that this vector can access
------------------------------------------------------------------------------------------
I have attached the minimal example and also the stack trace. This error is
probably being triggered somewhere in the estimate function. I am wondering if
I am missing anything here. It runs without any error message on 1 processor.
The vector you pass to KellyErrorEstimator::estimate() appears to only have
locally owned elements, but seems to not store ghost elements. But the
estimate() function obviously needs to access all "locally active" DoF values
on the locally owned cells, and so the vector you pass in needs to have ghost
elements.
You initialize this vector using
MatrixFree<3,double> matrix_free_data;
...
dealii::parallel::distributed::Vector<double> vec;
matrix_free_data.initialize_dof_vector(vec);
I don't know whether this function is supposed to initialize ghost entries in
the vector -- maybe one of the MatrixFree folks could help out with the answer?
Cheers
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
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.
