2nd fluid simulator in supercollider

clean-up: #46

here's another supercollider port of msafluid by Memo Akten. it is very similar to the MSAFluidSolver2D i ported from java, but this one is from the c++ version (openframeworks ofxMSAFluid) and has more features and a gui interface.

unfortunately it runs slower than the java port.

classes, helpfiles and an example are attached as a zip archive below.

Package icon MSAFluidSolver.zip17.34 KB

tweets 0220-0225

clean-up: #45

today i recorded and uploaded my latest six supercollider sctweets to my soundcloud account. they were all recorded using sc3.7a0 and audacity for trimming, normalising, adding a quick fadeout and exporting.
recording short soundfiles of these one-liners help non sc-users get an idea about what is going on and it's also good for archiving. they make it possible to, in the future when sc3.7a0 is long gone, go back and listen to what the code did produce.

a=SinOsc;play{RecordBuf.ar(c=InFeedback.ar,b=Buffer.alloc(s,9e4));HPF.ar(a.ar(99,c*6)/9+TGrains.ar(2,a ar:3,b,c+3,2,12),9)}// #SuperCollider


a=SinOsc;play{tanh((c=InFeedback.ar(0,2))+HPF.ar(a.ar(b=1/[5,4],a.ar(a.ar(b*1.1,a.ar(b*2))+a.ar(b*1.4,c,5,4).ceil*99)),9))}// #SuperCollider


a=SinOscFB;play{((c=InFeedback.ar(0,2).lag(b=1/67))+DelayL.ar(HPF.ar(a.ar([99,98]*50.666*a.ar(c+b*b,c),c%2),50),b,b)).tanh}// #SuperCollider


a=LFSaw;play{Splay.ar(BPF.ar(a.ar(f=Duty.ar(a.ar(a.ar(c=3/d=(2..6)))*a.ar(d)/c,0,Dseq(ceil(a ar:d)+d*99,inf))+a.ar(c)),f))}// #SuperCollider


a=SinOsc;play{d=BufRd.ar(2,b=LocalBuf(c=2e5,2).clear,a.ar([2,3]*9)*c,0);BufWr.ar(a.ar(3/[2,3])/3,b,a.ar([99,145]).abs*c);d}// #SuperCollider


a=LFSaw;play{b=(1..8)*99;Splay.ar(CombN.ar(Blip.ar(b/2+a.kr(-8/b,1,99),b/4+a.kr(1/b,1,99))*SinOsc.ar(8/b,a.ar(99/b)))).sin}// #SuperCollider

fluid simulator in supercollider

clean-up: #44

this is a supercollider port of Memo Akten's processing library MSAFluid. i looked at his MSAFluidSolver2D.java version 1.3.0 and rewrote it for sc. i also optimised it as much as i could but still it is slow and the framerate is very low. i blame the slow sc graphics. the fluid system itself should perform ok - it's just the drawing that doesn't live up to being useful.
the good news is that one doesn't need use the fancy graphics. sc is anyway made for sound synthesis and the fluid simulator can run and particles can be mapped to sound without displaying anything.

classes, helpfiles and examples are attached as a zip archive below.

msafluid2d from redFrik on Vimeo.

Package icon MSAFluidSolver2D.zip15.74 KB

failed tweet

clean-up: #43

during today's clean-up-of-old-projects i got stuck shrinking an old piece of supercollider code down to under 140 characters (sctweet). the shortest i could get it down to was 165 characters. still i thought i'd post the results of my efforts here - both if someone have suggestions or if someone want to study the process. it's the same code, just rewritten over and over to be shorter - using every trick i could think of.

of course one could imagine other goals than briefness - better sounding or less cpu taxing are two (probably better and more) common goals when optimising code in this manner.

n=Pbind(\tempo,2.5,\freq,Pseq([Pseq([103.8],14),Pseq([82.4],14)],inf),\amp,Pseq([2,1,1,2,2,1,1]/2,inf),\dur, 0.5);Ptpar([0,n,5,Pmul(\dur,2,Pmul(\freq,2,n)),15,Pmul(\freq,3,n),30,Pmul(\freq,3.5,n),80,Pmul(\dur,3,Pmul(\freq,5.035, n)),140,Pmul(\dur,1,Pmul(\freq,8,n))]).play

