MetaPhysicL/test/ct_types_unit.C at master · libMesh/MetaPhysicL · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
#include <metaphysicl/ct_types.h>
#include <metaphysicl/dualnumber.h>
#include <metaphysicl/dynamicsparsenumberarray.h>
#include <metaphysicl/dynamicsparsenumbervector.h>
#include "metaphysicl/metaphysicl_exceptions.h"
#include <metaphysicl/semidynamicsparsenumberarray.h>
#include <metaphysicl/numberarray.h>
#include <metaphysicl/numbervector.h>

#define METAPHYSICL_UNIT_ASSERT(expr)                                                              \
  if (!(expr))                                                                                     \
  metaphysicl_error()

#ifndef METAPHYSICL_COMMA
#define METAPHYSICL_COMMA ,
#endif

using namespace MetaPhysicL;

int
main()
{
  MetaPhysicL::enableFPE(true);

  NumberArray<1, float> f_na;
  DynamicSparseNumberArray<float, unsigned int> f_dsna;
  SemiDynamicSparseNumberArray<float, unsigned int, NWrapper<1>> f_sdsna;
  std::array<float, 1> f_a;
  std::vector<float> f_v;

  NumberArray<1, double> d_na;
  DynamicSparseNumberArray<double, unsigned int> d_dsna;
  SemiDynamicSparseNumberArray<double, unsigned int, NWrapper<1>> d_sdsna;
  std::array<double, 1> d_a;
  std::vector<double> d_v;

  // Test container replacement
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<decltype(f_na) METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA decltype(d_na)>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<decltype(f_dsna) METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA decltype(d_dsna)>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<decltype(f_sdsna) METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA decltype(d_sdsna)>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<decltype(f_a) METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA decltype(d_a)>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<decltype(f_v) METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA decltype(d_v)>::value);

  // Test atomic/nuclear algebraic type replacement
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<char METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<signed char METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<unsigned char METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<short int METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<unsigned short int METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<int METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<unsigned int METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<long METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<long long METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<unsigned long METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<unsigned long long METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<float METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<long double METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<DualNumber<double> METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(std::is_same<typename ReplaceAlgebraicType<
                              DynamicSparseNumberVector<double METAPHYSICL_COMMA unsigned int>
                                  METAPHYSICL_COMMA double>::type METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<
          NumberVector<1 METAPHYSICL_COMMA double> METAPHYSICL_COMMA double>::type
                       METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ReplaceAlgebraicType<double METAPHYSICL_COMMA float>::type
                       METAPHYSICL_COMMA float>::value);

  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ValueType<std::vector<double>>::type METAPHYSICL_COMMA double>::value);
  METAPHYSICL_UNIT_ASSERT(
      std::is_same<typename ValueType<std::array<double METAPHYSICL_COMMA 1>>::type
                       METAPHYSICL_COMMA double>::value);
}