moebiusband/NuSiF-Solver
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fix bug. Add output for compression.

Browse Source
This commit is contained in:
Jan Eitzinger 2024-03-06 06:02:40 +01:00
parent 6f3d9e73ef
commit 7ce84e2773
2 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
EnhancedSolver/2D-seq/src/grid.h
View File

@ -33,6 +33,6 @@ typedef struct {
static inline int gridIsFluid(Grid* g, int i, int j) static inline int gridIsFluid(Grid* g, int i, int j)
{ {
return g->s[j * (g->imax + 2) + i]; return g->s[j * (g->imax + 2) + i] == FLUID;
} }
#endif // __GRID_H_ #endif // __GRID_H_

3
EnhancedSolver/2D-seq/src/particletracing.c
View File

@ -110,6 +110,8 @@ static void compress(ParticleTracer* p)
Particle tempPool[p->totalParticles]; Particle tempPool[p->totalParticles];
int totalParticles = 0; int totalParticles = 0;
printf("Performing compression ...");
for (int i = 0; i < p->totalParticles; i++) { for (int i = 0; i < p->totalParticles; i++) {
if (memPool[i].flag == 1) { if (memPool[i].flag == 1) {
tempPool[totalParticles].x = memPool[i].x; tempPool[totalParticles].x = memPool[i].x;
@ -119,6 +121,7 @@ static void compress(ParticleTracer* p)
} }
} }
printf(" remove %d particles\n", p->totalParticles - totalParticles);
p->totalParticles = totalParticles; p->totalParticles = totalParticles;
p->removedParticles = 0; p->removedParticles = 0;
p->pointer = totalParticles + 1; p->pointer = totalParticles + 1;