studiorailgun/Renderer
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fix gravity test
Some checks failed
studiorailgun/Renderer/pipeline/head There was a failure building this commit
Details

Browse Source
This commit is contained in:
austin 2024-12-18 17:43:32 -05:00
parent dc79f1b0ba
commit efb37b1618
4 changed files with 13 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
src/main/c/includes/fluid/sim/pressurecell/solver_consts.h
View File

@ -47,6 +47,11 @@
*/ */
#define FLUID_PRESSURECELL_DIV_PRESSURE_CONST 5.0f #define FLUID_PRESSURECELL_DIV_PRESSURE_CONST 5.0f
/**
* Cutoff after which density is clamped to zero while diffusing
*/
#define FLUID_PRESSURECELL_MIN_DENSITY_CLAMP_CUTOFF 0.0001f
/** /**
* The number of times to relax most solvers * The number of times to relax most solvers
*/ */

3
src/main/c/src/fluid/sim/pressurecell/density.c
View File

@ -130,6 +130,9 @@ LIBRARY_API void pressurecell_diffuse_density(Environment * environment, Chunk *
) )
) * a ) * a
; ;
if(densityTemp[IX(x,y,z)] < FLUID_PRESSURECELL_MIN_DENSITY_CLAMP_CUTOFF){
densityTemp[IX(x,y,z)] = 0;
}
} }
} }
} }

4
src/test/c/fluid/sim/pressurecell/add_gravity_tests.c
View File

@ -65,10 +65,10 @@ int fluid_sim_pressurecell_add_gravity_test1(){
// //
expected = 0; expected = 0;
actual = currentChunk->v[CENTER_LOC][IX(1,1,DIM-2)]; actual = currentChunk->v[CENTER_LOC][IX(1,1,DIM-2)];
if(fabs(expected - actual) < 0.1f){ if(fabs(expected - actual) < FLUID_PRESSURE_CELL_ERROR_MARGIN){
rVal++; rVal++;
printf("Gravity not applied!\n"); printf("Gravity not applied!\n");
printf("%f \n", currentChunk->v[CENTER_LOC][IX(1,1,DIM-2)]); printf("grav(1,1,DIM-2): %f \n", currentChunk->v[CENTER_LOC][IX(1,1,DIM-2)]);
} }
return rVal; return rVal;

6
src/test/c/fluid/sim/pressurecell/sim_e2e_tests.c
View File

@ -95,10 +95,10 @@ int fluid_sim_pressurecell_sim_e2e_test2(){
// cell that originall had values // cell that originall had values
// //
expected = 0; expected = 0;
actual = currentChunk->d[CENTER_LOC][IX(1,2,DIM-2)]; actual = currentChunk->d[CENTER_LOC][IX(1,4,DIM-2)];
if(fabs(expected - actual) > FLUID_PRESSURE_CELL_ERROR_MARGIN){ if(fabs(expected - actual) > FLUID_PRESSURE_CELL_ERROR_MARGIN){
rVal += assertEqualsFloat(expected,actual,"Failed to recapture density into (1,2,DIM-2)! expected: %f actual: %f \n"); rVal += assertEqualsFloat(expected,actual,"Failed to recapture density into (1,4,DIM-2)! expected: %f actual: %f \n");
printf("%f \n", currentChunk->v[CENTER_LOC][IX(1,2,DIM-2)]); printf("%f \n", currentChunk->v[CENTER_LOC][IX(1,4,DIM-2)]);
} }
expected = MAX_FLUID_VALUE; expected = MAX_FLUID_VALUE;