summaryrefslogtreecommitdiff
path: root/src/powder.c
diff options
context:
space:
mode:
authorSimon Robertshaw <simon@hardwired.org.uk>2012-11-17 19:43:59 (GMT)
committer Simon Robertshaw <simon@hardwired.org.uk>2012-11-17 19:43:59 (GMT)
commite3594aba9e05c6865d396418c028049cda92c2f3 (patch)
tree4b825dc642cb6eb9a060e54bf8d69288fbee4904 /src/powder.c
parentfb43f7d23e99765ae093fc45608901cb5907d1d8 (diff)
downloadpowder-e3594aba9e05c6865d396418c028049cda92c2f3.zip
powder-e3594aba9e05c6865d396418c028049cda92c2f3.tar.gz
Remove old code
Diffstat (limited to 'src/powder.c')
-rw-r--r--src/powder.c3751
1 files changed, 0 insertions, 3751 deletions
diff --git a/src/powder.c b/src/powder.c
deleted file mode 100644
index 9f1970e..0000000
--- a/src/powder.c
+++ /dev/null
@@ -1,3751 +0,0 @@
-/**
- * Powder Toy - particle simulation
- *
- * This program is free software; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 3 of the License, or
- * (at your option) any later version.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
- * GNU General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program. If not, see <http://www.gnu.org/licenses/>.
- */
-
-#include <stdint.h>
-#include <math.h>
-#include <defines.h>
-#include <powder.h>
-#include <air.h>
-#include <misc.h>
-#include "gravity.h"
-#ifdef LUACONSOLE
-#include <luaconsole.h>
-#endif
-
-int wire_placed = 0;
-
-int lighting_recreate = 0;
-int force_stacking_check = 0;//whether to force a check for excessively stacked particles
-
-playerst player;
-playerst player2;
-
-playerst fighters[256]; //255 is the maximum number of fighters
-unsigned char fighcount = 0; //Contains the number of fighters
-
-particle *parts;
-particle *cb_parts;
-
-int airMode = 0;
-
-
-unsigned char bmap[YRES/CELL][XRES/CELL];
-unsigned char emap[YRES/CELL][XRES/CELL];
-
-unsigned char cb_bmap[YRES/CELL][XRES/CELL];
-unsigned char cb_emap[YRES/CELL][XRES/CELL];
-
-int pfree;
-
-unsigned pmap[YRES][XRES];
-int pmap_count[YRES][XRES];
-unsigned cb_pmap[YRES][XRES];
-unsigned photons[YRES][XRES];
-
-void get_gravity_field(int x, int y, float particleGrav, float newtonGrav, float *pGravX, float *pGravY)
-{
- *pGravX = newtonGrav*gravx[(y/CELL)*(XRES/CELL)+(x/CELL)];
- *pGravY = newtonGrav*gravy[(y/CELL)*(XRES/CELL)+(x/CELL)];
- switch (gravityMode)
- {
- default:
- case 0: //normal, vertical gravity
- *pGravY += particleGrav;
- break;
- case 1: //no gravity
- break;
- case 2: //radial gravity
- if (x-XCNTR != 0 || y-YCNTR != 0)
- {
- float pGravMult = particleGrav/sqrtf((x-XCNTR)*(x-XCNTR) + (y-YCNTR)*(y-YCNTR));
- *pGravX -= pGravMult * (float)(x - XCNTR);
- *pGravY -= pGravMult * (float)(y - YCNTR);
- }
- }
-}
-
-static int pn_junction_sprk(int x, int y, int pt)
-{
- unsigned r = pmap[y][x];
- if ((r & 0xFF) != pt)
- return 0;
- r >>= 8;
- if (parts[r].type != pt)
- return 0;
- if (parts[r].life != 0)
- return 0;
-
- parts[r].ctype = pt;
- part_change_type(r,x,y,PT_SPRK);
- parts[r].life = 4;
- return 1;
-}
-
-static void photoelectric_effect(int nx, int ny)//create sparks from PHOT when hitting PSCN and NSCN
-{
- unsigned r = pmap[ny][nx];
-
- if ((r&0xFF) == PT_PSCN) {
- if ((pmap[ny][nx-1] & 0xFF) == PT_NSCN ||
- (pmap[ny][nx+1] & 0xFF) == PT_NSCN ||
- (pmap[ny-1][nx] & 0xFF) == PT_NSCN ||
- (pmap[ny+1][nx] & 0xFF) == PT_NSCN)
- pn_junction_sprk(nx, ny, PT_PSCN);
- }
-}
-
-unsigned char can_move[PT_NUM][PT_NUM];
-
-void init_can_move()
-{
- // can_move[moving type][type at destination]
- // 0 = No move/Bounce
- // 1 = Swap
- // 2 = Both particles occupy the same space.
- // 3 = Varies, go run some extra checks
- int t, rt, stkm_move;
- for (rt=0;rt<PT_NUM;rt++)
- can_move[0][rt] = 0; // particles that don't exist shouldn't move...
- for (t=1;t<PT_NUM;t++)
- for (rt=0;rt<PT_NUM;rt++)
- can_move[t][rt] = 1;
- for (rt=1;rt<PT_NUM;rt++)
- {
- can_move[PT_PHOT][rt] = 2;
- }
- for (t=1;t<PT_NUM;t++)
- {
- for (rt=1;rt<PT_NUM;rt++)
- {
- // weight check, also prevents particles of same type displacing each other
- if (ptypes[t].weight <= ptypes[rt].weight || rt==PT_GEL) can_move[t][rt] = 0;
- if (t==PT_NEUT && ptypes[rt].properties&PROP_NEUTPASS)
- can_move[t][rt] = 2;
- if (t==PT_NEUT && ptypes[rt].properties&PROP_NEUTABSORB)
- can_move[t][rt] = 1;
- if (t==PT_NEUT && ptypes[rt].properties&PROP_NEUTPENETRATE)
- can_move[t][rt] = 1;
- if (ptypes[t].properties&PROP_NEUTPENETRATE && rt==PT_NEUT)
- can_move[t][rt] = 0;
- if (ptypes[t].properties&TYPE_ENERGY && ptypes[rt].properties&TYPE_ENERGY)
- can_move[t][rt] = 2;
- }
- }
- can_move[PT_DEST][PT_DMND] = 0;
- can_move[PT_DEST][PT_CLNE] = 0;
- can_move[PT_DEST][PT_PCLN] = 0;
- can_move[PT_DEST][PT_BCLN] = 0;
- can_move[PT_DEST][PT_PBCN] = 0;
- can_move[PT_BIZR][PT_FILT] = 2;
- can_move[PT_BIZRG][PT_FILT] = 2;
- for (t=0;t<PT_NUM;t++)
- {
- //spark shouldn't move
- can_move[PT_SPRK][t] = 0;
- stkm_move = 0;
- if (ptypes[t].properties & (TYPE_LIQUID | TYPE_GAS))
- stkm_move = 2;
- if (!t || t==PT_PRTO || t==PT_SPAWN || t==PT_SPAWN2)
- stkm_move = 2;
- can_move[PT_STKM][t] = stkm_move;
- can_move[PT_STKM2][t] = stkm_move;
- can_move[PT_FIGH][t] = stkm_move;
- }
- for (t=0;t<PT_NUM;t++)
- {
- // make them eat things
- can_move[t][PT_BHOL] = 1;
- can_move[t][PT_NBHL] = 1;
- can_move[t][PT_STKM] = 0;
- can_move[t][PT_STKM2] = 0;
- can_move[t][PT_FIGH] = 0;
- //INVIS behaviour varies with pressure
- can_move[t][PT_INVIS] = 3;
- //stop CNCT being displaced by other particles
- can_move[t][PT_CNCT] = 0;
- //void behaviour varies with powered state and ctype
- can_move[t][PT_PVOD] = 3;
- can_move[t][PT_VOID] = 3;
- }
- for (t=0;t<PT_NUM;t++)
- {
- if (t==PT_GLAS || t==PT_PHOT || t==PT_CLNE || t==PT_PCLN
- || t==PT_GLOW || t==PT_WATR || t==PT_DSTW || t==PT_SLTW
- || t==PT_ISOZ || t==PT_ISZS || t==PT_FILT || t==PT_INVIS
- || t==PT_QRTZ || t==PT_PQRT)
- can_move[PT_PHOT][t] = 2;
- }
- can_move[PT_ELEC][PT_LCRY] = 2;
- can_move[PT_ELEC][PT_EXOT] = 2;
- can_move[PT_NEUT][PT_EXOT] = 2;
- can_move[PT_PHOT][PT_LCRY] = 3;//varies according to LCRY life
-
- can_move[PT_PHOT][PT_BIZR] = 2;
- can_move[PT_ELEC][PT_BIZR] = 2;
- can_move[PT_PHOT][PT_BIZRG] = 2;
- can_move[PT_ELEC][PT_BIZRG] = 2;
- can_move[PT_PHOT][PT_BIZRS] = 2;
- can_move[PT_ELEC][PT_BIZRS] = 2;
-
- can_move[PT_NEUT][PT_INVIS] = 2;
- //whol eats anar
- can_move[PT_ANAR][PT_WHOL] = 1;
- can_move[PT_ANAR][PT_NWHL] = 1;
- can_move[PT_THDR][PT_THDR] = 2;
- can_move[PT_EMBR][PT_EMBR] = 2;
-}
-
-/*
- RETURN-value explenation
-1 = Swap
-0 = No move/Bounce
-2 = Both particles occupy the same space.
- */
-int eval_move(int pt, int nx, int ny, unsigned *rr)
-{
- unsigned r;
- int result;
-
- if (nx<0 || ny<0 || nx>=XRES || ny>=YRES)
- return 0;
-
- r = pmap[ny][nx];
- if (r)
- r = (r&~0xFF) | parts[r>>8].type;
- if (rr)
- *rr = r;
- if (pt>=PT_NUM || (r&0xFF)>=PT_NUM)
- return 0;
- result = can_move[pt][r&0xFF];
- if (result==3)
- {
- if ((pt==PT_PHOT || pt==PT_ELEC) && (r&0xFF)==PT_LCRY)
- result = (parts[r>>8].life > 5)? 2 : 0;
- if ((r&0xFF)==PT_INVIS)
- {
- if (pv[ny/CELL][nx/CELL]>4.0f || pv[ny/CELL][nx/CELL]<-4.0f) result = 2;
- else result = 0;
- }
- if ((r&0xFF)==PT_PVOD)
- {
- if (parts[r>>8].life == 10)
- {
- if(!parts[r>>8].ctype || (parts[r>>8].ctype==pt)!=(parts[r>>8].tmp&1))
- result = 1;
- else
- result = 0;
- }
- else result = 0;
- }
- if ((r&0xFF)==PT_VOID)
- {
- if(!parts[r>>8].ctype || (parts[r>>8].ctype==pt)!=(parts[r>>8].tmp&1))
- result = 1;
- else
- result = 0;
- }
- }
- if (bmap[ny/CELL][nx/CELL])
- {
- if (bmap[ny/CELL][nx/CELL]==WL_ALLOWGAS && !(ptypes[pt].properties&TYPE_GAS))// && ptypes[pt].falldown!=0 && pt!=PT_FIRE && pt!=PT_SMKE)
- return 0;
- if (bmap[ny/CELL][nx/CELL]==WL_ALLOWENERGY && !(ptypes[pt].properties&TYPE_ENERGY))// && ptypes[pt].falldown!=0 && pt!=PT_FIRE && pt!=PT_SMKE)
- return 0;
- if (bmap[ny/CELL][nx/CELL]==WL_ALLOWLIQUID && ptypes[pt].falldown!=2)
- return 0;
- if (bmap[ny/CELL][nx/CELL]==WL_ALLOWSOLID && ptypes[pt].falldown!=1)
- return 0;
- if (bmap[ny/CELL][nx/CELL]==WL_ALLOWAIR || bmap[ny/CELL][nx/CELL]==WL_WALL || bmap[ny/CELL][nx/CELL]==WL_WALLELEC)
- return 0;
- if (bmap[ny/CELL][nx/CELL]==WL_EWALL && !emap[ny/CELL][nx/CELL])
- return 0;
- if (bmap[ny/CELL][nx/CELL]==WL_EHOLE && !emap[ny/CELL][nx/CELL] && !(ptypes[pt].properties&TYPE_SOLID) && !(ptypes[r&0xFF].properties&TYPE_SOLID))
- return 2;
- }
- return result;
-}
-
-int try_move(int i, int x, int y, int nx, int ny)
-{
- unsigned r, e;
-
- if (x==nx && y==ny)
- return 1;
- if (nx<0 || ny<0 || nx>=XRES || ny>=YRES)
- return 1;
-
- e = eval_move(parts[i].type, nx, ny, &r);
-
- /* half-silvered mirror */
- if (!e && parts[i].type==PT_PHOT &&
- (((r&0xFF)==PT_BMTL && rand()<RAND_MAX/2) ||
- (pmap[y][x]&0xFF)==PT_BMTL))
- e = 2;
-
- if (!e) //if no movement
- {
- if (!(ptypes[parts[i].type].properties & TYPE_ENERGY))
- return 0;
- if (!legacy_enable && parts[i].type==PT_PHOT && r)//PHOT heat conduction
- {
- if ((r & 0xFF) == PT_COAL || (r & 0xFF) == PT_BCOL)
- parts[r>>8].temp = parts[i].temp;
-
- if ((r & 0xFF) < PT_NUM && ptypes[r&0xFF].hconduct && ((r&0xFF)!=PT_HSWC||parts[r>>8].life==10) && (r&0xFF)!=PT_FILT)
- parts[i].temp = parts[r>>8].temp = restrict_flt((parts[r>>8].temp+parts[i].temp)/2, MIN_TEMP, MAX_TEMP);
- }
- if ((parts[i].type==PT_NEUT || parts[i].type==PT_ELEC) && ((r&0xFF)==PT_CLNE || (r&0xFF)==PT_PCLN || (r&0xFF)==PT_BCLN || (r&0xFF)==PT_PBCN)) {
- if (!parts[r>>8].ctype)
- parts[r>>8].ctype = parts[i].type;
- }
- if ((r&0xFF)==PT_PRTI && (ptypes[parts[i].type].properties & TYPE_ENERGY))
- {
- int nnx, count;
- for (count=0; count<8; count++)
- {
- if (isign(x-nx)==isign(portal_rx[count]) && isign(y-ny)==isign(portal_ry[count]))
- break;
- }
- count = count%8;
- parts[r>>8].tmp = (int)((parts[r>>8].temp-73.15f)/100+1);
- if (parts[r>>8].tmp>=CHANNELS) parts[r>>8].tmp = CHANNELS-1;
- else if (parts[r>>8].tmp<0) parts[r>>8].tmp = 0;
- for ( nnx=0; nnx<80; nnx++)
- if (!portalp[parts[r>>8].tmp][count][nnx].type)
- {
- portalp[parts[r>>8].tmp][count][nnx] = parts[i];
- parts[i].type=PT_NONE;
- break;
- }
- }
- return 0;
- }
-
- if (e == 2) //if occupy same space
- {
- if (parts[i].type == PT_PHOT && (r&0xFF)==PT_GLOW && !parts[r>>8].life)
- if (rand() < RAND_MAX/30)
- {
- parts[r>>8].life = 120;
- create_gain_photon(i);
- }
- if (parts[i].type == PT_PHOT && (r&0xFF)==PT_FILT)
- {
- int temp_bin = (int)((parts[r>>8].temp-273.0f)*0.025f);
- if (temp_bin < 0) temp_bin = 0;
- if (temp_bin > 25) temp_bin = 25;
- if(!parts[r>>8].tmp){
- parts[i].ctype = 0x1F << temp_bin; //Assign Colour
- } else if(parts[r>>8].tmp==1){
- parts[i].ctype &= 0x1F << temp_bin; //Filter Colour
- } else if(parts[r>>8].tmp==2){
- parts[i].ctype |= 0x1F << temp_bin; //Add Colour
- } else if(parts[r>>8].tmp==3){
- parts[i].ctype &= ~(0x1F << temp_bin); //Subtract Colour
- }
- }
- if (parts[i].type == PT_NEUT && (r&0xFF)==PT_GLAS) {
- if (rand() < RAND_MAX/10)
- create_cherenkov_photon(i);
- }
- if (parts[i].type == PT_PHOT && (r&0xFF)==PT_INVIS && pv[ny/CELL][nx/CELL]<=4.0f && pv[ny/CELL][nx/CELL]>=-4.0f) {
- part_change_type(i,x,y,PT_NEUT);
- parts[i].ctype = 0;
- }
- if ((parts[i].type==PT_BIZR||parts[i].type==PT_BIZRG) && (r&0xFF)==PT_FILT)
- {
- int temp_bin = (int)((parts[r>>8].temp-273.0f)*0.025f);
- if (temp_bin < 0) temp_bin = 0;
- if (temp_bin > 25) temp_bin = 25;
- parts[i].ctype = 0x1F << temp_bin;
- }
- return 1;
- }
- //else e=1 , we are trying to swap the particles, return 0 no swap/move, 1 is still overlap/move, because the swap takes place later
-
- if (parts[i].type==PT_NEUT && (ptypes[r&0xFF].properties&PROP_NEUTABSORB))
- {
- kill_part(i);
- return 0;
- }
- if ((r&0xFF)==PT_VOID || (r&0xFF)==PT_PVOD) //this is where void eats particles
- {
- //void ctype already checked in eval_move
- kill_part(i);
- return 0;
- }
- if ((r&0xFF)==PT_BHOL || (r&0xFF)==PT_NBHL) //this is where blackhole eats particles
- {
- kill_part(i);
- if (!legacy_enable)
- {
- parts[r>>8].temp = restrict_flt(parts[r>>8].temp+parts[i].temp/2, MIN_TEMP, MAX_TEMP);//3.0f;
- }
-
- return 0;
- }
- if (((r&0xFF)==PT_WHOL||(r&0xFF)==PT_NWHL) && parts[i].type==PT_ANAR) //whitehole eats anar
- {
- kill_part(i);
- if (!legacy_enable)
- {
- parts[r>>8].temp = restrict_flt(parts[r>>8].temp- (MAX_TEMP-parts[i].temp)/2, MIN_TEMP, MAX_TEMP);
- }
-
- return 0;
- }
-
- if (parts[i].type==PT_CNCT && y<ny && (pmap[y+1][x]&0xFF)==PT_CNCT)//check below CNCT for another CNCT
- return 0;
-
- if ((bmap[y/CELL][x/CELL]==WL_EHOLE && !emap[y/CELL][x/CELL]) && !(bmap[ny/CELL][nx/CELL]==WL_EHOLE && !emap[ny/CELL][nx/CELL]))
- return 0;
-
- if(parts[i].type==PT_GBMB&&parts[i].life>0)
- return 0;
-
- e = r >> 8; //e is now the particle number at r (pmap[ny][nx])
- if (r)//the swap part, if we make it this far, swap
- {
- if (parts[i].type==PT_NEUT) {
- // target material is NEUTPENETRATE, meaning it gets moved around when neutron passes
- unsigned s = pmap[y][x];
- if (s && !(ptypes[s&0xFF].properties&PROP_NEUTPENETRATE))
- return 1; // if the element currently underneath neutron isn't NEUTPENETRATE, don't move anything except the neutron
- // if nothing is currently underneath neutron, only move target particle
- if (s)
- {
- pmap[ny][nx] = (s&~(0xFF))|parts[s>>8].type;
- parts[s>>8].x = nx;
- parts[s>>8].y = ny;
- }
- else pmap[ny][nx] = 0;
- parts[e].x = x;
- parts[e].y = y;
- pmap[y][x] = (e<<8)|parts[e].type;
- return 1;
- }
-
- if ((pmap[ny][nx]>>8)==e) pmap[ny][nx] = 0;
- parts[e].x += x-nx;
- parts[e].y += y-ny;
- pmap[(int)(parts[e].y+0.5f)][(int)(parts[e].x+0.5f)] = (e<<8)|parts[e].type;
- }
- return 1;
-}
-
-// try to move particle, and if successful update pmap and parts[i].x,y
-int do_move(int i, int x, int y, float nxf, float nyf)
-{
- int nx = (int)(nxf+0.5f), ny = (int)(nyf+0.5f), result;
- if (parts[i].type == PT_NONE)
- return 0;
- result = try_move(i, x, y, nx, ny);
- if (result)
- {
- int t = parts[i].type;
- parts[i].x = nxf;
- parts[i].y = nyf;
- if (ny!=y || nx!=x)
- {
- if ((pmap[y][x]>>8)==i) pmap[y][x] = 0;
- else if ((photons[y][x]>>8)==i) photons[y][x] = 0;
- if (nx<CELL || nx>=XRES-CELL || ny<CELL || ny>=YRES-CELL)//kill_part if particle is out of bounds
- {
- kill_part(i);
- return -1;
- }
- if (ptypes[t].properties & TYPE_ENERGY)
- photons[ny][nx] = t|(i<<8);
- else if (t)
- pmap[ny][nx] = t|(i<<8);
- }
- }
- return result;
-}
-
-
-static unsigned direction_to_map(float dx, float dy, int t)
-{
- // TODO:
- // Adding extra directions causes some inaccuracies.
- // Not adding them causes problems with some diagonal surfaces (photons absorbed instead of reflected).
- // For now, don't add them.
- // Solution may involve more intelligent setting of initial i0 value in find_next_boundary?
- // or rewriting normal/boundary finding code
-
- return (dx >= 0) |
- (((dx + dy) >= 0) << 1) | /* 567 */
- ((dy >= 0) << 2) | /* 4+0 */
- (((dy - dx) >= 0) << 3) | /* 321 */
- ((dx <= 0) << 4) |
- (((dx + dy) <= 0) << 5) |
- ((dy <= 0) << 6) |
- (((dy - dx) <= 0) << 7);
- /*
- return (dx >= -0.001) |
- (((dx + dy) >= -0.001) << 1) | // 567
- ((dy >= -0.001) << 2) | // 4+0
- (((dy - dx) >= -0.001) << 3) | // 321
- ((dx <= 0.001) << 4) |
- (((dx + dy) <= 0.001) << 5) |
- ((dy <= 0.001) << 6) |
- (((dy - dx) <= 0.001) << 7);
- }*/
-}
-
-static int is_blocking(int t, int x, int y)
-{
- if (t & REFRACT) {
- if (x<0 || y<0 || x>=XRES || y>=YRES)
- return 0;
- if ((pmap[y][x] & 0xFF) == PT_GLAS)
- return 1;
- return 0;
- }
-
- return !eval_move(t, x, y, NULL);
-}
-
-static int is_boundary(int pt, int x, int y)
-{
- if (!is_blocking(pt,x,y))
- return 0;
- if (is_blocking(pt,x,y-1) && is_blocking(pt,x,y+1) && is_blocking(pt,x-1,y) && is_blocking(pt,x+1,y))
- return 0;
- return 1;
-}
-
-static int find_next_boundary(int pt, int *x, int *y, int dm, int *em)
-{
- static int dx[8] = {1,1,0,-1,-1,-1,0,1};
- static int dy[8] = {0,1,1,1,0,-1,-1,-1};
- static int de[8] = {0x83,0x07,0x0E,0x1C,0x38,0x70,0xE0,0xC1};
- int i, ii, i0;
-
- if (*x <= 0 || *x >= XRES-1 || *y <= 0 || *y >= YRES-1)
- return 0;
-
- if (*em != -1) {
- i0 = *em;
- dm &= de[i0];
- } else
- i0 = 0;
-
- for (ii=0; ii<8; ii++) {
- i = (ii + i0) & 7;
- if ((dm & (1 << i)) && is_boundary(pt, *x+dx[i], *y+dy[i])) {
- *x += dx[i];
- *y += dy[i];
- *em = i;
- return 1;
- }
- }
-
- return 0;
-}
-
-int get_normal(int pt, int x, int y, float dx, float dy, float *nx, float *ny)
-{
- int ldm, rdm, lm, rm;
- int lx, ly, lv, rx, ry, rv;
- int i, j;
- float r, ex, ey;
-
- if (!dx && !dy)
- return 0;
-
- if (!is_boundary(pt, x, y))
- return 0;
-
- ldm = direction_to_map(-dy, dx, pt);
- rdm = direction_to_map(dy, -dx, pt);
- lx = rx = x;
- ly = ry = y;
- lv = rv = 1;
- lm = rm = -1;
-
- j = 0;
- for (i=0; i<SURF_RANGE; i++) {
- if (lv)
- lv = find_next_boundary(pt, &lx, &ly, ldm, &lm);
- if (rv)
- rv = find_next_boundary(pt, &rx, &ry, rdm, &rm);
- j += lv + rv;
- if (!lv && !rv)
- break;
- }
-
- if (j < NORMAL_MIN_EST)
- return 0;
-
- if ((lx == rx) && (ly == ry))
- return 0;
-
- ex = rx - lx;
- ey = ry - ly;
- r = 1.0f/hypot(ex, ey);
- *nx = ey * r;
- *ny = -ex * r;
-
- return 1;
-}
-
-int get_normal_interp(int pt, float x0, float y0, float dx, float dy, float *nx, float *ny)
-{
- int x, y, i;
-
- dx /= NORMAL_FRAC;
- dy /= NORMAL_FRAC;
-
- for (i=0; i<NORMAL_INTERP; i++) {
- x = (int)(x0 + 0.5f);
- y = (int)(y0 + 0.5f);
- if (is_boundary(pt, x, y))
- break;
- x0 += dx;
- y0 += dy;
- }
- if (i >= NORMAL_INTERP)
- return 0;
-
- if (pt == PT_PHOT)
- photoelectric_effect(x, y);
-
- return get_normal(pt, x, y, dx, dy, nx, ny);
-}
-
-//For soap only
-void detach(int i)
-{
- if ((parts[i].ctype&2) == 2)
- {
- if ((parts[parts[i].tmp].ctype&4) == 4)
- parts[parts[i].tmp].ctype ^= 4;
- }
-
- if ((parts[i].ctype&4) == 4)
- {
- if ((parts[parts[i].tmp2].ctype&2) == 2)
- parts[parts[i].tmp2].ctype ^= 2;
- }
-
- parts[i].ctype = 0;
-}
-
-void kill_part(int i)//kills particle number i
-{
- int x, y;
-
- // Remove from pmap even if type==0, otherwise infinite recursion occurs when flood fill deleting
- // a particle which sets type to 0 without calling kill_part (such as LIFE)
- x = (int)(parts[i].x+0.5f);
- y = (int)(parts[i].y+0.5f);
- if (x>=0 && y>=0 && x<XRES && y<YRES) {
- if ((pmap[y][x]>>8)==i)
- pmap[y][x] = 0;
- else if ((photons[y][x]>>8)==i)
- photons[y][x] = 0;
- }
-
- if (parts[i].type == PT_NONE) //This shouldn't happen anymore, but it's here just in case
- return;
-
- if (parts[i].type == PT_STKM)
- {
- player.spwn = 0;
- }
- if (parts[i].type == PT_STKM2)
- {
- player2.spwn = 0;
- }
- if (parts[i].type == PT_FIGH)
- {
- fighters[(unsigned char)parts[i].tmp].spwn = 0;
- fighcount--;
- }
- if (parts[i].type == PT_SPAWN)
- {
- ISSPAWN1 = 0;
- }
- if (parts[i].type == PT_SPAWN2)
- {
- ISSPAWN2 = 0;
- }
- if (parts[i].type == PT_SOAP)
- {
- detach(i);
- }
-
- parts[i].type = PT_NONE;
- parts[i].life = pfree;
- pfree = i;
-}
-
-#if defined(WIN32) && !defined(__GNUC__)
-_inline void part_change_type(int i, int x, int y, int t)
-#else
-inline void part_change_type(int i, int x, int y, int t)//changes the type of particle number i, to t. This also changes pmap at the same time.
-#endif
-{
- if (x<0 || y<0 || x>=XRES || y>=YRES || i>=NPART || t<0 || t>=PT_NUM)
- return;
- if (!ptypes[t].enabled)
- t = PT_NONE;
-
- if (parts[i].type == PT_STKM)
- player.spwn = 0;
-
- if (parts[i].type == PT_STKM2)
- player2.spwn = 0;
-
- if (parts[i].type == PT_FIGH)
- {
- fighters[(unsigned char)parts[i].tmp].spwn = 0;
- fighcount--;
- }
-
- parts[i].type = t;
- if (ptypes[t].properties & TYPE_ENERGY)
- {
- photons[y][x] = t|(i<<8);
- if ((pmap[y][x]>>8)==i)
- pmap[y][x] = 0;
- }
- else
- {
- pmap[y][x] = t|(i<<8);
- if ((photons[y][x]>>8)==i)
- photons[y][x] = 0;
- }
-}
-
-#if defined(WIN32) && !defined(__GNUC__)
-_inline int create_part(int p, int x, int y, int tv)
-#else
-inline int create_part(int p, int x, int y, int tv)//the function for creating a particle, use p=-1 for creating a new particle, -2 is from a brush, or a particle number to replace a particle.
-#endif
-{
- int i;
-
- int t = tv & 0xFF;
- int v = (tv >> 8) & 0xFF;
-
- if (x<0 || y<0 || x>=XRES || y>=YRES || ((t<=0 || t>=PT_NUM)&&t!=SPC_HEAT&&t!=SPC_COOL&&t!=SPC_AIR&&t!=SPC_VACUUM&&t!=SPC_PGRV&&t!=SPC_NGRV))
- return -1;
- if (t>=0 && t<PT_NUM && !ptypes[t].enabled)
- return -1;
- if(t==SPC_PROP) {
- return -1; //Prop tool works on a mouse click basic, make sure it doesn't do anything here
- }
-
- if (t==SPC_HEAT||t==SPC_COOL)
- {
- int r = pmap[y][x];
- if (!(r&0xFF))
- r = photons[y][x];
- if ((r&0xFF)!=PT_NONE&&(r&0xFF)<PT_NUM)
- {
- if (t==SPC_HEAT&&parts[r>>8].temp<MAX_TEMP)
- {
- float heatchange;
- int fast = ((sdl_mod & (KMOD_SHIFT)) && (sdl_mod & (KMOD_CTRL)));
- if ((r&0xFF)==PT_PUMP || (r&0xFF)==PT_GPMP)
- heatchange = fast?1.0f:.1f;
- else
- heatchange = fast?50.0f:4.0f;
-
- parts[r>>8].temp = restrict_flt(parts[r>>8].temp + heatchange, MIN_TEMP, MAX_TEMP);
- }
- if (t==SPC_COOL&&parts[r>>8].temp>MIN_TEMP)
- {
- float heatchange;
- int fast = ((sdl_mod & (KMOD_SHIFT)) && (sdl_mod & (KMOD_CTRL)));
- if ((r&0xFF)==PT_PUMP || (r&0xFF)==PT_GPMP)
- heatchange = fast?1.0f:.1f;
- else
- heatchange = fast?50.0f:4.0f;
-
- parts[r>>8].temp = restrict_flt(parts[r>>8].temp - heatchange, MIN_TEMP, MAX_TEMP);
- }
- return r>>8;
- }
- else
- {
- return -1;
- }
- }
- if (t==SPC_AIR)
- {
- pv[y/CELL][x/CELL] += 0.10f;
- return -1;
- }
- if (t==SPC_VACUUM)
- {
- pv[y/CELL][x/CELL] -= 0.10f;
- return -1;
- }
- if (t==SPC_PGRV)
- {
- gravmap[(y/CELL)*(XRES/CELL)+(x/CELL)] = 5;
- return -1;
- }
- if (t==SPC_NGRV)
- {
- gravmap[(y/CELL)*(XRES/CELL)+(x/CELL)] = -5;
- return -1;
- }
-
-
- if (t==PT_SPRK)
- {
- if((pmap[y][x]&0xFF)==PT_WIRE){
- parts[pmap[y][x]>>8].ctype=PT_DUST;
- }
- if (!((pmap[y][x]&0xFF)==PT_INST||(ptypes[pmap[y][x]&0xFF].properties&PROP_CONDUCTS)))
- return -1;
- if (parts[pmap[y][x]>>8].life!=0)
- return -1;
- if (p==-2 && (pmap[y][x]&0xFF)==PT_INST)
- {
- flood_INST(x, y, PT_SPRK, PT_INST);
- return pmap[y][x]>>8;
- }
- parts[pmap[y][x]>>8].type = PT_SPRK;
- parts[pmap[y][x]>>8].life = 4;
- parts[pmap[y][x]>>8].ctype = pmap[y][x]&0xFF;
- pmap[y][x] = (pmap[y][x]&~0xFF) | PT_SPRK;
- return pmap[y][x]>>8;
- }
- if (t==PT_SPAWN&&ISSPAWN1)
- return -1;
- if (t==PT_SPAWN2&&ISSPAWN2)
- return -1;
- if (p==-1)//creating from anything but brush
- {
- // If there is a particle, only allow creation if the new particle can occupy the same space as the existing particle
- // If there isn't a particle but there is a wall, check whether the new particle is allowed to be in it
- // (not "!=2" for wall check because eval_move returns 1 for moving into empty space)
- // If there's no particle and no wall, assume creation is allowed
- if (pmap[y][x] ? (eval_move(t, x, y, NULL)!=2) : (bmap[y/CELL][x/CELL] && eval_move(t, x, y, NULL)==0))
- {
- if ((pmap[y][x]&0xFF)!=PT_SPAWN&&(pmap[y][x]&0xFF)!=PT_SPAWN2)
- {
- if (t!=PT_STKM&&t!=PT_STKM2&&t!=PT_FIGH)
- {
- return -1;
- }
- }
- }
- if (pfree == -1)
- return -1;
- i = pfree;
- pfree = parts[i].life;
- }
- else if (p==-2)//creating from brush
- {
- if (pmap[y][x])
- {
- int pmaptype = (pmap[y][x]&0xFF);
- if ((
- (pmaptype==PT_STOR&&!(ptypes[t].properties&TYPE_SOLID))||
- pmaptype==PT_CLNE||
- pmaptype==PT_BCLN||
- pmaptype==PT_CONV||
- (pmaptype==PT_PCLN&&t!=PT_PSCN&&t!=PT_NSCN)||
- (pmaptype==PT_PBCN&&t!=PT_PSCN&&t!=PT_NSCN)
- )&&(
- t!=PT_CLNE&&t!=PT_PCLN&&
- t!=PT_BCLN&&t!=PT_STKM&&
- t!=PT_STKM2&&t!=PT_PBCN&&
- t!=PT_STOR&&t!=PT_FIGH)
- )
- {
- parts[pmap[y][x]>>8].ctype = t;
- if (t==PT_LIFE && v<NGOLALT && (pmap[y][x]&0xFF)!=PT_STOR) parts[pmap[y][x]>>8].tmp = v;
- }
- else if (pmaptype == PT_DTEC && pmaptype != t)
- {
- parts[pmap[y][x]>>8].ctype = t;
- if (t==PT_LIFE && v<NGOLALT) parts[pmap[y][x]>>8].tmp = v;
- }
- return -1;
- }
- if (photons[y][x] && (ptypes[t].properties & TYPE_ENERGY))
- return -1;
- if (pfree == -1)
- return -1;
- i = pfree;
- pfree = parts[i].life;
- }
- else if (p==-3)//skip pmap checks, e.g. for sing explosion
- {
- if (pfree == -1)
- return -1;
- i = pfree;
- pfree = parts[i].life;
- }
- else
- {
- int oldX = (int)(parts[p].x+0.5f);
- int oldY = (int)(parts[p].y+0.5f);
- if ((pmap[oldY][oldX]>>8)==p)
- pmap[oldY][oldX] = 0;
- if ((photons[oldY][oldX]>>8)==p)
- photons[oldY][oldX] = 0;
- i = p;
- }
-
- if (i>parts_lastActiveIndex) parts_lastActiveIndex = i;
-
- parts[i].dcolour = 0;
- parts[i].flags = 0;
- if (t==PT_GLAS)
- {
- parts[i].pavg[1] = pv[y/CELL][x/CELL];
- }
- else if (t==PT_QRTZ)
- {
- parts[i].pavg[1] = pv[y/CELL][x/CELL];
- }
- else
- {
- parts[i].pavg[0] = 0.0f;
- parts[i].pavg[1] = 0.0f;
- }
- if (t!=PT_STKM&&t!=PT_STKM2&&t!=PT_FIGH)//set everything to default values first, except for stickman.
- {
- parts[i].x = (float)x;
- parts[i].y = (float)y;
-#ifdef OGLR
- parts[i].lastX = (float)x;
- parts[i].lastY = (float)y;
-#endif
- parts[i].type = t;
- parts[i].vx = 0;
- parts[i].vy = 0;
- parts[i].life = 0;
- parts[i].ctype = 0;
- parts[i].temp = ptypes[t].heat;
- parts[i].tmp = 0;
- parts[i].tmp2 = 0;
- }
- //now set various properties that we want at spawn.
- /*if (ptypes[t].properties&PROP_LIFE) {
- int r;
- for (r = 0; r<NGOL; r++)
- if (t==goltype[r])
- parts[i].tmp = grule[r+1][9] - 1;
- }*/
- switch (t)
- {
- case PT_SOAP:
- parts[i].tmp = -1;
- parts[i].tmp2 = -1;
- break;
- case PT_ACID: case PT_CAUS:
- parts[i].life = 75;
- break;
- /*Testing
- case PT_WOOD:
- parts[i].life = 150;
- break;
- End Testing*/
- case PT_WARP:
- parts[i].life = rand()%95+70;
- break;
- case PT_FUSE:
- parts[i].life = 50;
- parts[i].tmp = 50;
- break;
- case PT_LIFE:
- if (v<NGOLALT)
- {
- parts[i].tmp = grule[v+1][9] - 1;
- parts[i].ctype = v;
- }
- break;
- case PT_DEUT:
- parts[i].life = 10;
- break;
- case PT_MERC:
- parts[i].tmp = 10;
- break;
- case PT_BRAY:
- parts[i].life = 30;
- break;
- case PT_GPMP: case PT_PUMP:
- parts[i].life = 10;
- break;
- case PT_SING:
- parts[i].life = rand()%50+60;
- break;
- case PT_QRTZ:
- parts[i].tmp = (rand()%11);
- break;
- case PT_PQRT:
- parts[i].tmp = (rand()%11);
- break;
- case PT_CLST:
- parts[i].tmp = (rand()%7);
- break;
- case PT_FSEP:
- parts[i].life = 50;
- break;
- case PT_COAL:
- parts[i].life = 110;
- parts[i].tmp = 50;
- break;
- case PT_IGNT:
- parts[i].life = 3;
- break;
- case PT_FRZW:
- parts[i].life = 100;
- break;
- case PT_PIPE:
- case PT_PPIP:
- parts[i].life = 60;
- break;
- case PT_BCOL:
- parts[i].life = 110;
- break;
- case PT_FIRE:
- parts[i].life = rand()%50+120;
- break;
- case PT_PLSM:
- parts[i].life = rand()%150+50;
- break;
- case PT_HFLM:
- parts[i].life = rand()%150+50;
- break;
- case PT_LAVA:
- parts[i].life = rand()%120+240;
- break;
- case PT_NBLE:
- parts[i].life = 0;
- break;
- case PT_ICEI:
- parts[i].ctype = PT_WATR;
- break;
- case PT_MORT:
- parts[i].vx = 2;
- break;
- case PT_EXOT:
- parts[i].life = 1000;
- parts[i].tmp = 244;
- break;
- case PT_EMBR:
- parts[i].life = 50;
- break;
- case PT_STKM:
- if (player.spwn==0)
- {
- parts[i].x = (float)x;
- parts[i].y = (float)y;
-#ifdef OGLR
- parts[i].lastX = (float)x;
- parts[i].lastY = (float)y;
-#endif
- parts[i].type = PT_STKM;
- parts[i].vx = 0;
- parts[i].vy = 0;
- parts[i].life = 100;
- parts[i].ctype = 0;
- parts[i].temp = ptypes[t].heat;
- STKM_init_legs(&player, i);
- player.spwn = 1;
- }
- else
- {
- return -1;
- }
- create_part(-3,x,y,PT_SPAWN);
- ISSPAWN1 = 1;
- break;
- case PT_STKM2:
- if (player2.spwn==0)
- {
- parts[i].x = (float)x;
- parts[i].y = (float)y;
-#ifdef OGLR
- parts[i].lastX = (float)x;
- parts[i].lastY = (float)y;
-#endif
- parts[i].type = PT_STKM2;
- parts[i].vx = 0;
- parts[i].vy = 0;
- parts[i].life = 100;
- parts[i].ctype = 0;
- parts[i].temp = ptypes[t].heat;
- STKM_init_legs(&player2, i);
- player2.spwn = 1;
- }
- else
- {
- return -1;
- }
- create_part(-3,x,y,PT_SPAWN2);
- ISSPAWN2 = 1;
- break;
- case PT_BIZR: case PT_BIZRG: case PT_BIZRS:
- parts[i].ctype = 0x47FFFF;
- break;
- case PT_DTEC:
- parts[i].tmp2 = 2;
- break;
- default:
- if (t==PT_FIGH)
- {
- unsigned char fcount = 0;
- while (fcount < 100 && fcount < (fighcount+1) && fighters[fcount].spwn==1) fcount++;
- if (fcount < 100 && fighters[fcount].spwn==0)
- {
- parts[i].x = (float)x;
- parts[i].y = (float)y;
-#ifdef OGLR
- parts[i].lastX = (float)x;
- parts[i].lastY = (float)y;
-#endif
- parts[i].type = PT_FIGH;
- parts[i].vx = 0;
- parts[i].vy = 0;
- parts[i].life = 100;
- parts[i].ctype = 0;
- parts[i].tmp = fcount;
- parts[i].temp = ptypes[t].heat;
- STKM_init_legs(&fighters[fcount], i);
- fighters[fcount].spwn = 1;
- fighters[fcount].elem = PT_DUST;
- fighcount++;
-
- return i;
- }
- return -1;
- }
- if (t==PT_PHOT)
- {
- float a = (rand()%8) * 0.78540f;
- parts[i].life = 680;
- parts[i].ctype = 0x3FFFFFFF;
- parts[i].vx = 3.0f*cosf(a);
- parts[i].vy = 3.0f*sinf(a);
- }
- if (t==PT_ELEC)
- {
- float a = (rand()%360)*3.14159f/180.0f;
- parts[i].life = 680;
- parts[i].vx = 2.0f*cosf(a);
- parts[i].vy = 2.0f*sinf(a);
- }
- if (t==PT_NEUT)
- {
- float r = (rand()%128+128)/127.0f;
- float a = (rand()%360)*3.14159f/180.0f;
- parts[i].life = rand()%480+480;
- parts[i].vx = r*cosf(a);
- parts[i].vy = r*sinf(a);
- }
- if (t==PT_TRON)
- {
- int randhue = rand()%360;
- int randomdir = rand()%4;
- parts[i].tmp = 1|(randomdir<<5)|(randhue<<7);//set as a head and a direction
- parts[i].tmp2 = 4;//tail
- parts[i].life = 5;
- }
- if (t==PT_LIGH)
- {
- float gx, gy, gsize;
- if (p!=-2)
- {
- parts[i].life=30;
- parts[i].temp=parts[i].life*150.0f; // temperature of the lighting shows the power of the lighting
- }
- get_gravity_field(x, y, 1.0f, 1.0f, &gx, &gy);
- gsize = gx*gx+gy*gy;
- if (gsize<0.0016f)
- {
- float angle = (rand()%6284)*0.001f;//(in radians, between 0 and 2*pi)
- gsize = sqrtf(gsize);
- // randomness in weak gravity fields (more randomness with weaker fields)
- gx += cosf(angle)*(0.04f-gsize);
- gy += sinf(angle)*(0.04f-gsize);
- }
- parts[i].tmp = (((int)(atan2f(-gy, gx)*(180.0f/M_PI)))+rand()%40-20+360)%360;
- parts[i].tmp2 = 4;
- }
- }
- //and finally set the pmap/photon maps to the newly created particle
- if (ptypes[t].properties & TYPE_ENERGY)
- photons[y][x] = t|(i<<8);
- else if (t!=PT_STKM && t!=PT_STKM2 && t!=PT_FIGH)
- pmap[y][x] = t|(i<<8);
-
- //Fancy dust effects for powder types
- if((ptypes[t].properties & TYPE_PART) && pretty_powder)
- {
- int colr, colg, colb, randa;
- randa = (rand()%30)-15;
- colr = (PIXR(ptypes[t].pcolors)+sandcolour_r+(rand()%20)-10+randa);
- colg = (PIXG(ptypes[t].pcolors)+sandcolour_g+(rand()%20)-10+randa);
- colb = (PIXB(ptypes[t].pcolors)+sandcolour_b+(rand()%20)-10+randa);
- colr = colr>255 ? 255 : (colr<0 ? 0 : colr);
- colg = colg>255 ? 255 : (colg<0 ? 0 : colg);
- colb = colb>255 ? 255 : (colb<0 ? 0 : colb);
- parts[i].dcolour = 0xFF000000 | (colr<<16) | (colg<<8) | colb;
- }
-
- return i;
-}
-
-static void create_gain_photon(int pp)//photons from PHOT going through GLOW
-{
- float xx, yy;
- int i, lr, temp_bin, nx, ny;
-
- if (pfree == -1)
- return;
- i = pfree;
-
- lr = rand() % 2;
-
- if (lr) {
- xx = parts[pp].x - 0.3*parts[pp].vy;
- yy = parts[pp].y + 0.3*parts[pp].vx;
- } else {
- xx = parts[pp].x + 0.3*parts[pp].vy;
- yy = parts[pp].y - 0.3*parts[pp].vx;
- }
-
- nx = (int)(xx + 0.5f);
- ny = (int)(yy + 0.5f);
-
- if (nx<0 || ny<0 || nx>=XRES || ny>=YRES)
- return;
-
- if ((pmap[ny][nx] & 0xFF) != PT_GLOW)
- return;
-
- pfree = parts[i].life;
- if (i>parts_lastActiveIndex) parts_lastActiveIndex = i;
-
- parts[i].type = PT_PHOT;
- parts[i].life = 680;
- parts[i].x = xx;
- parts[i].y = yy;
- parts[i].vx = parts[pp].vx;
- parts[i].vy = parts[pp].vy;
- parts[i].temp = parts[pmap[ny][nx] >> 8].temp;
- parts[i].tmp = 0;
- parts[i].pavg[0] = parts[i].pavg[1] = 0.0f;
- photons[ny][nx] = PT_PHOT|(i<<8);
-
- temp_bin = (int)((parts[i].temp-273.0f)*0.25f);
- if (temp_bin < 0) temp_bin = 0;
- if (temp_bin > 25) temp_bin = 25;
- parts[i].ctype = 0x1F << temp_bin;
-}
-
-static void create_cherenkov_photon(int pp)//photons from NEUT going through GLAS
-{
- int i, lr, nx, ny;
- float r, eff_ior;
-
- if (pfree == -1)
- return;
- i = pfree;
-
- nx = (int)(parts[pp].x + 0.5f);
- ny = (int)(parts[pp].y + 0.5f);
- if ((pmap[ny][nx] & 0xFF) != PT_GLAS)
- return;
-
- if (hypotf(parts[pp].vx, parts[pp].vy) < 1.44f)
- return;
-
- pfree = parts[i].life;
- if (i>parts_lastActiveIndex) parts_lastActiveIndex = i;
-
- lr = rand() % 2;
-
- parts[i].type = PT_PHOT;
- parts[i].ctype = 0x00000F80;
- parts[i].life = 680;
- parts[i].x = parts[pp].x;
- parts[i].y = parts[pp].y;
- parts[i].temp = parts[pmap[ny][nx] >> 8].temp;
- parts[i].tmp = 0;
- parts[i].pavg[0] = parts[i].pavg[1] = 0.0f;
- photons[ny][nx] = PT_PHOT|(i<<8);
-
- if (lr) {
- parts[i].vx = parts[pp].vx - 2.5f*parts[pp].vy;
- parts[i].vy = parts[pp].vy + 2.5f*parts[pp].vx;
- } else {
- parts[i].vx = parts[pp].vx + 2.5f*parts[pp].vy;
- parts[i].vy = parts[pp].vy - 2.5f*parts[pp].vx;
- }
-
- /* photons have speed of light. no discussion. */
- r = 1.269 / hypotf(parts[i].vx, parts[i].vy);
- parts[i].vx *= r;
- parts[i].vy *= r;
-}
-
-#if defined(WIN32) && !defined(__GNUC__)
-_inline void delete_part(int x, int y, int flags)//calls kill_part with the particle located at x,y
-#else
-inline void delete_part(int x, int y, int flags)//calls kill_part with the particle located at x,y
-#endif
-{
- unsigned i;
-
- if (x<0 || y<0 || x>=XRES || y>=YRES)
- return;
- if (photons[y][x]) {
- i = photons[y][x];
- } else {
- i = pmap[y][x];
- }
-
- if (!i)
- return;
- if (!(flags&BRUSH_SPECIFIC_DELETE) || parts[i>>8].type==SLALT || SLALT==0)//specific deletiom
- {
- kill_part(i>>8);
- }
- else if (ptypes[parts[i>>8].type].menusection==SEC)//specific menu deletion
- {
- kill_part(i>>8);
- }
- else
- return;
-}
-
-#if defined(WIN32) && !defined(__GNUC__)
-_inline int is_wire(int x, int y)
-#else
-inline int is_wire(int x, int y)
-#endif
-{
- return bmap[y][x]==WL_DETECT || bmap[y][x]==WL_EWALL || bmap[y][x]==WL_ALLOWLIQUID || bmap[y][x]==WL_WALLELEC || bmap[y][x]==WL_ALLOWALLELEC || bmap[y][x]==WL_EHOLE;
-}
-
-#if defined(WIN32) && !defined(__GNUC__)
-_inline int is_wire_off(int x, int y)
-#else
-inline int is_wire_off(int x, int y)
-#endif
-{
- return (bmap[y][x]==WL_DETECT || bmap[y][x]==WL_EWALL || bmap[y][x]==WL_ALLOWLIQUID || bmap[y][x]==WL_WALLELEC || bmap[y][x]==WL_ALLOWALLELEC || bmap[y][x]==WL_EHOLE) && emap[y][x]<8;
-}
-
-int get_wavelength_bin(int *wm)
-{
- int i, w0=30, wM=0;
-
- if (!*wm)
- return -1;
-
- for (i=0; i<30; i++)
- if (*wm & (1<<i)) {
- if (i < w0)
- w0 = i;
- if (i > wM)
- wM = i;
- }
-
- if (wM-w0 < 5)
- return (wM+w0)/2;
-
- i = rand() % (wM-w0-3);
- i += w0;
-
- *wm &= 0x1F << i;
- return i + 2;
-}
-
-void set_emap(int x, int y)
-{
- int x1, x2;
-
- if (!is_wire_off(x, y))
- return;
-
- // go left as far as possible
- x1 = x2 = x;
- while (x1>0)
- {
- if (!is_wire_off(x1-1, y))
- break;
- x1--;
- }
- while (x2<XRES/CELL-1)
- {
- if (!is_wire_off(x2+1, y))
- break;
- x2++;
- }
-
- // fill span
- for (x=x1; x<=x2; x++)
- emap[y][x] = 16;
-
- // fill children
-
- if (y>1 && x1==x2 &&
- is_wire(x1-1, y-1) && is_wire(x1, y-1) && is_wire(x1+1, y-1) &&
- !is_wire(x1-1, y-2) && is_wire(x1, y-2) && !is_wire(x1+1, y-2))
- set_emap(x1, y-2);
- else if (y>0)
- for (x=x1; x<=x2; x++)
- if (is_wire_off(x, y-1))
- {
- if (x==x1 || x==x2 || y>=YRES/CELL-1 ||
- is_wire(x-1, y-1) || is_wire(x+1, y-1) ||
- is_wire(x-1, y+1) || !is_wire(x, y+1) || is_wire(x+1, y+1))
- set_emap(x, y-1);
- }
-
- if (y<YRES/CELL-2 && x1==x2 &&
- is_wire(x1-1, y+1) && is_wire(x1, y+1) && is_wire(x1+1, y+1) &&
- !is_wire(x1-1, y+2) && is_wire(x1, y+2) && !is_wire(x1+1, y+2))
- set_emap(x1, y+2);
- else if (y<YRES/CELL-1)
- for (x=x1; x<=x2; x++)
- if (is_wire_off(x, y+1))
- {
- if (x==x1 || x==x2 || y<0 ||
- is_wire(x-1, y+1) || is_wire(x+1, y+1) ||
- is_wire(x-1, y-1) || !is_wire(x, y-1) || is_wire(x+1, y-1))
- set_emap(x, y+1);
- }
-}
-
-#if defined(WIN32) && !defined(__GNUC__)
-int parts_avg(int ci, int ni,int t)//t is the particle you are looking for, returns the particle between two particles
-#else
-inline int parts_avg(int ci, int ni,int t)
-#endif
-{
- if (t==PT_INSL)//to keep electronics working
- {
- int pmr = pmap[((int)(parts[ci].y+0.5f) + (int)(parts[ni].y+0.5f))/2][((int)(parts[ci].x+0.5f) + (int)(parts[ni].x+0.5f))/2];
- if (pmr)
- return parts[pmr>>8].type;
- else
- return PT_NONE;
- }
- else
- {
- int pmr2 = pmap[(int)((parts[ci].y + parts[ni].y)/2+0.5f)][(int)((parts[ci].x + parts[ni].x)/2+0.5f)];//seems to be more accurate.
- if (pmr2)
- {
- if (parts[pmr2>>8].type==t)
- return t;
- }
- else
- return PT_NONE;
- }
- return PT_NONE;
-}
-
-
-int nearest_part(int ci, int t, int max_d)
-{
- int distance = (max_d!=-1)?max_d:MAX_DISTANCE;
- int ndistance = 0;
- int id = -1;
- int i = 0;
- int cx = (int)parts[ci].x;
- int cy = (int)parts[ci].y;
- for (i=0; i<=parts_lastActiveIndex; i++)
- {
- if ((parts[i].type==t||(t==-1&&parts[i].type))&&!parts[i].life&&i!=ci)
- {
- ndistance = abs(cx-parts[i].x)+abs(cy-parts[i].y);// Faster but less accurate Older: sqrt(pow(cx-parts[i].x, 2)+pow(cy-parts[i].y, 2));
- if (ndistance<distance)
- {
- distance = ndistance;
- id = i;
- }
- }
- }
- return id;
-}
-
-void create_arc(int sx, int sy, int dx, int dy, int midpoints, int variance, int type, int flags)
-{
- int i;
- float xint, yint;
- int *xmid, *ymid;
- int voffset = variance/2;
- xmid = calloc(midpoints + 2, sizeof(int));
- ymid = calloc(midpoints + 2, sizeof(int));
- xint = (float)(dx-sx)/(float)(midpoints+1.0f);
- yint = (float)(dy-sy)/(float)(midpoints+1.0f);
- xmid[0] = sx;
- xmid[midpoints+1] = dx;
- ymid[0] = sy;
- ymid[midpoints+1] = dy;
-
- for(i = 1; i <= midpoints; i++)
- {
- ymid[i] = ymid[i-1]+yint;
- xmid[i] = xmid[i-1]+xint;
- }
-
- for(i = 0; i <= midpoints; i++)
- {
- if(i!=midpoints)
- {
- xmid[i+1] += (rand()%variance)-voffset;
- ymid[i+1] += (rand()%variance)-voffset;
- }
- create_line(xmid[i], ymid[i], xmid[i+1], ymid[i+1], 0, 0, type, flags);
- }
- free(xmid);
- free(ymid);
-}
-
-//the main function for updating particles
-void update_particles_i(pixel *vid, int start, int inc)
-{
- int i, j, x, y, t, nx, ny, r, surround_space, s, lt, rt, nt, nnx, nny, q, golnum, goldelete, z, neighbors, createdsomething;
- float mv, dx, dy, ix, iy, lx, ly, nrx, nry, dp, ctemph, ctempl, gravtot, gel_scale;
- int fin_x, fin_y, clear_x, clear_y, stagnant;
- float fin_xf, fin_yf, clear_xf, clear_yf;
- float nn, ct1, ct2, swappage;
- float pt = R_TEMP;
- float c_heat = 0.0f;
- int h_count = 0;
- int starti = (start*-1);
- int surround[8];
- int surround_hconduct[8];
- int lighting_ok=1;
- unsigned int elem_properties;
- float pGravX, pGravY, pGravD;
- int excessive_stacking_found = 0;
-
- if (sys_pause&&lighting_recreate>0)
- {
- for (i=0; i<=parts_lastActiveIndex; i++)
- {
- if (parts[i].type==PT_LIGH && parts[i].tmp2>0)
- {
- lighting_ok=0;
- break;
- }
- }
- }
-
- if (lighting_ok)
- lighting_recreate--;
-
- if (lighting_recreate<0)
- lighting_recreate=1;
-
- if (lighting_recreate>21)
- lighting_recreate=21;
-
- if (sys_pause&&!framerender)//do nothing if paused
- return;
-
- if (force_stacking_check || (rand()%10)==0)
- {
- force_stacking_check = 0;
- excessive_stacking_found = 0;
- for (y=0; y<YRES; y++)
- {
- for (x=0; x<XRES; x++)
- {
- // Use a threshold, since some particle stacking can be normal (e.g. BIZR + FILT)
- // Setting pmap_count[y][x] > NPART means BHOL will form in that spot
- if (pmap_count[y][x]>5)
- {
- if (bmap[y/CELL][x/CELL]==WL_EHOLE)
- {
- // Allow more stacking in E-hole
- if (pmap_count[y][x]>1500)
- {
- pmap_count[y][x] = pmap_count[y][x] + NPART;
- excessive_stacking_found = 1;
- }
- }
- // Random chance to turn into BHOL that increases with the amount of stacking, up to a threshold where it is certain to turn into BHOL
- else if (pmap_count[y][x]>1500 || (rand()%1600)<=(pmap_count[y][x]+100))
- {
- pmap_count[y][x] = pmap_count[y][x] + NPART;
- excessive_stacking_found = 1;
- }
- }
- }
- }
- if (excessive_stacking_found)
- {
- for (i=0; i<=parts_lastActiveIndex; i++)
- {
- if (parts[i].type)
- {
- t = parts[i].type;
- x = (int)(parts[i].x+0.5f);
- y = (int)(parts[i].y+0.5f);
- if (x>=0 && y>=0 && x<XRES && y<YRES && !(ptypes[t].properties&TYPE_ENERGY))
- {
- if (pmap_count[y][x]>=NPART)
- {
- if (pmap_count[y][x]>NPART)
- {
- create_part(i, x, y, PT_NBHL);
- parts[i].temp = MAX_TEMP;
- parts[i].tmp = pmap_count[y][x]-NPART;//strength of grav field
- if (parts[i].tmp>51200) parts[i].tmp = 51200;
- pmap_count[y][x] = NPART;
- }
- else
- {
- kill_part(i);
- }
- }
- }
- }
- }
- }
- }
-
- if (ISGRAV==1)//crappy grav color handling, i will change this someday
- {
- ISGRAV = 0;
- GRAV ++;
- GRAV_R = 60;
- GRAV_G = 0;
- GRAV_B = 0;
- GRAV_R2 = 30;
- GRAV_G2 = 30;
- GRAV_B2 = 0;
- for ( q = 0; q <= GRAV; q++)
- {
- if (GRAV_R >0 && GRAV_G==0)
- {
- GRAV_R--;
- GRAV_B++;
- }
- if (GRAV_B >0 && GRAV_R==0)
- {
- GRAV_B--;
- GRAV_G++;
- }
- if (GRAV_G >0 && GRAV_B==0)
- {
- GRAV_G--;
- GRAV_R++;
- }
- if (GRAV_R2 >0 && GRAV_G2==0)
- {
- GRAV_R2--;
- GRAV_B2++;
- }
- if (GRAV_B2 >0 && GRAV_R2==0)
- {
- GRAV_B2--;
- GRAV_G2++;
- }
- if (GRAV_G2 >0 && GRAV_B2==0)
- {
- GRAV_G2--;
- GRAV_R2++;
- }
- }
- if (GRAV>180) GRAV = 0;
-
- }
- if (ISLOVE==1)//LOVE element handling
- {
- ISLOVE = 0;
- for (ny=0; ny<YRES-4; ny++)
- {
- for (nx=0; nx<XRES-4; nx++)
- {
- r=pmap[ny][nx];
- if (!r)
- {
- continue;
- }
- else if ((ny<9||nx<9||ny>YRES-7||nx>XRES-10)&&parts[r>>8].type==PT_LOVE)
- kill_part(r>>8);
- else if (parts[r>>8].type==PT_LOVE)
- {
- love[nx/9][ny/9] = 1;
- }
-
- }
- }
- for (nx=9; nx<=XRES-18; nx++)
- {
- for (ny=9; ny<=YRES-7; ny++)
- {
- if (love[nx/9][ny/9]==1)
- {
- for ( nnx=0; nnx<9; nnx++)
- for ( nny=0; nny<9; nny++)
- {
- if (ny+nny>0&&ny+nny<YRES&&nx+nnx>=0&&nx+nnx<XRES)
- {
- rt=pmap[ny+nny][nx+nnx];
- if (!rt&&loverule[nnx][nny]==1)
- create_part(-1,nx+nnx,ny+nny,PT_LOVE);
- else if (!rt)
- continue;
- else if (parts[rt>>8].type==PT_LOVE&&loverule[nnx][nny]==0)
- kill_part(rt>>8);
- }
- }
- }
- love[nx/9][ny/9]=0;
- }
- }
- }
- if (ISLOLZ==1)//LOLZ element handling
- {
- ISLOLZ = 0;
- for (ny=0; ny<YRES-4; ny++)
- {
- for (nx=0; nx<XRES-4; nx++)
- {
- r=pmap[ny][nx];
- if (!r)
- {
- continue;
- }
- else if ((ny<9||nx<9||ny>YRES-7||nx>XRES-10)&&parts[r>>8].type==PT_LOLZ)
- kill_part(r>>8);
- else if (parts[r>>8].type==PT_LOLZ)
- {
- lolz[nx/9][ny/9] = 1;
- }
-
- }
- }
- for (nx=9; nx<=XRES-18; nx++)
- {
- for (ny=9; ny<=YRES-7; ny++)
- {
- if (lolz[nx/9][ny/9]==1)
- {
- for ( nnx=0; nnx<9; nnx++)
- for ( nny=0; nny<9; nny++)
- {
- if (ny+nny>0&&ny+nny<YRES&&nx+nnx>=0&&nx+nnx<XRES)
- {
- rt=pmap[ny+nny][nx+nnx];
- if (!rt&&lolzrule[nny][nnx]==1)
- create_part(-1,nx+nnx,ny+nny,PT_LOLZ);
- else if (!rt)
- continue;
- else if (parts[rt>>8].type==PT_LOLZ&&lolzrule[nny][nnx]==0)
- kill_part(rt>>8);
-
- }
- }
- }
- lolz[nx/9][ny/9]=0;
- }
- }
- }
- //wire!
- if(wire_placed == 1)
- {
- wire_placed = 0;
- for (nx=0; nx<XRES; nx++)
- {
- for (ny=0; ny<YRES; ny++)
- {
- r = pmap[ny][nx];
- if (!r)
- continue;
- if(parts[r>>8].type==PT_WIRE)
- parts[r>>8].tmp=parts[r>>8].ctype;
- }
- }
- }
-
- if (ppip_changed)
- {
- for (i=0; i<=parts_lastActiveIndex; i++)
- {
- if (parts[i].type==PT_PPIP)
- {
- parts[i].tmp |= (parts[i].tmp&0xE0000000)>>3;
- parts[i].tmp &= ~0xE0000000;
- }
- }
- ppip_changed = 0;
- }
-
- //game of life!
- if (ISGOL==1&&++CGOL>=GSPEED)//GSPEED is frames per generation
- {
- int createdsomething = 0;
- CGOL=0;
- ISGOL=0;
- for (ny=CELL; ny<YRES-CELL; ny++)
- {//go through every particle and set neighbor map
- for (nx=CELL; nx<XRES-CELL; nx++)
- {
- r = pmap[ny][nx];
- if (!r)
- {
- gol[ny][nx] = 0;
- continue;
- }
- else
- {
- //for ( golnum=1; golnum<=NGOL; golnum++) //This shouldn't be necessary any more.
- //{
- if (parts[r>>8].type==PT_LIFE/* && parts[r>>8].ctype==golnum-1*/)
- {
- golnum = parts[r>>8].ctype+1;
- if (golnum<=0 || golnum>NGOLALT) {
- parts[r>>8].type = PT_NONE;
- continue;
- }
- if (parts[r>>8].tmp == grule[golnum][9]-1) {
- gol[ny][nx] = golnum;
- for ( nnx=-1; nnx<2; nnx++)
- {
- for ( nny=-1; nny<2; nny++)//it will count itself as its own neighbor, which is needed, but will have 1 extra for delete check
- {
- rt = pmap[((ny+nny+YRES-3*CELL)%(YRES-2*CELL))+CELL][((nx+nnx+XRES-3*CELL)%(XRES-2*CELL))+CELL];
- if (!rt || (rt&0xFF)==PT_LIFE)
- {
- gol2[((ny+nny+YRES-3*CELL)%(YRES-2*CELL))+CELL][((nx+nnx+XRES-3*CELL)%(XRES-2*CELL))+CELL][golnum] ++;
- gol2[((ny+nny+YRES-3*CELL)%(YRES-2*CELL))+CELL][((nx+nnx+XRES-3*CELL)%(XRES-2*CELL))+CELL][0] ++;
- }
- }
- }
- } else {
- parts[r>>8].tmp --;
- if (parts[r>>8].tmp<=0)
- parts[r>>8].type = PT_NONE;//using kill_part makes it not work
- }
- }
- //}
- }
- }
- }
- for (ny=CELL; ny<YRES-CELL; ny++)
- { //go through every particle again, but check neighbor map, then update particles
- for (nx=CELL; nx<XRES-CELL; nx++)
- {
- r = pmap[ny][nx];
- neighbors = gol2[ny][nx][0];
- if (neighbors==0 || !((r&0xFF)==PT_LIFE || !(r&0xFF)))
- continue;
- for ( golnum = 1; golnum<=NGOL; golnum++)
- {
- goldelete = neighbors;
- if (gol[ny][nx]==0&&grule[golnum][goldelete]>=2&&gol2[ny][nx][golnum]>=(goldelete%2)+goldelete/2)
- {
- if (create_part(-1, nx, ny, PT_LIFE|((golnum-1)<<8)))
- createdsomething = 1;
- }
- else if (gol[ny][nx]==golnum&&(grule[golnum][goldelete-1]==0||grule[golnum][goldelete-1]==2))//subtract 1 because it counted itself
- {
- if (parts[r>>8].tmp==grule[golnum][9]-1)
- parts[r>>8].tmp --;
- }
- if (r && parts[r>>8].tmp<=0)
- parts[r>>8].type = PT_NONE;//using kill_part makes it not work
- }
- for ( z = 0; z<=NGOL; z++)
- gol2[ny][nx][z] = 0;//this improves performance A LOT compared to the memset, i was getting ~23 more fps with this.
- }
- }
- if (createdsomething)
- GENERATION ++;
- //memset(gol2, 0, sizeof(gol2));
- }
- if (ISWIRE>0)//wifi channel reseting
- {
- for ( q = 0; q<(int)(MAX_TEMP-73.15f)/100+2; q++)
- {
- wireless[q][0] = wireless[q][1];
- wireless[q][1] = 0;
- }
- ISWIRE--;
- }
- for (i=0; i<=parts_lastActiveIndex; i++)
- if (parts[i].type)
- {
- t = parts[i].type;
-#ifdef OGLR
- parts[i].lastX = parts[i].x;
- parts[i].lastY = parts[i].y;
-#endif
- if (t<0 || t>=PT_NUM)
- {
- kill_part(i);
- continue;
- }
- elem_properties = ptypes[t].properties;
- if (parts[i].life>0 && (elem_properties&PROP_LIFE_DEC))
- {
- // automatically decrease life
- parts[i].life--;
- if (parts[i].life<=0 && (elem_properties&(PROP_LIFE_KILL_DEC|PROP_LIFE_KILL)))
- {
- // kill on change to no life
- kill_part(i);
- continue;
- }
- }
- else if (parts[i].life<=0 && (elem_properties&PROP_LIFE_KILL))
- {
- // kill if no life
- kill_part(i);
- continue;
- }
- }
- //the main particle loop function, goes over all particles.
- for (i=0; i<=parts_lastActiveIndex; i++)
- if (parts[i].type)
- {
- t = parts[i].type;
- x = (int)(parts[i].x+0.5f);
- y = (int)(parts[i].y+0.5f);
-
- //this kills any particle out of the screen, or in a wall where it isn't supposed to go
- if (x<CELL || y<CELL || x>=XRES-CELL || y>=YRES-CELL ||
- (bmap[y/CELL][x/CELL] &&
- (bmap[y/CELL][x/CELL]==WL_WALL ||
- bmap[y/CELL][x/CELL]==WL_WALLELEC ||
- bmap[y/CELL][x/CELL]==WL_ALLOWAIR ||
- (bmap[y/CELL][x/CELL]==WL_DESTROYALL) ||
- (bmap[y/CELL][x/CELL]==WL_ALLOWLIQUID && ptypes[t].falldown!=2) ||
- (bmap[y/CELL][x/CELL]==WL_ALLOWSOLID && ptypes[t].falldown!=1) ||
- (bmap[y/CELL][x/CELL]==WL_ALLOWGAS && !(ptypes[t].properties&TYPE_GAS)) || //&& ptypes[t].falldown!=0 && parts[i].type!=PT_FIRE && parts[i].type!=PT_SMKE && parts[i].type!=PT_HFLM) ||
- (bmap[y/CELL][x/CELL]==WL_ALLOWENERGY && !(ptypes[t].properties&TYPE_ENERGY)) ||
- (bmap[y/CELL][x/CELL]==WL_DETECT && (t==PT_METL || t==PT_SPRK)) ||
- (bmap[y/CELL][x/CELL]==WL_EWALL && !emap[y/CELL][x/CELL])) && (t!=PT_STKM) && (t!=PT_STKM2) && (t!=PT_FIGH)))
- {
- kill_part(i);
- continue;
- }
- if (bmap[y/CELL][x/CELL]==WL_DETECT && emap[y/CELL][x/CELL]<8)
- set_emap(x/CELL, y/CELL);
-
- //adding to velocity from the particle's velocity
- vx[y/CELL][x/CELL] = vx[y/CELL][x/CELL]*ptypes[t].airloss + ptypes[t].airdrag*parts[i].vx;
- vy[y/CELL][x/CELL] = vy[y/CELL][x/CELL]*ptypes[t].airloss + ptypes[t].airdrag*parts[i].vy;
-
- if (t==PT_GAS||t==PT_NBLE)
- {
- if (pv[y/CELL][x/CELL]<3.5f)
- pv[y/CELL][x/CELL] += ptypes[t].hotair*(3.5f-pv[y/CELL][x/CELL]);
- if (y+CELL<YRES && pv[y/CELL+1][x/CELL]<3.5f)
- pv[y/CELL+1][x/CELL] += ptypes[t].hotair*(3.5f-pv[y/CELL+1][x/CELL]);
- if (x+CELL<XRES)
- {
- if (pv[y/CELL][x/CELL+1]<3.5f)
- pv[y/CELL][x/CELL+1] += ptypes[t].hotair*(3.5f-pv[y/CELL][x/CELL+1]);
- if (y+CELL<YRES && pv[y/CELL+1][x/CELL+1]<3.5f)
- pv[y/CELL+1][x/CELL+1] += ptypes[t].hotair*(3.5f-pv[y/CELL+1][x/CELL+1]);
- }
- }
- else//add the hotair variable to the pressure map, like black hole, or white hole.
- {
- pv[y/CELL][x/CELL] += ptypes[t].hotair;
- if (y+CELL<YRES)
- pv[y/CELL+1][x/CELL] += ptypes[t].hotair;
- if (x+CELL<XRES)
- {
- pv[y/CELL][x/CELL+1] += ptypes[t].hotair;
- if (y+CELL<YRES)
- pv[y/CELL+1][x/CELL+1] += ptypes[t].hotair;
- }
- }
-
- //Gravity mode by Moach
- switch (gravityMode)
- {
- default:
- case 0:
- pGravX = 0.0f;
- pGravY = ptypes[t].gravity;
- break;
- case 1:
- pGravX = pGravY = 0.0f;
- break;
- case 2:
- pGravD = 0.01f - hypotf((x - XCNTR), (y - YCNTR));
- pGravX = ptypes[t].gravity * ((float)(x - XCNTR) / pGravD);
- pGravY = ptypes[t].gravity * ((float)(y - YCNTR) / pGravD);
- }
- //Get some gravity from the gravity map
- if (t==PT_ANAR)
- {
- // perhaps we should have a ptypes variable for this
- pGravX -= gravx[(y/CELL)*(XRES/CELL)+(x/CELL)];
- pGravY -= gravy[(y/CELL)*(XRES/CELL)+(x/CELL)];
- }
- else if(t!=PT_STKM && t!=PT_STKM2 && t!=PT_FIGH && !(ptypes[t].properties & TYPE_SOLID))
- {
- pGravX += gravx[(y/CELL)*(XRES/CELL)+(x/CELL)];
- pGravY += gravy[(y/CELL)*(XRES/CELL)+(x/CELL)];
- }
- //velocity updates for the particle
- if (!(parts[i].flags&FLAG_MOVABLE))
- {
- parts[i].vx *= ptypes[t].loss;
- parts[i].vy *= ptypes[t].loss;
- }
- //particle gets velocity from the vx and vy maps
- parts[i].vx += ptypes[t].advection*vx[y/CELL][x/CELL] + pGravX;
- parts[i].vy += ptypes[t].advection*vy[y/CELL][x/CELL] + pGravY;
-
-
- if (ptypes[t].diffusion)//the random diffusion that gases have
- {
-#ifdef REALISTIC
- //The magic number controlls diffusion speed
- parts[i].vx += 0.05f*sqrtf(parts[i].temp)*ptypes[t].diffusion*(rand()/(0.5f*RAND_MAX)-1.0f);
- parts[i].vy += 0.05f*sqrtf(parts[i].temp)*ptypes[t].diffusion*(rand()/(0.5f*RAND_MAX)-1.0f);
-#else
- parts[i].vx += ptypes[t].diffusion*(rand()/(0.5f*RAND_MAX)-1.0f);
- parts[i].vy += ptypes[t].diffusion*(rand()/(0.5f*RAND_MAX)-1.0f);
-#endif
- }
-
- j = surround_space = nt = 0;//if nt is greater than 1 after this, then there is a particle around the current particle, that is NOT the current particle's type, for water movement.
- for (nx=-1; nx<2; nx++)
- for (ny=-1; ny<2; ny++) {
- if (nx||ny) {
- surround[j] = r = pmap[y+ny][x+nx];
- j++;
- if (!(r&0xFF))
- surround_space++;//there is empty space
- if ((r&0xFF)!=t)
- nt++;//there is nothing or a different particle
- }
- }
-
- gel_scale = 1.0f;
- if (t==PT_GEL)
- gel_scale = parts[i].tmp*2.55f;
-
- if (!legacy_enable)
- {
- if (y-2 >= 0 && y-2 < YRES && (ptypes[t].properties&TYPE_LIQUID) && (t!=PT_GEL || gel_scale>(1+rand()%255))) {//some heat convection for liquids
- r = pmap[y-2][x];
- if (!(!r || parts[i].type != (r&0xFF))) {
- if (parts[i].temp>parts[r>>8].temp) {
- swappage = parts[i].temp;
- parts[i].temp = parts[r>>8].temp;
- parts[r>>8].temp = swappage;
- }
- }
- }
-
- //heat transfer code
- h_count = 0;
-#ifdef REALISTIC
- if (t&&(t!=PT_HSWC||parts[i].life==10)&&(ptypes[t].hconduct*gel_scale))
- {
- float c_Cm = 0.0f;
-#else
- if (t&&(t!=PT_HSWC||parts[i].life==10)&&(ptypes[t].hconduct*gel_scale)>(rand()%250))
- {
- float c_Cm = 0.0f;
-#endif
- if (aheat_enable && !(ptypes[t].properties&PROP_NOAMBHEAT))
- {
-#ifdef REALISTIC
- c_heat = parts[i].temp*96.645/ptypes[t].hconduct*gel_scale*fabs(ptypes[t].weight)
- + hv[y/CELL][x/CELL]*100*(pv[y/CELL][x/CELL]+273.15f)/256;
- c_Cm = 96.645/ptypes[t].hconduct*gel_scale*fabs(ptypes[t].weight)
- + 100*(pv[y/CELL][x/CELL]+273.15f)/256;
- pt = c_heat/c_Cm;
- pt = restrict_flt(pt, -MAX_TEMP+MIN_TEMP, MAX_TEMP-MIN_TEMP);
- parts[i].temp = pt;
- //Pressure increase from heat (temporary)
- pv[y/CELL][x/CELL] += (pt-hv[y/CELL][x/CELL])*0.004;
- hv[y/CELL][x/CELL] = pt;
-#else
- c_heat = (hv[y/CELL][x/CELL]-parts[i].temp)*0.04;
- c_heat = restrict_flt(c_heat, -MAX_TEMP+MIN_TEMP, MAX_TEMP-MIN_TEMP);
- parts[i].temp += c_heat;
- hv[y/CELL][x/CELL] -= c_heat;
-#endif
- }
- c_heat = 0.0f; c_Cm = 0.0f;
- for (j=0; j<8; j++)
- {
- surround_hconduct[j] = i;
- r = surround[j];
- if (!r)
- continue;
- rt = r&0xFF;
-
- if (rt&&ptypes[rt].hconduct&&(rt!=PT_HSWC||parts[r>>8].life==10)
- &&(t!=PT_FILT||(rt!=PT_BRAY&&rt!=PT_BIZR&&rt!=PT_BIZRG))
- &&(rt!=PT_FILT||(t!=PT_BRAY&&t!=PT_PHOT&&t!=PT_BIZR&&t!=PT_BIZRG)))
- {
- surround_hconduct[j] = r>>8;
-#ifdef REALISTIC
- if (rt==PT_GEL)
- gel_scale = parts[r>>8].tmp*2.55f;
- else gel_scale = 1.0f;
-
- c_heat += parts[r>>8].temp*96.645/ptypes[rt].hconduct*gel_scale*fabs(ptypes[rt].weight);
- c_Cm += 96.645/ptypes[rt].hconduct*gel_scale*fabs(ptypes[rt].weight);
-#else
- c_heat += parts[r>>8].temp;
-#endif
- h_count++;
- }
- }
-#ifdef REALISTIC
- if (t==PT_GEL)
- gel_scale = parts[i].tmp*2.55f;
- else gel_scale = 1.0f;
-
- if (t == PT_PHOT)
- pt = (c_heat+parts[i].temp*96.645)/(c_Cm+96.645);
- else
- pt = (c_heat+parts[i].temp*96.645/ptypes[t].hconduct*gel_scale*fabs(ptypes[t].weight))/(c_Cm+96.645/ptypes[t].hconduct*gel_scale*fabs(ptypes[t].weight));
-
- c_heat += parts[i].temp*96.645/ptypes[t].hconduct*gel_scale*fabs(ptypes[t].weight);
- c_Cm += 96.645/ptypes[t].hconduct*gel_scale*fabs(ptypes[t].weight);
- parts[i].temp = restrict_flt(pt, MIN_TEMP, MAX_TEMP);
-#else
- pt = (c_heat+parts[i].temp)/(h_count+1);
- pt = parts[i].temp = restrict_flt(pt, MIN_TEMP, MAX_TEMP);
- for (j=0; j<8; j++)
- {
- parts[surround_hconduct[j]].temp = pt;
- }
-#endif
-
- ctemph = ctempl = pt;
- // change boiling point with pressure
- if ((ptypes[t].state==ST_LIQUID && ptransitions[t].tht>-1 && ptransitions[t].tht<PT_NUM && ptypes[ptransitions[t].tht].state==ST_GAS)
- || t==PT_LNTG || t==PT_SLTW)
- ctemph -= 2.0f*pv[y/CELL][x/CELL];
- else if ((ptypes[t].state==ST_GAS && ptransitions[t].tlt>-1 && ptransitions[t].tlt<PT_NUM && ptypes[ptransitions[t].tlt].state==ST_LIQUID)
- || t==PT_WTRV)
- ctempl -= 2.0f*pv[y/CELL][x/CELL];
- s = 1;
-
- //A fix for ice with ctype = 0
- if ((t==PT_ICEI || t==PT_SNOW) && (parts[i].ctype==0 || parts[i].ctype>=PT_NUM || parts[i].ctype==PT_ICEI || parts[i].ctype==PT_SNOW))
- parts[i].ctype = PT_WATR;
-
- if (ctemph>ptransitions[t].thv&&ptransitions[t].tht>-1) {
- // particle type change due to high temperature
-#ifdef REALISTIC
- float dbt = ctempl - pt;
- if (ptransitions[t].tht!=PT_NUM)
- {
- if (platent[t] <= (c_heat - (ptransitions[t].thv - dbt)*c_Cm))
- {
- pt = (c_heat - platent[t])/c_Cm;
- t = ptransitions[t].tht;
- }
- else
- {
- parts[i].temp = restrict_flt(ptransitions[t].thv - dbt, MIN_TEMP, MAX_TEMP);
- s = 0;
- }
- }
-#else
- if (ptransitions[t].tht!=PT_NUM)
- t = ptransitions[t].tht;
-#endif
- else if (t==PT_ICEI || t==PT_SNOW) {
- if (parts[i].ctype<PT_NUM&&parts[i].ctype!=t) {
- if (ptransitions[parts[i].ctype].tlt==t&&pt<=ptransitions[parts[i].ctype].tlv) s = 0;
- else {
-#ifdef REALISTIC
- //One ice table value for all it's kinds
- if (platent[t] <= (c_heat - (ptransitions[parts[i].ctype].tlv - dbt)*c_Cm))
- {
- pt = (c_heat - platent[t])/c_Cm;
- t = parts[i].ctype;
- parts[i].ctype = PT_NONE;
- parts[i].life = 0;
- }
- else
- {
- parts[i].temp = restrict_flt(ptransitions[parts[i].ctype].tlv - dbt, MIN_TEMP, MAX_TEMP);
- s = 0;
- }
-#else
- t = parts[i].ctype;
- parts[i].ctype = PT_NONE;
- parts[i].life = 0;
-#endif
- }
- }
- else s = 0;
- }
- else if (t==PT_SLTW) {
-#ifdef REALISTIC
- if (platent[t] <= (c_heat - (ptransitions[t].thv - dbt)*c_Cm))
- {
- pt = (c_heat - platent[t])/c_Cm;
-
- if (rand()%4==0) t = PT_SALT;
- else t = PT_WTRV;
- }
- else
- {
- parts[i].temp = restrict_flt(ptransitions[t].thv - dbt, MIN_TEMP, MAX_TEMP);
- s = 0;
- }
-#else
- if (rand()%4==0) t = PT_SALT;
- else t = PT_WTRV;
-#endif
- }
- else s = 0;
- } else if (ctempl<ptransitions[t].tlv&&ptransitions[t].tlt>-1) {
- // particle type change due to low temperature
-#ifdef REALISTIC
- float dbt = ctempl - pt;
- if (ptransitions[t].tlt!=PT_NUM)
- {
- if (platent[ptransitions[t].tlt] >= (c_heat - (ptransitions[t].tlv - dbt)*c_Cm))
- {
- pt = (c_heat + platent[ptransitions[t].tlt])/c_Cm;
- t = ptransitions[t].tlt;
- }
- else
- {
- parts[i].temp = restrict_flt(ptransitions[t].tlv - dbt, MIN_TEMP, MAX_TEMP);
- s = 0;
- }
- }
-#else
- if (ptransitions[t].tlt!=PT_NUM)
- t = ptransitions[t].tlt;
-#endif
- else if (t==PT_WTRV) {
- if (pt<273.0f) t = PT_RIME;
- else t = PT_DSTW;
- }
- else if (t==PT_LAVA) {
- if (parts[i].ctype>0 && parts[i].ctype<PT_NUM && parts[i].ctype!=PT_LAVA) {
- if (parts[i].ctype==PT_THRM&&pt>=ptransitions[PT_BMTL].thv) s = 0;
- else if (ptransitions[parts[i].ctype].tht==PT_LAVA) {
- if (pt>=ptransitions[parts[i].ctype].thv) s = 0;
- }
- else if (pt>=973.0f) s = 0; // freezing point for lava with any other (not listed in ptransitions as turning into lava) ctype
- if (s) {
- t = parts[i].ctype;
- parts[i].ctype = PT_NONE;
- if (t==PT_THRM) {
- parts[i].tmp = 0;
- t = PT_BMTL;
- }
- if (t==PT_PLUT)
- {
- parts[i].tmp = 0;
- t = PT_LAVA;
- }
- }
- }
- else if (pt<973.0f) t = PT_STNE;
- else s = 0;
- }
- else s = 0;
- }
- else s = 0;
-#ifdef REALISTIC
- pt = restrict_flt(pt, MIN_TEMP, MAX_TEMP);
- for (j=0; j<8; j++)
- {
- parts[surround_hconduct[j]].temp = pt;
- }
-#endif
- if (s) { // particle type change occurred
- if (t==PT_ICEI||t==PT_LAVA||t==PT_SNOW)
- parts[i].ctype = parts[i].type;
- if (!(t==PT_ICEI&&parts[i].ctype==PT_FRZW)) parts[i].life = 0;
- if (ptypes[t].state==ST_GAS&&ptypes[parts[i].type].state!=ST_GAS)
- pv[y/CELL][x/CELL] += 0.50f;
- part_change_type(i,x,y,t);
- if (t==PT_FIRE||t==PT_PLSM||t==PT_HFLM)
- parts[i].life = rand()%50+120;
- if (t==PT_LAVA) {
- if (parts[i].ctype==PT_BRMT) parts[i].ctype = PT_BMTL;
- else if (parts[i].ctype==PT_SAND) parts[i].ctype = PT_GLAS;
- else if (parts[i].ctype==PT_BGLA) parts[i].ctype = PT_GLAS;
- else if (parts[i].ctype==PT_PQRT) parts[i].ctype = PT_QRTZ;
- parts[i].life = rand()%120+240;
- }
- if (t==PT_NONE) {
- kill_part(i);
- goto killed;
- }
- }
-
- pt = parts[i].temp = restrict_flt(parts[i].temp, MIN_TEMP, MAX_TEMP);
- if (t==PT_LAVA) {
- parts[i].life = restrict_flt((parts[i].temp-700)/7, 0.0f, 400.0f);
- if (parts[i].ctype==PT_THRM&&parts[i].tmp>0)
- {
- parts[i].tmp--;
- parts[i].temp = 3500;
- }
- if (parts[i].ctype==PT_PLUT&&parts[i].tmp>0)
- {
- parts[i].tmp--;
- parts[i].temp = MAX_TEMP;
- }
- }
- }
- else parts[i].temp = restrict_flt(parts[i].temp, MIN_TEMP, MAX_TEMP);
- }
-
- if (t==PT_LIFE)
- {
- parts[i].temp = restrict_flt(parts[i].temp-50.0f, MIN_TEMP, MAX_TEMP);
- ISGOL=1;//means there is a life particle on screen
- }
- if (t==PT_WIRE)
- {
- wire_placed = 1;
- }
- //spark updates from walls
- if ((ptypes[t].properties&PROP_CONDUCTS) || t==PT_SPRK)
- {
- nx = x % CELL;
- if (nx == 0)
- nx = x/CELL - 1;
- else if (nx == CELL-1)
- nx = x/CELL + 1;
- else
- nx = x/CELL;
- ny = y % CELL;
- if (ny == 0)
- ny = y/CELL - 1;
- else if (ny == CELL-1)
- ny = y/CELL + 1;
- else
- ny = y/CELL;
- if (nx>=0 && ny>=0 && nx<XRES/CELL && ny<YRES/CELL)
- {
- if (t!=PT_SPRK)
- {
- if (emap[ny][nx]==12 && !parts[i].life)
- {
- part_change_type(i,x,y,PT_SPRK);
- parts[i].life = 4;
- parts[i].ctype = t;
- t = PT_SPRK;
- }
- }
- else if (bmap[ny][nx]==WL_DETECT || bmap[ny][nx]==WL_EWALL || bmap[ny][nx]==WL_ALLOWLIQUID || bmap[ny][nx]==WL_WALLELEC || bmap[ny][nx]==WL_ALLOWALLELEC || bmap[ny][nx]==WL_EHOLE)
- set_emap(nx, ny);
- }
- }
-
- //the basic explosion, from the .explosive variable
- if ((ptypes[t].explosive&2) && pv[y/CELL][x/CELL]>2.5f)
- {
- parts[i].life = rand()%80+180;
- parts[i].temp = restrict_flt(ptypes[PT_FIRE].heat + (ptypes[t].flammable/2), MIN_TEMP, MAX_TEMP);
- t = PT_FIRE;
- part_change_type(i,x,y,t);
- pv[y/CELL][x/CELL] += 0.25f * CFDS;
- }
-
-
- s = 1;
- gravtot = fabs(gravy[(y/CELL)*(XRES/CELL)+(x/CELL)])+fabs(gravx[(y/CELL)*(XRES/CELL)+(x/CELL)]);
- if (pv[y/CELL][x/CELL]>ptransitions[t].phv&&ptransitions[t].pht>-1) {
- // particle type change due to high pressure
- if (ptransitions[t].pht!=PT_NUM)
- t = ptransitions[t].pht;
- else if (t==PT_BMTL) {
- if (pv[y/CELL][x/CELL]>2.5f)
- t = PT_BRMT;
- else if (pv[y/CELL][x/CELL]>1.0f && parts[i].tmp==1)
- t = PT_BRMT;
- else s = 0;
- }
- else s = 0;
- } else if (pv[y/CELL][x/CELL]<ptransitions[t].plv&&ptransitions[t].plt>-1) {
- // particle type change due to low pressure
- if (ptransitions[t].plt!=PT_NUM)
- t = ptransitions[t].plt;
- else s = 0;
- } else if (gravtot>(ptransitions[t].phv/4.0f)&&ptransitions[t].pht>-1) {
- // particle type change due to high gravity
- if (ptransitions[t].pht!=PT_NUM)
- t = ptransitions[t].pht;
- else if (t==PT_BMTL) {
- if (gravtot>0.625f)
- t = PT_BRMT;
- else if (gravtot>0.25f && parts[i].tmp==1)
- t = PT_BRMT;
- else s = 0;
- }
- else s = 0;
- } else s = 0;
- if (s) { // particle type change occurred
- parts[i].life = 0;
- part_change_type(i,x,y,t);
- if (t==PT_FIRE)
- parts[i].life = rand()%50+120;
- if (t==PT_NONE) {
- kill_part(i);
- goto killed;
- }
- }
-
- //call the particle update function, if there is one
-#ifdef LUACONSOLE
- if (ptypes[t].update_func && lua_el_mode[t] != 2)
-#else
- if (ptypes[t].update_func)
-#endif
- {
- if ((*(ptypes[t].update_func))(i,x,y,surround_space,nt))
- continue;
- else if (t==PT_WARP)
- {
- // Warp does some movement in its update func, update variables to avoid incorrect data in pmap
- x = (int)(parts[i].x+0.5f);
- y = (int)(parts[i].y+0.5f);
- }
- }
-#ifdef LUACONSOLE
- if(lua_el_mode[t])
- {
- if(luacon_part_update(t,i,x,y,surround_space,nt))
- continue;
- // Need to update variables, in case they've been changed by Lua
- x = (int)(parts[i].x+0.5f);
- y = (int)(parts[i].y+0.5f);
- }
-#endif
- if (legacy_enable)//if heat sim is off
- update_legacy_all(i,x,y,surround_space,nt);
-
-killed:
- if (parts[i].type == PT_NONE)//if its dead, skip to next particle
- continue;
-
- if (!parts[i].vx&&!parts[i].vy)//if its not moving, skip to next particle, movement code it next
- continue;
-
-#if defined(WIN32) && !defined(__GNUC__)
- mv = max(fabsf(parts[i].vx), fabsf(parts[i].vy));
-#else
- mv = fmaxf(fabsf(parts[i].vx), fabsf(parts[i].vy));
-#endif
- if (mv < ISTP)
- {
- clear_x = x;
- clear_y = y;
- clear_xf = parts[i].x;
- clear_yf = parts[i].y;
- fin_xf = clear_xf + parts[i].vx;
- fin_yf = clear_yf + parts[i].vy;
- fin_x = (int)(fin_xf+0.5f);
- fin_y = (int)(fin_yf+0.5f);
- }
- else
- {
- // interpolate to see if there is anything in the way
- dx = parts[i].vx*ISTP/mv;
- dy = parts[i].vy*ISTP/mv;
- fin_xf = parts[i].x;
- fin_yf = parts[i].y;
- while (1)
- {
- mv -= ISTP;
- fin_xf += dx;
- fin_yf += dy;
- fin_x = (int)(fin_xf+0.5f);
- fin_y = (int)(fin_yf+0.5f);
- if (mv <= 0.0f)
- {
- // nothing found
- fin_xf = parts[i].x + parts[i].vx;
- fin_yf = parts[i].y + parts[i].vy;
- fin_x = (int)(fin_xf+0.5f);
- fin_y = (int)(fin_yf+0.5f);
- clear_xf = fin_xf-dx;
- clear_yf = fin_yf-dy;
- clear_x = (int)(clear_xf+0.5f);
- clear_y = (int)(clear_yf+0.5f);
- break;
- }
- if (fin_x<CELL || fin_y<CELL || fin_x>=XRES-CELL || fin_y>=YRES-CELL || pmap[fin_y][fin_x] || (bmap[fin_y/CELL][fin_x/CELL] && (bmap[fin_y/CELL][fin_x/CELL]==WL_DESTROYALL || !eval_move(t,fin_x,fin_y,NULL))))
- {
- // found an obstacle
- clear_xf = fin_xf-dx;
- clear_yf = fin_yf-dy;
- clear_x = (int)(clear_xf+0.5f);
- clear_y = (int)(clear_yf+0.5f);
- break;
- }
- if (bmap[fin_y/CELL][fin_x/CELL]==WL_DETECT && emap[fin_y/CELL][fin_x/CELL]<8)
- set_emap(fin_x/CELL, fin_y/CELL);
- }
- }
-
- stagnant = parts[i].flags & FLAG_STAGNANT;
- parts[i].flags &= ~FLAG_STAGNANT;
-
- if (t==PT_STKM || t==PT_STKM2 || t==PT_FIGH)
- {
- int nx, ny;
- //head movement, let head pass through anything
- parts[i].x += parts[i].vx;
- parts[i].y += parts[i].vy;
- nx = (int)((float)parts[i].x+0.5f);
- ny = (int)((float)parts[i].y+0.5f);
- if (ny!=y || nx!=x)
- {
- if ((pmap[y][x]>>8)==i) pmap[y][x] = 0;
- else if ((photons[y][x]>>8)==i) photons[y][x] = 0;
- if (nx<CELL || nx>=XRES-CELL || ny<CELL || ny>=YRES-CELL)
- {
- kill_part(i);
- continue;
- }
- if (ptypes[t].properties & TYPE_ENERGY)
- photons[ny][nx] = t|(i<<8);
- else if (t)
- pmap[ny][nx] = t|(i<<8);
- }
- }
- else if (ptypes[t].properties & TYPE_ENERGY)
- {
- if (t == PT_PHOT) {
- if (parts[i].flags&FLAG_SKIPMOVE)
- {
- parts[i].flags &= ~FLAG_SKIPMOVE;
- continue;
- }
-
- rt = pmap[fin_y][fin_x] & 0xFF;
- lt = pmap[y][x] & 0xFF;
-
- r = eval_move(PT_PHOT, fin_x, fin_y, NULL);
- if (((rt==PT_GLAS && lt!=PT_GLAS) || (rt!=PT_GLAS && lt==PT_GLAS)) && r) {
- if (!get_normal_interp(REFRACT|t, parts[i].x, parts[i].y, parts[i].vx, parts[i].vy, &nrx, &nry)) {
- kill_part(i);
- continue;
- }
-
- r = get_wavelength_bin(&parts[i].ctype);
- if (r == -1) {
- kill_part(i);
- continue;
- }
- nn = GLASS_IOR - GLASS_DISP*(r-15)/15.0f;
- nn *= nn;
- nrx = -nrx;
- nry = -nry;
- if (rt==PT_GLAS && lt!=PT_GLAS)
- nn = 1.0f/nn;
- ct1 = parts[i].vx*nrx + parts[i].vy*nry;
- ct2 = 1.0f - (nn*nn)*(1.0f-(ct1*ct1));
- if (ct2 < 0.0f) {
- // total internal reflection
- parts[i].vx -= 2.0f*ct1*nrx;
- parts[i].vy -= 2.0f*ct1*nry;
- fin_xf = parts[i].x;
- fin_yf = parts[i].y;
- fin_x = x;
- fin_y = y;
- } else {
- // refraction
- ct2 = sqrtf(ct2);
- ct2 = ct2 - nn*ct1;
- parts[i].vx = nn*parts[i].vx + ct2*nrx;
- parts[i].vy = nn*parts[i].vy + ct2*nry;
- }
- }
- }
- if (stagnant)//FLAG_STAGNANT set, was reflected on previous frame
- {
- // cast coords as int then back to float for compatibility with existing saves
- if (!do_move(i, x, y, (float)fin_x, (float)fin_y) && parts[i].type) {
- kill_part(i);
- continue;
- }
- }
- else if (!do_move(i, x, y, fin_xf, fin_yf))
- {
- if (parts[i].type == PT_NONE)
- continue;
- // reflection
- parts[i].flags |= FLAG_STAGNANT;
- if (t==PT_NEUT && 100>(rand()%1000))
- {
- kill_part(i);
- continue;
- }
- r = pmap[fin_y][fin_x];
-
- if (((r&0xFF)==PT_PIPE || (r&0xFF) == PT_PPIP) && !(parts[r>>8].tmp&0xFF))
- {
- parts[r>>8].tmp = (parts[r>>8].tmp&~0xFF) | parts[i].type;
- parts[r>>8].temp = parts[i].temp;
- parts[r>>8].tmp2 = parts[i].life;
- parts[r>>8].pavg[0] = parts[i].tmp;
- parts[r>>8].pavg[1] = parts[i].ctype;
- kill_part(i);
- continue;
- }
-
- // this should be replaced with a particle type attribute ("photwl" or something)
- if ((r & 0xFF) == PT_PSCN) parts[i].ctype = 0x00000000;
- if ((r & 0xFF) == PT_NSCN) parts[i].ctype = 0x00000000;
- if ((r & 0xFF) == PT_SPRK) parts[i].ctype = 0x00000000;
- if ((r & 0xFF) == PT_COAL) parts[i].ctype = 0x00000000;
- if ((r & 0xFF) == PT_BCOL) parts[i].ctype = 0x00000000;
- if ((r & 0xFF) == PT_PLEX) parts[i].ctype &= 0x1F00003E;
- if ((r & 0xFF) == PT_NITR) parts[i].ctype &= 0x0007C000;
- if ((r & 0xFF) == PT_NBLE) parts[i].ctype &= 0x3FFF8000;
- if ((r & 0xFF) == PT_LAVA) parts[i].ctype &= 0x3FF00000;
- if ((r & 0xFF) == PT_ACID) parts[i].ctype &= 0x1FE001FE;
- if ((r & 0xFF) == PT_DUST) parts[i].ctype &= 0x3FFFFFC0;
- if ((r & 0xFF) == PT_SNOW) parts[i].ctype &= 0x03FFFFFF;
- if ((r & 0xFF) == PT_GOO) parts[i].ctype &= 0x3FFAAA00;
- if ((r & 0xFF) == PT_PLNT) parts[i].ctype &= 0x0007C000;
- if ((r & 0xFF) == PT_PLUT) parts[i].ctype &= 0x001FCE00;
- if ((r & 0xFF) == PT_URAN) parts[i].ctype &= 0x003FC000;
-
- if (get_normal_interp(t, parts[i].x, parts[i].y, parts[i].vx, parts[i].vy, &nrx, &nry)) {
- dp = nrx*parts[i].vx + nry*parts[i].vy;
- parts[i].vx -= 2.0f*dp*nrx;
- parts[i].vy -= 2.0f*dp*nry;
- // leave the actual movement until next frame so that reflection of fast particles and refraction happen correctly
- } else {
- if (t!=PT_NEUT)
- kill_part(i);
- continue;
- }
- if (!(parts[i].ctype&0x3FFFFFFF)&&t!=PT_NEUT&&t!=PT_ELEC) {
- kill_part(i);
- continue;
- }
- }
- }
- else if (ptypes[t].falldown==0)
- {
- // gases and solids (but not powders)
- if (!do_move(i, x, y, fin_xf, fin_yf))
- {
- if (parts[i].type == PT_NONE)
- continue;
- // can't move there, so bounce off
- // TODO
- if (fin_x>x+ISTP) fin_x=x+ISTP;
- if (fin_x<x-ISTP) fin_x=x-ISTP;
- if (fin_y>y+ISTP) fin_y=y+ISTP;
- if (fin_y<y-ISTP) fin_y=y-ISTP;
- if (do_move(i, x, y, 0.25f+(float)(2*x-fin_x), 0.25f+fin_y))
- {
- parts[i].vx *= ptypes[t].collision;
- }
- else if (do_move(i, x, y, 0.25f+fin_x, 0.25f+(float)(2*y-fin_y)))
- {
- parts[i].vy *= ptypes[t].collision;
- }
- else
- {
- parts[i].vx *= ptypes[t].collision;
- parts[i].vy *= ptypes[t].collision;
- }
- }
- }
- else
- {
- if (water_equal_test && ptypes[t].falldown == 2 && 1>= rand()%400)//checking stagnant is cool, but then it doesn't update when you change it later.
- {
- if (!flood_water(x,y,i,y, parts[i].tmp2))
- goto movedone;
- }
- // liquids and powders
- if (!do_move(i, x, y, fin_xf, fin_yf))
- {
- if (parts[i].type == PT_NONE)
- continue;
- if (fin_x!=x && do_move(i, x, y, fin_xf, clear_yf))
- {
- parts[i].vx *= ptypes[t].collision;
- parts[i].vy *= ptypes[t].collision;
- }
- else if (fin_y!=y && do_move(i, x, y, clear_xf, fin_yf))
- {
- parts[i].vx *= ptypes[t].collision;
- parts[i].vy *= ptypes[t].collision;
- }
- else
- {
- s = 1;
- r = (rand()%2)*2-1;
- if ((clear_x!=x || clear_y!=y || nt || surround_space) &&
- (fabsf(parts[i].vx)>0.01f || fabsf(parts[i].vy)>0.01f))
- {
- // allow diagonal movement if target position is blocked
- // but no point trying this if particle is stuck in a block of identical particles
- dx = parts[i].vx - parts[i].vy*r;
- dy = parts[i].vy + parts[i].vx*r;
- if (fabsf(dy)>fabsf(dx))
- mv = fabsf(dy);
- else
- mv = fabsf(dx);
- dx /= mv;
- dy /= mv;
- if (do_move(i, x, y, clear_xf+dx, clear_yf+dy))
- {
- parts[i].vx *= ptypes[t].collision;
- parts[i].vy *= ptypes[t].collision;
- goto movedone;
- }
- swappage = dx;
- dx = dy*r;
- dy = -swappage*r;
- if (do_move(i, x, y, clear_xf+dx, clear_yf+dy))
- {
- parts[i].vx *= ptypes[t].collision;
- parts[i].vy *= ptypes[t].collision;
- goto movedone;
- }
- }
- if (ptypes[t].falldown>1 && !ngrav_enable && gravityMode==0 && parts[i].vy>fabsf(parts[i].vx))
- {
- s = 0;
- // stagnant is true if FLAG_STAGNANT was set for this particle in previous frame
- if (!stagnant || nt) //nt is if there is an something else besides the current particle type, around the particle
- rt = 30;//slight less water lag, although it changes how it moves a lot
- else
- rt = 10;
-
- if (t==PT_GEL)
- rt = parts[i].tmp*0.20f+5.0f;
-
- for (j=clear_x+r; j>=0 && j>=clear_x-rt && j<clear_x+rt && j<XRES; j+=r)
- {
- if (((pmap[fin_y][j]&0xFF)!=t || bmap[fin_y/CELL][j/CELL])
- && (s=do_move(i, x, y, (float)j, fin_yf)))
- {
- nx = (int)(parts[i].x+0.5f);
- ny = (int)(parts[i].y+0.5f);
- break;
- }
- if (fin_y!=clear_y && ((pmap[clear_y][j]&0xFF)!=t || bmap[clear_y/CELL][j/CELL])
- && (s=do_move(i, x, y, (float)j, clear_yf)))
- {
- nx = (int)(parts[i].x+0.5f);
- ny = (int)(parts[i].y+0.5f);
- break;
- }
- if ((pmap[clear_y][j]&0xFF)!=t || (bmap[clear_y/CELL][j/CELL] && bmap[clear_y/CELL][j/CELL]!=WL_STREAM))
- break;
- }
- if (parts[i].vy>0)
- r = 1;
- else
- r = -1;
- if (s==1)
- for (j=ny+r; j>=0 && j<YRES && j>=ny-rt && j<ny+rt; j+=r)
- {
- if (((pmap[j][nx]&0xFF)!=t || bmap[j/CELL][nx/CELL]) && do_move(i, nx, ny, (float)nx, (float)j))
- break;
- if ((pmap[j][nx]&255)!=t || (bmap[j/CELL][nx/CELL] && bmap[j/CELL][nx/CELL]!=WL_STREAM))
- break;
- }
- else if (s==-1) {} // particle is out of bounds
- else if ((clear_x!=x||clear_y!=y) && do_move(i, x, y, clear_xf, clear_yf)) {}
- else parts[i].flags |= FLAG_STAGNANT;
- parts[i].vx *= ptypes[t].collision;
- parts[i].vy *= ptypes[t].collision;
- }
- else if (ptypes[t].falldown>1 && fabsf(pGravX*parts[i].vx+pGravY*parts[i].vy)>fabsf(pGravY*parts[i].vx-pGravX*parts[i].vy))
- {
- float nxf, nyf, prev_pGravX, prev_pGravY, ptGrav = ptypes[t].gravity;
- s = 0;
- // stagnant is true if FLAG_STAGNANT was set for this particle in previous frame
- if (!stagnant || nt) //nt is if there is an something else besides the current particle type, around the particle
- rt = 30;//slight less water lag, although it changes how it moves a lot
- else
- rt = 10;
- nxf = clear_xf;
- nyf = clear_yf;
- for (j=0;j<rt;j++)
- {
- switch (gravityMode)
- {
- default:
- case 0:
- pGravX = 0.0f;
- pGravY = ptGrav;
- break;
- case 1:
- pGravX = pGravY = 0.0f;
- break;
- case 2:
- pGravD = 0.01f - hypotf((nx - XCNTR), (ny - YCNTR));
- pGravX = ptGrav * ((float)(nx - XCNTR) / pGravD);
- pGravY = ptGrav * ((float)(ny - YCNTR) / pGravD);
- }
- pGravX += gravx[(ny/CELL)*(XRES/CELL)+(nx/CELL)];
- pGravY += gravy[(ny/CELL)*(XRES/CELL)+(nx/CELL)];
- if (fabsf(pGravY)>fabsf(pGravX))
- mv = fabsf(pGravY);
- else
- mv = fabsf(pGravX);
- if (mv<0.0001f) break;
- pGravX /= mv;
- pGravY /= mv;
- if (j)
- {
- nxf += r*(pGravY*2.0f-prev_pGravY);
- nyf += -r*(pGravX*2.0f-prev_pGravX);
- }
- else
- {
- nxf += r*pGravY;
- nyf += -r*pGravX;
- }
- prev_pGravX = pGravX;
- prev_pGravY = pGravY;
- nx = (int)(nxf+0.5f);
- ny = (int)(nyf+0.5f);
- if (nx<0 || ny<0 || nx>=XRES || ny >=YRES)
- break;
- if ((pmap[ny][nx]&0xFF)!=t || bmap[ny/CELL][nx/CELL])
- {
- s = do_move(i, x, y, nxf, nyf);
- if (s)
- {
- nx = (int)(parts[i].x+0.5f);
- ny = (int)(parts[i].y+0.5f);
- break;
- }
- if (bmap[ny/CELL][nx/CELL]!=WL_STREAM)
- break;
- }
- }
- if (s==1)
- {
- clear_x = nx;
- clear_y = ny;
- for (j=0;j<rt;j++)
- {
- switch (gravityMode)
- {
- default:
- case 0:
- pGravX = 0.0f;
- pGravY = ptGrav;
- break;
- case 1:
- pGravX = pGravY = 0.0f;
- break;
- case 2:
- pGravD = 0.01f - hypotf((nx - XCNTR), (ny - YCNTR));
- pGravX = ptGrav * ((float)(nx - XCNTR) / pGravD);
- pGravY = ptGrav * ((float)(ny - YCNTR) / pGravD);
- }
- pGravX += gravx[(ny/CELL)*(XRES/CELL)+(nx/CELL)];
- pGravY += gravy[(ny/CELL)*(XRES/CELL)+(nx/CELL)];
- if (fabsf(pGravY)>fabsf(pGravX))
- mv = fabsf(pGravY);
- else
- mv = fabsf(pGravX);
- if (mv<0.0001f) break;
- pGravX /= mv;
- pGravY /= mv;
- nxf += pGravX;
- nyf += pGravY;
- nx = (int)(nxf+0.5f);
- ny = (int)(nyf+0.5f);
- if (nx<0 || ny<0 || nx>=XRES || ny>=YRES)
- break;
- if ((pmap[ny][nx]&0xFF)!=t || bmap[ny/CELL][nx/CELL])
- {
- s = do_move(i, clear_x, clear_y, nxf, nyf);
- if (s || bmap[ny/CELL][nx/CELL]!=WL_STREAM)
- break;
- }
- }
- }
- else if (s==-1) {} // particle is out of bounds
- else if ((clear_x!=x||clear_y!=y) && do_move(i, x, y, clear_xf, clear_yf)) {}
- else parts[i].flags |= FLAG_STAGNANT;
- parts[i].vx *= ptypes[t].collision;
- parts[i].vy *= ptypes[t].collision;
- }
- else
- {
- // if interpolation was done, try moving to last clear position
- if ((clear_x!=x||clear_y!=y) && do_move(i, x, y, clear_xf, clear_yf)) {}
- else parts[i].flags |= FLAG_STAGNANT;
- parts[i].vx *= ptypes[t].collision;
- parts[i].vy *= ptypes[t].collision;
- }
- }
- }
- }
-movedone:
- continue;
- }
-}
-
-int parts_lastActiveIndex = NPART-1;
-void update_particles(pixel *vid)//doesn't update the particles themselves, but some other things
-{
- int i, x, y, t;
- int lastPartUsed = 0;
- int lastPartUnused = -1;
-#ifdef MT
- int pt = 0, pc = 0;
- pthread_t *InterThreads;
-#endif
-
- memset(pmap, 0, sizeof(pmap));
- memset(pmap_count, 0, sizeof(pmap_count));
- memset(photons, 0, sizeof(photons));
- NUM_PARTS = 0;
- for (i=0; i<=parts_lastActiveIndex; i++)//the particle loop that resets the pmap/photon maps every frame, to update them.
- {
- if (parts[i].type)
- {
- t = parts[i].type;
- x = (int)(parts[i].x+0.5f);
- y = (int)(parts[i].y+0.5f);
- if (x>=0 && y>=0 && x<XRES && y<YRES)
- {
- if (ptypes[t].properties & TYPE_ENERGY)
- photons[y][x] = t|(i<<8);
- else
- {
- // Particles are sometimes allowed to go inside INVS and FILT
- // To make particles collide correctly when inside these elements, these elements must not overwrite an existing pmap entry from particles inside them
- if (!pmap[y][x] || (t!=PT_INVIS && t!= PT_FILT))
- pmap[y][x] = t|(i<<8);
- // Count number of particles at each location, for excess stacking check
- // (there are a few exceptions, including energy particles - currently no limit on stacking those)
- if (t!=PT_THDR && t!=PT_EMBR && t!=PT_FIGH && t!=PT_PLSM)
- pmap_count[y][x]++;
- }
- }
- lastPartUsed = i;
- NUM_PARTS ++;
- }
- else
- {
- if (lastPartUnused<0) pfree = i;
- else parts[lastPartUnused].life = i;
- lastPartUnused = i;
- }
- }
- if (lastPartUnused==-1)
- {
- if (parts_lastActiveIndex>=NPART-1) pfree = -1;
- else pfree = parts_lastActiveIndex+1;
- }
- else
- {
- if (parts_lastActiveIndex>=NPART-1) parts[lastPartUnused].life = -1;
- else parts[lastPartUnused].life = parts_lastActiveIndex+1;
- }
- parts_lastActiveIndex = lastPartUsed;
- if (!sys_pause||framerender)
- {
- for (y=0; y<YRES/CELL; y++)
- {
- for (x=0; x<XRES/CELL; x++)
- {
- if (emap[y][x])
- emap[y][x] --;
- bmap_blockair[y][x] = (bmap[y][x]==WL_WALL || bmap[y][x]==WL_WALLELEC || (bmap[y][x]==WL_EWALL && !emap[y][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]));
- }
- }
- }
-
- update_particles_i(vid, 0, 1);
-}
-
-void clear_area(int area_x, int area_y, int area_w, int area_h)
-{
- int cx = 0;
- int cy = 0;
- int i;
- for (cy=0; cy<area_h; cy++)
- {
- for (cx=0; cx<area_w; cx++)
- {
- bmap[(cy+area_y)/CELL][(cx+area_x)/CELL] = 0;
- delete_part(cx+area_x, cy+area_y, 0);
- }
- }
- for (i=0; i<MAXSIGNS; i++)
- {
- if (signs[i].x>=area_x && signs[i].x<area_x+area_w && signs[i].y>=area_y && signs[i].y<area_y+area_h)
- {
- signs[i].text[0] = 0;
- }
- }
-}
-
-void create_box(int x1, int y1, int x2, int y2, int c, int flags)
-{
- int i, j;
- if (c==SPC_PROP)
- return;
- if (x1>x2)
- {
- i = x2;
- x2 = x1;
- x1 = i;
- }
- if (y1>y2)
- {
- j = y2;
- y2 = y1;
- y1 = j;
- }
- for (j=y1; j<=y2; j++)
- for (i=x1; i<=x2; i++)
- create_parts(i, j, 0, 0, c, flags, 1);
-}
-
-int flood_prop_2(int x, int y, size_t propoffset, void * propvalue, int proptype, int parttype, char * bitmap)
-{
- int x1, x2, i, dy = 1;
- x1 = x2 = x;
- while (x1>=CELL)
- {
- if ((pmap[y][x1-1]&0xFF)!=parttype || bitmap[(y*XRES)+x1-1])
- {
- break;
- }
- x1--;
- }
- while (x2<XRES-CELL)
- {
- if ((pmap[y][x2+1]&0xFF)!=parttype || bitmap[(y*XRES)+x2+1])
- {
- break;
- }
- x2++;
- }
- for (x=x1; x<=x2; x++)
- {
- i = pmap[y][x]>>8;
- if(proptype==2){
- *((float*)(((char*)&parts[i])+propoffset)) = *((float*)propvalue);
- } else if(proptype==0) {
- *((int*)(((char*)&parts[i])+propoffset)) = *((int*)propvalue);
- } else if(proptype==1) {
- *((char*)(((char*)&parts[i])+propoffset)) = *((char*)propvalue);
- }
- bitmap[(y*XRES)+x] = 1;
- }
- if (y>=CELL+dy)
- for (x=x1; x<=x2; x++)
- if ((pmap[y-dy][x]&0xFF)==parttype && !bitmap[((y-dy)*XRES)+x])
- if (!flood_prop_2(x, y-dy, propoffset, propvalue, proptype, parttype, bitmap))
- return 0;
- if (y<YRES-CELL-dy)
- for (x=x1; x<=x2; x++)
- if ((pmap[y+dy][x]&0xFF)==parttype && !bitmap[((y+dy)*XRES)+x])
- if (!flood_prop_2(x, y+dy, propoffset, propvalue, proptype, parttype, bitmap))
- return 0;
- return 1;
-}
-
-int flood_prop(int x, int y, size_t propoffset, void * propvalue, int proptype)
-{
- int r = 0;
- char * bitmap = malloc(XRES*YRES); //Bitmap for checking
- memset(bitmap, 0, XRES*YRES);
- r = pmap[y][x];
- flood_prop_2(x, y, propoffset, propvalue, proptype, r&0xFF, bitmap);
- free(bitmap);
-}
-
-#define PMAP_CMP_CONDUCTIVE(pmap, t) (((pmap)&0xFF)==(t) || (((pmap)&0xFF)==PT_SPRK && parts[(pmap)>>8].ctype==(t)))
-
-int flood_INST(int x, int y, int fullc, int cm)
-{
- int c = fullc&0xFF;
- int x1, x2, dy = (c<PT_NUM)?1:CELL;
- int co = c;
- int coord_stack_limit = XRES*YRES;
- unsigned short (*coord_stack)[2];
- int coord_stack_size = 0;
- int created_something = 0;
-
- if (c>=PT_NUM)
- return 0;
-
- if (cm==-1)
- {
- if (c==0)
- {
- cm = pmap[y][x]&0xFF;
- if (!cm)
- return 0;
- }
- else
- cm = 0;
- }
-
- if ((pmap[y][x]&0xFF)!=cm || parts[pmap[y][x]>>8].life!=0)
- return 1;
-
- coord_stack = malloc(sizeof(unsigned short)*2*coord_stack_limit);
- coord_stack[coord_stack_size][0] = x;
- coord_stack[coord_stack_size][1] = y;
- coord_stack_size++;
-
- do
- {
- coord_stack_size--;
- x = coord_stack[coord_stack_size][0];
- y = coord_stack[coord_stack_size][1];
- x1 = x2 = x;
- // go left as far as possible
- while (x1>=CELL)
- {
- if ((pmap[y][x1-1]&0xFF)!=cm || parts[pmap[y][x1-1]>>8].life!=0)
- {
- break;
- }
- x1--;
- }
- // go right as far as possible
- while (x2<XRES-CELL)
- {
- if ((pmap[y][x2+1]&0xFF)!=cm || parts[pmap[y][x2+1]>>8].life!=0)
- {
- break;
- }
- x2++;
- }
- // fill span
- for (x=x1; x<=x2; x++)
- {
- if (create_part(-1, x, y, fullc)>=0)
- created_something = 1;
- }
-
- // add vertically adjacent pixels to stack
- // (wire crossing for INST)
- if (y>=CELL+1 && x1==x2 &&
- PMAP_CMP_CONDUCTIVE(pmap[y-1][x1-1], cm) && PMAP_CMP_CONDUCTIVE(pmap[y-1][x1], cm) && PMAP_CMP_CONDUCTIVE(pmap[y-1][x1+1], cm) &&
- !PMAP_CMP_CONDUCTIVE(pmap[y-2][x1-1], cm) && PMAP_CMP_CONDUCTIVE(pmap[y-2][x1], cm) && !PMAP_CMP_CONDUCTIVE(pmap[y-2][x1+1], cm))
- {
- // travelling vertically up, skipping a horizontal line
- if ((pmap[y-2][x1]&0xFF)==cm && !parts[pmap[y-2][x1]>>8].life)
- {
- coord_stack[coord_stack_size][0] = x1;
- coord_stack[coord_stack_size][1] = y-2;
- coord_stack_size++;
- if (coord_stack_size>=coord_stack_limit)
- {
- free(coord_stack);
- return -1;
- }
- }
- }
- else if (y>=CELL+1)
- {
- for (x=x1; x<=x2; x++)
- {
- if ((pmap[y-1][x]&0xFF)==cm && !parts[pmap[y-1][x]>>8].life)
- {
- if (x==x1 || x==x2 || y>=YRES-CELL-1 || !PMAP_CMP_CONDUCTIVE(pmap[y+1][x], cm))
- {
- // if at the end of a horizontal section, or if it's a T junction
- coord_stack[coord_stack_size][0] = x;
- coord_stack[coord_stack_size][1] = y-1;
- coord_stack_size++;
- if (coord_stack_size>=coord_stack_limit)
- {
- free(coord_stack);
- return -1;
- }
- }
- }
- }
- }
-
- if (y<YRES-CELL-1 && x1==x2 &&
- PMAP_CMP_CONDUCTIVE(pmap[y+1][x1-1], cm) && PMAP_CMP_CONDUCTIVE(pmap[y+1][x1], cm) && PMAP_CMP_CONDUCTIVE(pmap[y+1][x1+1], cm) &&
- !PMAP_CMP_CONDUCTIVE(pmap[y+2][x1-1], cm) && PMAP_CMP_CONDUCTIVE(pmap[y+2][x1], cm) && !PMAP_CMP_CONDUCTIVE(pmap[y+2][x1+1], cm))
- {
- // travelling vertically down, skipping a horizontal line
- if ((pmap[y+2][x1]&0xFF)==cm && !parts[pmap[y+2][x1]>>8].life)
- {
- coord_stack[coord_stack_size][0] = x1;
- coord_stack[coord_stack_size][1] = y+2;
- coord_stack_size++;
- if (coord_stack_size>=coord_stack_limit)
- {
- free(coord_stack);
- return -1;
- }
- }
- }
- else if (y<YRES-CELL-1)
- {
- for (x=x1; x<=x2; x++)
- {
- if ((pmap[y+1][x]&0xFF)==cm && !parts[pmap[y+1][x]>>8].life)
- {
- if (x==x1 || x==x2 || y<0 || !PMAP_CMP_CONDUCTIVE(pmap[y-1][x], cm))
- {
- // if at the end of a horizontal section, or if it's a T junction
- coord_stack[coord_stack_size][0] = x;
- coord_stack[coord_stack_size][1] = y+1;
- coord_stack_size++;
- if (coord_stack_size>=coord_stack_limit)
- {
- free(coord_stack);
- return -1;
- }
- }
-
- }
- }
- }
- } while (coord_stack_size>0);
- free(coord_stack);
- return created_something;
-}
-
-
-int flood_parts(int x, int y, int fullc, int cm, int bm, int flags)
-{
- int c = fullc&0xFF;
- int x1, x2, dy = (c<PT_NUM)?1:CELL;
- int co = c;
- int coord_stack_limit = XRES*YRES;
- unsigned short (*coord_stack)[2];
- int coord_stack_size = 0;
- int created_something = 0;
-
- if (c==SPC_PROP)
- return 0;
- if (cm==-1)
- {
- if (c==0)
- {
- cm = pmap[y][x]&0xFF;
- if (!cm)
- return 0;
- if ((flags&BRUSH_REPLACEMODE) && cm!=SLALT)
- return 0;
- }
- else
- cm = 0;
- }
- if (bm==-1)
- {
- if (c-UI_WALLSTART+UI_ACTUALSTART==WL_ERASE)
- {
- bm = bmap[y/CELL][x/CELL];
- if (!bm)
- return 0;
- if (bm==WL_WALL)
- cm = 0xFF;
- }
- else
- bm = 0;
- }
-
- if (((pmap[y][x]&0xFF)!=cm || bmap[y/CELL][x/CELL]!=bm )||( (flags&BRUSH_SPECIFIC_DELETE) && cm!=SLALT))
- return 1;
-
- coord_stack = malloc(sizeof(unsigned short)*2*coord_stack_limit);
- coord_stack[coord_stack_size][0] = x;
- coord_stack[coord_stack_size][1] = y;
- coord_stack_size++;
-
- do
- {
- coord_stack_size--;
- x = coord_stack[coord_stack_size][0];
- y = coord_stack[coord_stack_size][1];
- x1 = x2 = x;
- // go left as far as possible
- while (x1>=CELL)
- {
- if ((pmap[y][x1-1]&0xFF)!=cm || bmap[y/CELL][(x1-1)/CELL]!=bm)
- {
- break;
- }
- x1--;
- }
- // go right as far as possible
- while (x2<XRES-CELL)
- {
- if ((pmap[y][x2+1]&0xFF)!=cm || bmap[y/CELL][(x2+1)/CELL]!=bm)
- {
- break;
- }
- x2++;
- }
- // fill span
- for (x=x1; x<=x2; x++)
- {
- if (create_parts(x, y, 0, 0, fullc, flags, 1))
- created_something = 1;
- }
-
- // add vertically adjacent pixels to stack
- if (y>=CELL+dy)
- for (x=x1; x<=x2; x++)
- if ((pmap[y-dy][x]&0xFF)==cm && bmap[(y-dy)/CELL][x/CELL]==bm)
- {
- coord_stack[coord_stack_size][0] = x;
- coord_stack[coord_stack_size][1] = y-dy;
- coord_stack_size++;
- if (coord_stack_size>=coord_stack_limit)
- {
- free(coord_stack);
- return -1;
- }
- }
- if (y<YRES-CELL-dy)
- for (x=x1; x<=x2; x++)
- if ((pmap[y+dy][x]&0xFF)==cm && bmap[(y+dy)/CELL][x/CELL]==bm)
- {
- coord_stack[coord_stack_size][0] = x;
- coord_stack[coord_stack_size][1] = y+dy;
- coord_stack_size++;
- if (coord_stack_size>=coord_stack_limit)
- {
- free(coord_stack);
- return -1;
- }
- }
- } while (coord_stack_size>0);
- free(coord_stack);
- return created_something;
-}
-
-int flood_water(int x, int y, int i, int originaly, int check)
-{
- int x1 = 0,x2 = 0;
- // go left as far as possible
- x1 = x2 = x;
- if (!pmap[y][x])
- return 1;
-
- while (x1>=CELL)
- {
- if ((ptypes[(pmap[y][x1-1]&0xFF)].falldown)!=2)
- {
- break;
- }
- x1--;
- }
- while (x2<XRES-CELL)
- {
- if ((ptypes[(pmap[y][x2+1]&0xFF)].falldown)!=2)
- {
- break;
- }
- x2++;
- }
-
- // fill span
- for (x=x1; x<=x2; x++)
- {
- parts[pmap[y][x]>>8].tmp2 = !check;//flag it as checked, maybe shouldn't use .tmp2
- //check above, maybe around other sides too?
- if ( ((y-1) > originaly) && !pmap[y-1][x] && eval_move(parts[i].type, x, y-1, NULL))
- {
- int oldx = (int)(parts[i].x + 0.5f);
- int oldy = (int)(parts[i].y + 0.5f);
- pmap[y-1][x] = pmap[oldy][oldx];
- pmap[oldy][oldx] = 0;
- parts[i].x = x;
- parts[i].y = y-1;
- return 0;
- }
- }
- // fill children
-
- if (y>=CELL+1)
- for (x=x1; x<=x2; x++)
- if ((ptypes[(pmap[y-1][x]&0xFF)].falldown)==2 && parts[pmap[y-1][x]>>8].tmp2 == check)
- if (!flood_water(x, y-1, i, originaly, check))
- return 0;
- if (y<YRES-CELL-1)
- for (x=x1; x<=x2; x++)
- if ((ptypes[(pmap[y+1][x]&0xFF)].falldown)==2 && parts[pmap[y+1][x]>>8].tmp2 == check)
- if (!flood_water(x, y+1, i, originaly, check))
- return 0;
- return 1;
-}
-
-//wrapper around create_part to create TESC with correct tmp value
-int create_part_add_props(int p, int x, int y, int tv, int rx, int ry)
-{
- p=create_part(p, x, y, tv);
- if (tv==PT_TESC)
- {
- parts[p].tmp=rx*4+ry*4+7;
- if (parts[p].tmp>300)
- parts[p].tmp=300;
- }
- return p;
-}
-
-//this creates particles from a brush, don't use if you want to create one particle
-int create_parts(int x, int y, int rx, int ry, int c, int flags, int fill)
-{
- int i, j, r, f = 0, u, v, oy, ox, b = 0, dw = 0, stemp = 0, p, fn;
-
- int wall = c - 100;
- if (c==SPC_WIND || c==PT_FIGH)
- return 0;
-
- if(c==SPC_PROP){
- prop_edit_ui(vid_buf, x, y);
- return 0;
- }
- if (c == SPC_AIR || c == SPC_HEAT || c == SPC_COOL || c == SPC_VACUUM || c == SPC_PGRV || c == SPC_NGRV)
- fill = 1;
- for (r=UI_ACTUALSTART; r<=UI_ACTUALSTART+UI_WALLCOUNT; r++)
- {
- if (wall==r)
- {
- if (c == SPC_AIR || c == SPC_HEAT || c == SPC_COOL || c == SPC_VACUUM || c == SPC_PGRV || c == SPC_NGRV || wall == WL_SIGN)
- break;
- if (wall == WL_ERASE)
- b = 0;
- else
- b = wall;
- dw = 1;
- }
- }
- if (c == WL_FANHELPER)
- {
- b = WL_FANHELPER;
- dw = 1;
- }
- if (wall == WL_GRAV)
- {
- gravwl_timeout = 60;
- }
- if (c==PT_LIGH)
- {
- if (lighting_recreate>0 && rx+ry>0)
- return 0;
- p=create_part(-2, x, y, c);
- if (p!=-1)
- {
- parts[p].life=rx+ry;
- if (parts[p].life>55)
- parts[p].life=55;
- parts[p].temp=parts[p].life*150; // temperature of the lighting shows the power of the lighting
- lighting_recreate+=parts[p].life/2+1;
- return 1;
- }
- else return 0;
- }
- if (c == PT_STKM || c == PT_STKM2 || c == PT_FIGH)
- rx = ry = 0;
-
- if (dw==1)
- {
- ry = ry/CELL;
- rx = rx/CELL;
- x = x/CELL;
- y = y/CELL;
- x -= rx/2;
- y -= ry/2;
- for (ox=x; ox<=x+rx; ox++)
- {
- for (oy=y; oy<=y+rx; oy++)
- {
- if (ox>=0&&ox<XRES/CELL&&oy>=0&&oy<YRES/CELL)
- {
- i = ox;
- j = oy;
- if ((flags&BRUSH_SPECIFIC_DELETE) && b!=WL_FANHELPER)
- {
- if (bmap[j][i]==SLALT-100)
- {
- b = 0;
- if (SLALT==WL_GRAV) gravwl_timeout = 60;
- }
- else
- continue;
- }
- if (b==WL_FAN)
- {
- fvx[j][i] = 0.0f;
- fvy[j][i] = 0.0f;
- }
- if (b==WL_STREAM)
- {
- i = x + rx/2;
- j = y + ry/2;
- for (v=-1; v<2; v++)
- for (u=-1; u<2; u++)
- if (i+u>=0 && i+u<XRES/CELL &&
- j+v>=0 && j+v<YRES/CELL &&
- bmap[j+v][i+u] == WL_STREAM)
- return 1;
- bmap[j][i] = WL_STREAM;
- continue;
- }
- if (b==0 && bmap[j][i]==WL_GRAV) gravwl_timeout = 60;
- bmap[j][i] = b;
- }
- }
- }
- return 1;
- }
-
- if (c == SPC_AIR || c == SPC_HEAT || c == SPC_COOL || c == SPC_VACUUM || c == SPC_PGRV || c == SPC_NGRV)
- fn = 3;
- else if (c == 0 && !(flags&BRUSH_REPLACEMODE)) // delete
- fn = 0;
- else if ((flags&BRUSH_SPECIFIC_DELETE) && !(flags&BRUSH_REPLACEMODE)) // specific delete
- fn = 1;
- else if (flags&BRUSH_REPLACEMODE) // replace
- fn = 2;
- else // normal draw
- fn = 3;
-
- if (rx<=0) //workaround for rx == 0 crashing. todo: find a better fix later.
- {
- for (j = y - ry; j <= y + ry; j++)
- if (create_parts2(fn,x,j,c,rx,ry,flags))
- f = 1;
- }
- else
- {
- int tempy = y, i, j, jmax, oldy;
- // tempy is the smallest y value that is still inside the brush
- // jmax is the largest y value that is still inside the brush
- if (CURRENT_BRUSH == TRI_BRUSH)
- tempy = y + ry;
- for (i = x - rx; i <= x; i++) {
- oldy = tempy;
- // Fix a problem with the triangle brush which occurs if the bottom corner (the first point tested) isn't recognised as being inside the brush
- if (!InCurrentBrush(i-x,tempy-y,rx,ry))
- continue;
- while (InCurrentBrush(i-x,tempy-y,rx,ry))
- tempy = tempy - 1;
- tempy = tempy + 1;
- if (fill)
- {
- jmax = 2*y - tempy;
- if (CURRENT_BRUSH == TRI_BRUSH)
- jmax = y + ry;
- for (j = tempy; j <= jmax; j++) {
- if (create_parts2(fn,i,j,c,rx,ry,flags))
- f = 1;
- if (i!=x && create_parts2(fn,2*x-i,j,c,rx,ry,flags))
- f = 1;
- }
- }
- else
- {
- if ((oldy != tempy && CURRENT_BRUSH != SQUARE_BRUSH) || i == x-rx)
- oldy--;
- //if (CURRENT_BRUSH == TRI_BRUSH)
- // oldy = tempy;
- for (j = tempy; j <= oldy+1; j++) {
- int i2 = 2*x-i, j2 = 2*y-j;
- if (CURRENT_BRUSH == TRI_BRUSH)
- j2 = y+ry;
- if (create_parts2(fn,i,j,c,rx,ry,flags))
- f = 1;
- if (i2 != i && create_parts2(fn,i2,j,c,rx,ry,flags))
- f = 1;
- if (j2 != j && create_parts2(fn,i,j2,c,rx,ry,flags))
- f = 1;
- if (i2 != i && j2 != j && create_parts2(fn,i2,j2,c,rx,ry,flags))
- f = 1;
- }
- }
- }
- }
- return !f;
-}
-
-int create_parts2(int f, int x, int y, int c, int rx, int ry, int flags)
-{
- if (f == 0) //delete
- delete_part(x, y, 0);
- else if (f == 1) //specific delete
- delete_part(x, y, flags);
- else if (f == 2) //replace mode
- {
- if (x<0 || y<0 || x>=XRES || y>=YRES)
- return 0;
- if ((pmap[y][x]&0xFF)!=SLALT&&SLALT!=0)
- return 0;
- if ((pmap[y][x]))
- {
- delete_part(x, y, 0);
- if (c!=0)
- create_part_add_props(-2, x, y, c, rx, ry);
- }
- }
- else if (f == 3) //normal draw
- if (create_part_add_props(-2, x, y, c, rx, ry)==-1)
- return 1;
- return 0;
-}
-
-int InCurrentBrush(int i, int j, int rx, int ry)
-{
- switch(CURRENT_BRUSH)
- {
- case CIRCLE_BRUSH:
- return (pow((double)i,2)*pow((double)ry,2)+pow((double)j,2)*pow((double)rx,2)<=pow((double)rx,2)*pow((double)ry,2));
- break;
- case SQUARE_BRUSH:
- return (abs(i) <= rx && abs(j) <= ry);
- break;
- case TRI_BRUSH:
- return ((abs((rx+2*i)*ry+rx*j) + abs(2*rx*(j-ry)) + abs((rx-2*i)*ry+rx*j))<=(4*rx*ry));
- break;
- default:
- return 0;
- break;
- }
-}
-int get_brush_flags()
-{
- int flags = 0;
- if (REPLACE_MODE)
- flags |= BRUSH_REPLACEMODE;
- if (sdl_mod & (KMOD_CAPS))
- flags |= BRUSH_SPECIFIC_DELETE;
- if (sdl_mod & (KMOD_LALT) && sdl_mod & (KMOD_CTRL))
- flags |= BRUSH_SPECIFIC_DELETE;
- return flags;
-}
-void create_line(int x1, int y1, int x2, int y2, int rx, int ry, int c, int flags)
-{
- int cp=abs(y2-y1)>abs(x2-x1), x, y, dx, dy, sy, fill = 1;
- float e, de;
- if (c==SPC_PROP)
- return;
- if (cp)
- {
- y = x1;
- x1 = y1;
- y1 = y;
- y = x2;
- x2 = y2;
- y2 = y;
- }
- if (x1 > x2)
- {
- y = x1;
- x1 = x2;
- x2 = y;
- y = y1;
- y1 = y2;
- y2 = y;
- }
- dx = x2 - x1;
- dy = abs(y2 - y1);
- e = 0.0f;
- if (dx)
- de = dy/(float)dx;
- else
- de = 0.0f;
- y = y1;
- sy = (y1<y2) ? 1 : -1;
- for (x=x1; x<=x2; x++)
- {
- if (cp)
- create_parts(y, x, rx, ry, c, flags, fill);
- else
- create_parts(x, y, rx, ry, c, flags, fill);
- fill = 0;
- e += de;
- if (e >= 0.5f)
- {
- y += sy;
- if ((c==WL_EHOLE+100 || c==WL_ALLOWGAS+100 || c==WL_ALLOWENERGY+100 || c==WL_ALLOWALLELEC+100 || c==WL_ALLOWSOLID+100 || c==WL_ALLOWAIR+100 || c==WL_WALL+100 || c==WL_DESTROYALL+100 || c==WL_ALLOWLIQUID+100 || c==WL_FAN+100 || c==WL_STREAM+100 || c==WL_DETECT+100 || c==WL_EWALL+100 || c==WL_WALLELEC+100 || !(rx+ry))
- && ((y1<y2) ? (y<=y2) : (y>=y2)))
- {
- if (cp)
- create_parts(y, x, rx, ry, c, flags, fill);
- else
- create_parts(x, y, rx, ry, c, flags, fill);
- }
- e -= 1.0f;
- }
- }
-}
-
-#if defined(WIN32) && !defined(__GNUC__)
-void orbitalparts_get(int block1, int block2, int resblock1[], int resblock2[])
-#else
-inline void orbitalparts_get(int block1, int block2, int resblock1[], int resblock2[])
-#endif
-{
- resblock1[0] = (block1&0x000000FF);
- resblock1[1] = (block1&0x0000FF00)>>8;
- resblock1[2] = (block1&0x00FF0000)>>16;
- resblock1[3] = (block1&0xFF000000)>>24;
-
- resblock2[0] = (block2&0x000000FF);
- resblock2[1] = (block2&0x0000FF00)>>8;
- resblock2[2] = (block2&0x00FF0000)>>16;
- resblock2[3] = (block2&0xFF000000)>>24;
-}
-
-#if defined(WIN32) && !defined(__GNUC__)
-void orbitalparts_set(int *block1, int *block2, int resblock1[], int resblock2[])
-#else
-inline void orbitalparts_set(int *block1, int *block2, int resblock1[], int resblock2[])
-#endif
-{
- int block1tmp = 0;
- int block2tmp = 0;
-
- block1tmp = (resblock1[0]&0xFF);
- block1tmp |= (resblock1[1]&0xFF)<<8;
- block1tmp |= (resblock1[2]&0xFF)<<16;
- block1tmp |= (resblock1[3]&0xFF)<<24;
-
- block2tmp = (resblock2[0]&0xFF);
- block2tmp |= (resblock2[1]&0xFF)<<8;
- block2tmp |= (resblock2[2]&0xFF)<<16;
- block2tmp |= (resblock2[3]&0xFF)<<24;
-
- *block1 = block1tmp;
- *block2 = block2tmp;
-}