From: Jean-Paul Pelteret Date: Fri, 12 Oct 2018 13:18:08 +0000 (+0200) Subject: Add tests to check compilation with AD numbers: Physics NS (part 1) X-Git-Tag: v9.1.0-rc1~626^2~2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c252ec054aa92085d7a81c021655ee9d80679753;p=dealii.git Add tests to check compilation with AD numbers: Physics NS (part 1) Physics::Elasticity::Kinematics --- diff --git a/tests/ad_common_tests/physics_functions_01.h b/tests/ad_common_tests/physics_functions_01.h new file mode 100644 index 0000000000..a13a268fa2 --- /dev/null +++ b/tests/ad_common_tests/physics_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 functions in the Physics namespace compile +// with the various auto-differentiable number types: +// Kinematics + +#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: Kinematics, " + << "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 Tensor<2, dim, ADNumberType> F_inv = + invert(F); // Test, since this can be a problem too. + + const Tensor<2, dim, ADNumberType> F_iso = + Physics::Elasticity::Kinematics::F_iso(F); + const SymmetricTensor<2, dim, ADNumberType> F_vol = + Physics::Elasticity::Kinematics::F_vol(F); + const SymmetricTensor<2, dim, ADNumberType> C = + Physics::Elasticity::Kinematics::C(F); + const SymmetricTensor<2, dim, ADNumberType> b = + Physics::Elasticity::Kinematics::b(F); + + const SymmetricTensor<2, dim, ADNumberType> E = + Physics::Elasticity::Kinematics::E(F); + const SymmetricTensor<2, dim, ADNumberType> epsilon = + Physics::Elasticity::Kinematics::epsilon(grad_u); + const SymmetricTensor<2, dim, ADNumberType> e = + Physics::Elasticity::Kinematics::e(F); + + Tensor<2, dim, ADNumberType> dF_dt; + const Tensor<2, dim, ADNumberType> l = + Physics::Elasticity::Kinematics::l(F, dF_dt); + const SymmetricTensor<2, dim, ADNumberType> d = + Physics::Elasticity::Kinematics::d(F, dF_dt); + const Tensor<2, dim, ADNumberType> w = + Physics::Elasticity::Kinematics::w(F, dF_dt); +} diff --git a/tests/adolc/physics_functions_01_1.cc b/tests/adolc/physics_functions_01_1.cc new file mode 100644 index 0000000000..6c715454e1 --- /dev/null +++ b/tests/adolc/physics_functions_01_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: +// Kinematics +// +// AD number type: ADOL-C taped + +#include "../ad_common_tests/physics_functions_01.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_01_1.output b/tests/adolc/physics_functions_01_1.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/adolc/physics_functions_01_1.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/adolc/physics_functions_01_2.cc b/tests/adolc/physics_functions_01_2.cc new file mode 100644 index 0000000000..ff60b10adf --- /dev/null +++ b/tests/adolc/physics_functions_01_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: +// Kinematics +// +// AD number type: ADOL-C tapeless + +#include "../ad_common_tests/physics_functions_01.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_01_2.output b/tests/adolc/physics_functions_01_2.output new file mode 100644 index 0000000000..3c6dd257b1 --- /dev/null +++ b/tests/adolc/physics_functions_01_2.output @@ -0,0 +1,3 @@ + +DEAL::Double:OK +DEAL::OK diff --git a/tests/sacado/physics_functions_01_1.cc b/tests/sacado/physics_functions_01_1.cc new file mode 100644 index 0000000000..eb8044132e --- /dev/null +++ b/tests/sacado/physics_functions_01_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: +// Kinematics +// +// AD number type: Sacado DFad + +#include "../ad_common_tests/physics_functions_01.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_01_1.output b/tests/sacado/physics_functions_01_1.output new file mode 100644 index 0000000000..ea4ec510fe --- /dev/null +++ b/tests/sacado/physics_functions_01_1.output @@ -0,0 +1,4 @@ + +DEAL::Double:OK +DEAL::Float:OK +DEAL::OK diff --git a/tests/sacado/physics_functions_01_2.cc b/tests/sacado/physics_functions_01_2.cc new file mode 100644 index 0000000000..c7b0f5739b --- /dev/null +++ b/tests/sacado/physics_functions_01_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: +// Kinematics +// +// AD number type: Sacado DFad-DFad + +#include "../ad_common_tests/physics_functions_01.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_01_2.output b/tests/sacado/physics_functions_01_2.output new file mode 100644 index 0000000000..ea4ec510fe --- /dev/null +++ b/tests/sacado/physics_functions_01_2.output @@ -0,0 +1,4 @@ + +DEAL::Double:OK +DEAL::Float:OK +DEAL::OK diff --git a/tests/sacado/physics_functions_01_3.cc b/tests/sacado/physics_functions_01_3.cc new file mode 100644 index 0000000000..d87071d16f --- /dev/null +++ b/tests/sacado/physics_functions_01_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: +// Kinematics +// +// AD number type: Sacado Rad + +#include "../ad_common_tests/physics_functions_01.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_01_3.output b/tests/sacado/physics_functions_01_3.output new file mode 100644 index 0000000000..ea4ec510fe --- /dev/null +++ b/tests/sacado/physics_functions_01_3.output @@ -0,0 +1,4 @@ + +DEAL::Double:OK +DEAL::Float:OK +DEAL::OK diff --git a/tests/sacado/physics_functions_01_4.cc b/tests/sacado/physics_functions_01_4.cc new file mode 100644 index 0000000000..33841408e7 --- /dev/null +++ b/tests/sacado/physics_functions_01_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: +// Kinematics +// +// AD number type: Sacado Rad-DFad + +#include "../ad_common_tests/physics_functions_01.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_01_4.output b/tests/sacado/physics_functions_01_4.output new file mode 100644 index 0000000000..ea4ec510fe --- /dev/null +++ b/tests/sacado/physics_functions_01_4.output @@ -0,0 +1,4 @@ + +DEAL::Double:OK +DEAL::Float:OK +DEAL::OK