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
|
// cpm_map.cc -- "integer +/-" mapping i --> j = const +/- i
// $Header$
//
// jtutil::cpm_map::cpm_map # mirror map, specified by fixed point
// jtutil::cpm_map::cpm_map # generic map specified by sample point & sign
// jtutil::cpm_map::cpm_map # generic map specified by *fp* sample point & sign
//
// ***** template instantiations *****
//
#include <assert.h>
#include <stdio.h>
#include "stdc.h"
#include "util.hh"
#include "cpm_map.hh"
using jtutil::error_exit;
//******************************************************************************
//******************************************************************************
//******************************************************************************
//
// This function constructs a cpm_map object with a "-" sign and a
// specified fixed point (must be integer or half-integer) and domain.
// The sample point need not be in the map's domain/range.
//
namespace jtutil
{
template <typename fp>
cpm_map<fp>::cpm_map(int min_i_in, int max_i_in,
fp fixed_point)
: min_i_(min_i_in), max_i_(max_i_in),
map_is_plus_(false)
{
const fp d_offset = 2.0 * fixed_point;
if (! fuzzy<fp>::is_integer(d_offset))
then error_exit(ERROR_EXIT,
"***** cpm_map::cpm_map (mirror):\n"
" fixed_point=%g isn't (fuzzily) integral or half-integral!\n"
,
double(fixed_point)); /*NOTREACHED*/
offset_ = round<fp>::to_integer(d_offset);
// verify that we have setup correct
assert(
map_unchecked(fuzzy<fp>::floor (fixed_point))
==
fuzzy<fp>::ceiling(fixed_point)
);
}
} // namespace jtutil::
//******************************************************************************
//
// This function constructs a generic cpm_map object, with the mapping
// specified by a sample point sample_i --> sample_j and by sign.
// The sample point need not be in the map's domain/range.
//
namespace jtutil
{
template <typename fp>
cpm_map<fp>::cpm_map(int min_i_in, int max_i_in,
int sample_i, int sample_j,
bool map_is_plus_in)
: min_i_(min_i_in), max_i_(max_i_in),
offset_(map_is_plus_in ? sample_j - sample_i
: sample_j + sample_i),
map_is_plus_(map_is_plus_in)
{
// verify that we have setup correct
assert( map_unchecked(sample_i) == sample_j );
}
} // namespace jtutil::
//******************************************************************************
//
// This function constructs a generic cpm_map object, with the mapping
// specified by a *fp* sample point sample_i --> sample_j (which
// must specify an integer --> integer mapping, i.e. 4.2 --> 4.2 is
// ok for a + map, and 4.5 --> 4.5 is ok for a minus map, but 4.2 --> 4.7
// is never ok) and by sign. The sample point need not be in the map's
// domain/range.
//
namespace jtutil
{
template <typename fp>
cpm_map<fp>::cpm_map(int min_i_in, int max_i_in,
fp sample_i, fp sample_j,
bool map_is_plus_in)
: min_i_(min_i_in), max_i_(max_i_in),
map_is_plus_(map_is_plus_in)
{
const fp fp_offset = map_is_plus_in ? sample_j - sample_i
: sample_j + sample_i;
if (! fuzzy<fp>::is_integer(fp_offset))
then error_exit(ERROR_EXIT,
"***** cpm_map::cpm_map (generic via fp sample point):\n"
" fp_offset=%g isn't fuzzily integral!\n"
" ==> sample_i=%g --> sample_j=%g\n"
" doesn't fuzzily specify an integer --> integer mapping!\n"
,
double(fp_offset),
double(sample_i), double(sample_j)); /*NOTREACHED*/
offset_ = round<fp>::to_integer(fp_offset);
// verify that we have setup correct
assert(
map_unchecked( fuzzy<fp>::floor(sample_i))
==
(map_is_plus_in ? fuzzy<fp>::floor (sample_j)
: fuzzy<fp>::ceiling(sample_j)) );
}
} // namespace jtutil::
//******************************************************************************
//******************************************************************************
//******************************************************************************
//
// ***** template instantiations *****
//
// instantiation for <double>
template class jtutil::cpm_map<double>;
|
