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
|
// SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
// SPDX-FileCopyrightText: Bradley M. Bell <bradbell@seanet.com>
// SPDX-FileContributor: 2023-23 Bradley M. Bell
# include <cppad/local/val_graph/tape.hpp>
# include <limits>
//
bool test_nan(void)
{ bool ok = true;
//
// tape_t, Vector, addr_t, con_op_enum
using CppAD::local::val_graph::tape_t;
using CppAD::local::val_graph::Vector;
using CppAD::local::val_graph::addr_t;
using CppAD::local::val_graph::op_enum_t;
//
// nan
double nan = std::numeric_limits<double>::quiet_NaN();
//
// tape
tape_t<double> tape;
addr_t n_ind = 1;
addr_t index_of_nan = tape.set_ind(n_ind);
ok &= index_of_nan == n_ind;
//
// op_arg
Vector<addr_t> op_arg(2);
//
// dep_vec
Vector<addr_t> dep_vec(1);
//
// tape
tape.record_con_op(5.0); // not used
dep_vec[0] = tape.record_con_op(nan); // nan
//
// set_dep
tape.set_dep( dep_vec );
//
ok &= tape.n_op() == 2;
ok &= tape.con_vec().size() == 2;
ok &= tape.arg_vec().size() == 2;
//
bool trace = false;
//
// val_vec
Vector<double> val_vec( tape.n_val() );
for(addr_t i = 0; i < n_ind; ++i)
val_vec[i] = 0.0;
tape.eval(trace, val_vec);
//
// renumber
tape.renumber();
//
// dead_code
tape.dead_code();
//
ok &= tape.n_op() == 1;
ok &= tape.con_vec().size() == 1;
ok &= tape.arg_vec().size() == 1;
//
// val_vec
val_vec.resize( tape.n_val() );
for(addr_t i = 0; i < n_ind; ++i)
val_vec[i] = 0.0;
tape.eval(trace, val_vec);
//
// y
Vector<double> y(1);
dep_vec = tape.dep_vec();
y[0] = val_vec[ dep_vec[0] ];
//
// ok
ok &= std::isnan( y[0] );
//
return ok;
}
// END_C++
|
