Ejemplo n.º 1
0
        public static ComplexD Div(double a, ComplexD b)
        {
            var d = (b.real * b.real) + (b.imag * b.imag);
            var s = a / d;

            return(new ComplexD(s * b.real, -s * b.imag));
        }
Ejemplo n.º 2
0
        public static ComplexD Pow(ComplexD a, double b)
        {
            var p = Math.Atan2(a.imag, a.real);
            var m = Math.Pow(a.Mag2(), b * 0.5f);

            return(new ComplexD(m * Math.Cos(p * b), m * Math.Sin(p * b)));
        }
Ejemplo n.º 3
0
        public static ComplexD Div(ComplexD a, ComplexD b)
        {
            var d = (b.real * b.real) + (b.imag * b.imag);
            var s = 1.0 / d;

            return(new ComplexD(
                       ((a.real * b.real) + (a.imag * b.imag)) * s,
                       ((a.imag * b.real) - (a.real * b.imag)) * s
                       ));
        }
Ejemplo n.º 4
0
        private void DrawFilterCurve(
            Rect r,
            float[] coeffs,
            float lowGain,
            float midGain,
            float highGain,
            Color color,
            bool filled,
            double samplerate,
            double magScale)
        {
            var wm = (-2.0 * 3.1415926) / samplerate;

            AudioCurveRendering.AudioCurveEvaluator d = delegate(float x)
            {
                var w   = ComplexD.Exp(wm * GUIHelpers.MapNormalizedFrequency(x, samplerate, useLogScale, true));
                var lpf = ComplexD.Pow(
                    ((w * ((w * coeffs[0]) + coeffs[1])) + coeffs[2]) / ((w * ((w * coeffs[3]) + coeffs[4])) + 1.0f),
                    filterOrder
                    );
                var bpf1 = ComplexD.Pow(
                    ((w * ((w * coeffs[5]) + coeffs[6])) + coeffs[7]) / ((w * ((w * coeffs[8]) + coeffs[9])) + 1.0f),
                    filterOrder
                    );
                var bpf2 = ComplexD.Pow(
                    ((w * ((w * coeffs[10]) + coeffs[11])) + coeffs[12]) / ((w * ((w * coeffs[13]) + coeffs[14])) + 1.0f),
                    filterOrder
                    );
                var hpf = ComplexD.Pow(
                    ((w * ((w * coeffs[15]) + coeffs[16])) + coeffs[17]) / ((w * ((w * coeffs[18]) + coeffs[19])) + 1.0f),
                    filterOrder
                    );
                var h   = (lpf * lowGain).Mag2() + (bpf1 * bpf2 * midGain).Mag2() + (hpf * highGain).Mag2();
                var mag = masterGain + (10.0 * Math.Log10(h));
                return((float)(mag * magScale));
            };

            if (filled)
            {
                AudioCurveRendering.DrawFilledCurve(r, d, color);
            }
            else
            {
                AudioCurveRendering.DrawCurve(r, d, color);
            }
        }
        public void DrawFilterCurve(
            Rect r,
            float[] coeffs,
            bool lowGain,
            bool midGain,
            bool highGain,
            Color color,
            bool useLogScale,
            bool filled,
            double masterGain,
            double samplerate,
            double magScale)
        {
            var wm = (-2.0f * 3.1415926) / samplerate;

            var one = new ComplexD(1.0f, 0.0f);

            AudioCurveRendering.AudioCurveEvaluator d = delegate(float x)
            {
                var w  = ComplexD.Exp(wm * GUIHelpers.MapNormalizedFrequency(x, samplerate, useLogScale, true));
                var hl = !lowGain
                    ? one
                    : ((w * ((w * coeffs[0]) + coeffs[1])) + coeffs[2]) / ((w * ((w * coeffs[3]) + coeffs[4])) + 1.0f);
                var hp = !midGain
                    ? one
                    : ((w * ((w * coeffs[5]) + coeffs[6])) + coeffs[7]) / ((w * ((w * coeffs[8]) + coeffs[9])) + 1.0f);
                var hh = !highGain
                    ? one
                    : ((w * ((w * coeffs[10]) + coeffs[11])) + coeffs[12]) / ((w * ((w * coeffs[13]) + coeffs[14])) + 1.0f);
                var h   = hh * hp * hl;
                var mag = masterGain + (10.0 * Math.Log10(h.Mag2()));
                return((float)(mag * magScale));
            };

            if (filled)
            {
                AudioCurveRendering.DrawFilledCurve(r, d, color);
            }
            else
            {
                AudioCurveRendering.DrawCurve(r, d, color);
            }
        }
        public void DrawFilterCurve(
            Rect r,
            float[] coeffs,
            Color color,
            int numModes,
            bool useLogScale,
            bool filled,
            double samplerate,
            double magScale)
        {
            var wm = (-2.0f * 3.1415926) / samplerate;

            var zero = new ComplexD(0.0f, 0.0f);
            var one  = new ComplexD(1.0f, 0.0f);

            AudioCurveRendering.AudioCurveEvaluator d = delegate(float x)
            {
                var w   = ComplexD.Exp(wm * GUIHelpers.MapNormalizedFrequency(x, samplerate, useLogScale, true));
                var h   = zero;
                var num = numModes * 3;
                for (var n = 0; n < num; n += 3)
                {
                    h += (coeffs[n] * (one - (w * w))) / ((w * ((w * coeffs[n + 2]) + coeffs[n + 1])) + 1.0);
                }

                var mag = 10.0 * Math.Log10(h.Mag2());
                return((float)(mag * magScale));
            };

            if (filled)
            {
                AudioCurveRendering.DrawFilledCurve(r, d, color);
            }
            else
            {
                AudioCurveRendering.DrawCurve(r, d, color);
            }
        }
Ejemplo n.º 7
0
        public static ComplexD Div(ComplexD a, double b)
        {
            var s = 1.0 / b;

            return(new ComplexD(a.real * s, a.imag * s));
        }
Ejemplo n.º 8
0
 public static ComplexD Mul(double a, ComplexD b)
 {
     return(new ComplexD(a * b.real, a * b.imag));
 }
Ejemplo n.º 9
0
 public static ComplexD Mul(ComplexD a, double b)
 {
     return(new ComplexD(a.real * b, a.imag * b));
 }
Ejemplo n.º 10
0
 public static ComplexD Mul(ComplexD a, ComplexD b)
 {
     return(new ComplexD((a.real * b.real) - (a.imag * b.imag), (a.real * b.imag) + (a.imag * b.real)));
 }
Ejemplo n.º 11
0
 public static ComplexD Sub(double a, ComplexD b)
 {
     return(new ComplexD(a - b.real, -b.imag));
 }
Ejemplo n.º 12
0
 public static ComplexD Sub(ComplexD a, double b)
 {
     return(new ComplexD(a.real - b, a.imag));
 }
Ejemplo n.º 13
0
 public static ComplexD Sub(ComplexD a, ComplexD b)
 {
     return(new ComplexD(a.real - b.real, a.imag - b.imag));
 }
Ejemplo n.º 14
0
 public static ComplexD Add(double a, ComplexD b)
 {
     return(new ComplexD(a + b.real, b.imag));
 }
Ejemplo n.º 15
0
 public static ComplexD Add(ComplexD a, double b)
 {
     return(new ComplexD(a.real + b, a.imag));
 }
Ejemplo n.º 16
0
 public static ComplexD Add(ComplexD a, ComplexD b)
 {
     return(new ComplexD(a.real + b.real, a.imag + b.imag));
 }