Mate,
> ============================ Linking step-18
> m...@r12:/home/1/mate/software/deal.II.source/examples/step-18$ mpiexec -n
> 2 ./step-18
> Timestep 1 at time 1
> Cycle 0:
> Number of active cells: 3712 (by partition: 1856+1856)
> Number of degrees of freedom: 17226 (by partition: 8488+8738)
> Assembling system... norm of rhs is 2.31518e+10
> Solver converged in 133 iterations.
> Updating quadrature point data...
> Cycle 1:
> --------------------------------------------------------
> An error occurred in line <670> of file
> </home/1/mate/software/deal.II.source/deal.II/include/grid/tria_objects.h>
> in function
> unsigned int&
> dealii::internal::Triangulation::TriaObjects<G>::user_index(unsigned int)
> [with G = dealii::internal::Triangulation::TriaObject<3>]
> The violated condition was:
> user_data_type == data_unknown || user_data_type == data_index
> The name and call sequence of the exception was:
> ExcPointerIndexClash()
> Additional Information:
> (none)
I believe you need the following patch for a problem we have already fixed in
the current development sources.
Best
W.
Index: step-18.cc
===================================================================
--- step-18.cc (revision 18901)
+++ step-18.cc (revision 18902)
@@ -2509,11 +2509,10 @@
for (typename Triangulation<dim>::active_cell_iterator
cell = triangulation.begin_active();
cell != triangulation.end(); ++cell)
- {
- cell->clear_user_pointer();
- if (cell->subdomain_id() == this_mpi_process)
- ++our_cells;
- }
+ if (cell->subdomain_id() == this_mpi_process)
+ ++our_cells;
+
+ triangulation.clear_user_data();
// Next, allocate as many quadrature
// objects as we need. Since the
-------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
