Method sqr2

src/extended/PolyBLEPNode.cpp:246–271

    }

    double sqr2() const
    {
        double t1 = t + 0.875 + 0.25 * (pulseWidth - 0.5);
        t1 -= bitwise_or_zero(t1);

        double t2 = t + 0.375 + 0.25 * (pulseWidth - 0.5);
        t2 -= bitwise_or_zero(t2);

        // Square #1
        double y = t1 < 0.5 ? 1 : -1;

        y += blep(t1, freqInSecondsPerSample) - blep(t2, freqInSecondsPerSample);

        t1 += 0.5 * (1 - pulseWidth);
        t1 -= bitwise_or_zero(t1);

        t2 += 0.5 * (1 - pulseWidth);
        t2 -= bitwise_or_zero(t2);

        // Square #2
        y += t1 < 0.5 ? 1 : -1;

        y += blep(t1, freqInSecondsPerSample) - blep(t2, freqInSecondsPerSample);

        return amplitude * 0.5 * y;
    }

    double rect() const
    {