]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add a test.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 10 Jul 2016 18:52:30 +0000 (13:52 -0500)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 11 Jul 2016 11:56:38 +0000 (06:56 -0500)
tests/base/derivative_forms_02.cc [new file with mode: 0644]
tests/base/derivative_forms_02.output [new file with mode: 0644]

diff --git a/tests/base/derivative_forms_02.cc b/tests/base/derivative_forms_02.cc
new file mode 100644 (file)
index 0000000..edd6cd5
--- /dev/null
@@ -0,0 +1,62 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2000 - 2016 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.
+//
+// ---------------------------------------------------------------------
+
+// exercise DerivativeForm::norm()
+
+
+#include "../tests.h"
+#include <deal.II/base/derivative_form.h>
+
+
+template<int dim, int spacedim>
+void test()
+{
+  DerivativeForm<1,dim,spacedim> dF;
+  double dF_norm_sqr = 0;
+  for (unsigned int i=0; i<spacedim; ++i)
+    for (unsigned int j=0; j<dim; ++j)
+      {
+        dF[i][j] = (i+2*j+1);
+        dF_norm_sqr += (i+2*j+1)*(i+2*j+1);
+      }
+
+  DerivativeForm<2,dim,spacedim> ddF;
+  double ddF_norm_sqr = 0;
+  for (unsigned int i=0; i<spacedim; ++i)
+    for (unsigned int j=0; j<dim; ++j)
+      for (unsigned int k=0; k<dim; ++k)
+        {
+          ddF[i][j][k] = (i+2*j+3*k+1);
+          ddF_norm_sqr += (i+2*j+3*k+1)*(i+2*j+3*k+1);
+        }
+
+  // output the norms of these objects
+  deallog << "||dF||: " << dF.norm() << std::endl;
+  deallog << "||ddF||: " << ddF.norm() << std::endl;
+
+  Assert (std::fabs(dF.norm()-std::sqrt(dF_norm_sqr)) < 1e-12, ExcInternalError());
+  Assert (std::fabs(ddF.norm()-std::sqrt(ddF_norm_sqr)) < 1e-12, ExcInternalError());
+}
+
+int main()
+{
+  initlog();
+  test<1,1>();
+  test<1,2>();
+  test<1,3>();
+  test<2,2>();
+  test<2,3>();
+  test<3,3>();
+}
diff --git a/tests/base/derivative_forms_02.output b/tests/base/derivative_forms_02.output
new file mode 100644 (file)
index 0000000..d217408
--- /dev/null
@@ -0,0 +1,13 @@
+
+DEAL::||dF||: 1.00000
+DEAL::||ddF||: 1.00000
+DEAL::||dF||: 2.23607
+DEAL::||ddF||: 2.23607
+DEAL::||dF||: 3.74166
+DEAL::||ddF||: 3.74166
+DEAL::||dF||: 5.47723
+DEAL::||ddF||: 12.4900
+DEAL::||dF||: 8.00000
+DEAL::||ddF||: 17.0294
+DEAL::||dF||: 13.1909
+DEAL::||ddF||: 39.6863

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.