fluid_sim.js

elem = document.getElementById("container");
setCanvas(elem);

w = WIDTH;
h = HEIGHT;

t = 0;
i = 0;


ii = 5 ;


function gauss_seidel(A,x=[0,0,0],B){
    n=A.length;
    for(i=0;i<n;i++){
        d=B[i];
        for(j=0;j<n;j++){
            if(i!=j){
                d-=A[i][j]*x[j];
            }
        }
        x[i]=d/A[i][i];
    }
    return x;
}
/*
a=[[4, 1, 2],[3, 5, 1],[1, 1, 3]]
x=[0,0,0];
b=[4,7,3];

for(m=0;m<25;m++){
    x=gauss_seidel(a,x,b);
    console.log(x);
}
*/


    



function draw() {

   clearCanvas();
   simulate(0.1);
   for(i=1;i<N;i++){
         for(j=1;j<N;j++){
             val=sqrt(u[IX(i,j)]*u[IX(i,j)]+v[IX(i,j)]*v[IX(i,j)]);
             c=`hsl(${mapRange(val,0,0.0001,0,360)},100%,50%)`;
             new rect(i,j,ii,ii,c,1,"#fff",0);
         }
   }
  
   requestAnimationFrame(draw);


}

///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////

draw();