RE: [deal.II] Computing the solution gradient at the quadrature point on a face

2021-08-04 Thread Michael Li
Hi Jean-Paul, Yes, I set up the boundary check when taking care of the traction. It worked as expected using Physics::Transformations::nansons_formula(face_normal, F). To make a quick test, I did not linearize the load stiffness but just used the load at last Newton-Raphson step, so it is like a

Re: [deal.II] Re: Grid deformation after load-balancing

2021-08-04 Thread Wolfgang Bangerth
On 8/4/21 8:16 AM, blais...@gmail.com wrote: I would not know how to do this either. It would require that every manifold has some sort of "center of mass" variable that would be used to define the centroid of the manifold, but that would be some pretty heavy stuff to refactor and I am not

Re: [deal.II] Re: Grid deformation after load-balancing

2021-08-04 Thread shahab.g...@gmail.com
Dear Wolfgang, Bruno and Peter, Thank you for your helps. It was my fault for not reading the entire function's documentation before asking the question here. Best Shahab On Wednesday, August 4, 2021 at 10:16:44 AM UTC-4 blais...@gmail.com wrote: > Yeah I expected this... > I would not know

Re: [deal.II] Re: Grid deformation after load-balancing

2021-08-04 Thread blais...@gmail.com
Yeah I expected this... I would not know how to do this either. It would require that every manifold has some sort of "center of mass" variable that would be used to define the centroid of the manifold, but that would be some pretty heavy stuff to refactor and I am not sure it would lead to a

Re: [deal.II] Output file

2021-08-04 Thread Wolfgang Bangerth
On 8/4/21 7:49 AM, Hermes Sampedro wrote: I am working with step-29 and I was wondering what is the best way to import the solution vectors (real, imaginary and intensity solutions) into Matlab. Is it possible to export one point of the solution into a txt file? Or what would be the easiest

Re: [deal.II] Re: Grid deformation after load-balancing

2021-08-04 Thread Wolfgang Bangerth
On 8/4/21 6:57 AM, blais...@gmail.com wrote: Is that not a bit of an ill-defined behavior for the triangulations? Shouldnt the manifolds also translate themselves when the triangulation is translated? This is a discussion we had at some point when we looked at the GridTools::transform()

Re: [deal.II] step-35 in modified

2021-08-04 Thread Wolfgang Bangerth
On 8/4/21 4:20 AM, SebG wrote: a colleague and I have also worked with step-35 and we programmed a parallelized version whose layout is similar to that of step-32. Currently, the code is not suitable to be directly used as a tutorial, but, if you are interested, we could refactor it and send

[deal.II] Re: Grid deformation after load-balancing

2021-08-04 Thread 'peterrum' via deal.II User Group
I must admit I don't know how to transform manifolds, let alone if it is possible! Maybe someone else has an idea? Separate if this, could you try out what I have described and see if it resolves the problem. At least than, we have an idea what the cause is what is at least a better stage we

[deal.II] Output file

2021-08-04 Thread Hermes Sampedro
Hello, I am working with step-29 and I was wondering what is the best way to import the solution vectors (real, imaginary and intensity solutions) into Matlab. Is it possible to export one point of the solution into a txt file? Or what would be the easiest way using Dealii for that? Thank

[deal.II] Re: Grid deformation after load-balancing

2021-08-04 Thread blais...@gmail.com
Peter this is what I suspected :). Is that not a bit of an ill-defined behavior for the triangulations? Shouldnt the manifolds also translate themselves when the triangulation is translated? Clearly, it is possible to reset the manifold every time the triangulation is moved, but this becomes a

Re: [deal.II] step-35 in modified

2021-08-04 Thread SebG
Dear all, a colleague and I have also worked with step-35 and we programmed a parallelized version whose layout is similar to that of step-32. Currently, the code is not suitable to be directly used as a tutorial, but, if you are interested, we could refactor it and send it through github as a

[deal.II] Interpretation of convergence results

2021-08-04 Thread SebG
Dear all, I would like to address a problem related to the interpretation of the results of a convergence test. I am performing a temporal convergence test to test a solver for advection-diffusion problems. The test is based on a manufactured solution on a unit square domain. The Peclet number

Re: [deal.II] profiling and parallel performance of deal.II user codes

2021-08-04 Thread richard....@gmail.com
Hi all, thanks for the many tips and suggestions, I really appreciate you spending your time and effort helping me out! I set up valgrind and kcachegrind, which I found exceptionally easy and can get started now - perfect! In case anyone reads this in the future, I had to use "mpirun -n 4

Re: [deal.II] Use a coarse grid solution as initial condition for a finer grid

2021-08-04 Thread vachanpo...@gmail.com
That brings a big simplification, at least when the aim is to just restart a simulation and not to map the solution onto some other grid. Thanks for all the help! On Tuesday, August 3, 2021 at 9:18:45 PM UTC+5:30 Wolfgang Bangerth wrote: > On 8/3/21 8:38 AM, vachan potluri wrote: > > > > I

[deal.II] Re: Grid deformation after load-balancing

2021-08-04 Thread 'peterrum' via deal.II User Group
Hi Shahab, My best guess is that the attached manifold is not consistent to what you need after the grid transformations. You have to be aware that the attached manifold is used for generating new vertices (as it is needed during AMR/load-balancing in regions that is newly owned by a rank).

Re: [deal.II] profiling and parallel performance of deal.II user codes

2021-08-04 Thread heena patel
Dear Richard, I recently attend summer school in which I learn about paraver but that is the visual tool for performance from Barcelona SuperComputing center. But it comes with one more pre- processor and a post-processor. Check the link below.

Re: [deal.II] Able to use umfpack in debug configuration, but not in release configuration

2021-08-04 Thread Jean-Paul Pelteret
Hi Lucas, I’m glad that you managed to solve this issue yourself. > Given that the problem was a silly mistake on my end, ought I delete the > forum post? Given that you’ve explained what you did to fix the problem, this is useful information for all. So we’ll leave the thread in place just

Re: [deal.II] Computing the solution gradient at the quadrature point on a face

2021-08-04 Thread Jean-Paul Pelteret
Hi Michael, The one thing that stands out about the snippet of code that you shared is that there is no check to if (cell->face(face)->at_boundary() == true && )... Do you still have something like that? Otherwise, after a cursory glance, everything appears to be alright with this. It mimics

Re: [deal.II] profiling and parallel performance of deal.II user codes

2021-08-04 Thread Jean-Paul Pelteret
To add to this great list of tools (which we should document on our Wiki), there’s also the LIKWID performance monitoring and benchmarking suite which is able to provide some very low-level metrics and tools to help benchmarking. https://github.com/RRZE-HPC/likwid