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
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
|
#! /bin/bash -e
# SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
# SPDX-FileCopyrightText: Bradley M. Bell <bradbell@seanet.com>
# SPDX-FileContributor: 2003-22 Bradley M. Bell
# ----------------------------------------------------------------------------
#!/bin/bash -e
name=`echo $0 | sed -e 's|^bug/||' -e 's|\.sh$||'`
if [ "$0" != "bug/$name.sh" ]
then
echo 'usage: bug/alloc_global.sh'
exit 1
fi
# -----------------------------------------------------------------------------
if [ -e build/bug ]
then
rm -r build/bug
fi
mkdir -p build/bug
cd build/bug
cmake ../..
# -----------------------------------------------------------------------------
cat << EOF
This test results in a very long error message with boost-1.55.0 and
gcc-4.9.2. Here is the gist of the message:
/usr/include/boost/numeric/ublas/detail/matrix_assign.hpp:33:35:
error: no match for ‘operator<’
... snip ...
return norm_inf (e1 - e2) < epsilon *
^
The following change to
/usr/include/boost/numeric/ublas/detail/matrix_assign.hpp
seems to fix the problem:
Old Text:
return norm_inf (e1 - e2) < epsilon *
std::max<S> (std::max<S> (norm_inf (e1), norm_inf (e2)), min_norm);
New Text:
S norm_1 = norm_inf(e1);
S norm_2 = norm_inf(e2);
S norm_diff = norm_inf(e1 - e2);
return norm_diff < epsilon * std::max( std::max<S>(norm_1, norm_2) , min_norm );
EOF
cat << EOF > boost_lu.cpp
#include <cppad/cppad.hpp>
# define NUMERIC_LIMITS_FUN(name) \
static CppAD::AD<double> name(void) \
{ return static_cast< CppAD::AD<double> > ( \
std::numeric_limits<double>::name() \
); \
}
# define NUMERIC_LIMITS_BOOL(name) \
static const bool name = \
std::numeric_limits<double>::name;
# define NUMERIC_LIMITS_INT(name) \
static const int name = \
std::numeric_limits<double>::name;
namespace std {
/// Specialization of numeric_limits< CppAD::AD<double> >
template <>
class numeric_limits< CppAD::AD<double> > {
public:
// has_denorm
static const float_denorm_style has_denorm =
std::numeric_limits<double>::has_denorm;
// round_style
static const float_round_style round_style =
std::numeric_limits<double>::round_style;
// bool
NUMERIC_LIMITS_BOOL(is_specialized);
NUMERIC_LIMITS_BOOL(is_signed);
NUMERIC_LIMITS_BOOL(is_integer);
NUMERIC_LIMITS_BOOL(is_exact);
NUMERIC_LIMITS_BOOL(has_infinity);
NUMERIC_LIMITS_BOOL(has_quiet_NaN);
NUMERIC_LIMITS_BOOL(has_signaling_NaN);
NUMERIC_LIMITS_BOOL(has_denorm_loss);
NUMERIC_LIMITS_BOOL(is_iec559);
NUMERIC_LIMITS_BOOL(is_bounded);
NUMERIC_LIMITS_BOOL(is_modulo);
NUMERIC_LIMITS_BOOL(traps);
NUMERIC_LIMITS_BOOL(tinyness_before);
// int
NUMERIC_LIMITS_INT(digits);
NUMERIC_LIMITS_INT(digits10);
NUMERIC_LIMITS_INT(radix);
NUMERIC_LIMITS_INT(min_exponent);
NUMERIC_LIMITS_INT(min_exponent10);
NUMERIC_LIMITS_INT(max_exponent);
NUMERIC_LIMITS_INT(max_exponent10);
/// functions
NUMERIC_LIMITS_FUN( epsilon)
NUMERIC_LIMITS_FUN( min )
NUMERIC_LIMITS_FUN( max )
};
}
#include <boost/numeric/ublas/lu.hpp>
int main() {
typedef CppAD::AD<double> T;
boost::numeric::ublas::matrix<T> a(5,5);
boost::numeric::ublas::permutation_matrix<std::size_t> pert(5);
// lu decomposition
const std::size_t s = lu_factorize(a, pert);
return 0;
}
EOF
echo "g++ -I../../include --std=c++11 -g $name.cpp -o $name >& $name.log"
if ! g++ -I../../include --std=c++11 -g $name.cpp -o $name >& $name.log
then
cat << EOF > $name.sed
s|\\[|&\\n|g
s|\\]|&\\n|g
s|[‘{}]|&\\n|g
s|[a-zA-Z0-9_]* *=|\\n&|g
#
s|boost::numeric::ublas::||g
s|/usr/include/boost/numeric/ublas/||g
#
s|CppAD::AD<double> *|AD|g
s|<matrix<AD>, basic_unit_lower<> >|<matrix_AD_ulower>|g
s|<matrix<AD>, basic_upper<> >|<matrix_AD_upper>|g
s|scalar_minus<AD, AD>|scalar_minus_AD|g
s|triangular_adaptor<\\([^<>]*\\)>|triangular_\\1|g
s|triangular_matrix_AD_ulower|AD_ulower|g
s|triangular_matrix_AD_upper|AD_upper|g
s|matrix_matrix_prod<AD_ulower, AD_upper, AD> *|AD_prod_ulower_upper|g
s|<AD_ulower, AD_upper, AD_prod_ulower_upper> *|<AD_prod_ulower_upper>|g
s|matrix_matrix_binary<AD_prod_ulower_upper>|AD_prod_ulower_upper|g
EOF
echo "sed -f $name.sed $name.log > ../$name.log"
sed -f $name.sed $name.log > ../$name.log
echo "$name.sh: Compilation Error: see build/bug/$name.log"
exit 1
fi
#
echo "./$name"
if ! ./$name
then
echo
echo "$name.sh: Execution Error"
exit 1
fi
echo
# ------------------------------------------------------------------------------
echo "$name.sh: OK"
exit 0
|
