matthew/The-Powder-Toy
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Prevent ambient heat going through walls

Browse Source 
Still doesn't conserve energy, it might disappear when blown towards a wall,
but it will no longer go through the wall.
This commit is contained in:
jacksonmj 2012-02-05 00:30:38 +08:00 committed by Simon Robertshaw
parent 06cb98f705
commit ccb88b6b3b
1 changed files with 14 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
src/air.c
View File

@ -10,6 +10,7 @@ float vx[YRES/CELL][XRES/CELL], ovx[YRES/CELL][XRES/CELL];
float vy[YRES/CELL][XRES/CELL], ovy[YRES/CELL][XRES/CELL];
float pv[YRES/CELL][XRES/CELL], opv[YRES/CELL][XRES/CELL];
unsigned char bmap_blockair[YRES/CELL][XRES/CELL];
unsigned char bmap_blockairh[YRES/CELL][XRES/CELL];
float cb_vx[YRES/CELL][XRES/CELL];
float cb_vy[YRES/CELL][XRES/CELL];
@ -38,7 +39,12 @@ void make_kernel(void) //used for velocity
void update_airh(void)
{
int x, y, i, j;
float dh, dx, dy, f, tx, ty;
float odh, dh, dx, dy, f, tx, ty;
for (y=0; y<YRES/CELL; y++)
for (x=0; x<XRES/CELL; x++)
{
bmap_blockairh[y][x] = (bmap[y][x]==WL_WALL || bmap[y][x]==WL_WALLELEC || bmap[y][x]==WL_GRAV || (bmap[y][x]==WL_EWALL && !emap[y][x]));
}
for (i=0; i<YRES/CELL; i++) //reduces pressure/velocity on the edges every frame
{
hv[i][0] = 295.15f;
@ -68,10 +74,7 @@ void update_airh(void)
{
if (y+j>0 && y+j<YRES/CELL-2 &&
x+i>0 && x+i<XRES/CELL-2 &&
bmap[y+j][x+i]!=WL_WALL &&
bmap[y+j][x+i]!=WL_WALLELEC &&
bmap[y+j][x+i]!=WL_GRAV &&
(bmap[y+j][x+i]!=WL_EWALL || emap[y+j][x+i]))
!bmap_blockairh[y+j][x+i])
{
f = kernel[i+1+(j+1)*3];
dh += hv[y+j][x+i]*f;
@ -95,11 +98,12 @@ void update_airh(void)
ty -= j;
if (i>=2 && i<XRES/CELL-3 && j>=2 && j<YRES/CELL-3)
{
odh = dh;
dh *= 1.0f - AIR_VADV;
dh += AIR_VADV*(1.0f-tx)*(1.0f-ty)*hv[j][i];
dh += AIR_VADV*tx*(1.0f-ty)*hv[j][i+1];
dh += AIR_VADV*(1.0f-tx)*ty*hv[j+1][i];
dh += AIR_VADV*tx*ty*hv[j+1][i+1];
dh += AIR_VADV*(1.0f-tx)*(1.0f-ty)*(bmap_blockairh[j][i] ? odh : hv[j][i]);
dh += AIR_VADV*tx*(1.0f-ty)*(bmap_blockairh[j][i+1] ? odh : hv[j][i+1]);
dh += AIR_VADV*(1.0f-tx)*ty*(bmap_blockairh[j+1][i] ? odh : hv[j+1][i]);
dh += AIR_VADV*tx*ty*(bmap_blockairh[j+1][i+1] ? odh : hv[j+1][i+1]);
}
if(!gravityMode){ //Vertical gravity only for the time being
float airdiff = dh-hv[y][x];