[deal.II] Re: Questions on parallel::distributed::GridRefinement::refine_and_coarsen_fixed_number in deal.II 9.0.1 released

2019-10-10 Thread Wolfgang Bangerth
> Dear Wolfgang Bangerth, > >   Recently, I use your deal.II (Version 9.0.1) It's really not *my* program any more -- more than 250 people have contributed to it by now :-) > to solve a group of PDEs > (haveing a wave-like solution whose amplitude is modulated by another > smoother > fie

Re: [deal.II] Large Memory Consumption making BlockDynamicSparsityPattern

2019-10-10 Thread Wolfgang Bangerth
On 10/9/19 10:26 AM, Matteo Frigo wrote: > > I'm saying that making a DynamicSparsityPattern by using the procedure > described above is more expensive (from the memory point of view) respect > using TrilinosWrappers::BlockSparsityPattern. > I noted this trouble trying to run some test cases wit

Re: [deal.II] Arpack solver reports 4294967295 iterations before converging

2019-10-10 Thread Wolfgang Bangerth
On 10/9/19 7:47 AM, 'Maxi Miller' via deal.II User Group wrote: > This number already comes from solver_control, given as > > solver_control.last_step (); Hm, that's really quite unfortunate. It took me a couple minutes of looking into the code, but I *think* what is happening is that the Arpack

[deal.II] Re: Query regarding DoFTools::dof_indices_with_subdomain_association()

2019-10-10 Thread Doug Shi-Dong
Hello Vachan, Sounds like you're implementing nodal DG, hence why you only need values and gradients at the *quadrature* points from the neighbor. Something you might want to consider rather than communicating them at solution points, in case you ever decide to overintegrate. You could still a

Re: [deal.II] Convergence rate

2019-10-10 Thread luca.heltai
Dear Felix, by any chance, did you take a look at step 10? https://www.dealii.org/current/doxygen/deal.II/step_10.html This step explains a little bit what to do when you want to solve on curved domains with high order finite elements. In particular, you need to ensure that the mapping you ar

[deal.II] Re: Convergence rate

2019-10-10 Thread Félix Bunel
Thanks for your answer ! I'm working in 2D on a disk ! This is my mesh on the 4 cycle of refinement (which I never use in practice because it's not refined enough for what I want to do). I don't think I have a mapping object... What i do to initalize my system on phi is : //On génère le dofhan

[deal.II] Re: Convergence rate

2019-10-10 Thread Bruno Blais
A quick question, since you are working on a sphere, are you specifying a mapping of the same order as your phi? On Wednesday, 9 October 2019 08:57:45 UTC-4, Félix Bunel wrote: > > Hello everyone. > > I'm having some trouble to understand the convergence rate i'm observing > in my code. > > Here