]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Merge from mainline.
authorBruno Turcksin <bruno.turcksin@gmail.com>
Thu, 2 May 2013 14:24:19 +0000 (14:24 +0000)
committerBruno Turcksin <bruno.turcksin@gmail.com>
Thu, 2 May 2013 14:24:19 +0000 (14:24 +0000)
git-svn-id: https://svn.dealii.org/branches/branch_bigger_global_dof_indices_4@29430 0785d39b-7218-0410-832d-ea1e28bc413d

1  2 
deal.II/include/deal.II/lac/arpack_solver.h
deal.II/include/deal.II/lac/matrix_lib.h
deal.II/include/deal.II/numerics/vector_tools.h
deal.II/source/distributed/tria.cc
deal.II/source/lac/block_sparsity_pattern.cc
tests/distributed_grids/2d_refinement_07.cc
tests/distributed_grids/2d_refinement_07/cmp/generic
tests/mpi/matrix_free_01.cc
tests/mpi/matrix_free_03.cc
tests/mpi/matrix_free_04.cc

Simple merge
index 3df0ecb92ce0b940236c3b916e163b065219446a,b09aeb3df9797193831cafca53d9c94977f098d2..cf8e5f3ba82df80d67f5c353a616e7feed127da9
@@@ -296,18 -296,18 +296,18 @@@ template <class SparsityPatternBase
  void
  BlockSparsityPatternBase<SparsityPatternBase>::print_gnuplot(std::ostream &out) const
  {
 -  unsigned int k=0;
 -  for (unsigned int ib=0; ib<n_block_rows(); ++ib)
 +  size_type k=0;
 +  for (size_type ib=0; ib<n_block_rows(); ++ib)
      {
 -      for (unsigned int i=0; i<block(ib,0).n_rows(); ++i)
 +      for (size_type i=0; i<block(ib,0).n_rows(); ++i)
          {
 -          unsigned int l=0;
 -          for (unsigned int jb=0; jb<n_block_cols(); ++jb)
 +          size_type l=0;
 +          for (size_type jb=0; jb<n_block_cols(); ++jb)
              {
                const SparsityPatternBase &b = block(ib,jb);
 -              for (unsigned int j=0; j<b.n_cols(); ++j)
 +              for (size_type j=0; j<b.n_cols(); ++j)
                  if (b.exists(i,j))
-                   out << l+j << " " << -static_cast<signed int>(i+k) << std::endl;;
+                   out << l+j << " " << -static_cast<signed int>(i+k) << std::endl;
                l += b.n_cols();
              }
          }
index 0000000000000000000000000000000000000000,61292d4e6fd2d5b5b940a41097ce41fe3c9dfc30..6e2fc3afee5c61bb6fcbb9b8918c732afa0081a4
mode 000000,100644..100644
--- /dev/null
@@@ -1,0 -1,74 +1,74 @@@
 -  for (unsigned int i=0;i<35;++i)
+ //---------------------------------------------------------------------------
+ //    $Id: 3d_refinement_01.cc 29075 2013-03-27 15:07:57Z heister $
+ //    Version: $Name$
+ //
+ //    Copyright (C) 2008, 2009 by the deal.II authors
+ //
+ //    This file is subject to QPL and may not be  distributed
+ //    without copyright and license information. Please refer
+ //    to the file deal.II/doc/license.html for the  text  and
+ //    further information on this license.
+ //
+ //---------------------------------------------------------------------------
+ // Test that we abort if you refine further than p4est::max_level
+ #include "../tests.h"
+ #include "coarse_grid_common.h"
+ #include <deal.II/base/logstream.h>
+ #include <deal.II/base/tensor.h>
+ #include <deal.II/grid/tria.h>
+ #include <deal.II/distributed/tria.h>
+ #include <deal.II/grid/tria_accessor.h>
+ #include <deal.II/grid/tria_iterator.h>
+ #include <deal.II/grid/grid_generator.h>
+ #include <deal.II/grid/grid_out.h>
+ #include <fstream>
+ template<int dim>
+ void test(std::ostream& /*out*/)
+ {
+   parallel::distributed::Triangulation<dim> tr(MPI_COMM_WORLD);
+   GridGenerator::hyper_cube(tr);
 -      std::cout << "cells: " << tr.n_active_cells() << " level:" << tr.n_levels() << std::endl;
++  for (unsigned int i=0;i<30;++i)
+     {
++      deallog << "cells: " << tr.n_active_cells() << " level:" << tr.n_levels() << std::endl;
+       typename parallel::distributed::Triangulation<dim>::cell_iterator it;
+       it=tr.begin_active();
+       while (it->level()<static_cast<int>(tr.n_levels()-1))
+       ++it;
+       
+       it->set_refine_flag();
+       try
+       {
+         tr.execute_coarsening_and_refinement ();
+       }
+       catch (ExceptionBase &e)
+       {
+         deallog << e.get_exc_name() << std::endl;
+       }
+     }
+   
+ }
+ int main(int argc, char *argv[])
+ {
+   deal_II_exceptions::disable_abort_on_exception();
+   Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1);
+   std::ofstream logfile("2d_refinement_07/output");
+   deallog.attach(logfile);
+   deallog.depth_console(0);
+   deallog.threshold_double(1.e-10);
+   deallog.push("2d");
+   test<2>(logfile);
+   deallog.pop();
+ }
index 0000000000000000000000000000000000000000,eb291e12252c33667c68494a10849db3726c85e2..909caf44e777e124c984ea1c75f5be72d9d7ebf6
mode 000000,100644..100644
--- /dev/null
@@@ -1,0 -1,21 +1,32 @@@
 -DEAL:3d::cells: 1 level:1
 -DEAL:3d::cells: 8 level:2
 -DEAL:3d::cells: 15 level:3
 -DEAL:3d::cells: 22 level:4
 -DEAL:3d::cells: 29 level:5
 -DEAL:3d::cells: 36 level:6
 -DEAL:3d::cells: 43 level:7
 -DEAL:3d::cells: 50 level:8
 -DEAL:3d::cells: 57 level:9
 -DEAL:3d::cells: 64 level:10
 -DEAL:3d::cells: 71 level:11
 -DEAL:3d::cells: 78 level:12
 -DEAL:3d::cells: 85 level:13
 -DEAL:3d::cells: 92 level:14
 -DEAL:3d::cells: 99 level:15
 -DEAL:3d::cells: 106 level:16
 -DEAL:3d::cells: 113 level:17
 -DEAL:3d::cells: 120 level:18
 -DEAL:3d::cells: 127 level:19
 -DEAL:3d::ExcMessage("Fatal Error: maximum refinement level of p4est reached.")
++DEAL:2d::cells: 1 level:1
++DEAL:2d::cells: 4 level:2
++DEAL:2d::cells: 7 level:3
++DEAL:2d::cells: 10 level:4
++DEAL:2d::cells: 13 level:5
++DEAL:2d::cells: 16 level:6
++DEAL:2d::cells: 19 level:7
++DEAL:2d::cells: 22 level:8
++DEAL:2d::cells: 25 level:9
++DEAL:2d::cells: 28 level:10
++DEAL:2d::cells: 31 level:11
++DEAL:2d::cells: 34 level:12
++DEAL:2d::cells: 37 level:13
++DEAL:2d::cells: 40 level:14
++DEAL:2d::cells: 43 level:15
++DEAL:2d::cells: 46 level:16
++DEAL:2d::cells: 49 level:17
++DEAL:2d::cells: 52 level:18
++DEAL:2d::cells: 55 level:19
++DEAL:2d::cells: 58 level:20
++DEAL:2d::cells: 61 level:21
++DEAL:2d::cells: 64 level:22
++DEAL:2d::cells: 67 level:23
++DEAL:2d::cells: 70 level:24
++DEAL:2d::cells: 73 level:25
++DEAL:2d::cells: 76 level:26
++DEAL:2d::cells: 79 level:27
++DEAL:2d::cells: 82 level:28
++DEAL:2d::cells: 85 level:29
++DEAL:2d::cells: 88 level:30
++DEAL:2d::ExcMessage("Fatal Error: maximum refinement level of p4est reached.")
Simple merge
Simple merge
Simple merge

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.