]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add new no_normal_flux_constraints test that used to throw a FPE.
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Sat, 2 Dec 2017 10:04:47 +0000 (11:04 +0100)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Tue, 5 Dec 2017 10:00:02 +0000 (11:00 +0100)
tests/numerics/no_flux_15.cc [new file with mode: 0644]
tests/numerics/no_flux_15.output [new file with mode: 0644]

diff --git a/tests/numerics/no_flux_15.cc b/tests/numerics/no_flux_15.cc
new file mode 100644 (file)
index 0000000..263b8d2
--- /dev/null
@@ -0,0 +1,70 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2017 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 at
+// the top level of the deal.II distribution.
+//
+// ---------------------------------------------------------------------
+
+
+
+// Similar to no_flux_14 for a higher order element including a high-order
+// mapping. The purpose here is to test that the result of
+// Manifold::normal_vector function inside of
+// VectorTools::compute_no_normal_flux_constraints is accurate for higher
+// order mappings.
+
+
+#include "../tests.h"
+
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/lac/constraint_matrix.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/mapping_q.h>
+#include <deal.II/numerics/vector_tools.h>
+#include <deal.II/grid/manifold_lib.h>
+
+template <int dim>
+void
+check ()
+{
+  SphericalManifold<dim> spherical;
+  Triangulation<dim> tria;
+  GridGenerator::hyper_shell(tria, Point<dim>(), 0.5, 1., 96, true);
+  tria.set_all_manifold_ids(0);
+  tria.set_manifold(0, spherical);
+
+  ConstraintMatrix cm;
+  MappingQ<dim> mapping(4);
+
+  FESystem<dim> fe(FE_Q<dim>(2),dim);
+  DoFHandler<dim> dofh(tria);
+
+  dofh.distribute_dofs (fe);
+
+  std::set<types::boundary_id> no_normal_flux_boundaries;
+  no_normal_flux_boundaries.insert (0);
+  no_normal_flux_boundaries.insert (1);
+  VectorTools::compute_no_normal_flux_constraints (dofh, 0, no_normal_flux_boundaries, cm, mapping);
+
+  cm.print (deallog.get_file_stream ());
+}
+
+
+
+int main ()
+{
+  initlog();
+  deallog.get_file_stream().precision(8);
+
+  check<3> ();
+}
diff --git a/tests/numerics/no_flux_15.output b/tests/numerics/no_flux_15.output
new file mode 100644 (file)
index 0000000..c0088ff
--- /dev/null
@@ -0,0 +1,693 @@
+
+    1 = 0
+    4 3:  -0.36602540
+    4 5:  0.36602540
+    7 6:  -0.36602540
+    7 8:  -0.36602540
+    10 9:  -1.0000000
+    25 24:  -0.17216360
+    25 26:  -0.17216360
+    28 27:  -0.64705172
+    28 29:  0.20377478
+    31 30:  -0.17216360
+    31 32:  0.17216360
+    34 33:  -0.64705172
+    34 35:  -0.20377478
+    73 72:  -0.39037473
+    83 81:  1.0000000
+    83 82:  1.0000000
+    84 85:  -0.36602540
+    84 86:  0.36602540
+    93 94:  -0.61579979
+    93 95:  0.61579979
+    97 96:  -0.61579979
+    97 98:  0.61579979
+    99 100:  -0.64705172
+    99 101:  0.20377478
+    127 126:  -1.0000000
+    127 128:  0.49102848
+    137 135:  -1.0000000
+    137 136:  -1.0000000
+    138 139:  -0.36602540
+    138 140:  -0.36602540
+    148 147:  -0.61579979
+    148 149:  -0.61579979
+    150 151:  -0.64705172
+    150 152:  -0.20377478
+    153 154:  -0.61579979
+    153 155:  -0.61579979
+    181 180:  -1.0000000
+    181 182:  -0.49102848
+    189 = 0
+    195 196:  -0.17216360
+    195 197:  0.17216360
+    198 199:  -0.17216360
+    198 200:  -0.17216360
+    216 217:  -0.39037473
+    226 = 0
+    229 228:  -0.36602540
+    229 230:  0.36602540
+    232 231:  -0.36602540
+    232 233:  -0.36602540
+    235 234:  -1.0000000
+    238 237:  -0.17216360
+    238 239:  -0.17216360
+    241 240:  -0.64705172
+    241 242:  0.20377478
+    244 243:  -0.17216360
+    244 245:  0.17216360
+    247 246:  -0.64705172
+    247 248:  -0.20377478
+    274 273:  -0.39037473
+    281 279:  1.0000000
+    281 280:  1.0000000
+    282 283:  -0.36602540
+    282 284:  0.36602540
+    285 286:  -0.61579979
+    285 287:  0.61579979
+    289 288:  -0.61579979
+    289 290:  0.61579979
+    291 292:  -0.64705172
+    291 293:  0.20377478
+    310 309:  -1.0000000
+    310 311:  0.49102848
+    317 315:  -1.0000000
+    317 316:  -1.0000000
+    318 319:  -0.36602540
+    318 320:  -0.36602540
+    322 321:  -0.61579979
+    322 323:  -0.61579979
+    324 325:  -0.64705172
+    324 326:  -0.20377478
+    327 328:  -0.61579979
+    327 329:  -0.61579979
+    346 345:  -1.0000000
+    346 347:  -0.49102848
+    351 = 0
+    354 355:  -0.17216360
+    354 356:  0.17216360
+    357 358:  -0.17216360
+    357 359:  -0.17216360
+    369 370:  -0.39037473
+    377 375:  0.36602540
+    377 376:  0.36602540
+    380 378:  1.0000000
+    389 387:  0.64705172
+    389 388:  0.20377478
+    392 390:  0.61579979
+    392 391:  0.61579979
+    393 394:  -0.20377478
+    393 395:  0.64705172
+    422 420:  1.0000000
+    422 421:  0.49102848
+    431 = 0
+    434 432:  0.36602540
+    434 433:  -0.36602540
+    443 441:  0.17216360
+    443 442:  -0.17216360
+    446 444:  0.17216360
+    446 445:  0.17216360
+    449 447:  0.64705172
+    449 448:  -0.20377478
+    476 474:  0.39037473
+    483 484:  0.36602540
+    483 485:  0.36602540
+    489 490:  0.20377478
+    489 491:  0.64705172
+    492 493:  0.17216360
+    492 494:  0.17216360
+    510 512:  0.39037473
+    521 519:  1.0000000
+    521 520:  -1.0000000
+    527 525:  0.61579979
+    527 526:  -0.61579979
+    528 529:  0.61579979
+    528 530:  0.61579979
+    548 546:  1.0000000
+    548 547:  -0.49102848
+    557 555:  0.36602540
+    557 556:  0.36602540
+    560 558:  1.0000000
+    563 561:  0.64705172
+    563 562:  0.20377478
+    566 564:  0.61579979
+    566 565:  0.61579979
+    567 568:  -0.20377478
+    567 569:  0.64705172
+    587 585:  1.0000000
+    587 586:  0.49102848
+    593 = 0
+    596 594:  0.36602540
+    596 595:  -0.36602540
+    599 597:  0.17216360
+    599 598:  -0.17216360
+    602 600:  0.17216360
+    602 601:  0.17216360
+    605 603:  0.64705172
+    605 604:  -0.20377478
+    623 621:  0.39037473
+    627 628:  0.36602540
+    627 629:  0.36602540
+    630 631:  0.20377478
+    630 632:  0.64705172
+    633 634:  0.17216360
+    633 635:  0.17216360
+    645 647:  0.39037473
+    653 651:  1.0000000
+    653 652:  -1.0000000
+    656 654:  0.61579979
+    656 655:  -0.61579979
+    657 658:  0.61579979
+    657 659:  0.61579979
+    671 669:  1.0000000
+    671 670:  -0.49102848
+    676 675:  0.36602540
+    676 677:  0.36602540
+    680 679:  1.0000000
+    688 687:  0.20377478
+    688 689:  0.64705172
+    691 690:  0.17216360
+    691 692:  0.17216360
+    694 693:  -0.20377478
+    694 695:  0.64705172
+    721 722:  0.39037473
+    731 729:  -1.0000000
+    731 730:  1.0000000
+    734 732:  -0.36602540
+    734 733:  0.36602540
+    743 741:  -0.61579979
+    743 742:  0.61579979
+    745 744:  0.61579979
+    745 746:  0.61579979
+    749 747:  -0.20377478
+    749 748:  0.64705172
+    776 774:  -0.49102848
+    776 775:  1.0000000
+    785 783:  0.20377478
+    785 784:  0.64705172
+    794 792:  0.49102848
+    794 793:  1.0000000
+    803 801:  -0.17216360
+    803 802:  0.17216360
+    812 811:  0.39037473
+    820 819:  0.36602540
+    820 821:  0.36602540
+    824 823:  1.0000000
+    826 825:  0.20377478
+    826 827:  0.64705172
+    829 828:  0.17216360
+    829 830:  0.17216360
+    832 831:  -0.20377478
+    832 833:  0.64705172
+    850 851:  0.39037473
+    857 855:  -1.0000000
+    857 856:  1.0000000
+    860 858:  -0.36602540
+    860 859:  0.36602540
+    863 861:  -0.61579979
+    863 862:  0.61579979
+    865 864:  0.61579979
+    865 866:  0.61579979
+    869 867:  -0.20377478
+    869 868:  0.64705172
+    887 885:  -0.49102848
+    887 886:  1.0000000
+    893 891:  0.20377478
+    893 892:  0.64705172
+    899 897:  0.49102848
+    899 898:  1.0000000
+    905 903:  -0.17216360
+    905 904:  0.17216360
+    911 910:  0.39037473
+    917 915:  1.0000000
+    917 916:  -1.0000000
+    918 919:  0.36602540
+    918 920:  0.36602540
+    922 921:  0.36602540
+    922 923:  -0.36602540
+    925 924:  1.0000000
+    940 939:  0.61579979
+    940 941:  -0.61579979
+    942 943:  0.64705172
+    942 944:  0.20377478
+    945 946:  0.61579979
+    945 947:  0.61579979
+    949 948:  0.64705172
+    949 950:  -0.20377478
+    988 987:  1.0000000
+    988 989:  -0.49102848
+    996 = 0
+    999 1000:  0.36602540
+    999 1001:  -0.36602540
+    1008 1009:  0.17216360
+    1008 1010:  -0.17216360
+    1011 1012:  0.17216360
+    1011 1013:  0.17216360
+    1014 1015:  0.64705172
+    1014 1016:  -0.20377478
+    1041 1042:  0.39037473
+    1051 1050:  0.17216360
+    1051 1052:  -0.17216360
+    1054 1053:  0.64705172
+    1054 1055:  0.20377478
+    1069 1068:  0.39037473
+    1077 1078:  0.61579979
+    1077 1079:  -0.61579979
+    1087 1086:  1.0000000
+    1087 1088:  0.49102848
+    1097 1095:  1.0000000
+    1097 1096:  -1.0000000
+    1098 1099:  0.36602540
+    1098 1100:  0.36602540
+    1102 1101:  0.36602540
+    1102 1103:  -0.36602540
+    1105 1104:  1.0000000
+    1108 1107:  0.61579979
+    1108 1109:  -0.61579979
+    1110 1111:  0.64705172
+    1110 1112:  0.20377478
+    1113 1114:  0.61579979
+    1113 1115:  0.61579979
+    1117 1116:  0.64705172
+    1117 1118:  -0.20377478
+    1144 1143:  1.0000000
+    1144 1145:  -0.49102848
+    1149 = 0
+    1152 1153:  0.36602540
+    1152 1154:  -0.36602540
+    1155 1156:  0.17216360
+    1155 1157:  -0.17216360
+    1158 1159:  0.17216360
+    1158 1160:  0.17216360
+    1161 1162:  0.64705172
+    1161 1163:  -0.20377478
+    1179 1180:  0.39037473
+    1186 1185:  0.17216360
+    1186 1187:  -0.17216360
+    1189 1188:  0.64705172
+    1189 1190:  0.20377478
+    1198 1197:  0.39037473
+    1203 1204:  0.61579979
+    1203 1205:  -0.61579979
+    1210 1209:  1.0000000
+    1210 1211:  0.49102848
+    1215 1216:  -0.36602540
+    1215 1217:  -0.36602540
+    1220 1218:  -1.0000000
+    1227 1228:  -0.20377478
+    1227 1229:  -0.64705172
+    1230 1231:  -0.17216360
+    1230 1232:  -0.17216360
+    1233 1234:  0.20377478
+    1233 1235:  -0.64705172
+    1260 1262:  -0.39037473
+    1271 1269:  -1.0000000
+    1271 1270:  -1.0000000
+    1274 1272:  -0.36602540
+    1274 1273:  -0.36602540
+    1283 1281:  -0.61579979
+    1283 1282:  -0.61579979
+    1284 1285:  -0.61579979
+    1284 1286:  -0.61579979
+    1289 1287:  -0.64705172
+    1289 1288:  -0.20377478
+    1316 1314:  -1.0000000
+    1316 1315:  -0.49102848
+    1325 1323:  -0.64705172
+    1325 1324:  0.20377478
+    1334 1332:  -1.0000000
+    1334 1333:  0.49102848
+    1343 1341:  -0.17216360
+    1343 1342:  -0.17216360
+    1352 1350:  -0.39037473
+    1359 1360:  -0.36602540
+    1359 1361:  -0.36602540
+    1364 1362:  -1.0000000
+    1365 1366:  -0.20377478
+    1365 1367:  -0.64705172
+    1368 1369:  -0.17216360
+    1368 1370:  -0.17216360
+    1371 1372:  0.20377478
+    1371 1373:  -0.64705172
+    1389 1391:  -0.39037473
+    1397 1395:  -1.0000000
+    1397 1396:  -1.0000000
+    1400 1398:  -0.36602540
+    1400 1399:  -0.36602540
+    1403 1401:  -0.61579979
+    1403 1402:  -0.61579979
+    1404 1405:  -0.61579979
+    1404 1406:  -0.61579979
+    1409 1407:  -0.64705172
+    1409 1408:  -0.20377478
+    1427 1425:  -1.0000000
+    1427 1426:  -0.49102848
+    1433 1431:  -0.64705172
+    1433 1432:  0.20377478
+    1439 1437:  -1.0000000
+    1439 1438:  0.49102848
+    1445 1443:  -0.17216360
+    1445 1444:  -0.17216360
+    1451 1449:  -0.39037473
+    1456 1455:  -0.36602540
+    1456 1457:  -0.36602540
+    1460 1459:  -1.0000000
+    1468 1467:  -0.20377478
+    1468 1469:  -0.64705172
+    1471 1470:  -0.61579979
+    1471 1472:  -0.61579979
+    1475 1473:  -0.20377478
+    1475 1474:  -0.64705172
+    1502 1500:  -0.49102848
+    1502 1501:  -1.0000000
+    1510 = 0
+    1513 1512:  0.36602540
+    1513 1514:  -0.36602540
+    1522 1521:  0.17216360
+    1522 1523:  -0.17216360
+    1525 1524:  -0.17216360
+    1525 1526:  -0.17216360
+    1528 1527:  0.20377478
+    1528 1529:  -0.64705172
+    1555 1556:  -0.39037473
+    1565 1563:  0.20377478
+    1565 1564:  -0.64705172
+    1574 1573:  -0.39037473
+    1582 1581:  0.61579979
+    1582 1583:  -0.61579979
+    1592 1590:  0.49102848
+    1592 1591:  -1.0000000
+    1600 1599:  -0.36602540
+    1600 1601:  -0.36602540
+    1604 1603:  -1.0000000
+    1606 1605:  -0.20377478
+    1606 1607:  -0.64705172
+    1609 1608:  -0.61579979
+    1609 1610:  -0.61579979
+    1613 1611:  -0.20377478
+    1613 1612:  -0.64705172
+    1631 1629:  -0.49102848
+    1631 1630:  -1.0000000
+    1636 = 0
+    1639 1638:  0.36602540
+    1639 1640:  -0.36602540
+    1642 1641:  0.17216360
+    1642 1643:  -0.17216360
+    1645 1644:  -0.17216360
+    1645 1646:  -0.17216360
+    1648 1647:  0.20377478
+    1648 1649:  -0.64705172
+    1666 1667:  -0.39037473
+    1673 1671:  0.20377478
+    1673 1672:  -0.64705172
+    1679 1678:  -0.39037473
+    1684 1683:  0.61579979
+    1684 1685:  -0.61579979
+    1691 1689:  0.49102848
+    1691 1690:  -1.0000000
+    1695 1696:  -0.36602540
+    1695 1697:  0.36602540
+    1699 1698:  -1.0000000
+    1707 1708:  -0.64705172
+    1707 1709:  0.20377478
+    1710 1711:  -0.17216360
+    1710 1712:  0.17216360
+    1713 1714:  -0.64705172
+    1713 1715:  -0.20377478
+    1740 1741:  -0.39037473
+    1751 1749:  1.0000000
+    1751 1750:  1.0000000
+    1753 1752:  -0.36602540
+    1753 1754:  0.36602540
+    1762 1761:  -0.61579979
+    1762 1763:  0.61579979
+    1764 1765:  -0.61579979
+    1764 1766:  0.61579979
+    1768 1767:  -0.64705172
+    1768 1769:  0.20377478
+    1795 1794:  -1.0000000
+    1795 1796:  0.49102848
+    1804 1803:  -0.64705172
+    1804 1805:  -0.20377478
+    1813 1812:  -1.0000000
+    1813 1814:  -0.49102848
+    1822 1821:  -0.17216360
+    1822 1823:  0.17216360
+    1831 1830:  -0.39037473
+    1839 1840:  -0.36602540
+    1839 1841:  0.36602540
+    1843 1842:  -1.0000000
+    1845 1846:  -0.64705172
+    1845 1847:  0.20377478
+    1848 1849:  -0.17216360
+    1848 1850:  0.17216360
+    1851 1852:  -0.64705172
+    1851 1853:  -0.20377478
+    1869 1870:  -0.39037473
+    1877 1875:  1.0000000
+    1877 1876:  1.0000000
+    1879 1878:  -0.36602540
+    1879 1880:  0.36602540
+    1882 1881:  -0.61579979
+    1882 1883:  0.61579979
+    1884 1885:  -0.61579979
+    1884 1886:  0.61579979
+    1888 1887:  -0.64705172
+    1888 1889:  0.20377478
+    1906 1905:  -1.0000000
+    1906 1907:  0.49102848
+    1912 1911:  -0.64705172
+    1912 1913:  -0.20377478
+    1918 1917:  -1.0000000
+    1918 1919:  -0.49102848
+    1924 1923:  -0.17216360
+    1924 1925:  0.17216360
+    1930 1929:  -0.39037473
+    1937 1935:  0.36602540
+    1937 1936:  -0.36602540
+    1940 1938:  1.0000000
+    1949 1947:  0.64705172
+    1949 1948:  -0.20377478
+    1952 1950:  0.61579979
+    1952 1951:  -0.61579979
+    1953 1954:  0.20377478
+    1953 1955:  0.64705172
+    1982 1980:  1.0000000
+    1982 1981:  -0.49102848
+    1991 = 0
+    1994 1992:  0.36602540
+    1994 1993:  0.36602540
+    2003 2001:  0.17216360
+    2003 2002:  0.17216360
+    2006 2004:  0.17216360
+    2006 2005:  -0.17216360
+    2009 2007:  0.64705172
+    2009 2008:  0.20377478
+    2036 2034:  0.39037473
+    2043 2044:  -0.20377478
+    2043 2045:  0.64705172
+    2052 2054:  0.39037473
+    2063 2061:  0.61579979
+    2063 2062:  0.61579979
+    2072 2070:  1.0000000
+    2072 2071:  0.49102848
+    2081 2079:  0.36602540
+    2081 2080:  -0.36602540
+    2084 2082:  1.0000000
+    2087 2085:  0.64705172
+    2087 2086:  -0.20377478
+    2090 2088:  0.61579979
+    2090 2089:  -0.61579979
+    2091 2092:  0.20377478
+    2091 2093:  0.64705172
+    2111 2109:  1.0000000
+    2111 2110:  -0.49102848
+    2117 = 0
+    2120 2118:  0.36602540
+    2120 2119:  0.36602540
+    2123 2121:  0.17216360
+    2123 2122:  0.17216360
+    2126 2124:  0.17216360
+    2126 2125:  -0.17216360
+    2129 2127:  0.64705172
+    2129 2128:  0.20377478
+    2147 2145:  0.39037473
+    2151 2152:  -0.20377478
+    2151 2153:  0.64705172
+    2157 2159:  0.39037473
+    2165 2163:  0.61579979
+    2165 2164:  0.61579979
+    2171 2169:  1.0000000
+    2171 2170:  0.49102848
+    2177 2175:  -0.36602540
+    2177 2176:  0.36602540
+    2180 2179:  1.0000000
+    2189 2187:  -0.20377478
+    2189 2188:  0.64705172
+    2192 2190:  -0.17216360
+    2192 2191:  0.17216360
+    2195 2193:  0.20377478
+    2195 2194:  0.64705172
+    2222 2221:  0.39037473
+    2231 2229:  -1.0000000
+    2231 2230:  1.0000000
+    2233 2232:  0.36602540
+    2233 2234:  0.36602540
+    2242 2241:  0.61579979
+    2242 2243:  0.61579979
+    2246 2244:  -0.61579979
+    2246 2245:  0.61579979
+    2248 2247:  0.20377478
+    2248 2249:  0.64705172
+    2276 2274:  -0.49102848
+    2276 2275:  1.0000000
+    2284 2283:  -0.20377478
+    2284 2285:  0.64705172
+    2294 2292:  0.49102848
+    2294 2293:  1.0000000
+    2302 2301:  0.17216360
+    2302 2303:  0.17216360
+    2311 2312:  0.39037473
+    2321 2319:  -0.36602540
+    2321 2320:  0.36602540
+    2324 2323:  1.0000000
+    2327 2325:  -0.20377478
+    2327 2326:  0.64705172
+    2330 2328:  -0.17216360
+    2330 2329:  0.17216360
+    2333 2331:  0.20377478
+    2333 2332:  0.64705172
+    2351 2350:  0.39037473
+    2357 2355:  -1.0000000
+    2357 2356:  1.0000000
+    2359 2358:  0.36602540
+    2359 2360:  0.36602540
+    2362 2361:  0.61579979
+    2362 2363:  0.61579979
+    2366 2364:  -0.61579979
+    2366 2365:  0.61579979
+    2368 2367:  0.20377478
+    2368 2369:  0.64705172
+    2387 2385:  -0.49102848
+    2387 2386:  1.0000000
+    2392 2391:  -0.20377478
+    2392 2393:  0.64705172
+    2399 2397:  0.49102848
+    2399 2398:  1.0000000
+    2404 2403:  0.17216360
+    2404 2405:  0.17216360
+    2410 2411:  0.39037473
+    2415 2416:  0.36602540
+    2415 2417:  -0.36602540
+    2419 2418:  1.0000000
+    2427 2428:  0.64705172
+    2427 2429:  -0.20377478
+    2430 2431:  0.61579979
+    2430 2432:  -0.61579979
+    2434 2433:  0.64705172
+    2434 2435:  0.20377478
+    2461 2460:  1.0000000
+    2461 2462:  0.49102848
+    2469 2470:  0.17216360
+    2469 2471:  -0.17216360
+    2472 2473:  0.64705172
+    2472 2474:  0.20377478
+    2487 2488:  0.39037473
+    2497 2496:  0.64705172
+    2497 2498:  -0.20377478
+    2506 2505:  0.39037473
+    2515 2514:  1.0000000
+    2515 2516:  -0.49102848
+    2523 2524:  0.36602540
+    2523 2525:  -0.36602540
+    2527 2526:  1.0000000
+    2529 2530:  0.64705172
+    2529 2531:  -0.20377478
+    2532 2533:  0.61579979
+    2532 2534:  -0.61579979
+    2536 2535:  0.64705172
+    2536 2537:  0.20377478
+    2554 2553:  1.0000000
+    2554 2555:  0.49102848
+    2559 2560:  0.17216360
+    2559 2561:  -0.17216360
+    2562 2563:  0.64705172
+    2562 2564:  0.20377478
+    2571 2572:  0.39037473
+    2578 2577:  0.64705172
+    2578 2579:  -0.20377478
+    2584 2583:  0.39037473
+    2590 2589:  1.0000000
+    2590 2591:  -0.49102848
+    2597 2595:  -0.36602540
+    2597 2596:  -0.36602540
+    2600 2598:  -1.0000000
+    2609 2607:  -0.64705172
+    2609 2608:  -0.20377478
+    2612 2610:  -0.17216360
+    2612 2611:  -0.17216360
+    2615 2613:  -0.64705172
+    2615 2614:  0.20377478
+    2642 2640:  -0.39037473
+    2651 2649:  -0.61579979
+    2651 2650:  -0.61579979
+    2652 2653:  -0.20377478
+    2652 2654:  -0.64705172
+    2669 2667:  -1.0000000
+    2669 2668:  -0.49102848
+    2676 2677:  0.20377478
+    2676 2678:  -0.64705172
+    2687 2685:  -1.0000000
+    2687 2686:  0.49102848
+    2694 2696:  -0.39037473
+    2705 2703:  -0.36602540
+    2705 2704:  -0.36602540
+    2708 2706:  -1.0000000
+    2711 2709:  -0.64705172
+    2711 2710:  -0.20377478
+    2714 2712:  -0.17216360
+    2714 2713:  -0.17216360
+    2717 2715:  -0.64705172
+    2717 2716:  0.20377478
+    2735 2733:  -0.39037473
+    2741 2739:  -0.61579979
+    2741 2740:  -0.61579979
+    2742 2743:  -0.20377478
+    2742 2744:  -0.64705172
+    2753 2751:  -1.0000000
+    2753 2752:  -0.49102848
+    2757 2758:  0.20377478
+    2757 2759:  -0.64705172
+    2765 2763:  -1.0000000
+    2765 2764:  0.49102848
+    2769 2771:  -0.39037473
+    2777 2776:  -1.0000000
+    2782 2781:  0.20377478
+    2782 2783:  -0.64705172
+    2786 2784:  0.20377478
+    2786 2785:  -0.64705172
+    2804 2802:  0.49102848
+    2804 2803:  -1.0000000
+    2812 2811:  -0.20377478
+    2812 2813:  -0.64705172
+    2821 2822:  -0.39037473
+    2831 2829:  -0.20377478
+    2831 2830:  -0.64705172
+    2840 2839:  -0.39037473
+    2849 2847:  -0.49102848
+    2849 2848:  -1.0000000
+    2858 2857:  -1.0000000
+    2860 2859:  0.20377478
+    2860 2861:  -0.64705172
+    2864 2862:  0.20377478
+    2864 2863:  -0.64705172
+    2876 2874:  0.49102848
+    2876 2875:  -1.0000000
+    2881 2880:  -0.20377478
+    2881 2882:  -0.64705172
+    2887 2888:  -0.39037473
+    2894 2892:  -0.20377478
+    2894 2893:  -0.64705172
+    2900 2899:  -0.39037473
+    2906 2904:  -0.49102848
+    2906 2905:  -1.0000000

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.