Re: [deal.II] Re: Epetra_FEVector not being created as a LinearMap

2018-05-14 Thread josh . hanophy
Pull request submitted. On Saturday, May 12, 2018 at 2:33:36 AM UTC-7, Daniel Arndt wrote: > > Josh, > > The key issue for my application is whether the make_trilinos_map method >> of the IndexSet class is called with overlapping equal to true or false. At >> the moment, when make_trilinos_map i

Re: [deal.II] Re: Epetra_FEVector not being created as a LinearMap

2018-05-11 Thread josh . hanophy
The key issue for my application is whether the make_trilinos_map method of the IndexSet class is called with overlapping equal to true or false. At the moment, when make_trilinos_map is called from the reinit method for Vector in the TrilinosWrappers∷MPI namespace, overlapping is always set t

Re: [deal.II] Re: Epetra_FEVector not being created as a LinearMap

2018-05-11 Thread josh . hanophy
The key issue for my application is whether the make_trilinos_map method of the IndexSet class is called with overlapping equal to true or false. At the moment, when make_trilinos_map is called from the reinit method for Vector in the TrilinosWrappers∷MPI namespace, overlapping is always set t

[deal.II] Re: Epetra_FEVector not being created as a LinearMap

2018-05-10 Thread josh . hanophy
I solved my problem. In the code I included below, I was using the Vector object from the namespace dealii::LinearAlgebraTrilinos::MPI. When I switch to using Vector objects from the namespace dealii::LinearAlgebra::EpetraWrappers (and add the trilinos_epetra_vector.h include file), the Epetr

Re: [deal.II] Epetra_FEVector not being created as a LinearMap

2018-05-08 Thread josh . hanophy
Yes I apologize for not including code. A program that runs is included below. Note I tried to remove all the non-deal.ii include files, but I might have missed a few. Principally, I'm not sure whether deal.ii means to create an Epetra_FEVector with LinearMap() equal to true or not. In the sour

Re: [deal.II] Epetra_FEVector not being created as a LinearMap

2018-05-08 Thread josh . hanophy
> > Yes I apologize for not including code. A program that runs is included > below. Principally, I'm not sure whether deal.ii means to create an > Epetra_FEVector with LinearMap() equal to true or not. In the source code, > the word "linear" is used when constructing Epetra_FEVector object.

[deal.II] Epetra_FEVector not being created as a LinearMap

2018-05-08 Thread josh . hanophy
I am using Dealii with the Ifpack_Hypre package in Trilinos. Everything works well except that I need an Epetra_FEVector vector for which vector.Map().LinearMap() evaluates to true. However, it appears I am getting a discontinuous mapping, but I don't understand why based on what I see in the d