Hi,
I've tried running step 17 after installing PETSc and it runs fine. Now I am
trying to code up poisson type equation using Discontinuous Galerkin
Interior Penalty method using step 17 . It compiles without any errors but
when I try to run the program I get the following error :
terminate called after throwing an instance of
'dealii::PETScWrappers::internal::VectorReference::ExcAccessToNonlocalElement'
what(): --------------------------------------------------------
An error occurred in line <140> of file
</gpfs/runtime/opt/dealii/7.0.0/source/lac/petsc_vector_base.cc> in function
dealii::PETScWrappers::internal::VectorReference::operator PetscScalar()
const
The violated condition was:
(index >= static_cast<unsigned int>(begin)) && (index <
static_cast<unsigned int>(end))
The name and call sequence of the exception was:
ExcAccessToNonlocalElement (index, begin, end-1)
Additional Information:
You tried to access element 27 of a distributed vector, but only elements 9
through 17 are stored locally and can be accessed.
After some debugging , I got to know that the error is because of the
following statement:
fe_v_neighbor.get_function_values(old_f,
f_values_neighbor);
I comment this out and everything works fine. Here old_f is a variable I
calculate from advection equation and substitute in the poisson equation.
I get what the error says but I do need to access the neighboring element in
IP methods. How can I do this with PETSc? The program runs fine in serial
computing even when I use
mpiexec -np 1
--
Thanks,
gk
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
