From bef55e11535598acbccc7fe1503fa837c55aebd8 Mon Sep 17 00:00:00 2001 From: Jean-Paul Pelteret Date: Fri, 12 Oct 2018 15:16:35 +0200 Subject: [PATCH] Add tests to check compilation with AD numbers: Tensor class --- tests/ad_common_tests/tensor_functions_01.h | 76 +++++++++++++++++++++ tests/adolc/tensor_functions_01_1.cc | 39 +++++++++++ tests/adolc/tensor_functions_01_1.output | 3 + tests/adolc/tensor_functions_01_2.cc | 39 +++++++++++ tests/adolc/tensor_functions_01_2.output | 3 + tests/sacado/tensor_functions_01_1.cc | 39 +++++++++++ tests/sacado/tensor_functions_01_1.output | 3 + tests/sacado/tensor_functions_01_2.cc | 39 +++++++++++ tests/sacado/tensor_functions_01_2.output | 3 + tests/sacado/tensor_functions_01_3.cc | 39 +++++++++++ tests/sacado/tensor_functions_01_3.output | 3 + tests/sacado/tensor_functions_01_4.cc | 39 +++++++++++ tests/sacado/tensor_functions_01_4.output | 3 + 13 files changed, 328 insertions(+) create mode 100644 tests/ad_common_tests/tensor_functions_01.h create mode 100644 tests/adolc/tensor_functions_01_1.cc create mode 100644 tests/adolc/tensor_functions_01_1.output create mode 100644 tests/adolc/tensor_functions_01_2.cc create mode 100644 tests/adolc/tensor_functions_01_2.output create mode 100644 tests/sacado/tensor_functions_01_1.cc create mode 100644 tests/sacado/tensor_functions_01_1.output create mode 100644 tests/sacado/tensor_functions_01_2.cc create mode 100644 tests/sacado/tensor_functions_01_2.output create mode 100644 tests/sacado/tensor_functions_01_3.cc create mode 100644 tests/sacado/tensor_functions_01_3.output create mode 100644 tests/sacado/tensor_functions_01_4.cc create mode 100644 tests/sacado/tensor_functions_01_4.output diff --git a/tests/ad_common_tests/tensor_functions_01.h b/tests/ad_common_tests/tensor_functions_01.h new file mode 100644 index 0000000000..4b4c819dca --- /dev/null +++ b/tests/ad_common_tests/tensor_functions_01.h @@ -0,0 +1,76 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 - 2018 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. +// +// --------------------------------------------------------------------- + + +// Header file: +// Test to check that the various Tensor functions compile +// with the various auto-differentiable number types + +#include + +#include + +#include + +#include "../tests.h" + +using namespace dealii; +namespace AD = dealii::Differentiation::AD; + +template +void +test_tensor() +{ + typedef + typename AD::NumberTraits::ad_type ADNumberType; + + std::cout << "*** Test Tensor functions, " + << "dim = " << Utilities::to_string(dim) << ", " + << "Type code: " << static_cast(ad_type_code) << std::endl; + + const ADNumberType a = 1.0; + const Tensor<1, dim, ADNumberType> v1; + const Tensor<1, dim, ADNumberType> v2; + const Tensor<2, dim, ADNumberType> A( + unit_symmetric_tensor()); + const Tensor<2, dim, ADNumberType> B( + unit_symmetric_tensor()); + + const Tensor<2, dim, ADNumberType> C1 = A + B; + const Tensor<2, dim, ADNumberType> C2 = A - B; + const Tensor<2, dim, ADNumberType> C3 = A * B; + const Tensor<2, dim, ADNumberType> C4 = a * A; + const Tensor<2, dim, ADNumberType> C5 = A * a; + const Tensor<2, dim, ADNumberType> C6 = A / a; + + const ADNumberType det_A = determinant(A); + const ADNumberType tr_A = trace(A); + const Tensor<2, dim, ADNumberType> A_inv = invert(A); + const Tensor<2, dim, ADNumberType> A_T = transpose(A); + const Tensor<2, dim, ADNumberType> A_adj = adjugate(A); + const Tensor<2, dim, ADNumberType> A_cof = cofactor(A); + const ADNumberType A_l1_norm = l1_norm(A); + const ADNumberType A_linf_norm = linfty_norm(A); + + const ADNumberType A_ddot_B = double_contract(A, B); + const Tensor<2, dim, ADNumberType> A_dot_B = contract<1, 0>(A, B); + const ADNumberType sp_A_B = scalar_product(A, B); + const Tensor<4, dim, ADNumberType> op_A_B = outer_product(A, B); + + if (dim == 2) + const Tensor<1, dim, ADNumberType> v3 = cross_product_2d(v1); + if (dim == 3) + const Tensor<1, dim, ADNumberType> v3 = cross_product_3d(v1, v2); +} diff --git a/tests/adolc/tensor_functions_01_1.cc b/tests/adolc/tensor_functions_01_1.cc new file mode 100644 index 0000000000..57925c53c6 --- /dev/null +++ b/tests/adolc/tensor_functions_01_1.cc @@ -0,0 +1,39 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 - 2018 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. +// +// --------------------------------------------------------------------- + + +// Test to check that the various Tensor functions compile +// with the various auto-differentiable number types +// +// AD number type: ADOL-C taped + +#include "../ad_common_tests/tensor_functions_01.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_tensor<2, double, AD::NumberTypes::adolc_taped>(); + test_tensor<3, double, AD::NumberTypes::adolc_taped>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/adolc/tensor_functions_01_1.output b/tests/adolc/tensor_functions_01_1.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/adolc/tensor_functions_01_1.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/adolc/tensor_functions_01_2.cc b/tests/adolc/tensor_functions_01_2.cc new file mode 100644 index 0000000000..274fb9582b --- /dev/null +++ b/tests/adolc/tensor_functions_01_2.cc @@ -0,0 +1,39 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 - 2018 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. +// +// --------------------------------------------------------------------- + + +// Test to check that the various Tensor functions compile +// with the various auto-differentiable number types +// +// AD number type: ADOL-C tapeless + +#include "../ad_common_tests/tensor_functions_01.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_tensor<2, double, AD::NumberTypes::adolc_tapeless>(); + test_tensor<3, double, AD::NumberTypes::adolc_tapeless>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/adolc/tensor_functions_01_2.output b/tests/adolc/tensor_functions_01_2.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/adolc/tensor_functions_01_2.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/sacado/tensor_functions_01_1.cc b/tests/sacado/tensor_functions_01_1.cc new file mode 100644 index 0000000000..610aa78d30 --- /dev/null +++ b/tests/sacado/tensor_functions_01_1.cc @@ -0,0 +1,39 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 - 2018 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. +// +// --------------------------------------------------------------------- + + +// Test to check that the various Tensor functions compile +// with the various auto-differentiable number types +// +// AD number type: Sacado DFad + +#include "../ad_common_tests/tensor_functions_01.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_tensor<2, double, AD::NumberTypes::sacado_dfad>(); + test_tensor<3, double, AD::NumberTypes::sacado_dfad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/sacado/tensor_functions_01_1.output b/tests/sacado/tensor_functions_01_1.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/sacado/tensor_functions_01_1.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/sacado/tensor_functions_01_2.cc b/tests/sacado/tensor_functions_01_2.cc new file mode 100644 index 0000000000..caff61dbf6 --- /dev/null +++ b/tests/sacado/tensor_functions_01_2.cc @@ -0,0 +1,39 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 - 2018 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. +// +// --------------------------------------------------------------------- + + +// Test to check that the various Tensor functions compile +// with the various auto-differentiable number types +// +// AD number type: Sacado DFad-DFad + +#include "../ad_common_tests/tensor_functions_01.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_tensor<2, double, AD::NumberTypes::sacado_dfad_dfad>(); + test_tensor<3, double, AD::NumberTypes::sacado_dfad_dfad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/sacado/tensor_functions_01_2.output b/tests/sacado/tensor_functions_01_2.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/sacado/tensor_functions_01_2.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/sacado/tensor_functions_01_3.cc b/tests/sacado/tensor_functions_01_3.cc new file mode 100644 index 0000000000..2389c40b57 --- /dev/null +++ b/tests/sacado/tensor_functions_01_3.cc @@ -0,0 +1,39 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 - 2018 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. +// +// --------------------------------------------------------------------- + + +// Test to check that the various Tensor functions compile +// with the various auto-differentiable number types +// +// AD number type: Sacado Rad + +#include "../ad_common_tests/tensor_functions_01.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_tensor<2, double, AD::NumberTypes::sacado_rad>(); + test_tensor<3, double, AD::NumberTypes::sacado_rad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/sacado/tensor_functions_01_3.output b/tests/sacado/tensor_functions_01_3.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/sacado/tensor_functions_01_3.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/sacado/tensor_functions_01_4.cc b/tests/sacado/tensor_functions_01_4.cc new file mode 100644 index 0000000000..b57f559ec5 --- /dev/null +++ b/tests/sacado/tensor_functions_01_4.cc @@ -0,0 +1,39 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2016 - 2018 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. +// +// --------------------------------------------------------------------- + + +// Test to check that the various Tensor functions compile +// with the various auto-differentiable number types +// +// AD number type: Sacado Rad-DFad + +#include "../ad_common_tests/tensor_functions_01.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_tensor<2, double, AD::NumberTypes::sacado_rad_dfad>(); + test_tensor<3, double, AD::NumberTypes::sacado_rad_dfad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/sacado/tensor_functions_01_4.output b/tests/sacado/tensor_functions_01_4.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/sacado/tensor_functions_01_4.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK -- 2.39.5