int main(int argc, char *argv[]){
int nbParts=20000;
float minMax[3][2]={{-1,1},{-1,1},{-1,1}};
float cellSize=2*0.02*1.001;// 2*particle radius + 0.1%
float dt=0.002,simTime=5;
float dtFrame=0.02; //0.04 <=> 25 frames per seconds
int nbDtPerFrame,i;
char filename[]="toto.par";
float mass;
System *sys;
clock_t time;
/*
FILE *file;
float r[3];
file=fopen("filter.plt","w");
fprintf(file,"VARIABLES = x w\n");
r[1]=r[2]=0;
for(i=0;i<=100;i++){
r[0]=i/100.;
fprintf(file,"%e\t%e\n",i/100.,kernelPoly6(r,1));
}
fclose(file);
return 0;*/
nbDtPerFrame=dtFrame/dt;
sys=createSystem(nbParts,minMax);
tmpAddWalls_1(sys);
setGrid(sys, cellSize);
//printGrid(sys); //test nico
//return 0;
mass=1000*(minMax[0][1]-minMax[0][0])*(minMax[1][1]-minMax[1][0])*(minMax[2][1]-minMax[2][0])/nbParts;
setParticles(sys,0,nbParts-1,mass/4); //for water, ro0=1000kg/m3, so 1000 parts/Liter => mi=1e-3kg
setForces(sys);
generateParFile(sys,filename,CREATE);
i=0;
time=clock();
while(sys->time<=simTime){
printf("time = %e (i=%d)\n",sys->time,i);fflush(stdout);
updateSys(sys,dt);
if((++i)%nbDtPerFrame==0)
generateParFile(sys,filename,UPDATE);
}
time=clock()-time;
generateParFile(sys,filename,CLOSE);
deleteSystem(sys);
updateSys(NULL,0);
#ifdef USE_MULTI_THREADING
cleanThreads();
#endif
printf("done in %f sec!\n", (double)(time)/CLOCKS_PER_SEC);
return 0;
}
static void __exit rwlockSimpleCleanup(void) {
printk(KERN_INFO "[PP] Exit routine called.\n");
mainTask = current;
set_current_state(TASK_INTERRUPTIBLE);
cleanThreads();
printk(KERN_INFO "[PP] Called cleanThreads, will now sleep\n");
schedule();
printk(KERN_INFO "[PP] WAKE\n");
kfree(readers);
kfree(writers);
kfree(readCount);
kfree(writeCount);
kfree(foo);
kfree(baseThreadName);
printk(KERN_INFO "[PP] Cleaning up module.\n");
}
