#include <math.h>
#include <stddef.h>
#include "ts.h"

bool ts_test_ulcer_index( void )
{
    fp_t const x[] = { 10, 12, 11, 13, 12, 14, 11, 13, 12, 15 };
    fp_t const * const inp[1] = { x };

    /* k=3, cumulative HHV (running max from index 0) */
    /* HHV: 10, 12, 12, 13, 13, 14, 14, 14, 14, 15 */
    /* d%:  0, 0, 8.333, 0, 7.692, 0, 21.429, 7.143, 14.286, 0 */
    /* d^2: 0, 0, 69.444, 0, 59.172, 0, 459.184, 51.020, 204.082, 0 */
    /* SMA3 of d^2, then sqrt: */
    /* i=2: window [0,1,2], d^2=[0,0,69.444], mean=23.148, UI=4.811 */
    /* i=3: window [1,2,3], d^2=[0,69.444,0], mean=23.148, UI=4.811 */
    /* i=4: window [2,3,4], d^2=[69.444,0,59.172], mean=42.872, UI=6.548 */
    /* i=5: window [3,4,5], d^2=[0,59.172,0], mean=19.724, UI=4.441 */
    /* i=6: window [4,5,6], d^2=[59.172,0,459.184], mean=172.785, UI=13.145 */
    /* i=7: window [5,6,7], d^2=[0,459.184,51.020], mean=170.069, UI=13.041 */
    /* i=8: window [6,7,8], d^2=[459.184,51.020,204.082], mean=238.095, UI=15.430 */
    /* i=9: window [7,8,9], d^2=[51.020,204.082,0], mean=85.034, UI=9.221 */
    fp_t const e[] = { 4.811252, 4.811252, 6.547672, 4.441156, 13.144774, 13.041013, 15.430335, 9.221389 };
    fp_t const * const exp[1] = { e };

    uint_t k = 3;
    return ts_test(ts_ulcer_index, ts_ulcer_index_offset, inp, exp, sizeof( x ) / sizeof( *x ), (ts_constant_t *)&k, &ts_ulcer_index_p);
}