Fork Code: Walsh-Hadamart transformations in supercollider

name

code content

// more explanation at https://technogems.blogspot.com/2019/05/walsh-hadamart-transformations-in.html
(
var walshfun;
var walshspectrum = Array.newClear(256).fill(0);
var total_result = [];
var resultbuffer = Buffer.new(s, 256*100*256*2);

s.waitForBoot({
    s.freeAllBuffers;
    s.sync;

    ~walsh_transform = {
        | values, rescale=true|
        var number = values.size;
        var output = Array.newClear(number);
        var temp = Array.newClear(number);
        var scale = 1.0;
        if ((number.nextPowerOfTwo != number), {
            "Error: length of values must be power of two.".postln;
            nil ;
        }, /*else*/ {
            /* Mainumber loop. Iteratively compute the hadamard transform */
            var level = 1;
            var stage = 2;
            var work = values;
            number.do({
                |i|
                if ((i.even), {
                    work[i] = work[i] + work[i+1];
                    work[i+1] = work[i] - (2*work[i+1]);
                });
            });

            2.for(number.log2.ceil, {
                |stage|
                var m = 2.pow(level);
                var jj = 0;
                var k = 0;

                while ({k < (number-1)}, {
                    jj.forBy(jj+m-2, 2, {
                        |j|
                        output[k] = work[j] + work[j+m];
                        output[k+1] = work[j] - work[j+m];
                        output[k+2] = work[j+1] - work[j+m+1];
                        output[k+3] = work[j+1] + work[j+m+1];
                        k = k + 4;
                    });

                    jj = jj + (2*m);
                });

                // swap to prepare for next iteration
                temp = work;
                work = output;
                output = temp;

                level = level + 1;
            });

            if (rescale, {
                work/number;
            }, {
                work;
            });
        });
    };

    // listen to the base functions (nothing spectacular here! they sound like PWM as expected)
    total_result = 256.collect({
        | walshidx |
        var hundred_copies = [];
        var localspectrum = walshspectrum.copy();
        localspectrum[walshidx] = 1;
        localspectrum = ~walsh_transform.value(values:localspectrum, rescale:false); // inverse transform
        hundred_copies = 100.collect({
            localspectrum;
        }).flatten;
        (hundred_copies*0.01);
    });
    total_result = total_result.flatten();
    // now turn it into a stereo signal
    resultbuffer = Buffer.loadCollection(s, total_result.stutter(2)/2, 2);

    s.sync;

    resultbuffer.play;

});
)

code description use markdown for formating

category tags comma separated, i.g. "wild, siren" (do not enter default SC class names, please)

ancestor(s) comma separated identificators, i.g. "1-C,1-1,1-4M,1-x"

Private? the code will be accessible by direct url and not visible in public activity

comment of change