]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
update tests/bits/find_cell_6
authorheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 30 May 2012 20:42:36 +0000 (20:42 +0000)
committerheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 30 May 2012 20:42:36 +0000 (20:42 +0000)
git-svn-id: https://svn.dealii.org/trunk@25577 0785d39b-7218-0410-832d-ea1e28bc413d

tests/bits/find_cell_6.cc
tests/bits/find_cell_6/cmp/generic

index ab1fb06b8788480f2f0d7e5e394b29e689ac91cd..ea104818285a26a28380395b540e88fef8c4cc74 100644 (file)
 //
 //----------------------------  find_cell_5.cc  ---------------------------
 
-// find_active_cell_around_point() stops working if some cells are refined.
-// this is most likely a rounding error, but it is frustrating.
-// check2 triggers the following error:
-/*
-An error occurred in line <838> of file </scratch/deal-trunk/deal.II/source/grid/grid_tools.cc> in function
-    std::pair<typename Container<dim, spacedim>::active_cell_iterator, dealii::Point<dim> > dealii::GridTools::find_active_cell_around_point(const dealii::Mapping<dim, spacedim>&, const Container<dim, spacedim>&, const dealii::Point<spacedim>&) [with int dim = 3, Container = dealii::Triangulation, int spacedim = 3, typename Container<dim, spacedim>::active_cell_iterator = dealii::TriaActiveIterator<dealii::CellAccessor<3, 3> >]
-The violated condition was: 
-    best_cell.first.state() == IteratorState::valid
-The name and call sequence of the exception was:
-    ExcPointNotFound<spacedim>(p)
-Additional Information: 
-The point <0.789000 0.313000 0.170000> could not be found inside any of the subcells of a coarse grid cell.*/
+// find_active_cell_around_point() stops working if some cells are refined in 3d.
+// this is caused by a bug in GridTools::find_cells_adjacent_to_vertex that got fixed in r 25562.
 
 #include "../tests.h"
 #include <deal.II/base/logstream.h>
@@ -67,10 +57,11 @@ void check2 ()
       if (idx==21)
        cell->set_refine_flag();      
     }
-  
   tria.execute_coarsening_and_refinement ();
+
   deallog << inside(tria, p2) << std::endl;
-  GridTools::find_active_cell_around_point(tria, p2); //triggers exception
+
+  GridTools::find_active_cell_around_point(tria, p2); //triggered exception
 }
 
 void check1 ()
@@ -79,9 +70,9 @@ void check1 ()
   GridGenerator::hyper_cube (tria);
   tria.refine_global (3);
 
-  for (int i=0;i<5;++i)
+  for (int i=0;i<3;++i)
     {
-      for (int j=0;j<10000;++j)
+      for (int j=0;j<1000;++j)
        {
          Point<3> p ((rand()%1000)/1000.0,(rand()%1000)/1000.0,(rand()%1000)/1000.0);
          if (!inside(tria, p))
index 1b77c2663db464b69004a61c1a1e636d751ffbd9..e46c0ac61e74783bbee0bfd4c4830ca7575e8a49 100644 (file)
@@ -1,3 +1,3 @@
 
-DEAL::1.2
-DEAL::<0.00000 0.500000> <0.500000 0.500000> <0.00000 1.00000> <0.500000 1.00000> 
+DEAL::1
+DEAL::1

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.