Dear Wolfgang, I have not been able to find the issue. I am considering looking for an alternative. I am using 2 components to describe a complex solution (for real and imaginary parts). Assuming that I keep the ordering by default and that I export my solution vector (to work in an external software), what I need in the end is to know for each value of the solution vector, if it corresponds to the real or imaginary part. I realized by debugging that the order is not consistent and it changes depending on the number of dof. I would like to ask if there is any function or if you could suggest a way to get the vector that contains information about the component type for each node (for example a mask type, where 0 is real and 1 imaginary component). The way I implemented and managed the components, assembly, etc is identical to the one described in step-29.
I would really appreciate your help on this matter. Thank you very much. Regards Hermes El martes, 5 de abril de 2022 a las 16:41:49 UTC+2, Wolfgang Bangerth escribió: > On 4/5/22 04:47, Hermes Sampedro wrote: > > > > Now I am trying to order the dof according to the components. However, I > am > > experiencing that when using /DoFRenumbering::component_wise/ and > running with > > more than 1 MPI rank, my solution vector has inf values. Which, does not > > happen when using only 1 process, or using the default dof ordering. I > am not > > sure where the problem could be. May I ask for some suggestions to see > if I > > can identify where the issue can be? > > This is an excellent problem to check your debugging skills. It *may* of > course be related in some way to the renumbering, but there appears to be > no > logical connection between the two. On the other hand, it is easy to check > for > the presence of inf values (for example, by outputting the norm of a > vector). > There are two ways for inf values to appear in the solution vector: The > right > hand side vector has inf values, or the iterative solver does not > converge. > Both are easy enough to test. If it's the solver, you need to figure out > why > that is. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > www: http://www.math.colostate.edu/~bangerth/ > > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/39dd3a68-2785-4076-aa94-3f372e9c84a2n%40googlegroups.com.
