Changeset 4879:3d39c745c8f7 in roaraudio

Timestamp:

04/26/11 14:36:32 (12 years ago)

Author:

phi

Branch:

default

Phase:

public

Message:

added roar_math_diffpoly(), roar_math_intpoly() and roar_math_numintpoly()

Files:

• include/libroardsp/poly.h (modified) (1 diff)
• libroardsp/poly.c (modified) (1 diff)

include/libroardsp/poly.h

@@ -51 +51 @@
 float roar_math_cvpoly_4x4 (float * poly, float t);
 
+int roar_math_diffpoly(float * poly, int len);
+int roar_math_intpoly(float * poly, int len, float c);
+
+float roar_math_numintpoly(float * poly, int len, float st, float et);
+
 #endif
 

libroardsp/poly.c

@@ -141 +141 @@
 }
 
+int roar_math_diffpoly(float * poly, int len) {
+ int i;
+
+ for (i = 1; i < len; i++) {
+  poly[i-1] = poly[i]*(i+1);
+ }
+
+ poly[len-1] = 0;
+
+ return 0;
+}
+
+int roar_math_intpoly(float * poly, int len, float c) {
+ int i;
+
+ for (i = len - 1; i > 0; i--) {
+  poly[i] = poly[i-1]/i;
+ }
+
+ poly[0] = c;
+
+ return 0;
+}
+
+float roar_math_numintpoly(float * poly, int len, float st, float et) {
+ float ipoly[8];
+ float upper, lower;
+
+ if (len > 8) {
+  roar_err_set(ROAR_ERROR_RANGE);
+  return -1.;
+ }
+
+ memcpy(ipoly, poly, sizeof(float)*len);
+ roar_math_intpoly(ipoly, len, 0);
+
+ lower = roar_math_cvpoly(ipoly, st, len+1);
+ upper = roar_math_cvpoly(ipoly, et, len+1);
+
+ //printf("lower(%f)=%f, upper(%f)=%f\n", st, lower, et, upper);
+
+ return upper-lower;
+}
+
 //ll