diff options
Diffstat (limited to 'lib/starc_bands.c')
| -rw-r--r-- | lib/starc_bands.c | 15 |
1 files changed, 12 insertions, 3 deletions
diff --git a/lib/starc_bands.c b/lib/starc_bands.c index 0c545f8..4975952 100644 --- a/lib/starc_bands.c +++ b/lib/starc_bands.c @@ -26,7 +26,7 @@ ts_parameter_t const ts_starc_bands_p = { .n_inputs = 4, .n_outputs = 3, .n_cons uint_t ts_starc_bands_offset( ts_constant_t const *c ) { ( void )c; - return ( c ? ((ts_starc_bands_params_t *)c)->k : 0 ); + return ( c ? ((ts_starc_bands_params_t const *)c)->k : 0 ); } int ts_starc_bands( @@ -37,11 +37,20 @@ int ts_starc_bands( { TS_PROLOGUE( ts_starc_bands_p ); if ( !c ) return -1; - uint_t k = ((ts_starc_bands_params_t *)c)->k; + TS_CHECK_INPUT(x, x_n, n); + uint_t k = ((ts_starc_bands_params_t const *)c)->k; if (k == 0) return -1; if (n < k + 1) return -1; - fp_t mult = ((ts_starc_bands_params_t *)c)->mult; + fp_t mult = ((ts_starc_bands_params_t const *)c)->mult; + if (mult < 0 || isnan(mult) || isinf(mult)) return -1; + + for( uint_t i = 0; i < k; ++i ) + { + y[0][i] = NAN; + y[1][i] = NAN; + y[2][i] = NAN; + } for( uint_t i = k; i < n; ++i ) { |