On Tue, Mar 8, 2016 at 11:35 AM, Harshad Sahasrabudhe <hsaha...@purdue.edu> wrote:
> Hi, > > I'm running into a problem with ghost vectors created with PetscVector in > LibMesh 0.9.5 PETSc 3.5.4. When I call the operator(), the call fails at > the assert local_index < _local_size (petsc_vector.h:1214). The _local_size > variable in PetscVector is 0 (after _get_array is called) even if the > *ghost > *indices vector is non-zero while creating the vector using the > constructor: > > PetscVector (const Parallel::Communicator &comm_in, const > numeric_index_type N, const numeric_index_type n_local, const std::vector< > numeric_index_type > &*ghost*, const ParallelType type=AUTOMATIC) > > I checked that the ghost indices are out of the scope of first_local_dof > and last_local_dof in DofMap. I am localizing the vector using the function > > system->solution->localize(ghost_vector, ghost_indices) > > I don't get any error when I run in debug mode. I assume you mean "opt" mode, otherwise that doesn't make sense because the error you are referring to above is an assert, which is only active in debug mode. > Could this be a compilation > issue? I have no clue what this could be. Any help is appreciated! > I think it's more likely that your application code is doing something incorrectly... -- John ------------------------------------------------------------------------------ Transform Data into Opportunity. Accelerate data analysis in your applications with Intel Data Analytics Acceleration Library. Click to learn more. http://makebettercode.com/inteldaal-eval _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
