From 4f1ee28d687d583601d41ff58e1618b381d2675f Mon Sep 17 00:00:00 2001
From: Katja <katjav@users.sourceforge.net>
Date: Sun, 6 Nov 2011 14:41:44 +0000
Subject: made creb compliant with double precision

- changed float to t_float
- adapted subnormal detection

svn path=/trunk/externals/creb/; revision=15706
---
 modules/bdiag~.c | 26 +++++++++++++-------------
 1 file changed, 13 insertions(+), 13 deletions(-)

(limited to 'modules/bdiag~.c')

diff --git a/modules/bdiag~.c b/modules/bdiag~.c
index ac0a363..be7df79 100644
--- a/modules/bdiag~.c
+++ b/modules/bdiag~.c
@@ -43,8 +43,8 @@ typedef struct bdiag
 } t_bdiag;
 
 
-static float randfloat(void){
-  float r = rand ();
+static t_float randfloat(void){
+  t_float r = rand ();
   r /= (RAND_MAX/2);
   r -= 1;
   return r;
@@ -94,22 +94,22 @@ static void bdiag_eigen(t_bdiag *x, t_floatarg index,
 static void bdiag_timefreq(t_bdiag *x, t_floatarg index, 
 			   t_floatarg time, t_floatarg freq)
 {
-  float r,a,b,n;
-  float sr = sys_getsr() / (float)x->x_ctl.c_order;
+  t_float r,a,b,n;
+  t_float sr = sys_getsr() / (t_float)x->x_ctl.c_order;
 
   /* time in ms */
   time *= 0.001;
 
-  if (time < 0.0f) time = 0.0f;
-  r = pow(0.001f, 1.0f / (time * sr));
-  if (r < 0.0f) r = 0.0f;
-  if (r > 1.0f) r = 1.0f;
+  if (time < 0.0) time = 0.0;
+  r = pow(0.001, 1.0 / (time * sr));
+  if (r < 0.0) r = 0.0;
+  if (r > 1.0) r = 1.0;
 
   a = cos(2*M_PI*freq/sr);
   b = sin(2*M_PI*freq/sr);
  
   /* normalize to be sure */
-  n = 1.0f / sqrt(a*a + b*b);
+  n = 1.0 / sqrt(a*a + b*b);
   a *= n;
   b *= n;
 
@@ -120,14 +120,14 @@ static void bdiag_preset(t_bdiag *x, t_floatarg preset)
 {
   int p = preset;
   int i;
-  float a, b, w, r;
+  t_float a, b, w, r;
 
   switch(p){
   case 0:
     post("preset 0");
     for (i=0; i<x->x_ctl.c_order/2; i++){
       w = randfloat() * .001;
-      r = 1. - (((float)i + 1.)/1000.);
+      r = 1. - (((t_float)i + 1.)/1000.);
       a = cos(w) * r;
       b = sin(w) * r;
       post("%f %f %f %f", w, r, a, b);
@@ -145,8 +145,8 @@ static t_int *bdiag_perform(t_int *w)
 {
 
 
-  t_float *in     = (float *)(w[3]);
-  t_float *out    = (float *)(w[4]);
+  t_float *in     = (t_float *)(w[3]);
+  t_float *out    = (t_float *)(w[4]);
   t_bdiagctl *ctl  = (t_bdiagctl *)(w[1]);
 
   t_float *eigen = ctl->c_eigen;
-- 
cgit v1.2.1