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:
-------------------------------------------------------- 
An error occurred in line <565> of file </home/roland/Downloads/dealii/
source/lac/trilinos_vector.cc> in function 
    void dealii::TrilinosWrappers::MPI::Vector::compress(dealii::
VectorOperation::values) 
The violated condition was:  
    result.max-result.min<1e-5 
Additional information:  
    Not all processors agree whether the last operation on this vector was 
an addition or a set operation. This will prevent the compress() operation 
from succeeding. 
 
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::Vector local_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? The program 
runs (more or less fine) when compiling in release mode.
Thanks!

-- 
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.

Reply via email to