diff options
Diffstat (limited to 'tran/resonvc.alg')
-rw-r--r-- | tran/resonvc.alg | 33 |
1 files changed, 33 insertions, 0 deletions
diff --git a/tran/resonvc.alg b/tran/resonvc.alg new file mode 100644 index 0000000..74945b8 --- /dev/null +++ b/tran/resonvc.alg @@ -0,0 +1,33 @@ +(RESONVC-ALG +(NAME "resonvc") +(ARGUMENTS ("sound_type" "s1") ("sound_type" "hz") ("double" "bw") + ("int" "normalization")) +(INLINE-INTERPOLATION T) +(INTERNAL-SCALING s1) +(ALWAYS-SCALE hz) +(START (MAX s1 hz)) +(TERMINATE (MIN s1 hz)) +(LOGICAL-STOP (MIN s1)) +(SAMPLE-RATE (MAX s1)) +(STEP-FUNCTION hz) +(STATE ("double" "scale1" "s1->scale") + ("double" "c3co" "exp(bw * -PI2 / s1->sr)") + ("double" "c3p1" "susp->c3co + 1.0") + ("double" "c3t4" "susp->c3co * 4.0") + ("double" "omc3" "1.0 - susp->c3co") + ("double" "c2" "0.0") + ("double" "c1" "0.0") + ("int" "normalization" "normalization") + ("double" "y1" "0.0") + ("double" "y2" "0.0; + hz->scale = (sample_type) (hz->scale * (PI2 / s1->sr))")) +(DEPENDS ("c2" "hz" "c3t4 * cos(hz) / c3p1") + ("c1" "hz" "(normalization == 0 ? scale1 : + (normalization == 1 ? omc3 * sqrt(1.0 - c2 * c2 / c3t4) : + sqrt(c3p1 * c3p1 - c2 * c2) * omc3 / c3p1)) * scale1")) +(CONSTANT "c1" "c2" "c3co" "c3p1" "c3t4" "omc3" "normalization" "scale1") +(FORCE-INTO-REGISTER c3t4 c3p1 normalization omc3 scale1) +(INNER-LOOP "{ double y0 = c1 * s1 + c2 * y1 - c3co * y2; + output = (sample_type) y0; + y2 = y1; y1 = y0; }") +) |