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 <bange...@colostate.edu> 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: bange...@colostate.edu > 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 dealii+unsubscr...@googlegroups.com. > 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