n=Pbind(\tempo,2.5,f=\freq,Pseq([Pseq([103.8],14),Pseq([82.4],14)],inf),\amp,Pseq([2,1,1,2,2,1,1]/2,inf),d=\dur, 0.5);Ptpar([0,n,5,Pmul(d,2,Pmul(\freq,2,n)),15,Pmul(f,3,n),30,Pmul(f,3.5,n),80,Pmul(d,3,Pmul(f,5.035, n)),140,Pmul(d,1,Pmul(f,8,n))]).play

n=Pbind(\tempo,2.5,f=\freq,Pseq([Pseq([103.8],14),Pseq([82.4],14)],inf),\amp,Pseq([2,1,1,2,2,1,1]/2,inf),d=\dur, 0.5);Ptpar([0,n,5,Pmul(d,2,Pmul(f,2,n)),15,Pmul(f,3,n),30,Pmul(f,3.5,n),80,Pmul(d,3,Pmul(f,5.035, n)),140,Pmul(d,1,Pmul(f,8,n))]).play

n=Pbind(\tempo,2.5,f=\freq,Pseq([Pseq([103.8],14),Pseq([82.4],14)],inf),\amp,Pseq([2,1,1,2,2,1,1]/2,inf),d=\dur, 0.5);Ptpar([0,n,5,Pmul(d,2,Pmul(f,2,n)),15,Pmul(f,3,n),30,Pmul(f,3.5,n),80,Pmul(d,3,Pmul(f,5.035, n)),140,Pmul(f,8,n)]).play







TempoClock.tempo= 1 //need to reset default clock tempo





n={|x,y=1|Pbind(\freq,Pseq([104,82.5] stutter:14,99)*x,\amp,Pseq([2,1,1,2,2,1,1]/2,99),\dur,y/5)};Ptpar([0,n.(1),2,n.(2,2),6,n.(3),12,n.(3.5),32,n.(5.035,3),56,n.(8)]).play

n={|x,y=1|Pbind(\freq,Pseq([104,82.5] stutter:14,99)*x,\amp,Pseq([2,1,1,2,2,1,1]/2,99),\dur,y/5)};Ptpar([0,n.(1),2,n.(2,2),6,n.(3),12,n.(3.5),32,n.(5.04,3),56,n.(8)]).play

a=[0,1,1,2,2,2,6,3,1,12,3.5,1,32,5.04,3,56,8,1];{|i|Pbind(\lag,a[b=i*3],\freq,Pseq([104,82.5] stutter:14,99)*a[b+1],\amp,Pseq([2,1,1,2,2,1,1]/2,99),\dur,a[b+2]/5).play}!6

[0,1,1,2,2,2,6,3,1,12,3.5,1,32,5.04,3,56,8,1].clump(3).do{|x|Pbind(\freq,Pseq([104,82.5] stutter:14,99)*x[1],\amp,Pseq([2,1,1,2,2,1,1]/2,99),\dur,x[2]/5,\lag,x[0]).play}

deleted quarks

clean-up: #42

started deleting my red quarks

it failed.

update 141024:
i did delete all my red quarks on the 4th of october (because of my 3rd october transition to green phase) but failed with the svn commit. the error i got was Transmitting file data ...svn: E160028: Commit failed (details follow): svn: E160028: Directory '/redSys' is out of date. i saw it as a sign that they shouldn't go.
but today i figured out that i had checked out the http version and that i needed the https to be able to commit. so after the command... svn relocate https://svn.code.sf.net/p/quarks/code i can now commit changes to the quarks svn repository again.


clean-up: #41

yet another small piece of code i found in my supercollider extensions folder.

it's very basic and just generates a pwm signal (pulse-width modulation). you could do the same thing quicker with the standard LFPulse, but sometimes it's good to roll your own just to see another approach.

