From 09566d9dcc7fab7bb7231f0ed5f15900b5460d5e Mon Sep 17 00:00:00 2001 From: Jean-Paul Pelteret Date: Fri, 12 Oct 2018 15:18:44 +0200 Subject: [PATCH] Add tests to check compilation with AD numbers: Physics NS (part 2) Physics::Elasticity::StandardTensors class --- tests/ad_common_tests/physics_functions_02.h | 61 ++++++++++++++++++++ tests/adolc/physics_functions_02_1.cc | 40 +++++++++++++ tests/adolc/physics_functions_02_1.output | 3 + tests/adolc/physics_functions_02_2.cc | 40 +++++++++++++ tests/adolc/physics_functions_02_2.output | 3 + tests/sacado/physics_functions_02_1.cc | 48 +++++++++++++++ tests/sacado/physics_functions_02_1.output | 4 ++ tests/sacado/physics_functions_02_2.cc | 48 +++++++++++++++ tests/sacado/physics_functions_02_2.output | 4 ++ tests/sacado/physics_functions_02_3.cc | 48 +++++++++++++++ tests/sacado/physics_functions_02_3.output | 4 ++ tests/sacado/physics_functions_02_4.cc | 48 +++++++++++++++ tests/sacado/physics_functions_02_4.output | 4 ++ 13 files changed, 355 insertions(+) create mode 100644 tests/ad_common_tests/physics_functions_02.h create mode 100644 tests/adolc/physics_functions_02_1.cc create mode 100644 tests/adolc/physics_functions_02_1.output create mode 100644 tests/adolc/physics_functions_02_2.cc create mode 100644 tests/adolc/physics_functions_02_2.output create mode 100644 tests/sacado/physics_functions_02_1.cc create mode 100644 tests/sacado/physics_functions_02_1.output create mode 100644 tests/sacado/physics_functions_02_2.cc create mode 100644 tests/sacado/physics_functions_02_2.output create mode 100644 tests/sacado/physics_functions_02_3.cc create mode 100644 tests/sacado/physics_functions_02_3.output create mode 100644 tests/sacado/physics_functions_02_4.cc create mode 100644 tests/sacado/physics_functions_02_4.output diff --git a/tests/ad_common_tests/physics_functions_02.h b/tests/ad_common_tests/physics_functions_02.h new file mode 100644 index 0000000000..2ec3bdcdfe --- /dev/null +++ b/tests/ad_common_tests/physics_functions_02.h @@ -0,0 +1,61 @@ +// --------------------------------------------------------------------- +// +// 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 functions in the Physics namespace compile +// with the various auto-differentiable number types: +// StandardTensors + +#include +#include + +#include + +#include +#include + +#include + +#include "../tests.h" + +using namespace dealii; +namespace AD = dealii::Differentiation::AD; + +template +void +test_physics() +{ + typedef + typename AD::NumberTraits::ad_type ADNumberType; + + std::cout << "*** Test physics functions: Standard tensors, " + << "dim = " << Utilities::to_string(dim) << ", " + << "Type code: " << static_cast(ad_type_code) << std::endl; + + Tensor<2, dim, ADNumberType> grad_u; + const Tensor<2, dim, ADNumberType> F = + Physics::Elasticity::Kinematics::F(grad_u); + + const SymmetricTensor<2, dim, ADNumberType> ddet_F_dC = + Physics::Elasticity::StandardTensors::ddet_F_dC(F); + const SymmetricTensor<4, dim, ADNumberType> dC_inv_dC = + Physics::Elasticity::StandardTensors::dC_inv_dC(F); + + const SymmetricTensor<4, dim, ADNumberType> Dev_P = + Physics::Elasticity::StandardTensors::Dev_P(F); + const SymmetricTensor<4, dim, ADNumberType> Dev_P_T = + Physics::Elasticity::StandardTensors::Dev_P_T(F); +} diff --git a/tests/adolc/physics_functions_02_1.cc b/tests/adolc/physics_functions_02_1.cc new file mode 100644 index 0000000000..5d6bafb256 --- /dev/null +++ b/tests/adolc/physics_functions_02_1.cc @@ -0,0 +1,40 @@ +// --------------------------------------------------------------------- +// +// 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 functions in the Physics namespace compile +// with the various auto-differentiable number types: +// Standard tensors +// +// AD number type: ADOL-C taped + +#include "../ad_common_tests/physics_functions_02.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_physics<2, double, AD::NumberTypes::adolc_taped>(); + test_physics<3, double, AD::NumberTypes::adolc_taped>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/adolc/physics_functions_02_1.output b/tests/adolc/physics_functions_02_1.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/adolc/physics_functions_02_1.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/adolc/physics_functions_02_2.cc b/tests/adolc/physics_functions_02_2.cc new file mode 100644 index 0000000000..d17f96cc75 --- /dev/null +++ b/tests/adolc/physics_functions_02_2.cc @@ -0,0 +1,40 @@ +// --------------------------------------------------------------------- +// +// 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 functions in the Physics namespace compile +// with the various auto-differentiable number types: +// Standard tensors +// +// AD number type: ADOL-C tapeless + +#include "../ad_common_tests/physics_functions_02.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_physics<2, double, AD::NumberTypes::adolc_tapeless>(); + test_physics<3, double, AD::NumberTypes::adolc_tapeless>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/adolc/physics_functions_02_2.output b/tests/adolc/physics_functions_02_2.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/adolc/physics_functions_02_2.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/sacado/physics_functions_02_1.cc b/tests/sacado/physics_functions_02_1.cc new file mode 100644 index 0000000000..d071399260 --- /dev/null +++ b/tests/sacado/physics_functions_02_1.cc @@ -0,0 +1,48 @@ +// --------------------------------------------------------------------- +// +// 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 functions in the Physics namespace compile +// with the various auto-differentiable number types: +// Standard tensors +// +// AD number type: Sacado DFad + +#include "../ad_common_tests/physics_functions_02.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_physics<2, double, AD::NumberTypes::adolc_taped>(); + test_physics<3, double, AD::NumberTypes::adolc_taped>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog.push("Float"); + { + test_physics<2, float, AD::NumberTypes::adolc_taped>(); + test_physics<3, float, AD::NumberTypes::adolc_taped>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/sacado/physics_functions_02_1.output b/tests/sacado/physics_functions_02_1.output new file mode 100644 index 0000000000..ea4ec510fe --- /dev/null +++ b/tests/sacado/physics_functions_02_1.output @@ -0,0 +1,4 @@ + +DEAL::Double:OK +DEAL::Float:OK +DEAL::OK diff --git a/tests/sacado/physics_functions_02_2.cc b/tests/sacado/physics_functions_02_2.cc new file mode 100644 index 0000000000..7a07c37cff --- /dev/null +++ b/tests/sacado/physics_functions_02_2.cc @@ -0,0 +1,48 @@ +// --------------------------------------------------------------------- +// +// 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 functions in the Physics namespace compile +// with the various auto-differentiable number types: +// Standard tensors +// +// AD number type: Sacado DFad-DFad + +#include "../ad_common_tests/physics_functions_02.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_physics<2, double, AD::NumberTypes::sacado_dfad_dfad>(); + test_physics<3, double, AD::NumberTypes::sacado_dfad_dfad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog.push("Float"); + { + test_physics<2, float, AD::NumberTypes::sacado_dfad_dfad>(); + test_physics<3, float, AD::NumberTypes::sacado_dfad_dfad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/sacado/physics_functions_02_2.output b/tests/sacado/physics_functions_02_2.output new file mode 100644 index 0000000000..ea4ec510fe --- /dev/null +++ b/tests/sacado/physics_functions_02_2.output @@ -0,0 +1,4 @@ + +DEAL::Double:OK +DEAL::Float:OK +DEAL::OK diff --git a/tests/sacado/physics_functions_02_3.cc b/tests/sacado/physics_functions_02_3.cc new file mode 100644 index 0000000000..cbc75d322f --- /dev/null +++ b/tests/sacado/physics_functions_02_3.cc @@ -0,0 +1,48 @@ +// --------------------------------------------------------------------- +// +// 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 functions in the Physics namespace compile +// with the various auto-differentiable number types: +// Standard tensors +// +// AD number type: Sacado Rad + +#include "../ad_common_tests/physics_functions_02.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_physics<2, double, AD::NumberTypes::sacado_rad>(); + test_physics<3, double, AD::NumberTypes::sacado_rad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog.push("Float"); + { + test_physics<2, float, AD::NumberTypes::sacado_rad>(); + test_physics<3, float, AD::NumberTypes::sacado_rad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/sacado/physics_functions_02_3.output b/tests/sacado/physics_functions_02_3.output new file mode 100644 index 0000000000..ea4ec510fe --- /dev/null +++ b/tests/sacado/physics_functions_02_3.output @@ -0,0 +1,4 @@ + +DEAL::Double:OK +DEAL::Float:OK +DEAL::OK diff --git a/tests/sacado/physics_functions_02_4.cc b/tests/sacado/physics_functions_02_4.cc new file mode 100644 index 0000000000..777c9a875d --- /dev/null +++ b/tests/sacado/physics_functions_02_4.cc @@ -0,0 +1,48 @@ +// --------------------------------------------------------------------- +// +// 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 functions in the Physics namespace compile +// with the various auto-differentiable number types: +// Standard tensors +// +// AD number type: Sacado Rad-DFad + +#include "../ad_common_tests/physics_functions_02.h" +#include "../tests.h" + +int +main() +{ + initlog(); + + deallog.push("Double"); + { + test_physics<2, double, AD::NumberTypes::sacado_rad_dfad>(); + test_physics<3, double, AD::NumberTypes::sacado_rad_dfad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog.push("Float"); + { + test_physics<2, float, AD::NumberTypes::sacado_rad_dfad>(); + test_physics<3, float, AD::NumberTypes::sacado_rad_dfad>(); + deallog << "OK" << std::endl; + } + deallog.pop(); + + deallog << "OK" << std::endl; +} diff --git a/tests/sacado/physics_functions_02_4.output b/tests/sacado/physics_functions_02_4.output new file mode 100644 index 0000000000..ea4ec510fe --- /dev/null +++ b/tests/sacado/physics_functions_02_4.output @@ -0,0 +1,4 @@ + +DEAL::Double:OK +DEAL::Float:OK +DEAL::OK -- 2.39.5