add test

diff --git a/tests/base/index_set_64bit_01.cc b/tests/base/index_set_64bit_01.cc
new file mode 100644 (file)
index 0000000..1d4bccb
--- /dev/null
+++ b/tests/base/index_set_64bit_01.cc
@@ -0,0 +1,130 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2019 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE.md at
+// the top level directory of deal.II.
+//
+// ---------------------------------------------------------------------
+
+
+// test IndexSet with 64 bit indices
+
+#include <deal.II/base/index_set.h>
+
+#include "../tests.h"
+
+
+void
+print(const IndexSet &s)
+{
+  deallog << "IndexSet: size " << s.size() << "\n"
+          << "n_elements " << s.n_elements() << "\n"
+          << "n_intervals " << s.n_intervals() << "\n"
+          << "largest_range_starting_index " << s.largest_range_starting_index()
+          << "\n"
+          << "is_element(1) " << s.is_element(1) << "\n";
+
+  deallog << "contents:" << std::endl;
+  s.print(deallog.get_file_stream());
+}
+
+
+void
+test()
+{
+  const unsigned long long base   = 10000000000ull; // 10 billion
+  const unsigned long long N      = 5 * base;
+  const unsigned long long large  = 2 * base;
+  const unsigned long long large2 = 3 * base;
+  const unsigned long long large3 = 4 * base;
+
+  {
+    deallog << "* complete:" << std::endl;
+    IndexSet s = complete_index_set(N);
+    Assert(s.size() == N, ExcInternalError());
+    Assert(s.n_elements() == N, ExcInternalError());
+    print(s);
+
+    IndexSet::size_type x = s.index_within_set(large2);
+    deallog << "global: " << large2 << " index_within_set: " << x << std::endl;
+    IndexSet::size_type y = s.nth_index_in_set(large2);
+    deallog << "local: " << large2 << " nth_index_in_set: " << y << std::endl;
+  }
+  {
+    deallog << "* interval far in:" << std::endl;
+    IndexSet s(N);
+    s.add_range(large, large2);
+    Assert(s.size() == N, ExcInternalError());
+    print(s);
+
+    IndexSet::size_type x = s.index_within_set(large + 1);
+    deallog << "global: " << large + 1 << " index_within_set: " << x
+            << std::endl;
+    IndexSet::size_type y = s.nth_index_in_set(3);
+    deallog << "local: " << 3 << " nth_index_in_set: " << y << std::endl;
+  }
+  {
+    deallog << "get view:" << std::endl;
+    IndexSet s(N);
+    s.add_range(large, large2);
+    s = s.get_view(large - 1, large2 - 1);
+    print(s);
+  }
+  {
+    deallog << "* complete minus interval:" << std::endl;
+    IndexSet s = complete_index_set(N);
+    IndexSet other(N);
+    other.add_range(large, large2);
+    s.subtract_set(other);
+
+    Assert(s.size() == N, ExcInternalError());
+    print(s);
+
+    IndexSet::size_type x = s.index_within_set(large2 + 1);
+    deallog << "global: " << large2 + 1 << " index_within_set: " << x
+            << std::endl;
+    IndexSet::size_type y = s.nth_index_in_set(x);
+    deallog << "local: " << x << " nth_index_in_set: " << y << std::endl;
+  }
+
+  {
+    deallog << "* add_indices far in:" << std::endl;
+    IndexSet s(N);
+    IndexSet other(10);
+    other.add_index(1);
+    other.add_range(4, 7);
+    s.add_indices(other, large);
+    Assert(s.size() == N, ExcInternalError());
+    print(s);
+  }
+
+  {
+    deallog << "* two large ranges:" << std::endl;
+    IndexSet s(N);
+    s.add_range(1, large / 2);
+    s.add_range(large, large3);
+    print(s);
+
+    IndexSet::size_type x = s.index_within_set(large2);
+    deallog << "global: " << large2 << " index_within_set: " << x << std::endl;
+    IndexSet::size_type y = s.nth_index_in_set(x);
+    deallog << "local: " << x << " nth_index_in_set: " << y << std::endl;
+  }
+}
+
+
+
+int
+main()
+{
+  initlog();
+
+  test();
+}

diff --git a/tests/base/index_set_64bit_01.with_64bit_indices=on.output b/tests/base/index_set_64bit_01.with_64bit_indices=on.output
new file mode 100644 (file)
index 0000000..afbb6cf
--- /dev/null
+++ b/tests/base/index_set_64bit_01.with_64bit_indices=on.output
@@ -0,0 +1,57 @@
+
+DEAL::* complete:
+DEAL::IndexSet: size 50000000000
+n_elements 50000000000
+n_intervals 1
+largest_range_starting_index 0
+is_element(1) 1
+contents:
+{[0,49999999999]}
+DEAL::global: 30000000000 index_within_set: 30000000000
+DEAL::local: 30000000000 nth_index_in_set: 4230196224
+DEAL::* interval far in:
+DEAL::IndexSet: size 50000000000
+n_elements 10000000000
+n_intervals 1
+largest_range_starting_index 0
+is_element(1) 0
+contents:
+{[20000000000,29999999999]}
+DEAL::global: 20000000001 index_within_set: 1
+DEAL::local: 3 nth_index_in_set: 20000000003
+DEAL::get view:
+DEAL::IndexSet: size 10000000000
+n_elements 9999999999
+n_intervals 1
+largest_range_starting_index 0
+is_element(1) 1
+contents:
+{[1,9999999999]}
+DEAL::* complete minus interval:
+DEAL::IndexSet: size 50000000000
+n_elements 40000000000
+n_intervals 2
+largest_range_starting_index 0
+is_element(1) 1
+contents:
+{[0,19999999999], [30000000000,49999999999]}
+DEAL::global: 30000000001 index_within_set: 20000000001
+DEAL::local: 20000000001 nth_index_in_set: 2820130817
+DEAL::* add_indices far in:
+DEAL::IndexSet: size 50000000000
+n_elements 4
+n_intervals 2
+largest_range_starting_index 1
+is_element(1) 0
+contents:
+{20000000001, [20000000004,20000000006]}
+DEAL::* two large ranges:
+DEAL::IndexSet: size 50000000000
+n_elements 29999999999
+n_intervals 2
+largest_range_starting_index 1410065407
+is_element(1) 1
+contents:
+{[1,9999999999], [20000000000,39999999999]}
+DEAL::global: 30000000000 index_within_set: 19999999999
+DEAL::local: 19999999999 nth_index_in_set: 2820130816