[deal.II] Questions related to FE_Nothing

Hello all, I am using deal.ii to solve a system of Hamilton Jacobi equations. I need to use FE_Nothing element in a part of my domain. I read a couple of posts in this forum regarding this but I have some questions. I will appreciate if anyone can answer these. Since I need to use MPI, and I

[deal.II] Using DataOut with MappingCollection

Dear all, As described in the forum: https://groups.google.com/forum/?utm_source=digest&utm_medium=email#!topic/dealii/ijJFbOTZxkk I would like to reproduce step-27 but with curved boundaries with the use of MappingCollection. Everything seems to work fine, but I noticed that data_out does not

Re: [deal.II] Re: Errors when using MUMPS/PETSc LU

> It shows the same issues, as expected. I will follow your advice and try to > use a different solver. Then you have to simplify your problem as much as possible until we can reproduce it. > Still, would it be possible for you to comment a bit more on those > MPI_Op_free errors? This can happ

Re: [deal.II] installation process on Ubuntu

Hi Juan Carlos, Reading the error more closely it could be a compiler issue. But you’ve caught me literally at the end of my work day, so I can’t investigate this issue right now. I would say that commit babb6ae7b2ef589a0998383214087fd60466beb6

[deal.II] Re: installation process on Ubuntu

Dear Juan Carlos, Can you please tell me what the commit hash of the deal.II respository that you're working from is? I did briefly introduce a bug that showed similar symptoms but I thought that I fixed it. It might also have to do with your version of Trilinos, which appears to be older than

Re: [deal.II] Re: Errors when using MUMPS/PETSc LU

On 11/17/2017 07:37 AM, Lucas Campos wrote: Invalid MPI_Op, error stack: MPI_Op_free(111): MPI_Op_free(op=0x7ffcf6298dac) failed MPI_Op_free(75).: Null Op pointer If this error would lead the the issues I am having, is up to discussion. I tried using PETSc's SolverPr

Re: [deal.II] Re: Errors when using MUMPS/PETSc LU

Dear Timo, Also, do you have recommendations for Valgrind flags other than --track-origins=yes --leak-check=full? Lucas On 17 November 2017 at 15:37, Lucas Campos wrote: > Dear Timo, > > Thanks for you advice. I am running the program in three different > computers -- my notebook, my research

Re: [deal.II] Re: Errors when using MUMPS/PETSc LU

Dear Timo, Thanks for you advice. I am running the program in three different computers -- my notebook, my research group's server and the local cluster. In all of them I have this (small) change to suddenly find the nan. According to MUST, there is clearly a problem inside deal.II, MUMPS or PETS

Re: [deal.II] Re: Read in tethex grid

This code is somewhere else. There is already an exception in tethex when the data-size is not equal to sizeof(double). So I think it is ok. Best regards, Felix -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/de

Re: [deal.II] Re: Errors when using MUMPS/PETSc LU

Lucas, those kind of bugs are hard to find. Honestly, the bug could still be in your code, inside deal.II, inside PETSc, inside MUMPS, or related to the software/hardware you are running on. I know this won't be of much help, but I would suggest you try a different solver to see if MUMPS is the p

[deal.II] Re: Errors when using MUMPS/PETSc LU

Sorry, I forgot to include the link to MUST. Here it is: https://doc.itc.rwth-aachen.de/display/CCP/Project+MUST On Friday, 17 November 2017 14:12:18 UTC+1, Lucas Campos wrote: > > Dear all, > > First of all, a bit of context: > I am trying to debug an error in my application where randomly I sta

[deal.II] Errors when using MUMPS/PETSc LU

Dear all, First of all, a bit of context: I am trying to debug an error in my application where randomly I start seeing nan's. The probability of this increases with the number of MPI processors I use, so it looks like it is a data race of some sort. Any advice on the best way to find the error

Re: [deal.II] A cube rather than a ball with GridGenerator::hyper_ball

Dear Alberto, > GridGenerator::hyper_ball ( triangulation, center, cell_radius ); > > static SphericalManifold volume_boundary; > > > triangulation.refine_global ( 3 ); > > triangulation.set_all_manifold_ids_on_boundary(0); > > triangulation.set_manifold (0, volume_boundary); Y