Pwm {
        *ar {|freq= 100, width= 0.5|
                ^LFSaw.ar(freq)>(width* -2+1);
        *kr {|freq= 100, width= 0.5|
                ^LFSaw.kr(freq)>(width* -2+1);
//mouse controls pwm duty cycle
{Pwm.ar(300, MouseX.kr(0, 1))!2*0.2}.play
//simulating voltage
{Amplitude.ar(Pwm.ar(300, MouseX.kr(0, 1))).poll; DC.ar(0)}.play


clean-up: #40

found an old extension to Wouter Snoei's nice MasterEQ class (it's in his wslib quark). my extension just takes the current setting and post code that recreates that. it's handy for when you want to hardcode a equalizer into some piece of code and when you don't want to be dependent on wslib being installed.

save as extMasterEQ.sc in your extensions folder.

//redFrik 120901
//post synthdef code for current preset
//requires wslib from quarks

+MasterEQ {
        *manual {|strip= true|
                var params;
                if(eq.notNil, {
                        params= eq[\frdb][0][[0, 2, 1]];
                        if(params[2]!=0 or:{strip.not}, {
                                "input= BLowShelf.ar(input, %, %, %);".format(*params).postln;
                        params= eq[\frdb][1][[0, 2, 1]];
                        if(params[2]!=0 or:{strip.not}, {
                                "input= BPeakEQ.ar(input, %, %, %);".format(*params).postln;
                        params= eq[\frdb][2][[0, 2, 1]];
                        if(params[2]!=0 or:{strip.not}, {
                                "input= BPeakEQ.ar(input, %, %, %);".format(*params).postln;
                        params= eq[\frdb][3][[0, 2, 1]];
                        if(params[2]!=0 or:{strip.not}, {
                                "input= BPeakEQ.ar(input, %, %, %);".format(*params).postln;
                        params= eq[\frdb][4][[0, 2, 1]];
                        if(params[2]!=0 or:{strip.not}, {
                                "input= BHiShelf.ar(input, %, %, %);".format(*params).postln;
                }, {
                        "start MasterEQ first".warn;

MasterEQ.manual(false) //keep zero gain poles


clean-up: #39

since long i had this idea of sonifying what my laptop is doing. before, when laptops had spinning harddrives, one could at least hear the faint little noises when there was disk activity. but today with the sd disks, this insight is gone.

below is a draft of how it could work. it tracks activity in the system.log file. to try it start supercollider and run the code. then launch some apps, save files, disconnect network etc etc. it's quite fun to have it running in the background for a while.

one day i plan to do a more thorough sonification, tracking many more things. here is a good reference for what to tap into in the future.

sudo iosnoop
sudo execsnoop -v
sudo opensnoop -ve
sudo dtruss -n SuperCollider
sudo errinfo -c
sudo iotop -CP 1
man -k dtrace

//osx only
s.latency= 0.05;
        var num= 10;
        var delta= 0.1;
        var maxlen= 256;
        var rate= 0.5;
        var curr= List.new;
        var mySplit= {|str|
                var res= "";
                var arr= [];
                var separator= Char.nl.ascii;
                var tab= Char.tab.ascii;
                str.do{|chr, i|
                        if(chr.ascii!=separator or:{str.clipAt(i+1).ascii==tab}, {
                                res= res++chr;
                        }, {
                                arr= arr.add(res);
                                res= "";
        var read= Routine({
                var res, arr, last, new;
                        res= ("tail -n"+num+"/var/log/system.log").unixCmdGetStdOut;
                        if(res!=last, {
                                last= res;
                                arr= mySplit.value(res);
                                new= arr.select{|x| curr.indexOfEqual(x).isNil};
                        delta.wait;     //time between checking system log
        var play= Routine({
                var data;
                        if(curr.size>0, {
                                data= curr.pop;
                                ("buf"++curr.size++": ").post;
                                        Synth(\log, [\data, data.ascii.extend(maxlen, 0)]);
                        }, {
        SynthDef(\log, {|out= 0, amp= 0.5, minFreq= 300, maxFreq= 14400|
                var data= Control.names([\data]).ir(Array.fill(maxlen, 0));
                var freq= Duty.ar(1/maxlen*rate, 0, Dseq(data), 2);
                var src= SinOsc.ar(freq.linexp(0, 255, minFreq, maxFreq), 0, (freq>0).lag(0.02));
                OffsetOut.ar(out, Pan2.ar(src, 0, amp));


Subscribe to f0blog RSS