Re: [deal.II] Re: Particle contact detection

2019-10-25 Thread Rene Gassmoeller
Hi Shahab, I think this problem does not come from the vertex_to_cell_map but from the way you try to use it to find the neighbors. After calling vertex_to_cell_map you do not need to call GridTools::find_active_cell_around_point, because you already now the cells a vertex belongs to. Try the

[deal.II] Re: Particle contact detection

2019-10-25 Thread Shahab Golshan
Hi Rene, I have a problem in finding the list neighbor cells. This is the code I wrote: * int iter3=0; for (Triangulation<3>::active_cell_iterator cell = tr.begin_active(); cell != tr.end(); ++cell) { for (unsigned int v = 0; v < GeometryInfo<3>::vertices_per_cell;

Re: [deal.II] How to output more than one value at the location of a vertex of the mesh?

2019-10-25 Thread Wolfgang Bangerth
On 10/25/19 11:06 AM, 'Andreas Hegendörfer' via deal.II User Group wrote: > >             for (unsigned int i=0; iget_fe().dofs_per_cell; ++i) >              { >             strain=fe_values[displacement].symmetric_gradient (i, q_point) > *solution [local_dof_indices[i]]+st

Re: [deal.II] How to output more than one value at the location of a vertex of the mesh?

2019-10-25 Thread 'Andreas Hegendörfer' via deal . II User Group
Hello Wolfgang, thank you very much for your response. I am sorry to send you parts of my code, but I see no other way to explain what I do. This is how I calculate the strain: template void ElasticProblem::calculate_strain() { hp::QCollection q_collection_strain; Quadrature piezo_quadrature_f

Re: [deal.II] How to output more than one value at the location of a vertex of the mesh?

2019-10-25 Thread Wolfgang Bangerth
On 10/25/19 5:03 AM, 'Andreas Hegendörfer' via deal.II User Group wrote: > > I am only able to output only one value at a support point, what makes > my solution looks continous (especially in the middle). The objective is > to make my output of the strain to look like the ouput of > DataPostpr

[deal.II] Re: optimization, SolverBFGS

2019-10-25 Thread Bruno Turcksin
Juan, Take a look at https://stackoverflow.com/questions/17131768/how-to-directly-bind-a-member-function-to-an-stdfunction-in-visual-studio-11 I advise using a lambda like the second reply suggests Best, Bruno On Friday, October 25, 2019 at 10:53:08 AM UTC-4, Juan Carlos Araujo Cabarcas wr

[deal.II] optimization, SolverBFGS

2019-10-25 Thread Juan Carlos Araujo Cabarcas
Dear all, I would like to use SolverBFGS from deal.II/optimization/ in my project, and I try to follow the documentation in https://www.dealii.org/current/doxygen/deal.II/classSolverBFGS.html First, I would like to mention that the documentation would really improve by adding a minimal exam

Re: [deal.II] Docker MPI

2019-10-25 Thread luca.heltai
You can do export PATH=$(spack location -i mpich)/bin:$PATH to prepend the path of mpich to your PATH. This is not done automatically. Alternatively, you could also run, as root: cd /usr/local spack view add . mpich which aliases all mpi commands to /usr/local/bin I usually do this in a cus

Re: [deal.II] deal stopped working with latest macOS update (10.15)

2019-10-25 Thread luca.heltai
Dear Ester, please take a look at the release page: https://github.com/dealii/dealii/releases/tag/v9.1.1 There is a package for Mac OS Catalina. In Catalina, permissions are much more restrictive, and you’ll have to explicitly enable the application to run (the first time), by right clicking o