filtering inf/nan in pow~
static t_class *pow_class;
-typedef struct _pow
-{
+typedef struct _pow{
t_object x_obj;
t_inlet *x_inlet;
}t_pow;
-static t_int *pow_perform(t_int *w)
-{
+static t_int *pow_perform(t_int *w){
int nblock = (int)(w[1]);
t_float *in1 = (t_float *)(w[2]);
t_float *in2 = (t_float *)(w[3]);
t_float *out = (t_float *)(w[4]);
- while (nblock--)
- {
- t_float f1 = *in1++;
- t_float f2 = *in2++;
- *out++ = powf(f2, f1);
+ while (nblock--){
+ float f1 = *in1++;
+ float f2 = *in2++;
+ *out++ = (f2 == 0 && f1 < 0) ||
+ (f2 < 0 && (f1 - (int)f1) != 0) ?
+ 0 : pow(f2, f1);
}
return (w + 5);
}
-static void pow_dsp(t_pow *x, t_signal **sp)
-{
+static void pow_dsp(t_pow *x, t_signal **sp){
dsp_add(pow_perform, 4, sp[0]->s_n,
sp[0]->s_vec, sp[1]->s_vec, sp[2]->s_vec);
}
-static void *pow_free(t_pow *x)
-{
+static void *pow_free(t_pow *x){
inlet_free(x->x_inlet);
return (void *)x;
}
-static void *pow_new(t_floatarg f)
-{
+static void *pow_new(t_floatarg f){
t_pow *x = (t_pow *)pd_new(pow_class);
x->x_inlet = inlet_new((t_object *)x, (t_pd *)x, &s_signal, &s_signal);
pd_float((t_pd *)x->x_inlet, f);
@@ -49,8 +45,7 @@ static void *pow_new(t_floatarg f)
return (x);
}
-void pow_tilde_setup(void)
-{
+void pow_tilde_setup(void){
pow_class = class_new(gensym("Pow~"), (t_newmethod)pow_new,
(t_method)pow_free, sizeof(t_pow), CLASS_DEFAULT, A_DEFFLOAT, 0);
class_addcreator((t_newmethod)pow_new, gensym("cyclone/pow~"), A_GIMME, 0);
@@ -59,7 +54,6 @@ void pow_tilde_setup(void)
class_addmethod(pow_class, (t_method)pow_dsp, gensym("dsp"), A_CANT, 0);
}
-void Pow_tilde_setup(void)
-{
+void Pow_tilde_setup(void){
pow_tilde_setup();
}
