Re: [deal.II] Extracting a component from Block FESystem and DofHandler

2016-08-16 Thread Wolfgang Bangerth
On 08/16/2016 09:34 AM, Spencer Patty wrote: Supposing velocity_block represents the block number associated to the velocity components, I noticed that solution.block{velocity_block} would be a TrilinosWrapper::MPI::Vector of the velocity as desired and it is possible to extract the sub FiniteEl

Re: [deal.II] Re: Extracting a component from Block FESystem and DofHandler

2016-08-16 Thread Bruno Turcksin
Sorry, I misunderstood what you were trying to do. Like you found out, the problem is that you cannot "extract" the DoFHandler because the DoFHandler has no block structure and it will have to be rebuilt. Everything else can be extracted. If your interface returns the Triangulation and the FESystem

[deal.II] Re: Extracting a component from Block FESystem and DofHandler

2016-08-16 Thread Spencer Patty
On Tuesday, August 16, 2016 at 12:35:05 PM UTC-5, Spencer Patty wrote: > > Bruno, > > Thanks for responding! > > > > On Tuesday, August 16, 2016 at 11:16:41 AM UTC-5, Bruno Turcksin wrote: >> >> Spencer, >> >> I think that passing around the entire system is the way to go but what I >> would d

[deal.II] Re: Extracting a component from Block FESystem and DofHandler

2016-08-16 Thread Spencer Patty
Bruno, Thanks for responding! On Tuesday, August 16, 2016 at 11:16:41 AM UTC-5, Bruno Turcksin wrote: > > Spencer, > > I think that passing around the entire system is the way to go but what I > would do is to use component mask ( > http://dealii.org/developer/doxygen/deal.II/DEALGlossary.ht

[deal.II] Re: Extracting a component from Block FESystem and DofHandler

2016-08-16 Thread Bruno Turcksin
Spencer, I think that passing around the entire system is the way to go but what I would do is to use component mask (http://dealii.org/developer/doxygen/deal.II/DEALGlossary.html#GlossComponentMask) and block mask (http://dealii.org/developer/doxygen/deal.II/DEALGlossary.html#GlossBlockMask)

[deal.II] Extracting a component from Block FESystem and DofHandler

2016-08-16 Thread Spencer Patty
Dear All, I have a class which represents a velocity field (I am using a variant of the pressure correction navier stokes model with two phases) and a class which represents a transport model along a provided velocity field. I am adding in some terms into the velocity model (currently surface

Re: [deal.II] Question concerning BlockSparsityPattern.copy_from() member function

2016-08-16 Thread Dustin Kumor
Dear Martin, the update to the most recent development version solved the problem. Now time needed to finish the copy operation is in the expected range. Thanks a lot for your support. Best, Dustin -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

Re: [deal.II] Question concerning BlockSparsityPattern.copy_from() member function

2016-08-16 Thread Dustin Kumor
Dear Martin, I'm using 8.5.0-pre. But I haven't updated the code for a long time. I will do this first, then compile the library again and let you know about the results. Best, Dustin -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://grou

[deal.II] Re: Data exchanges in IO

2016-08-16 Thread Daniel Arndt
Bruce, [...] 3. Like in 2., but for adding values to individual elements. Use > VectorOperation::add. > > I guess 3) is obscure. Does this mean that you are adding a new element > where none existed before or that you are incrementing an existing element? > If the former, aren't the total numb