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
|
// Ultimate Oscillator
//
// ultimate_oscillator
//
// Ultimate Oscillator
//
//
//
// Inputs: 4
// Outputs: 1
// Constants: 0
#include <math.h>
#include "ts.h"
ts_parameter_t const ts_ultimate_oscillator_p = { .n_inputs = 4, .n_outputs = 1, .n_constants = 0, .params_size = 0 };
uint_t ts_ultimate_oscillator_offset( ts_constant_t const *c )
{
( void )c;
return ( 54 );
}
int ts_ultimate_oscillator(
fp_t * const *y, uint_t y_n,
fp_t const * const *x, uint_t x_n,
uint_t n, ts_constant_t const *c
)
{
TS_PROLOGUE( ts_ultimate_oscillator_p );
(void)c;
for( uint_t i = 54; i < n; ++i )
{
fp_t bp7 = 0, tr7 = 0;
fp_t bp14 = 0, tr14 = 0;
fp_t bp28 = 0, tr28 = 0;
for( uint_t j = 0; j < 7; ++j )
{
uint_t idx = i - 6 + j;
fp_t ll = ts_l( x, idx );
fp_t hh = ts_h( x, idx );
for( uint_t k = 1; k < 7; ++k )
{
ll = fmin( ll, ts_l( x, idx - k ) );
hh = fmax( hh, ts_h( x, idx - k ) );
}
bp7 += ts_c( x, idx ) - ll;
tr7 += hh - ll;
}
for( uint_t j = 0; j < 14; ++j )
{
uint_t idx = i - 13 + j;
fp_t ll = ts_l( x, idx );
fp_t hh = ts_h( x, idx );
for( uint_t k = 1; k < 14; ++k )
{
ll = fmin( ll, ts_l( x, idx - k ) );
hh = fmax( hh, ts_h( x, idx - k ) );
}
bp14 += ts_c( x, idx ) - ll;
tr14 += hh - ll;
}
for( uint_t j = 0; j < 28; ++j )
{
uint_t idx = i - 27 + j;
fp_t ll = ts_l( x, idx );
fp_t hh = ts_h( x, idx );
for( uint_t k = 1; k < 28; ++k )
{
ll = fmin( ll, ts_l( x, idx - k ) );
hh = fmax( hh, ts_h( x, idx - k ) );
}
bp28 += ts_c( x, idx ) - ll;
tr28 += hh - ll;
}
fp_t denom = 7.0 * tr7 + 3.0 * tr14 + tr28;
fp_t numer = 7.0 * bp7 + 3.0 * bp14 + bp28;
y[0][i] = denom != 0 ? 100.0 * numer / denom : 50.0;
}
return 0;
}
|