I am coming out of a loop where I am doing hanging_node_constraints.distribute_local_to_global(cell_residual, local_dof_indices, local_residual);
so I already was wondering why I get that error. Am Dienstag, 10. Oktober 2017 20:22:21 UTC+2 schrieb Timo Heister: > > It might also be that you are forgetting to compress() your vector > after you assemble it. The error would then trigger on the next > operation, which is the set_zero() in this case. What is the last > thing you have been doing with this vector? > > On Tue, Oct 10, 2017 at 9:18 AM, Wolfgang Bangerth > <[email protected] <javascript:>> wrote: > > On 10/10/2017 05:15 AM, 'Maxi Miller' via deal.II User Group wrote: > >> > >> > >> When trying to run my program on an i7-CPU with 4/8 threads, I get the > >> following error when calling set_zero() from the > ConstraintMatrix-class: > >> | > >> -------------------------------------------------------- > >> Anerror occurred inline <565>of file > >> </home/roland/Downloads/dealii/source/lac/trilinos_vector.cc>infunction > >> > >> > voiddealii::TrilinosWrappers::MPI::Vector::compress(dealii::VectorOperation::values) > > > >> Theviolated condition was: > >> result.max-result.min<1e-5 > >> Additionalinformation: > >> Notall processors agree whether the lastoperation on thisvector was an > >> addition ora setoperation.Thiswill prevent the compress()operation > >> fromsucceeding. > >> > >> Stacktrace: > >> ----------- > >> #0 /opt/dealII/lib/libdeal_II.g.so.9.0.0-pre: > >> > dealii::TrilinosWrappers::MPI::Vector::compress(dealii::VectorOperation::values) > > > >> #1 /opt/dealII/lib/libdeal_II.g.so.9.0.0-pre: > >> #2 /opt/dealII/lib/libdeal_II.g.so.9.0.0-pre: void > >> > dealii::ConstraintMatrix::set_zero<dealii::TrilinosWrappers::MPI::Vector>(dealii::TrilinosWrappers::MPI::Vector&) > > > >> const > >> #3 main: MinimalSurfaceProblem<2>::compute_residual(double, double, > >> double) > >> #4 main: MinimalSurfaceProblem<2>::run() > >> #5 main: main > >> -------------------------------------------------------- > >> | > >> > >> The code where this line arises is > >> | > >> > >> > LinearAlgebraTrilinos::MPI::Vectorlocal_residual(dof_handler.locally_owned_dofs(),mpi_communicator); > > > >> //Code in between for setting values > >> > >> hanging_node_constraints.set_zero(local_residual); > >> > >> | > >> > >> > >> It only happens if running with more than 4 threads, i.e. with mpirun > -n 5 > >> or larger. What is the problem here, and how can I fix that? > > > > That would sound like a bug somewhere. Can you create a small program > that > > demonstrates this and that you could share so that others could look at > it? > > It doesn't have to do anything useful, just demonstrate the problem. For > > example, you could start from your existing program and strip assembly > and > > solver from it. > > > > Best > > W. > > > > > > > > -- > > ------------------------------------------------------------------------ > > Wolfgang Bangerth email: [email protected] > <javascript:> > > www: > > > https://urldefense.proofpoint.com/v2/url?u=http-3A__www.math.colostate.edu_-7Ebangerth_&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=R5lvg9JC99XvuTgScgbY_QFS80R7PEA2q0EPwDy7VQw&m=FDYQy9r_-sUm-_WjnWoeIPlURA8ea5TpOVtCtg5BBE4&s=Kcs0rW4kEvR2VF170rElcgkdyMqv9obffaWQEm_ZAm4&e= > > > -- > > The deal.II project is located at > > > https://urldefense.proofpoint.com/v2/url?u=http-3A__www.dealii.org_&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=R5lvg9JC99XvuTgScgbY_QFS80R7PEA2q0EPwDy7VQw&m=FDYQy9r_-sUm-_WjnWoeIPlURA8ea5TpOVtCtg5BBE4&s=SKIvxONbceO3FmrAH9KojFbKCSd-anplulITEuNMEp0&e= > > > For mailing list/forum options, see > > > https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_forum_dealii-3Fhl-3Den&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=R5lvg9JC99XvuTgScgbY_QFS80R7PEA2q0EPwDy7VQw&m=FDYQy9r_-sUm-_WjnWoeIPlURA8ea5TpOVtCtg5BBE4&s=hKawPWSQ0WMEnRH1iQtU3I0vhGKZUx4O0WI95oK_5B0&e= > > > --- 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] <javascript:>. > > For more options, visit > > > https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_optout&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=R5lvg9JC99XvuTgScgbY_QFS80R7PEA2q0EPwDy7VQw&m=FDYQy9r_-sUm-_WjnWoeIPlURA8ea5TpOVtCtg5BBE4&s=SpbJ1nIWvdzEUmE8TiMiISu1pTjD0x7L1YWbiWAiCK8&e= > > > . > > > > -- > Timo Heister > http://www.math.clemson.edu/~heister/ > -- 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]. For more options, visit https://groups.google.com/d/optout.
