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]
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].
For more options, visit https://groups.google.com/d/optout.
