[deal.II] value_list function to compute_nonzero_normal_flux_constraints

2017-12-08 Thread Giovanna Bucci
Dear all, I am trying to use VectorTools::compute_nonzero_normal_flux_constraints to apply a normal displacement field to part of the boundary. The code works if I simply take a ConstantFunction to impose as normal flux. I have tried defining a more general function as a list of points and v

Re: [deal.II] Re: mapping_collection and step-27

2017-12-08 Thread Juan Carlos Araujo Cabarcas
Thanks again, Now it passes the *interpolate_boundary_values *call. However, it is conflicting now with: MatrixTools::apply_boundary_values (bval, system_matrix,

Re: [deal.II] getting norm of a ghosted vector in the distributed adaptive refinement code

2017-12-08 Thread Wolfgang Bangerth
On 12/08/2017 09:14 AM, Marek Čapek wrote: An error occurred in line <1099> of file in function     dealii::IndexSet dealii::TrilinosWrappers::VectorBase::locally_owned_elements() const The violated condition was:     owned_elements

[deal.II] Re: Moving vertices of parallel triangulation breaks periodic face pair match

2017-12-08 Thread Sambit Das
Hello Dr. Arndt, Thank you for your explanation. I followed your suggestions to create a periodic match check without calling GridTools::collect_periodic_faces(..). I use the PeriodicFacePairs

Re: [deal.II] Re: mapping_collection and step-27

2017-12-08 Thread Wolfgang Bangerth
Merge: fd97710 cceeb7c Author: Wolfgang Bangerth Date:   Tue Dec 5 07:24:15 2017 -0700     Merge pull request #5579 from bangerth/more-complex-instantiations     More complex instantiations However, after configuration, installation and compilation of my code I obtain the error: Scanni

[deal.II] deal.II Newsletter #14

2017-12-08 Thread Rene Gassmoeller
Hello everyone! This is deal.II newsletter #14. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #5602: Add another missing instantiation for complex-valued functions.

Re: [deal.II] getting norm of a ghosted vector in the distributed adaptive refinement code

2017-12-08 Thread gandalfhaha
I have noticed, that orderings of operations in https://www.dealii.org/8.5.0/doxygen/deal.II/step_15.html#MinimalSurfaceProblemrefine_mesh and in https://www.dealii.org/8.5.0/doxygen/deal.II/step_15.html#MinimalSurfaceProblemrefine_mesh is different, namely In the formerdistribute_dofs, in

[deal.II] getting norm of a ghosted vector in the distributed adaptive refinement code

2017-12-08 Thread Marek Čapek
Hello, I develop distributed code like in step-40 and with the adaptive mesh refinement like in step-42. I was using ghosted Trilinos MPI vectors like LA::MPI::Vector solution_vel_n; initialized by solution_vel_n.reinit (locally_owned_dofs_vel, locally_relevant_dofs_ve

Re: [deal.II] Seg fault using Extracting hanging node dofs in MPI

2017-12-08 Thread Wolfgang Bangerth
On 12/07/2017 11:55 PM, Phani Motamarri wrote: This works fine in serial but when running in parallel on 8 procs, the above function gives a seg fault at line 854 of dof_tools.cc // dof numbers on the centers of the lines bounding this // face for (unsigned int line=0; line<4; ++li

Re: [deal.II] PETSc Sparse LU Preallocation

2017-12-08 Thread Wolfgang Bangerth
On 12/08/2017 01:21 AM, Lucas Campos wrote: > Yes, you need good preconditioners! I am investigating this now. I have to admit I was overwhelmed with the plethora of options, but I am going through Barrett et al's Template Book now, and I hope it will help. That's pretty old. Take a look a

Re: [deal.II] PETSc Sparse LU Preallocation

2017-12-08 Thread Lucas Campos
> Yes, you need good preconditioners! I am investigating this now. I have to admit I was overwhelmed with the plethora of options, but I am going through Barrett et al's Template Book now, and I hope it will help. On 7 December 2017 at 18:26, Wolfgang Bangerth wrote: > > Thanks for your thoug