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
|
// horizon_sequence.cc -- describes sequence of horizons a processor works on
// $Header$
//
// horizon_sequence::horizon_sequence
// horizon_sequence::~horizon_sequence
// horizon_sequence::sequence_string
// horizon_sequence::append_hn
// horizon_sequence::next_posn
//
#include <stdio.h>
#include <assert.h>
#include "../include/stdc.h"
#include "../jtutil/util.hh"
#include "horizon_sequence.hh"
//******************************************************************************
//
// This function constructs an horizon sequence, which can hold a maximum
// of N_horizons horizon numbers.
//
horizon_sequence::horizon_sequence(int N_horizons_in)
: N_horizons_(N_horizons_in),
my_N_horizons_(0), // sequence starts out empty
posn_(-1),
my_hn_(new int[N_horizons_in])
{ }
//******************************************************************************
//
// This function destroys a horizon sequence.
//
horizon_sequence::~horizon_sequence()
{
delete[] my_hn_;
}
//******************************************************************************
//
// This function returns a (pointer to a) C-style string showing all
// the horizon numbers in the sequence, separated by the C-style string
// sep , eg. "2,3,5".
//
// Note the result points into a private static buffer.
//
char* horizon_sequence::sequence_string(const char sep[]) const
{
const int N_hn_buffer = 10;
char hn_buffer[N_hn_buffer];
const int N_buffer = 100;
static char buffer[N_buffer];
buffer[0] = '\0';
for (int pos = 0 ; pos < my_N_horizons_ ; ++pos)
{
if (pos > 0)
then AHFinderDirect_Strlcat(buffer, sep, N_buffer);
snprintf(hn_buffer, N_hn_buffer, "%d", my_hn_[pos]);
AHFinderDirect_Strlcat(buffer, hn_buffer, N_buffer);
}
return buffer;
}
//******************************************************************************
//
// This function appends hn to the sequence. It returns the new value
// of my_N_horizons().
//
int horizon_sequence::append_hn(int hn)
{
assert( hn > 0 ); // can only append genuine horizons
assert( my_N_horizons_ < N_horizons_ ); // make sure there's space for it
my_hn_[my_N_horizons_++] = hn;
posn_ = 0;
return my_N_horizons_;
}
//******************************************************************************
//
// This function computes the internal position immediately following
// a given internal position in the sequence.
//
// Arguments:
// p = (in) The current internal position, with posn_ semantics
//
// Results:
// This function returns the next internal position after p.
//
int horizon_sequence::next_posn(int pos)
const
{
return (pos < 0) ? pos-1
: (pos+1 < my_N_horizons_) ? pos+1
: -1;
}
|
