Re: [deal.II] Hashing / mapping Point locations to material ID

2021-02-06 Thread Corbin Foucart
Thank you for the fast and clear answer, Wolfgang. Just add such a hash function for Point. It's true that two nearby points > will > have different hash values, but that shouldn't stop you from using hashes > if > you query the same points over and over. This is definitely doable. > That

Re: [deal.II] Hashing / mapping Point locations to material ID

2021-02-06 Thread Wolfgang Bangerth
On 2/6/21 1:55 PM, Corbin Foucart wrote: The compiler is rightly complaining that there isn't a hash method for the Point object, and it seems like this could indeed be dangerous, because of floating point comparison later between Point objects during lookup. Just add such a hash function fo

[deal.II] Hashing / mapping Point locations to material ID

2021-02-06 Thread Corbin Foucart
Hello all, I'm attempting to set material ID based on x-location. I'd like to create std::unordered_map between x-locations of cell centers and material id, and then refer to that map to quickly set material_id later. The compiler is rightly complaining that there isn't a hash method for the

Re: [deal.II] Outputting auxiliary cell-level data along with solution data

2021-02-06 Thread Corbin Foucart
Hi Praveen, Thanks for the helpful answer. I'm new to Visit; would you mind specifying how to view the subset of data (e.g., faces with boundary ID 0)? Best, Corbin On Thursday, November 15, 2018 at 9:31:16 PM UTC-5 Praveen C wrote: > You can use this > > https://dealii.org/developer/doxygen/

[deal.II] Understanding loops in matrix-free

2021-02-06 Thread 'David' via deal.II User Group
Sorry for messing up the topic. I should be Understanding loops in matrix-free. I wanted to insert a figure of the source code rather than the google groups formatting and it didn't work for some reason. David schrieb am Samstag, 6. Februar 2021 um 19:15:07 UTC+1: > Hi there, > > I'm currently

[deal.II] Understanding

2021-02-06 Thread 'David' via deal.II User Group
Hi there, I'm currently trying to pack my cell operations into one of the matrix-free loop functions. In my first version, I implemented the loops manually by using (sorry for the odd formatting): ``` local_apply_cell(*matrix_free, system_rhs,

Re: [deal.II] Re: Multilevel DoFRenumbering for a distributed Triangulation

2021-02-06 Thread julius...@gmail.com
Perfect, thank you very much. Unfortunately, the make_sparsity_pattern() method I was looking for is displayed after make_flux_sparsity_pattern() in my online docu., so I overlooked them, sorry for that. That's a good starting point to add a method for each level and I need some flux sparsity f