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. Could this be a compilation issue? I have no clue what this could be. Any help is appreciated! Thanks, Harshad ------------------------------------------------------------------------------ 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
