X-Git-Url: https://deadsoftware.ru/gitweb?a=blobdiff_plain;f=src%2Fgame%2Fg_grid.pas;h=049868c7821e7784afa08dbe18e68ab97c0d6a31;hb=a0b4df67e96c6b5e55e95269a87975df67692e6b;hp=8ff01a629340b582ad58c0cb40046f8b80c251a3;hpb=829537fc7b2b29baaab7287e54fba2cdd8d00b5c;p=d2df-sdl.git diff --git a/src/game/g_grid.pas b/src/game/g_grid.pas index 8ff01a6..049868c 100644 --- a/src/game/g_grid.pas +++ b/src/game/g_grid.pas @@ -15,7 +15,12 @@ *) // universal spatial grid {$INCLUDE ../shared/a_modes.inc} -{$DEFINE grid_use_buckets} +{$IF DEFINED(D2F_DEBUG)} + {.$DEFINE D2F_DEBUG_RAYTRACE} + {.$DEFINE D2F_DEBUG_XXQ} + {.$DEFINE D2F_DEBUG_MOVER} +{$ENDIF} +{$DEFINE GRID_USE_ORTHO_ACCEL} unit g_grid; interface @@ -27,94 +32,385 @@ type generic TBodyGridBase = class(TObject) public type TGridQueryCB = function (obj: ITP; tag: Integer): Boolean is nested; // return `true` to stop + type TGridRayQueryCB = function (obj: ITP; tag: Integer; x, y, prevx, prevy: Integer): Boolean is nested; // return `true` to stop + + type TCellQueryCB = procedure (x, y: Integer) is nested; // top-left cell corner coords + + const TagDisabled = $40000000; + const TagFullMask = $3fffffff; private const - GridDefaultTileSize = 32; + GridDefaultTileSize = 32; // must be power of two! GridCellBucketSize = 8; // WARNING! can't be less than 2! - private + public type PBodyProxyRec = ^TBodyProxyRec; TBodyProxyRec = record private mX, mY, mWidth, mHeight: Integer; // aabb - mQueryMark: DWord; // was this object visited at this query? + mQueryMark: LongWord; // was this object visited at this query? mObj: ITP; - mTag: Integer; + mTag: Integer; // `TagDisabled` set: disabled ;-) nextLink: TBodyProxyId; // next free or nothing private procedure setup (aX, aY, aWidth, aHeight: Integer; aObj: ITP; aTag: Integer); + + function getTag (): Integer; inline; + procedure setTag (v: Integer); inline; + + function getEnabled (): Boolean; inline; + procedure setEnabled (v: Boolean); inline; + + function getX1 (): Integer; inline; + function getY1 (): Integer; inline; + + public + property x: Integer read mX; + property y: Integer read mY; + property width: Integer read mWidth; + property height: Integer read mHeight; + property tag: Integer read getTag write setTag; + property enabled: Boolean read getEnabled write setEnabled; + property obj: ITP read mObj; + + property x0: Integer read mX; + property y0: Integer read mY; + property x1: Integer read getX1; + property y1: Integer read getY1; end; + private + type PGridCell = ^TGridCell; TGridCell = record - {$IFDEF grid_use_buckets} bodies: array [0..GridCellBucketSize-1] of Integer; // -1: end of list - {$ELSE} - body: Integer; - {$ENDIF} next: Integer; // in this cell; index in mCells end; - TGridInternalCB = function (grida: Integer): Boolean of object; // return `true` to stop + TCellArray = array of TGridCell; + + TGridInternalCB = function (grida: Integer; bodyId: TBodyProxyId): Boolean of object; // return `true` to stop + + private + //mTileSize: Integer; + const mTileSize = GridDefaultTileSize; + type TGetProxyFn = function (pxidx: Integer): PBodyProxyRec of object; + + public + const tileSize = mTileSize; + + type + TAtPointEnumerator = record + private + mCells: TCellArray; + curidx, curbki: Integer; + getpx: TGetProxyFn; + public + constructor Create (acells: TCellArray; aidx: Integer; agetpx: TGetProxyFn); + function MoveNext (): Boolean; inline; + function getCurrent (): PBodyProxyRec; inline; + property Current: PBodyProxyRec read getCurrent; + end; private - mTileSize: Integer; mMinX, mMinY: Integer; // so grids can start at any origin mWidth, mHeight: Integer; // in tiles mGrid: array of Integer; // mWidth*mHeight, index in mCells - mCells: array of TGridCell; // cell pool + mCells: TCellArray; // cell pool mFreeCell: Integer; // first free cell index or -1 - mLastQuery: DWord; + mLastQuery: LongWord; mUsedCells: Integer; mProxies: array of TBodyProxyRec; mProxyFree: TBodyProxyId; // free mProxyCount: Integer; // currently used mProxyMaxCount: Integer; + mInQuery: Boolean; - mUData: TBodyProxyId; // for inserter/remover - mTagMask: Integer; // for iterator - mItCB: TGridQueryCB; // for iterator + public + dbgShowTraceLog: Boolean; + {$IF DEFINED(D2F_DEBUG)} + dbgRayTraceTileHitCB: TCellQueryCB; + {$ENDIF} private - function allocCell: Integer; + function allocCell (): Integer; procedure freeCell (idx: Integer); // `next` is simply overwritten function allocProxy (aX, aY, aWidth, aHeight: Integer; aObj: ITP; aTag: Integer): TBodyProxyId; procedure freeProxy (body: TBodyProxyId); - procedure insert (body: TBodyProxyId); - procedure remove (body: TBodyProxyId); + procedure insertInternal (body: TBodyProxyId); + procedure removeInternal (body: TBodyProxyId); + + function forGridRect (x, y, w, h: Integer; cb: TGridInternalCB; bodyId: TBodyProxyId): Boolean; - function forGridRect (x, y, w, h: Integer; cb: TGridInternalCB): Boolean; + function inserter (grida: Integer; bodyId: TBodyProxyId): Boolean; + function remover (grida: Integer; bodyId: TBodyProxyId): Boolean; - function inserter (grida: Integer): Boolean; - function remover (grida: Integer): Boolean; - function iterator (grida: Integer): Boolean; + function getProxyEnabled (pid: TBodyProxyId): Boolean; inline; + procedure setProxyEnabled (pid: TBodyProxyId; val: Boolean); inline; + + function getGridWidthPx (): Integer; inline; + function getGridHeightPx (): Integer; inline; + + function getProxyById (idx: TBodyProxyId): PBodyProxyRec; inline; public - constructor Create (aMinPixX, aMinPixY, aPixWidth, aPixHeight: Integer; aTileSize: Integer=GridDefaultTileSize); + constructor Create (aMinPixX, aMinPixY, aPixWidth, aPixHeight: Integer{; aTileSize: Integer=GridDefaultTileSize}); destructor Destroy (); override; - function insertBody (aObj: ITP; ax, ay, aWidth, aHeight: Integer; aTag: Integer=0): TBodyProxyId; - procedure removeBody (aObj: TBodyProxyId); // WARNING! this WILL destroy proxy! + function insertBody (aObj: ITP; ax, ay, aWidth, aHeight: Integer; aTag: Integer=-1): TBodyProxyId; + procedure removeBody (body: TBodyProxyId); // WARNING! this WILL destroy proxy! + + procedure moveBody (body: TBodyProxyId; nx, ny: Integer); + procedure resizeBody (body: TBodyProxyId; nw, nh: Integer); + procedure moveResizeBody (body: TBodyProxyId; nx, ny, nw, nh: Integer); + + function insideGrid (x, y: Integer): Boolean; inline; + + // `false` if `body` is surely invalid + function getBodyXY (body: TBodyProxyId; out rx, ry: Integer): Boolean; inline; + function getBodyWH (body: TBodyProxyId; out rw, rh: Integer): Boolean; inline; + function getBodyDims (body: TBodyProxyId; out rx, ry, rw, rh: Integer): Boolean; inline; + + //WARNING: don't modify grid while any query is in progress (no checks are made!) + // you can set enabled/disabled flag, tho (but iterator can still return objects disabled inside it) + // no callback: return `true` on the first hit + function forEachInAABB (x, y, w, h: Integer; cb: TGridQueryCB; tagmask: Integer=-1; allowDisabled: Boolean=false): ITP; + + //WARNING: don't modify grid while any query is in progress (no checks are made!) + // you can set enabled/disabled flag, tho (but iterator can still return objects disabled inside it) + // no callback: return object on the first hit or nil + function forEachAtPoint (x, y: Integer; cb: TGridQueryCB; tagmask: Integer=-1; exittag: PInteger=nil): ITP; + + function atCellInPoint (x, y: Integer): TAtPointEnumerator; + + //WARNING: don't modify grid while any query is in progress (no checks are made!) + // you can set enabled/disabled flag, tho (but iterator can still return objects disabled inside it) + // cb with `(nil)` will be called before processing new tile + // no callback: return object of the nearest hit or nil + // if `inverted` is true, trace will register bodies *exluding* tagmask + //WARNING: don't change tags in callbacks here! + function traceRay (const x0, y0, x1, y1: Integer; cb: TGridRayQueryCB; tagmask: Integer=-1): ITP; overload; + function traceRay (out ex, ey: Integer; const ax0, ay0, ax1, ay1: Integer; cb: TGridRayQueryCB; tagmask: Integer=-1): ITP; + + // return `false` if we're still inside at the end + // line should be either strict horizontal, or strict vertical, otherwise an exception will be thrown + function traceOrthoRayWhileIn (out ex, ey: Integer; ax0, ay0, ax1, ay1: Integer; tagmask: Integer=-1): Boolean; + + //WARNING: don't modify grid while any query is in progress (no checks are made!) + // you can set enabled/disabled flag, tho (but iterator can still return objects disabled inside it) + // trace line along the grid, calling `cb` for all objects in passed cells, in no particular order + //WARNING: don't change tags in callbacks here! + function forEachAlongLine (ax0, ay0, ax1, ay1: Integer; cb: TGridQueryCB; tagmask: Integer=-1; log: Boolean=false): ITP; + + // debug + procedure forEachBodyCell (body: TBodyProxyId; cb: TCellQueryCB); + function forEachInCell (x, y: Integer; cb: TGridQueryCB): ITP; + procedure dumpStats (); - procedure moveBody (body: TBodyProxyId; dx, dy: Integer); - procedure resizeBody (body: TBodyProxyId; sx, sy: Integer); - procedure moveResizeBody (body: TBodyProxyId; dx, dy, sx, sy: Integer); + public + //WARNING! no sanity checks! + property proxyEnabled[pid: TBodyProxyId]: Boolean read getProxyEnabled write setProxyEnabled; - function forEachInAABB (x, y, w, h: Integer; cb: TGridQueryCB; tagmask: Integer=-1): Boolean; + property gridX0: Integer read mMinX; + property gridY0: Integer read mMinY; + property gridWidth: Integer read getGridWidthPx; // in pixels + property gridHeight: Integer read getGridHeightPx; // in pixels - procedure dumpStats (); + property proxy[idx: TBodyProxyId]: PBodyProxyRec read getProxyById; end; +// you are not supposed to understand this +// returns `true` if there is an intersection, and enter coords +// enter coords will be equal to (x0, y0) if starting point is inside the box +// if result is `false`, `inx` and `iny` are undefined +function lineAABBIntersects (x0, y0, x1, y1: Integer; bx, by, bw, bh: Integer; out inx, iny: Integer): Boolean; + +function distanceSq (x0, y0, x1, y1: Integer): Integer; inline; + +procedure swapInt (var a: Integer; var b: Integer); inline; +function minInt (a, b: Integer): Integer; inline; +function maxInt (a, b: Integer): Integer; inline; + + implementation uses - SysUtils, e_log; + SysUtils, e_log, g_console, utils; + + +// ////////////////////////////////////////////////////////////////////////// // +procedure swapInt (var a: Integer; var b: Integer); inline; var t: Integer; begin t := a; a := b; b := t; end; +function minInt (a, b: Integer): Integer; inline; begin if (a < b) then result := a else result := b; end; +function maxInt (a, b: Integer): Integer; inline; begin if (a > b) then result := a else result := b; end; + +function distanceSq (x0, y0, x1, y1: Integer): Integer; inline; begin result := (x1-x0)*(x1-x0)+(y1-y0)*(y1-y0); end; + + +// ////////////////////////////////////////////////////////////////////////// // +// you are not supposed to understand this +// returns `true` if there is an intersection, and enter coords +// enter coords will be equal to (x0, y0) if starting point is inside the box +// if result is `false`, `inx` and `iny` are undefined +function lineAABBIntersects (x0, y0, x1, y1: Integer; bx, by, bw, bh: Integer; out inx, iny: Integer): Boolean; +var + wx0, wy0, wx1, wy1: Integer; // window coordinates + stx, sty: Integer; // "steps" for x and y axes + dsx, dsy: Integer; // "lengthes" for x and y axes + dx2, dy2: Integer; // "double lengthes" for x and y axes + xd, yd: Integer; // current coord + e: Integer; // "error" (as in bresenham algo) + rem: Integer; + //!term: Integer; + d0, d1: PInteger; + xfixed: Boolean; + temp: Integer; +begin + result := false; + // why not + inx := x0; + iny := y0; + if (bw < 1) or (bh < 1) then exit; // impossible box + + if (x0 = x1) and (y0 = y1) then + begin + // check this point + result := (x0 >= bx) and (y0 >= by) and (x0 < bx+bw) and (y0 < by+bh); + exit; + end; + + // check if staring point is inside the box + if (x0 >= bx) and (y0 >= by) and (x0 < bx+bw) and (y0 < by+bh) then begin result := true; exit; end; + + // clip rectange + wx0 := bx; + wy0 := by; + wx1 := bx+bw-1; + wy1 := by+bh-1; + + // horizontal setup + if (x0 < x1) then + begin + // from left to right + if (x0 > wx1) or (x1 < wx0) then exit; // out of screen + stx := 1; // going right + end + else + begin + // from right to left + if (x1 > wx1) or (x0 < wx0) then exit; // out of screen + stx := -1; // going left + x0 := -x0; + x1 := -x1; + wx0 := -wx0; + wx1 := -wx1; + swapInt(wx0, wx1); + end; + + // vertical setup + if (y0 < y1) then + begin + // from top to bottom + if (y0 > wy1) or (y1 < wy0) then exit; // out of screen + sty := 1; // going down + end + else + begin + // from bottom to top + if (y1 > wy1) or (y0 < wy0) then exit; // out of screen + sty := -1; // going up + y0 := -y0; + y1 := -y1; + wy0 := -wy0; + wy1 := -wy1; + swapInt(wy0, wy1); + end; + + dsx := x1-x0; + dsy := y1-y0; + + if (dsx < dsy) then + begin + d0 := @yd; + d1 := @xd; + swapInt(x0, y0); + swapInt(x1, y1); + swapInt(dsx, dsy); + swapInt(wx0, wy0); + swapInt(wx1, wy1); + swapInt(stx, sty); + end + else + begin + d0 := @xd; + d1 := @yd; + end; + + dx2 := 2*dsx; + dy2 := 2*dsy; + xd := x0; + yd := y0; + e := 2*dsy-dsx; + //!term := x1; + + xfixed := false; + if (y0 < wy0) then + begin + // clip at top + temp := dx2*(wy0-y0)-dsx; + xd += temp div dy2; + rem := temp mod dy2; + if (xd > wx1) then exit; // x is moved out of clipping rect, nothing to do + if (xd+1 >= wx0) then + begin + yd := wy0; + e -= rem+dsx; + if (rem > 0) then begin Inc(xd); e += dy2; end; + xfixed := true; + end; + end; + + if (not xfixed) and (x0 < wx0) then + begin + // clip at left + temp := dy2*(wx0-x0); + yd += temp div dx2; + rem := temp mod dx2; + if (yd > wy1) or (yd = wy1) and (rem >= dsx) then exit; + xd := wx0; + e += rem; + if (rem >= dsx) then begin Inc(yd); e -= dx2; end; + end; + + (* + if (y1 > wy1) then + begin + // clip at bottom + temp := dx2*(wy1-y0)+dsx; + term := x0+temp div dy2; + rem := temp mod dy2; + if (rem = 0) then Dec(term); + end; + + if (term > wx1) then term := wx1; // clip at right + + Inc(term); // draw last point + //if (term = xd) then exit; // this is the only point, get out of here + *) + + if (sty = -1) then yd := -yd; + if (stx = -1) then begin xd := -xd; {!term := -term;} end; + //!dx2 -= dy2; + + inx := d0^; + iny := d1^; + result := true; +end; // ////////////////////////////////////////////////////////////////////////// // @@ -131,20 +427,91 @@ begin end; +function TBodyGridBase.TBodyProxyRec.getTag (): Integer; inline; +begin + result := mTag and TagFullMask; +end; + +procedure TBodyGridBase.TBodyProxyRec.setTag (v: Integer); inline; +begin + mTag := (mTag and TagDisabled) or (v and TagFullMask); +end; + +function TBodyGridBase.TBodyProxyRec.getEnabled (): Boolean; inline; +begin + result := ((mTag and TagDisabled) = 0); +end; + +procedure TBodyGridBase.TBodyProxyRec.setEnabled (v: Boolean); inline; +begin + if v then mTag := mTag and (not TagDisabled) else mTag := mTag or TagDisabled; +end; + +function TBodyGridBase.TBodyProxyRec.getX1 (): Integer; inline; +begin + result := mX+mWidth-1; +end; + +function TBodyGridBase.TBodyProxyRec.getY1 (): Integer; inline; +begin + result := mY+mHeight-1; +end; + + +// ////////////////////////////////////////////////////////////////////////// // +constructor TBodyGridBase.TAtPointEnumerator.Create (acells: TCellArray; aidx: Integer; agetpx: TGetProxyFn); +begin + mCells := acells; + curidx := aidx; + curbki := -1; + getpx := agetpx; +end; + + +function TBodyGridBase.TAtPointEnumerator.MoveNext (): Boolean; inline; +begin + while (curidx <> -1) do + begin + while (curbki < GridCellBucketSize) do + begin + Inc(curbki); + if (mCells[curidx].bodies[curbki] = -1) then break; + result := true; + exit; + end; + curidx := mCells[curidx].next; + curbki := -1; + end; + result := false; +end; + + +function TBodyGridBase.TAtPointEnumerator.getCurrent (): PBodyProxyRec; inline; +begin + result := getpx(mCells[curidx].bodies[curbki]); +end; + + // ////////////////////////////////////////////////////////////////////////// // -constructor TBodyGridBase.Create (aMinPixX, aMinPixY, aPixWidth, aPixHeight: Integer; aTileSize: Integer=GridDefaultTileSize); +constructor TBodyGridBase.Create (aMinPixX, aMinPixY, aPixWidth, aPixHeight: Integer{; aTileSize: Integer=GridDefaultTileSize}); var idx: Integer; begin + dbgShowTraceLog := false; + {$IF DEFINED(D2F_DEBUG)} + dbgRayTraceTileHitCB := nil; + {$ENDIF} + { if aTileSize < 1 then aTileSize := 1; if aTileSize > 8192 then aTileSize := 8192; // arbitrary limit - if aPixWidth < aTileSize then aPixWidth := aTileSize; - if aPixHeight < aTileSize then aPixHeight := aTileSize; mTileSize := aTileSize; + } + if (aPixWidth < mTileSize) then aPixWidth := mTileSize; + if (aPixHeight < mTileSize) then aPixHeight := mTileSize; mMinX := aMinPixX; mMinY := aMinPixY; - mWidth := (aPixWidth+aTileSize-1) div aTileSize; - mHeight := (aPixHeight+aTileSize-1) div aTileSize; + mWidth := (aPixWidth+mTileSize-1) div mTileSize; + mHeight := (aPixHeight+mTileSize-1) div mTileSize; SetLength(mGrid, mWidth*mHeight); SetLength(mCells, mWidth*mHeight); SetLength(mProxies, 8192); @@ -152,11 +519,8 @@ begin // init free list for idx := 0 to High(mCells) do begin - {$IFDEF grid_use_buckets} mCells[idx].bodies[0] := -1; - {$ELSE} - mCells[idx].body := -1; - {$ENDIF} + mCells[idx].bodies[GridCellBucketSize-1] := -1; // "has free room" flag mCells[idx].next := idx+1; end; mCells[High(mCells)].next := -1; // last cell @@ -170,9 +534,6 @@ begin mProxyFree := 0; mProxyCount := 0; mProxyMaxCount := 0; - mUData := 0; - mTagMask := -1; - mItCB := nil; e_WriteLog(Format('created grid with size: %dx%d (tile size: %d); pix: %dx%d', [mWidth, mHeight, mTileSize, mWidth*mTileSize, mHeight*mTileSize]), MSG_NOTIFY); end; @@ -186,6 +547,7 @@ begin end; +// ////////////////////////////////////////////////////////////////////////// // procedure TBodyGridBase.dumpStats (); var idx, mcb, cidx, cnt: Integer; @@ -206,9 +568,154 @@ begin end; -function TBodyGridBase.allocCell: Integer; +procedure TBodyGridBase.forEachBodyCell (body: TBodyProxyId; cb: TCellQueryCB); +var + g, f, cidx: Integer; + cc: PGridCell; +begin + if (body < 0) or (body > High(mProxies)) or not assigned(cb) then exit; + for g := 0 to High(mGrid) do + begin + cidx := mGrid[g]; + while (cidx <> -1) do + begin + cc := @mCells[cidx]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + if (cc.bodies[f] = body) then cb((g mod mWidth)*mTileSize+mMinX, (g div mWidth)*mTileSize+mMinY); + end; + // next cell + cidx := cc.next; + end; + end; +end; + + +function TBodyGridBase.forEachInCell (x, y: Integer; cb: TGridQueryCB): ITP; +var + f, cidx: Integer; + cc: PGridCell; +begin + result := Default(ITP); + if not assigned(cb) then exit; + Dec(x, mMinX); + Dec(y, mMinY); + if (x < 0) or (y < 0) or (x >= mWidth*mTileSize) or (y > mHeight*mTileSize) then exit; + cidx := mGrid[(y div mTileSize)*mWidth+(x div mTileSize)]; + while (cidx <> -1) do + begin + cc := @mCells[cidx]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + if cb(mProxies[cc.bodies[f]].mObj, mProxies[cc.bodies[f]].mTag) then begin result := mProxies[cc.bodies[f]].mObj; exit; end; + end; + // next cell + cidx := cc.next; + end; +end; + + +// ////////////////////////////////////////////////////////////////////////// // +function TBodyGridBase.getGridWidthPx (): Integer; inline; begin result := mWidth*mTileSize; end; +function TBodyGridBase.getGridHeightPx (): Integer; inline; begin result := mHeight*mTileSize; end; + + +function TBodyGridBase.insideGrid (x, y: Integer): Boolean; inline; +begin + // fix coords + Dec(x, mMinX); + Dec(y, mMinY); + result := (x >= 0) and (y >= 0) and (x < mWidth*mTileSize) and (y < mHeight*mTileSize); +end; + + +function TBodyGridBase.getBodyXY (body: TBodyProxyId; out rx, ry: Integer): Boolean; inline; +begin + if (body >= 0) and (body < Length(mProxies)) then + begin + with mProxies[body] do begin rx := mX; ry := mY; end; + result := true; + end + else + begin + rx := 0; + ry := 0; + result := false; + end; +end; + + +function TBodyGridBase.getBodyWH (body: TBodyProxyId; out rw, rh: Integer): Boolean; inline; +begin + if (body >= 0) and (body < Length(mProxies)) then + begin + with mProxies[body] do begin rw := mWidth; rh := mHeight; end; + result := true; + end + else + begin + rw := 0; + rh := 0; + result := false; + end; +end; + + +function TBodyGridBase.getBodyDims (body: TBodyProxyId; out rx, ry, rw, rh: Integer): Boolean; inline; +begin + if (body >= 0) and (body < Length(mProxies)) then + begin + with mProxies[body] do begin rx := mX; ry := mY; rw := mWidth; rh := mHeight; end; + result := true; + end + else + begin + rx := 0; + ry := 0; + rw := 0; + rh := 0; + result := false; + end; +end; + + + +// ////////////////////////////////////////////////////////////////////////// // +function TBodyGridBase.getProxyEnabled (pid: TBodyProxyId): Boolean; inline; +begin + if (pid >= 0) then result := ((mProxies[pid].mTag and TagDisabled) = 0) else result := false; +end; + + +procedure TBodyGridBase.setProxyEnabled (pid: TBodyProxyId; val: Boolean); inline; +begin + if (pid >= 0) then + begin + if val then + begin + mProxies[pid].mTag := mProxies[pid].mTag and not TagDisabled; + end + else + begin + mProxies[pid].mTag := mProxies[pid].mTag or TagDisabled; + end; + end; +end; + + +function TBodyGridBase.getProxyById (idx: TBodyProxyId): PBodyProxyRec; inline; +begin + if (idx >= 0) and (idx < High(mProxies)) then result := @mProxies[idx] else result := nil; +end; + + +// ////////////////////////////////////////////////////////////////////////// // +function TBodyGridBase.allocCell (): Integer; var idx: Integer; + pc: PGridCell; begin if (mFreeCell < 0) then begin @@ -217,18 +724,16 @@ begin SetLength(mCells, mFreeCell+32768); // arbitrary number for idx := mFreeCell to High(mCells) do begin - {$IFDEF grid_use_buckets} mCells[idx].bodies[0] := -1; - {$ELSE} - mCells[idx].body := -1; - {$ENDIF} + mCells[idx].bodies[GridCellBucketSize-1] := -1; // 'has free room' flag mCells[idx].next := idx+1; end; mCells[High(mCells)].next := -1; // last cell end; result := mFreeCell; - mFreeCell := mCells[result].next; - mCells[result].next := -1; + pc := @mCells[result]; + mFreeCell := pc.next; + pc.next := -1; Inc(mUsedCells); //e_WriteLog(Format('grid: allocated new cell #%d (total: %d)', [result, mUsedCells]), MSG_NOTIFY); end; @@ -236,17 +741,21 @@ end; procedure TBodyGridBase.freeCell (idx: Integer); begin - if (idx >= 0) and (idx < High(mCells)) then + if (idx >= 0) and (idx < Length(mCells)) then begin - //if mCells[idx].body = -1 then exit; // the thing that should not be - //mCells[idx].body := -1; - mCells[idx].next := mFreeCell; + with mCells[idx] do + begin + bodies[0] := -1; + bodies[GridCellBucketSize-1] := -1; // 'has free room' flag + next := mFreeCell; + end; mFreeCell := idx; Dec(mUsedCells); end; end; +// ////////////////////////////////////////////////////////////////////////// // function TBodyGridBase.allocProxy (aX, aY, aWidth, aHeight: Integer; aObj: ITP; aTag: Integer): TBodyProxyId; var olen, idx: Integer; @@ -285,9 +794,13 @@ begin end; -function TBodyGridBase.forGridRect (x, y, w, h: Integer; cb: TGridInternalCB): Boolean; +// ////////////////////////////////////////////////////////////////////////// // +function TBodyGridBase.forGridRect (x, y, w, h: Integer; cb: TGridInternalCB; bodyId: TBodyProxyId): Boolean; +const + tsize = mTileSize; var gx, gy: Integer; + gw, gh: Integer; begin result := false; if (w < 1) or (h < 1) or not assigned(cb) then exit; @@ -296,243 +809,519 @@ begin Dec(y, mMinY); // go on if (x+w <= 0) or (y+h <= 0) then exit; - if (x >= mWidth*mTileSize) or (y >= mHeight*mTileSize) then exit; - for gy := y div mTileSize to (y+h-1) div mTileSize do + gw := mWidth; + gh := mHeight; + //tsize := mTileSize; + if (x >= gw*tsize) or (y >= gh*tsize) then exit; + for gy := y div tsize to (y+h-1) div tsize do begin if (gy < 0) then continue; - if (gy >= mHeight) then break; - for gx := x div mTileSize to (x+w-1) div mTileSize do + if (gy >= gh) then break; + for gx := x div tsize to (x+w-1) div tsize do begin if (gx < 0) then continue; - if (gx >= mWidth) then break; - if (cb(gy*mWidth+gx)) then begin result := true; exit; end; + if (gx >= gw) then break; + result := cb(gy*gw+gx, bodyId); + if result then exit; end; end; end; -function TBodyGridBase.inserter (grida: Integer): Boolean; +// ////////////////////////////////////////////////////////////////////////// // +function TBodyGridBase.inserter (grida: Integer; bodyId: TBodyProxyId): Boolean; var cidx: Integer; - pc: PInteger; - {$IFDEF grid_use_buckets} + pc: Integer; pi: PGridCell; f: Integer; - {$ENDIF} begin result := false; // never stop // add body to the given grid cell - pc := @mGrid[grida]; - {$IFDEF grid_use_buckets} - if (pc^ <> -1) then + pc := mGrid[grida]; + if (pc <> -1) then begin - pi := @mCells[pc^]; - f := 0; - for f := 0 to High(TGridCell.bodies) do + {$IF DEFINED(D2F_DEBUG)} + cidx := pc; + while (cidx <> -1) do + begin + pi := @mCells[cidx]; + for f := 0 to GridCellBucketSize-1 do + begin + if (pi.bodies[f] = -1) then break; + if (pi.bodies[f] = bodyId) then raise Exception.Create('trying to insert already inserted proxy'); + end; + cidx := pi.next; + end; + {$ENDIF} + cidx := pc; + while (cidx <> -1) do begin - if (pi.bodies[f] = -1) then + pi := @mCells[cidx]; + // check "has room" flag + if (pi.bodies[GridCellBucketSize-1] = -1) then begin // can add here - pi.bodies[f] := mUData; - if (f+1 < Length(TGridCell.bodies)) then pi.bodies[f+1] := -1; - exit; + for f := 0 to GridCellBucketSize-1 do + begin + if (pi.bodies[f] = -1) then + begin + pi.bodies[f] := bodyId; + if (f+1 < GridCellBucketSize) then pi.bodies[f+1] := -1; + exit; + end; + end; + raise Exception.Create('internal error in grid inserter'); end; + // no room, go to next cell in list (if there is any) + cidx := pi.next; end; + // no room in cells, add new cell to list end; // either no room, or no cell at all cidx := allocCell(); - mCells[cidx].bodies[0] := mUData; - mCells[cidx].bodies[1] := -1; - mCells[cidx].next := pc^; - pc^ := cidx; - {$ELSE} - cidx := allocCell(); - //e_WriteLog(Format(' 01: allocated cell for grid coords (%d,%d), body coords:(%d,%d): #%d', [gx, gy, dx, dy, cidx]), MSG_NOTIFY); - mCells[cidx].body := mUData; - mCells[cidx].next := pc^; - pc^ := cidx; - {$ENDIF} + pi := @mCells[cidx]; + pi.bodies[0] := bodyId; + pi.bodies[1] := -1; + pi.next := pc; + mGrid[grida] := cidx; end; - -procedure TBodyGridBase.insert (body: TBodyProxyId); +procedure TBodyGridBase.insertInternal (body: TBodyProxyId); var px: PBodyProxyRec; - oudata: Integer; begin if (body < 0) or (body > High(mProxies)) then exit; // just in case px := @mProxies[body]; - oudata := mUData; - mUData := body; - forGridRect(px.mX, px.mY, px.mWidth, px.mHeight, inserter); - mUData := oudata; + forGridRect(px.mX, px.mY, px.mWidth, px.mHeight, inserter, body); end; -function TBodyGridBase.remover (grida: Integer): Boolean; +// assume that we cannot have one object added to bucket twice +function TBodyGridBase.remover (grida: Integer; bodyId: TBodyProxyId): Boolean; var - pidx, idx, tmp, f: Integer; + f, c: Integer; + pidx, cidx: Integer; pc: PGridCell; begin result := false; // never stop // find and remove cell - pidx := -1; - idx := mGrid[grida]; - while (idx >= 0) do + pidx := -1; // previous cell index + cidx := mGrid[grida]; // current cell index + while (cidx <> -1) do begin - tmp := mCells[idx].next; - {$IFDEF grid_use_buckets} - pc := @mCells[idx]; - f := 0; - while (f < High(TGridCell.bodies)) do + pc := @mCells[cidx]; + for f := 0 to GridCellBucketSize-1 do begin - if (pc.bodies[f] = mUData) then + if (pc.bodies[f] = bodyId) then begin // i found her! if (f = 0) and (pc.bodies[1] = -1) then begin // this cell contains no elements, remove it - tmp := mCells[idx].next; - if (pidx = -1) then mGrid[grida] := tmp else mCells[pidx].next := tmp; - freeCell(idx); - end - else + if (pidx = -1) then mGrid[grida] := pc.next else mCells[pidx].next := pc.next; + freeCell(cidx); + exit; + end; + // remove element from bucket + for c := f to GridCellBucketSize-2 do begin - // remove element from bucket - Inc(f); - while (f < High(TGridCell.bodies)) do - begin - pc.bodies[f-1] := pc.bodies[f]; - if (pc.bodies[f] = -1) then break; - Inc(f); - end; - pc.bodies[High(TGridCell.bodies)] := -1; // just in case + pc.bodies[c] := pc.bodies[c+1]; + if (pc.bodies[c] = -1) then break; end; - exit; // assume that we cannot have one object added to bucket twice + pc.bodies[GridCellBucketSize-1] := -1; // "has free room" flag + exit; end; - Inc(f); end; - {$ELSE} - if (mCells[idx].body = mUData) then - begin - if (pidx = -1) then mGrid[grida] := tmp else mCells[pidx].next := tmp; - freeCell(idx); - exit; // assume that we cannot have one object added to bucket twice - end; - {$ENDIF} - pidx := idx; - idx := tmp; + pidx := cidx; + cidx := pc.next; end; end; - -// absolutely not tested -procedure TBodyGridBase.remove (body: TBodyProxyId); +procedure TBodyGridBase.removeInternal (body: TBodyProxyId); var px: PBodyProxyRec; - oudata: Integer; begin if (body < 0) or (body > High(mProxies)) then exit; // just in case px := @mProxies[body]; - oudata := mUData; - mUData := body; - forGridRect(px.mX, px.mY, px.mWidth, px.mHeight, remover); - mUData := oudata; + forGridRect(px.mX, px.mY, px.mWidth, px.mHeight, remover, body); end; -function TBodyGridBase.insertBody (aObj: ITP; aX, aY, aWidth, aHeight: Integer; aTag: Integer=0): TBodyProxyId; +// ////////////////////////////////////////////////////////////////////////// // +function TBodyGridBase.insertBody (aObj: ITP; aX, aY, aWidth, aHeight: Integer; aTag: Integer=-1): TBodyProxyId; begin + aTag := aTag and TagFullMask; result := allocProxy(aX, aY, aWidth, aHeight, aObj, aTag); - insert(result); + insertInternal(result); end; -procedure TBodyGridBase.removeBody (aObj: TBodyProxyId); +procedure TBodyGridBase.removeBody (body: TBodyProxyId); begin - if (aObj < 0) or (aObj > High(mProxies)) then exit; // just in case - remove(aObj); - freeProxy(aObj); + if (body < 0) or (body > High(mProxies)) then exit; // just in case + removeInternal(body); + freeProxy(body); end; -procedure TBodyGridBase.moveResizeBody (body: TBodyProxyId; dx, dy, sx, sy: Integer); +// ////////////////////////////////////////////////////////////////////////// // +procedure TBodyGridBase.moveResizeBody (body: TBodyProxyId; nx, ny, nw, nh: Integer); var px: PBodyProxyRec; + x0, y0, w, h: Integer; begin if (body < 0) or (body > High(mProxies)) then exit; // just in case - if ((dx = 0) and (dy = 0) and (sx = 0) and (sy = 0)) then exit; - remove(body); px := @mProxies[body]; - Inc(px.mX, dx); - Inc(px.mY, dy); - Inc(px.mWidth, sx); - Inc(px.mHeight, sy); - insert(body); -end; - -procedure TBodyGridBase.moveBody (body: TBodyProxyId; dx, dy: Integer); -begin - moveResizeBody(body, dx, dy, 0, 0); -end; - -procedure TBodyGridBase.resizeBody (body: TBodyProxyId; sx, sy: Integer); -begin - moveResizeBody(body, 0, 0, sx, sy); + x0 := px.mX; + y0 := px.mY; + w := px.mWidth; + h := px.mHeight; + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format('proxy #%d: MOVERESIZE: xg=%d;yg=%d;w=%d;h=%d;nx=%d;ny=%d;nw=%d;nh=%d', [body, x0-mMinX, y0-mMinY, w, h, nx-mMinX, ny-mMinY, nw, nh]), MSG_NOTIFY); + {$ENDIF} + if (nx = x0) and (ny = y0) and (nw = w) and (nh = h) then exit; + // map -> grid + Dec(x0, mMinX); + Dec(y0, mMinY); + Dec(nx, mMinX); + Dec(ny, mMinY); + // did any corner crossed tile boundary? + if (x0 div mTileSize <> nx div mTileSize) or + (y0 div mTileSize <> ny div mTileSize) or + ((x0+w) div mTileSize <> (nx+nw) div mTileSize) or + ((y0+h) div mTileSize <> (ny+nh) div mTileSize) then + begin + removeInternal(body); + px.mX := nx+mMinX; + px.mY := ny+mMinY; + px.mWidth := nw; + px.mHeight := nh; + insertInternal(body); + end + else + begin + px.mX := nx+mMinX; + px.mY := ny+mMinY; + px.mWidth := nw; + px.mHeight := nh; + end; end; - -function TBodyGridBase.iterator (grida: Integer): Boolean; +//TODO: optimize for horizontal/vertical moves +procedure TBodyGridBase.moveBody (body: TBodyProxyId; nx, ny: Integer); var - idx: Integer; px: PBodyProxyRec; - {$IFDEF grid_use_buckets} - pi: PGridCell; - f: Integer; - {$ENDIF} + x0, y0: Integer; + ogx0, ogx1, ogy0, ogy1: Integer; // old grid rect + ngx0, ngx1, ngy0, ngy1: Integer; // new grid rect + gx, gy: Integer; + gw, gh: Integer; + pw, ph: Integer; begin - result := false; - idx := mGrid[grida]; - while (idx >= 0) do + if (body < 0) or (body > High(mProxies)) then exit; // just in case + // check if tile coords was changed + px := @mProxies[body]; + x0 := px.mX; + y0 := px.mY; + if (nx = x0) and (ny = y0) then exit; + // map -> grid + Dec(x0, mMinX); + Dec(y0, mMinY); + Dec(nx, mMinX); + Dec(ny, mMinY); + // check for heavy work + pw := px.mWidth; + ph := px.mHeight; + ogx0 := x0 div mTileSize; + ogy0 := y0 div mTileSize; + ngx0 := nx div mTileSize; + ngy0 := ny div mTileSize; + ogx1 := (x0+pw-1) div mTileSize; + ogy1 := (y0+ph-1) div mTileSize; + ngx1 := (nx+pw-1) div mTileSize; + ngy1 := (ny+ph-1) div mTileSize; + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format('proxy #%d: checkmove: xg=%d;yg=%d;w=%d;h=%d;nx=%d;ny=%d og:(%d,%d)-(%d,%d); ng:(%d,%d)-(%d,%d)', [body, x0, y0, pw, ph, nx, ny, ogx0, ogy0, ogx1, ogy1, ngx0, ngy0, ngx1, ngy1]), MSG_NOTIFY); + {$ENDIF} + if (ogx0 <> ngx0) or (ogy0 <> ngy0) or (ogx1 <> ngx1) or (ogy1 <> ngy1) then begin - {$IFDEF grid_use_buckets} - pi := @mCells[idx]; - for f := 0 to High(TGridCell.bodies) do - begin - if (pi.bodies[f] = -1) then break; - px := @mProxies[pi.bodies[f]]; - if (px.mQueryMark <> mLastQuery) and ((mTagMask = -1) or ((px.mTag and mTagMask) <> 0)) then - begin - //e_WriteLog(Format(' query #%d body hit: (%d,%d)-(%dx%d) tag:%d', [mLastQuery, mCells[idx].body.mX, mCells[idx].body.mY, mCells[idx].body.mWidth, mCells[idx].body.mHeight, mCells[idx].body.mTag]), MSG_NOTIFY); - px.mQueryMark := mLastQuery; - if (mItCB(px.mObj, px.mTag)) then begin result := true; exit; end; - end; - end; - idx := pi.next; - {$ELSE} - if (mCells[idx].body <> -1) then + // crossed tile boundary, do heavy work + gw := mWidth; + gh := mHeight; + // cycle with old rect, remove body where it is necessary + // optimized for horizontal moves + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format('proxy #%d: xg=%d;yg=%d;w=%d;h=%d;nx=%d;ny=%d og:(%d,%d)-(%d,%d); ng:(%d,%d)-(%d,%d)', [body, x0, y0, pw, ph, nx, ny, ogx0, ogy0, ogx1, ogy1, ngx0, ngy0, ngx1, ngy1]), MSG_NOTIFY); + {$ENDIF} + // remove stale marks + if not ((ogy0 >= gh) or (ogy1 < 0)) and + not ((ogx0 >= gw) or (ogx1 < 0)) then begin - px := @mProxies[mCells[idx].body]; - if (px.mQueryMark <> mLastQuery) and ((mTagMask = -1) or ((px.mTag and mTagMask) <> 0)) then + if (ogx0 < 0) then ogx0 := 0; + if (ogy0 < 0) then ogy0 := 0; + if (ogx1 > gw-1) then ogx1 := gw-1; + if (ogy1 > gh-1) then ogy1 := gh-1; + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format(' norm og:(%d,%d)-(%d,%d)', [ogx0, ogy0, ogx1, ogy1]), MSG_NOTIFY); + {$ENDIF} + for gx := ogx0 to ogx1 do begin - //e_WriteLog(Format(' query #%d body hit: (%d,%d)-(%dx%d) tag:%d', [mLastQuery, mCells[idx].body.mX, mCells[idx].body.mY, mCells[idx].body.mWidth, mCells[idx].body.mHeight, mCells[idx].body.mTag]), MSG_NOTIFY); - px.mQueryMark := mLastQuery; - if (mItCB(px.mObj, px.mTag)) then begin result := true; exit; end; - end; - end; - idx := mCells[idx].next; - {$ENDIF} + if (gx < ngx0) or (gx > ngx1) then + begin + // this column is completely outside of new rect + for gy := ogy0 to ogy1 do + begin + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format(' remove0:(%d,%d)', [gx, gy]), MSG_NOTIFY); + {$ENDIF} + remover(gy*gw+gx, body); + end; + end + else + begin + // heavy checks + for gy := ogy0 to ogy1 do + begin + if (gy < ngy0) or (gy > ngy1) then + begin + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format(' remove1:(%d,%d)', [gx, gy]), MSG_NOTIFY); + {$ENDIF} + remover(gy*gw+gx, body); + end; + end; + end; + end; + end; + // cycle with new rect, add body where it is necessary + if not ((ngy0 >= gh) or (ngy1 < 0)) and + not ((ngx0 >= gw) or (ngx1 < 0)) then + begin + if (ngx0 < 0) then ngx0 := 0; + if (ngy0 < 0) then ngy0 := 0; + if (ngx1 > gw-1) then ngx1 := gw-1; + if (ngy1 > gh-1) then ngy1 := gh-1; + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format(' norm ng:(%d,%d)-(%d,%d)', [ngx0, ngy0, ngx1, ngy1]), MSG_NOTIFY); + {$ENDIF} + for gx := ngx0 to ngx1 do + begin + if (gx < ogx0) or (gx > ogx1) then + begin + // this column is completely outside of old rect + for gy := ngy0 to ngy1 do + begin + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format(' insert0:(%d,%d)', [gx, gy]), MSG_NOTIFY); + {$ENDIF} + inserter(gy*gw+gx, body); + end; + end + else + begin + // heavy checks + for gy := ngy0 to ngy1 do + begin + if (gy < ogy0) or (gy > ogy1) then + begin + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format(' insert1:(%d,%d)', [gx, gy]), MSG_NOTIFY); + {$ENDIF} + inserter(gy*gw+gx, body); + end; + end; + end; + end; + end; + // done + end + else + begin + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format('proxy #%d: GRID OK: xg=%d;yg=%d;w=%d;h=%d;nx=%d;ny=%d og:(%d,%d)-(%d,%d); ng:(%d,%d)-(%d,%d)', [body, x0, y0, pw, ph, nx, ny, ogx0, ogy0, ogx1, ogy1, ngx0, ngy0, ngx1, ngy1]), MSG_NOTIFY); + {$ENDIF} + end; + // update coordinates + px.mX := nx+mMinX; + px.mY := ny+mMinY; +end; + +procedure TBodyGridBase.resizeBody (body: TBodyProxyId; nw, nh: Integer); +var + px: PBodyProxyRec; + x0, y0, w, h: Integer; +begin + if (body < 0) or (body > High(mProxies)) then exit; // just in case + // check if tile coords was changed + px := @mProxies[body]; + x0 := px.mX-mMinX; + y0 := px.mY-mMinY; + w := px.mWidth; + h := px.mHeight; + {$IF DEFINED(D2F_DEBUG_MOVER)} + e_WriteLog(Format('proxy #%d: RESIZE: xg=%d;yg=%d;w=%d;h=%d;nw=%d;nh=%d', [body, x0, y0, w, h, nw, nh]), MSG_NOTIFY); + {$ENDIF} + if ((x0+w) div mTileSize <> (x0+nw) div mTileSize) or + ((y0+h) div mTileSize <> (y0+nh) div mTileSize) then + begin + // crossed tile boundary, do heavy work + removeInternal(body); + px.mWidth := nw; + px.mHeight := nh; + insertInternal(body); + end + else + begin + // nothing to do with the grid, just fix size + px.mWidth := nw; + px.mHeight := nh; + end; +end; + + +// ////////////////////////////////////////////////////////////////////////// // +function TBodyGridBase.atCellInPoint (x, y: Integer): TAtPointEnumerator; +var + cidx: Integer = -1; +begin + Dec(x, mMinX); + Dec(y, mMinY); + if (x >= 0) and (y >= 0) and (x < mWidth*mTileSize) and (y < mHeight*mTileSize) then cidx := mGrid[(y div mTileSize)*mWidth+(x div mTileSize)]; + result := TAtPointEnumerator.Create(mCells, cidx, getProxyById); +end; + + +// ////////////////////////////////////////////////////////////////////////// // +// no callback: return `true` on the first hit +function TBodyGridBase.forEachAtPoint (x, y: Integer; cb: TGridQueryCB; tagmask: Integer=-1; exittag: PInteger=nil): ITP; +var + f: Integer; + idx, curci: Integer; + cc: PGridCell = nil; + px: PBodyProxyRec; + lq: LongWord; + ptag: Integer; +begin + result := Default(ITP); + if (exittag <> nil) then exittag^ := 0; + tagmask := tagmask and TagFullMask; + if (tagmask = 0) then exit; + + {$IF DEFINED(D2F_DEBUG_XXQ)} + if (assigned(cb)) then e_WriteLog(Format('0: grid pointquery: (%d,%d)', [x, y]), MSG_NOTIFY); + {$ENDIF} + + // make coords (0,0)-based + Dec(x, mMinX); + Dec(y, mMinY); + if (x < 0) or (y < 0) or (x >= mWidth*mTileSize) or (y >= mHeight*mTileSize) then exit; + + curci := mGrid[(y div mTileSize)*mWidth+(x div mTileSize)]; + + {$IF DEFINED(D2F_DEBUG_XXQ)} + if (assigned(cb)) then e_WriteLog(Format('1: grid pointquery: (%d,%d) (%d,%d) %d', [x, y, (x div mTileSize), (y div mTileSize), curci]), MSG_NOTIFY); + {$ENDIF} + + // restore coords + Inc(x, mMinX); + Inc(y, mMinY); + + // increase query counter + Inc(mLastQuery); + if (mLastQuery = 0) then + begin + // just in case of overflow + mLastQuery := 1; + for idx := 0 to High(mProxies) do mProxies[idx].mQueryMark := 0; + end; + lq := mLastQuery; + + {$IF DEFINED(D2F_DEBUG_XXQ)} + if (assigned(cb)) then e_WriteLog(Format('2: grid pointquery: (%d,%d); lq=%u', [x, y, lq]), MSG_NOTIFY); + {$ENDIF} + + while (curci <> -1) do + begin + {$IF DEFINED(D2F_DEBUG_XXQ)} + if (assigned(cb)) then e_WriteLog(Format(' cell #%d', [curci]), MSG_NOTIFY); + {$ENDIF} + cc := @mCells[curci]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + px := @mProxies[cc.bodies[f]]; + {$IF DEFINED(D2F_DEBUG_XXQ)} + if (assigned(cb)) then e_WriteLog(Format(' proxy #%d; qm:%u; tag:%08x; tagflag:%d %u', [cc.bodies[f], px.mQueryMark, px.mTag, (px.mTag and tagmask), LongWord(px.mObj)]), MSG_NOTIFY); + {$ENDIF} + // shit. has to do it this way, so i can change tag in callback + if (px.mQueryMark <> lq) then + begin + px.mQueryMark := lq; + ptag := px.mTag; + if ((ptag and TagDisabled) = 0) and ((ptag and tagmask) <> 0) and + (x >= px.mX) and (y >= px.mY) and (x < px.mX+px.mWidth) and (y < px.mY+px.mHeight) then + begin + if assigned(cb) then + begin + if cb(px.mObj, ptag) then + begin + result := px.mObj; + if (exittag <> nil) then exittag^ := ptag; + exit; + end; + end + else + begin + result := px.mObj; + if (exittag <> nil) then exittag^ := ptag; + exit; + end; + end; + end; + end; + curci := cc.next; end; end; -function TBodyGridBase.forEachInAABB (x, y, w, h: Integer; cb: TGridQueryCB; tagmask: Integer=-1): Boolean; + +// ////////////////////////////////////////////////////////////////////////// // +// no callback: return `true` on the first hit +function TBodyGridBase.forEachInAABB (x, y, w, h: Integer; cb: TGridQueryCB; tagmask: Integer=-1; allowDisabled: Boolean=false): ITP; +const + tsize = mTileSize; var idx: Integer; - otagmask: Integer; - ocb: TGridQueryCB; + gx, gy: Integer; + curci: Integer; + f: Integer; + cc: PGridCell = nil; + px: PBodyProxyRec; + lq: LongWord; + gw: Integer; + x0, y0: Integer; + ptag: Integer; begin - result := false; - if not assigned(cb) then exit; + result := Default(ITP); + if (w < 1) or (h < 1) then exit; + tagmask := tagmask and TagFullMask; + if (tagmask = 0) then exit; + + x0 := x; + y0 := y; + + // fix coords + Dec(x, mMinX); + Dec(y, mMinY); + + gw := mWidth; + //tsize := mTileSize; + + if (x+w <= 0) or (y+h <= 0) then exit; + if (x >= gw*tsize) or (y >= mHeight*tsize) then exit; + + if mInQuery then raise Exception.Create('recursive queries aren''t supported'); + mInQuery := true; // increase query counter Inc(mLastQuery); @@ -543,14 +1332,1244 @@ begin for idx := 0 to High(mProxies) do mProxies[idx].mQueryMark := 0; end; //e_WriteLog(Format('grid: query #%d: (%d,%d)-(%dx%d)', [mLastQuery, minx, miny, maxx, maxy]), MSG_NOTIFY); + lq := mLastQuery; + + // go on + for gy := y div tsize to (y+h-1) div tsize do + begin + if (gy < 0) then continue; + if (gy >= mHeight) then break; + for gx := x div tsize to (x+w-1) div tsize do + begin + if (gx < 0) then continue; + if (gx >= gw) then break; + // process cells + curci := mGrid[gy*gw+gx]; + while (curci <> -1) do + begin + cc := @mCells[curci]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + px := @mProxies[cc.bodies[f]]; + // shit. has to do it this way, so i can change tag in callback + if (px.mQueryMark = lq) then continue; + px.mQueryMark := lq; + ptag := px.mTag; + if (not allowDisabled) and ((ptag and TagDisabled) <> 0) then continue; + if ((ptag and tagmask) = 0) then continue; + if (x0 >= px.mX+px.mWidth) or (y0 >= px.mY+px.mHeight) then continue; + if (x0+w <= px.mX) or (y0+h <= px.mY) then continue; + if assigned(cb) then + begin + if cb(px.mObj, ptag) then begin result := px.mObj; mInQuery := false; exit; end; + end + else + begin + result := px.mObj; + mInQuery := false; + exit; + end; + end; + curci := cc.next; + end; + end; + end; + + mInQuery := false; +end; + + +// ////////////////////////////////////////////////////////////////////////// // +// no callback: return `true` on the nearest hit +function TBodyGridBase.traceRay (const x0, y0, x1, y1: Integer; cb: TGridRayQueryCB; tagmask: Integer=-1): ITP; +var + ex, ey: Integer; +begin + result := traceRay(ex, ey, x0, y0, x1, y1, cb, tagmask); +end; + + +// no callback: return `true` on the nearest hit +// you are not supposed to understand this +function TBodyGridBase.traceRay (out ex, ey: Integer; const ax0, ay0, ax1, ay1: Integer; cb: TGridRayQueryCB; tagmask: Integer=-1): ITP; +const + tsize = mTileSize; +var + wx0, wy0, wx1, wy1: Integer; // window coordinates + stx, sty: Integer; // "steps" for x and y axes + dsx, dsy: Integer; // "lengthes" for x and y axes + dx2, dy2: Integer; // "double lengthes" for x and y axes + xd, yd: Integer; // current coord + e: Integer; // "error" (as in bresenham algo) + rem: Integer; + term: Integer; + xptr, yptr: PInteger; + xfixed: Boolean; + temp: Integer; + prevx, prevy: Integer; + lastDistSq: Integer; + ccidx, curci: Integer; + hasUntried: Boolean; + lastGA: Integer = -1; + ga, x, y: Integer; + lastObj: ITP; + wasHit: Boolean = false; + gw, gh, minx, miny, maxx, maxy: Integer; + cc: PGridCell; + px: PBodyProxyRec; + lq: LongWord; + f, ptag, distSq: Integer; + x0, y0, x1, y1: Integer; + //swapped: Boolean = false; // true: xd is yd, and vice versa + // horizontal walker + {$IFDEF GRID_USE_ORTHO_ACCEL} + wklen, wkstep: Integer; + //wksign: Integer; + hopt: Boolean; + {$ENDIF} + // skipper + xdist, ydist: Integer; +begin + result := Default(ITP); + lastObj := Default(ITP); + tagmask := tagmask and TagFullMask; + ex := ax1; // why not? + ey := ay1; // why not? + if (tagmask = 0) then exit; + + if (ax0 = ax1) and (ay0 = ay1) then + begin + result := forEachAtPoint(ax0, ay0, nil, tagmask, @ptag); + if (result <> nil) then + begin + if assigned(cb) and not cb(result, ptag, ax0, ay0, ax0, ay0) then result := Default(ITP); + end; + exit; + end; + + lastDistSq := distanceSq(ax0, ay0, ax1, ay1)+1; + + gw := mWidth; + gh := mHeight; + minx := mMinX; + miny := mMinY; + maxx := gw*tsize-1; + maxy := gh*tsize-1; + + {$IF DEFINED(D2F_DEBUG_RAYTRACE)} + if assigned(dbgRayTraceTileHitCB) then e_WriteLog(Format('TRACING: (%d,%d)-(%d,%d) [(%d,%d)-(%d,%d)]; maxdistsq=%d', [ax0, ay0, ax1, ay1, minx, miny, maxx, maxy, lastDistSq]), MSG_NOTIFY); + {$ENDIF} + + x0 := ax0; + y0 := ay0; + x1 := ax1; + y1 := ay1; + + // offset query coords to (0,0)-based + Dec(x0, minx); + Dec(y0, miny); + Dec(x1, minx); + Dec(y1, miny); + + // clip rectange + wx0 := 0; + wy0 := 0; + wx1 := maxx; + wy1 := maxy; + + // horizontal setup + if (x0 < x1) then + begin + // from left to right + if (x0 > wx1) or (x1 < wx0) then exit; // out of screen + stx := 1; // going right + end + else + begin + // from right to left + if (x1 > wx1) or (x0 < wx0) then exit; // out of screen + stx := -1; // going left + x0 := -x0; + x1 := -x1; + wx0 := -wx0; + wx1 := -wx1; + swapInt(wx0, wx1); + end; + + // vertical setup + if (y0 < y1) then + begin + // from top to bottom + if (y0 > wy1) or (y1 < wy0) then exit; // out of screen + sty := 1; // going down + end + else + begin + // from bottom to top + if (y1 > wy1) or (y0 < wy0) then exit; // out of screen + sty := -1; // going up + y0 := -y0; + y1 := -y1; + wy0 := -wy0; + wy1 := -wy1; + swapInt(wy0, wy1); + end; - otagmask := mTagMask; - mTagMask := tagmask; - ocb := mItCB; - mItCB := cb; - result := forGridRect(x, y, w, h, iterator); - mTagMask := otagmask; - mItCB := ocb; + dsx := x1-x0; + dsy := y1-y0; + + if (dsx < dsy) then + begin + //swapped := true; + xptr := @yd; + yptr := @xd; + swapInt(x0, y0); + swapInt(x1, y1); + swapInt(dsx, dsy); + swapInt(wx0, wy0); + swapInt(wx1, wy1); + swapInt(stx, sty); + end + else + begin + xptr := @xd; + yptr := @yd; + end; + + dx2 := 2*dsx; + dy2 := 2*dsy; + xd := x0; + yd := y0; + e := 2*dsy-dsx; + term := x1; + + xfixed := false; + if (y0 < wy0) then + begin + // clip at top + temp := dx2*(wy0-y0)-dsx; + xd += temp div dy2; + rem := temp mod dy2; + if (xd > wx1) then exit; // x is moved out of clipping rect, nothing to do + if (xd+1 >= wx0) then + begin + yd := wy0; + e -= rem+dsx; + if (rem > 0) then begin Inc(xd); e += dy2; end; + xfixed := true; + end; + end; + + if (not xfixed) and (x0 < wx0) then + begin + // clip at left + temp := dy2*(wx0-x0); + yd += temp div dx2; + rem := temp mod dx2; + if (yd > wy1) or (yd = wy1) and (rem >= dsx) then exit; + xd := wx0; + e += rem; + if (rem >= dsx) then begin Inc(yd); e -= dx2; end; + end; + + if (y1 > wy1) then + begin + // clip at bottom + temp := dx2*(wy1-y0)+dsx; + term := x0+temp div dy2; + rem := temp mod dy2; + if (rem = 0) then Dec(term); + end; + + if (term > wx1) then term := wx1; // clip at right + + Inc(term); // draw last point + //if (term = xd) then exit; // this is the only point, get out of here + + if (sty = -1) then yd := -yd; + if (stx = -1) then begin xd := -xd; term := -term; end; + dx2 -= dy2; + + // first move, to skip starting point + // DON'T DO THIS! loop will take care of that + if (xd = term) then + begin + //FIXME! + result := forEachAtPoint(ax0, ay0, nil, tagmask, @ptag); + if (result <> nil) then + begin + if assigned(cb) then + begin + if cb(result, ptag, ax0, ay0, ax0, ay0) then + begin + ex := ax0; + ey := ay0; + end + else + begin + result := nil; + end; + end + else + begin + ex := ax0; + ey := ay0; + end; + end; + exit; + end; + + prevx := xptr^+minx; + prevy := yptr^+miny; + (* + // move coords + if (e >= 0) then begin yd += sty; e -= dx2; end else e += dy2; + xd += stx; + // done? + if (xd = term) then exit; + *) + + {$IF DEFINED(D2F_DEBUG)} + if (xptr^ < 0) or (yptr^ < 0) or (xptr^ >= gw*tsize) and (yptr^ >= gh*tsize) then raise Exception.Create('raycaster internal error (0)'); + {$ENDIF} + // DON'T DO THIS! loop will take care of that + //lastGA := (yptr^ div tsize)*gw+(xptr^ div tsize); + //ccidx := mGrid[lastGA]; + + {$IF DEFINED(D2F_DEBUG_RAYTRACE)} + //if assigned(dbgRayTraceTileHitCB) then e_WriteLog('1:TRACING!', MSG_NOTIFY); + {$ENDIF} + + //if (dbgShowTraceLog) then e_WriteLog(Format('raycast start: (%d,%d)-(%d,%d); xptr^=%d; yptr^=%d', [ax0, ay0, ax1, ay1, xptr^, yptr^]), MSG_NOTIFY); + + if mInQuery then raise Exception.Create('recursive queries aren''t supported'); + mInQuery := true; + + // increase query counter + Inc(mLastQuery); + if (mLastQuery = 0) then + begin + // just in case of overflow + mLastQuery := 1; + for f := 0 to High(mProxies) do mProxies[f].mQueryMark := 0; + end; + lq := mLastQuery; + + {$IFDEF GRID_USE_ORTHO_ACCEL} + // if this is strict horizontal/vertical trace, use optimized codepath + if (ax0 = ax1) or (ay0 = ay1) then + begin + // horizontal trace: walk the whole tiles, calculating mindist once for each proxy in cell + // stx < 0: going left, otherwise `stx` is > 0, and we're going right + // vertical trace: walk the whole tiles, calculating mindist once for each proxy in cell + // stx < 0: going up, otherwise `stx` is > 0, and we're going down + hopt := (ay0 = ay1); // horizontal? + if (stx < 0) then begin {wksign := -1;} wklen := -(term-xd); end else begin {wksign := 1;} wklen := term-xd; end; + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln('optimized htrace; wklen=%d', [wklen]); + {$ENDIF} + ga := (yptr^ div tsize)*gw+(xptr^ div tsize); + // one of those will never change + x := xptr^+minx; + y := yptr^+miny; + while (wklen > 0) do + begin + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' htrace; ga=%d; x=%d, y=%d; y=%d; y=%d', [ga, xptr^+minx, yptr^+miny, y, ay0]); + {$ENDIF} + // new tile? + if (ga <> lastGA) then + begin + lastGA := ga; + ccidx := mGrid[lastGA]; + // convert coords to map (to avoid ajdusting coords inside the loop) + if hopt then x := xptr^+minx else y := yptr^+miny; + while (ccidx <> -1) do + begin + cc := @mCells[ccidx]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + px := @mProxies[cc.bodies[f]]; + ptag := px.mTag; + if ((ptag and TagDisabled) = 0) and ((ptag and tagmask) <> 0) and (px.mQueryMark <> lq) and + // constant coord should be inside + ((hopt and (y >= px.y0) and (y <= px.y1)) or + ((not hopt) and (x >= px.x0) and (x <= px.x1))) then + begin + px.mQueryMark := lq; // mark as processed + // inside the proxy? + if (hopt and (x > px.x0) and (x < px.x1)) or + ((not hopt) and (y > px.y0) and (y < px.y1)) then + begin + // setup prev[xy] + if assigned(cb) then + begin + if cb(px.mObj, ptag, x, y, x, y) then + begin + result := px.mObj; + ex := x; + ey := y; + mInQuery := false; + exit; + end; + end + else + begin + distSq := distanceSq(ax0, ay0, x, y); + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' EMBEDDED hhit(%d): a=(%d,%d), h=(%d,%d), distsq=%d; lastsq=%d', [cc.bodies[f], ax0, ay0, x, y, distSq, lastDistSq]); + {$ENDIF} + if (distSq < lastDistSq) then + begin + ex := x; + ey := y; + result := px.mObj; + mInQuery := false; + exit; + end; + end; + continue; + end; + // remember this hitpoint if it is nearer than an old one + // setup prev[xy] + if hopt then + begin + // horizontal trace + prevy := y; + y := yptr^+miny; + if (stx < 0) then + begin + // going left + if (x < px.x1) then continue; // not on the right edge + x := px.x1; + prevx := x+1; + end + else + begin + // going right + if (x > px.x0) then continue; // not on the left edge + x := px.x0; + prevx := x-1; + end; + end + else + begin + // vertical trace + prevx := x; + x := xptr^+minx; + if (stx < 0) then + begin + // going up + if (y < px.y1) then continue; // not on the bottom edge + y := px.y1; + prevy := x+1; + end + else + begin + // going down + if (y > px.y0) then continue; // not on the top edge + y := px.y0; + prevy := y-1; + end; + end; + if assigned(cb) then + begin + if cb(px.mObj, ptag, x, y, prevx, prevy) then + begin + result := px.mObj; + ex := prevx; + ey := prevy; + mInQuery := false; + exit; + end; + end + else + begin + distSq := distanceSq(ax0, ay0, prevx, prevy); + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' hhit(%d): a=(%d,%d), h=(%d,%d), p=(%d,%d), distsq=%d; lastsq=%d', [cc.bodies[f], ax0, ay0, x, y, prevx, prevy, distSq, lastDistSq]); + {$ENDIF} + if (distSq < lastDistSq) then + begin + wasHit := true; + lastDistSq := distSq; + ex := prevx; + ey := prevy; + lastObj := px.mObj; + end; + end; + end; + end; + // next cell + ccidx := cc.next; + end; + if wasHit and not assigned(cb) then begin result := lastObj; mInQuery := false; exit; end; + if assigned(cb) and cb(nil, 0, x, y, x, y) then begin result := lastObj; mInQuery := false; exit; end; + end; + // skip to next tile + if hopt then + begin + if (stx > 0) then + begin + // to the right + wkstep := ((xptr^ or (mTileSize-1))+1)-xptr^; + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' right step: wklen=%d; wkstep=%d', [wklen, wkstep]); + {$ENDIF} + if (wkstep >= wklen) then break; + Inc(xptr^, wkstep); + Inc(ga); + end + else + begin + // to the left + wkstep := xptr^-((xptr^ and (not (mTileSize-1)))-1); + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' left step: wklen=%d; wkstep=%d', [wklen, wkstep]); + {$ENDIF} + if (wkstep >= wklen) then break; + Dec(xptr^, wkstep); + Dec(ga); + end; + end + else + begin + if (stx > 0) then + begin + // to the down + wkstep := ((yptr^ or (mTileSize-1))+1)-yptr^; + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' down step: wklen=%d; wkstep=%d', [wklen, wkstep]); + {$ENDIF} + if (wkstep >= wklen) then break; + Inc(yptr^, wkstep); + Inc(ga, mWidth); + end + else + begin + // to the up + wkstep := yptr^-((yptr^ and (not (mTileSize-1)))-1); + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' up step: wklen=%d; wkstep=%d', [wklen, wkstep]); + {$ENDIF} + if (wkstep >= wklen) then break; + Dec(yptr^, wkstep); + Dec(ga, mWidth); + end; + end; + Dec(wklen, wkstep); + end; + // we can travel less than one cell + if wasHit and not assigned(cb) then result := lastObj else begin ex := ax1; ey := ay1; end; + mInQuery := false; + exit; + end; + {$ENDIF} + + {$IF DEFINED(D2F_DEBUG_RAYTRACE)} + if assigned(dbgRayTraceTileHitCB) then dbgRayTraceTileHitCB((xptr^ div tsize*tsize)+minx, (yptr^ div tsize*tsize)+miny); + {$ENDIF} + + //e_LogWritefln('*********************', []); + ccidx := -1; + // can omit checks + while (xd <> term) do + begin + // check cell(s) + {$IF DEFINED(D2F_DEBUG)} + if (xptr^ < 0) or (yptr^ < 0) or (xptr^ >= gw*tsize) and (yptr^ >= gh*tsize) then raise Exception.Create('raycaster internal error (0)'); + {$ENDIF} + // new tile? + ga := (yptr^ div tsize)*gw+(xptr^ div tsize); + {$IF DEFINED(D2F_DEBUG_RAYTRACE)} + if assigned(dbgRayTraceTileHitCB) then e_WriteLog(Format(' xd=%d; term=%d; gx=%d; gy=%d; ga=%d; lastga=%d', [xd, term, xptr^, yptr^, ga, lastGA]), MSG_NOTIFY); + {$ENDIF} + if (ga <> lastGA) then + begin + // yes + {$IF DEFINED(D2F_DEBUG)} + if assigned(dbgRayTraceTileHitCB) then dbgRayTraceTileHitCB((xptr^ div tsize*tsize)+minx, (yptr^ div tsize*tsize)+miny); + {$ENDIF} + if (ccidx <> -1) then + begin + // signal cell completion + if assigned(cb) then + begin + if cb(nil, 0, xptr^+minx, yptr^+miny, prevx, prevy) then begin result := lastObj; mInQuery := false; exit; end; + end + else if wasHit then + begin + result := lastObj; + mInQuery := false; + exit; + end; + end; + lastGA := ga; + ccidx := mGrid[lastGA]; + end; + // has something to process in this tile? + if (ccidx <> -1) then + begin + // process cell + curci := ccidx; + hasUntried := false; // this will be set to `true` if we have some proxies we still want to process at the next step + // convert coords to map (to avoid ajdusting coords inside the loop) + x := xptr^+minx; + y := yptr^+miny; + // process cell list + while (curci <> -1) do + begin + cc := @mCells[curci]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + px := @mProxies[cc.bodies[f]]; + ptag := px.mTag; + if ((ptag and TagDisabled) = 0) and ((ptag and tagmask) <> 0) and (px.mQueryMark <> lq) then + begin + // can we process this proxy? + if (x >= px.mX) and (y >= px.mY) and (x < px.mX+px.mWidth) and (y < px.mY+px.mHeight) then + begin + px.mQueryMark := lq; // mark as processed + if assigned(cb) then + begin + if cb(px.mObj, ptag, x, y, prevx, prevy) then + begin + result := px.mObj; + ex := prevx; + ey := prevy; + mInQuery := false; + exit; + end; + end + else + begin + // remember this hitpoint if it is nearer than an old one + distSq := distanceSq(ax0, ay0, prevx, prevy); + {$IF DEFINED(D2F_DEBUG_RAYTRACE)} + if assigned(dbgRayTraceTileHitCB) then e_WriteLog(Format(' hit(%d): a=(%d,%d), h=(%d,%d), p=(%d,%d); distsq=%d; lastsq=%d', [cc.bodies[f], ax0, ay0, x, y, prevx, prevy, distSq, lastDistSq]), MSG_NOTIFY); + {$ENDIF} + if (distSq < lastDistSq) then + begin + wasHit := true; + lastDistSq := distSq; + ex := prevx; + ey := prevy; + lastObj := px.mObj; + end; + end; + end + else + begin + // this is possibly interesting proxy, set "has more to check" flag + hasUntried := true; + end; + end; + end; + // next cell + curci := cc.next; + end; + // still has something interesting in this cell? + if not hasUntried then + begin + // nope, don't process this cell anymore; signal cell completion + ccidx := -1; + if assigned(cb) then + begin + if cb(nil, 0, x, y, prevx, prevy) then begin result := lastObj; mInQuery := false; exit; end; + end + else if wasHit then + begin + result := lastObj; + mInQuery := false; + exit; + end; + end; + end; + if (ccidx = -1) then + begin + // move to cell edge, as we have nothing to trace here anymore + if (stx < 0) then xdist := xd and (not (mTileSize-1)) else xdist := xd or (mTileSize-1); + if (sty < 0) then ydist := yd and (not (mTileSize-1)) else ydist := yd or (mTileSize-1); + //e_LogWritefln('0: swapped=%d; xd=%d; yd=%d; stx=%d; sty=%d; e=%d; dx2=%d; dy2=%d; term=%d; xdist=%d; ydist=%d', [swapped, xd, yd, stx, sty, e, dx2, dy2, term, xdist, ydist]); + while (xd <> xdist) and (yd <> ydist) do + begin + // step + xd += stx; + if (e >= 0) then begin yd += sty; e -= dx2; end else e += dy2; + //e_LogWritefln(' xd=%d; yd=%d', [xd, yd]); + if (xd = term) then break; + end; + //e_LogWritefln('1: swapped=%d; xd=%d; yd=%d; stx=%d; sty=%d; e=%d; dx2=%d; dy2=%d; term=%d; xdist=%d; ydist=%d', [swapped, xd, yd, stx, sty, e, dx2, dy2, term, xdist, ydist]); + if (xd = term) then break; + end; + //putPixel(xptr^, yptr^); + // move coords + prevx := xptr^+minx; + prevy := yptr^+miny; + if (e >= 0) then begin yd += sty; e -= dx2; end else e += dy2; + xd += stx; + end; + // we can travel less than one cell + if wasHit and not assigned(cb) then + begin + result := lastObj; + end + else + begin + ex := ax1; // why not? + ey := ay1; // why not? + end; + + mInQuery := false; +end; + + +// ////////////////////////////////////////////////////////////////////////// // +//FIXME! optimize this with real tile walking +function TBodyGridBase.forEachAlongLine (ax0, ay0, ax1, ay1: Integer; cb: TGridQueryCB; tagmask: Integer=-1; log: Boolean=false): ITP; +const + tsize = mTileSize; +var + wx0, wy0, wx1, wy1: Integer; // window coordinates + stx, sty: Integer; // "steps" for x and y axes + dsx, dsy: Integer; // "lengthes" for x and y axes + dx2, dy2: Integer; // "double lengthes" for x and y axes + xd, yd: Integer; // current coord + e: Integer; // "error" (as in bresenham algo) + rem: Integer; + term: Integer; + xptr, yptr: PInteger; + xfixed: Boolean; + temp: Integer; + ccidx, curci: Integer; + lastGA: Integer = -1; + ga: Integer; + gw, gh, minx, miny, maxx, maxy: Integer; + cc: PGridCell; + px: PBodyProxyRec; + lq: LongWord; + f, ptag: Integer; + x0, y0, x1, y1: Integer; + //swapped: Boolean = false; // true: xd is yd, and vice versa + // horizontal walker + {$IFDEF GRID_USE_ORTHO_ACCEL} + wklen, wkstep: Integer; + //wksign: Integer; + hopt: Boolean; + {$ENDIF} + // skipper + xdist, ydist: Integer; +begin + log := false; + result := Default(ITP); + tagmask := tagmask and TagFullMask; + if (tagmask = 0) or not assigned(cb) then exit; + + if (ax0 = ax1) and (ay0 = ay1) then + begin + result := forEachAtPoint(ax0, ay0, cb, tagmask, @ptag); + exit; + end; + + gw := mWidth; + gh := mHeight; + minx := mMinX; + miny := mMinY; + maxx := gw*tsize-1; + maxy := gh*tsize-1; + + x0 := ax0; + y0 := ay0; + x1 := ax1; + y1 := ay1; + + // offset query coords to (0,0)-based + Dec(x0, minx); + Dec(y0, miny); + Dec(x1, minx); + Dec(y1, miny); + + // clip rectange + wx0 := 0; + wy0 := 0; + wx1 := maxx; + wy1 := maxy; + + // horizontal setup + if (x0 < x1) then + begin + // from left to right + if (x0 > wx1) or (x1 < wx0) then exit; // out of screen + stx := 1; // going right + end + else + begin + // from right to left + if (x1 > wx1) or (x0 < wx0) then exit; // out of screen + stx := -1; // going left + x0 := -x0; + x1 := -x1; + wx0 := -wx0; + wx1 := -wx1; + swapInt(wx0, wx1); + end; + + // vertical setup + if (y0 < y1) then + begin + // from top to bottom + if (y0 > wy1) or (y1 < wy0) then exit; // out of screen + sty := 1; // going down + end + else + begin + // from bottom to top + if (y1 > wy1) or (y0 < wy0) then exit; // out of screen + sty := -1; // going up + y0 := -y0; + y1 := -y1; + wy0 := -wy0; + wy1 := -wy1; + swapInt(wy0, wy1); + end; + + dsx := x1-x0; + dsy := y1-y0; + + if (dsx < dsy) then + begin + //swapped := true; + xptr := @yd; + yptr := @xd; + swapInt(x0, y0); + swapInt(x1, y1); + swapInt(dsx, dsy); + swapInt(wx0, wy0); + swapInt(wx1, wy1); + swapInt(stx, sty); + end + else + begin + xptr := @xd; + yptr := @yd; + end; + + dx2 := 2*dsx; + dy2 := 2*dsy; + xd := x0; + yd := y0; + e := 2*dsy-dsx; + term := x1; + + xfixed := false; + if (y0 < wy0) then + begin + // clip at top + temp := dx2*(wy0-y0)-dsx; + xd += temp div dy2; + rem := temp mod dy2; + if (xd > wx1) then exit; // x is moved out of clipping rect, nothing to do + if (xd+1 >= wx0) then + begin + yd := wy0; + e -= rem+dsx; + if (rem > 0) then begin Inc(xd); e += dy2; end; + xfixed := true; + end; + end; + + if (not xfixed) and (x0 < wx0) then + begin + // clip at left + temp := dy2*(wx0-x0); + yd += temp div dx2; + rem := temp mod dx2; + if (yd > wy1) or (yd = wy1) and (rem >= dsx) then exit; + xd := wx0; + e += rem; + if (rem >= dsx) then begin Inc(yd); e -= dx2; end; + end; + + if (y1 > wy1) then + begin + // clip at bottom + temp := dx2*(wy1-y0)+dsx; + term := x0+temp div dy2; + rem := temp mod dy2; + if (rem = 0) then Dec(term); + end; + + if (term > wx1) then term := wx1; // clip at right + + Inc(term); // draw last point + //if (term = xd) then exit; // this is the only point, get out of here + + if (sty = -1) then yd := -yd; + if (stx = -1) then begin xd := -xd; term := -term; end; + dx2 -= dy2; + + // first move, to skip starting point + // DON'T DO THIS! loop will take care of that + if (xd = term) then + begin + result := forEachAtPoint(ax0, ay0, cb, tagmask, @ptag); + exit; + end; + + (* + // move coords + if (e >= 0) then begin yd += sty; e -= dx2; end else e += dy2; + xd += stx; + // done? + if (xd = term) then exit; + *) + + {$IF DEFINED(D2F_DEBUG)} + if (xptr^ < 0) or (yptr^ < 0) or (xptr^ >= gw*tsize) and (yptr^ >= gh*tsize) then raise Exception.Create('raycaster internal error (0)'); + {$ENDIF} + // DON'T DO THIS! loop will take care of that + //lastGA := (yptr^ div tsize)*gw+(xptr^ div tsize); + //ccidx := mGrid[lastGA]; + + if mInQuery then raise Exception.Create('recursive queries aren''t supported'); + mInQuery := true; + + // increase query counter + Inc(mLastQuery); + if (mLastQuery = 0) then + begin + // just in case of overflow + mLastQuery := 1; + for f := 0 to High(mProxies) do mProxies[f].mQueryMark := 0; + end; + lq := mLastQuery; + + {$IFDEF GRID_USE_ORTHO_ACCEL} + // if this is strict horizontal/vertical trace, use optimized codepath + if (ax0 = ax1) or (ay0 = ay1) then + begin + // horizontal trace: walk the whole tiles, calculating mindist once for each proxy in cell + // stx < 0: going left, otherwise `stx` is > 0, and we're going right + // vertical trace: walk the whole tiles, calculating mindist once for each proxy in cell + // stx < 0: going up, otherwise `stx` is > 0, and we're going down + hopt := (ay0 = ay1); // horizontal? + if (stx < 0) then begin {wksign := -1;} wklen := -(term-xd); end else begin {wksign := 1;} wklen := term-xd; end; + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln('optimized htrace; wklen=%d', [wklen]); + {$ENDIF} + ga := (yptr^ div tsize)*gw+(xptr^ div tsize); + while (wklen > 0) do + begin + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' htrace; ga=%d; x=%d, y=%d; ay0=%d', [ga, xptr^+minx, yptr^+miny, ay0]); + {$ENDIF} + // new tile? + if (ga <> lastGA) then + begin + lastGA := ga; + ccidx := mGrid[lastGA]; + // convert coords to map (to avoid ajdusting coords inside the loop) + while (ccidx <> -1) do + begin + cc := @mCells[ccidx]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + px := @mProxies[cc.bodies[f]]; + ptag := px.mTag; + if ((ptag and TagDisabled) = 0) and ((ptag and tagmask) <> 0) and (px.mQueryMark <> lq) then + begin + px.mQueryMark := lq; // mark as processed + if assigned(cb) then + begin + if cb(px.mObj, ptag) then begin result := px.mObj; mInQuery := false; exit; end; + end + else + begin + result := px.mObj; + mInQuery := false; + exit; + end; + end; + end; + // next cell + ccidx := cc.next; + end; + end; + // skip to next tile + if hopt then + begin + if (stx > 0) then + begin + // to the right + wkstep := ((xptr^ or (mTileSize-1))+1)-xptr^; + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' right step: wklen=%d; wkstep=%d', [wklen, wkstep]); + {$ENDIF} + if (wkstep >= wklen) then break; + Inc(xptr^, wkstep); + Inc(ga); + end + else + begin + // to the left + wkstep := xptr^-((xptr^ and (not (mTileSize-1)))-1); + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' left step: wklen=%d; wkstep=%d', [wklen, wkstep]); + {$ENDIF} + if (wkstep >= wklen) then break; + Dec(xptr^, wkstep); + Dec(ga); + end; + end + else + begin + if (stx > 0) then + begin + // to the down + wkstep := ((yptr^ or (mTileSize-1))+1)-yptr^; + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' down step: wklen=%d; wkstep=%d', [wklen, wkstep]); + {$ENDIF} + if (wkstep >= wklen) then break; + Inc(yptr^, wkstep); + Inc(ga, mWidth); + end + else + begin + // to the up + wkstep := yptr^-((yptr^ and (not (mTileSize-1)))-1); + {$IF DEFINED(D2F_DEBUG)} + if dbgShowTraceLog then e_LogWritefln(' up step: wklen=%d; wkstep=%d', [wklen, wkstep]); + {$ENDIF} + if (wkstep >= wklen) then break; + Dec(yptr^, wkstep); + Dec(ga, mWidth); + end; + end; + Dec(wklen, wkstep); + end; + mInQuery := false; + exit; + end; + {$ENDIF} + + {$IF DEFINED(D2F_DEBUG_RAYTRACE)} + if assigned(dbgRayTraceTileHitCB) then dbgRayTraceTileHitCB((xptr^ div tsize*tsize)+minx, (yptr^ div tsize*tsize)+miny); + {$ENDIF} + + ccidx := -1; + // can omit checks + while (xd <> term) do + begin + // check cell(s) + {$IF DEFINED(D2F_DEBUG)} + if (xptr^ < 0) or (yptr^ < 0) or (xptr^ >= gw*tsize) and (yptr^ >= gh*tsize) then raise Exception.Create('raycaster internal error (0)'); + {$ENDIF} + // new tile? + ga := (yptr^ div tsize)*gw+(xptr^ div tsize); + {$IF DEFINED(D2F_DEBUG_RAYTRACE)} + if assigned(dbgRayTraceTileHitCB) then e_WriteLog(Format(' xd=%d; term=%d; gx=%d; gy=%d; ga=%d; lastga=%d', [xd, term, xptr^, yptr^, ga, lastGA]), MSG_NOTIFY); + {$ENDIF} + if (ga <> lastGA) then + begin + // yes + {$IF DEFINED(D2F_DEBUG)} + if assigned(dbgRayTraceTileHitCB) then dbgRayTraceTileHitCB((xptr^ div tsize*tsize)+minx, (yptr^ div tsize*tsize)+miny); + {$ENDIF} + lastGA := ga; + ccidx := mGrid[lastGA]; + end; + // has something to process in this tile? + if (ccidx <> -1) then + begin + // process cell + curci := ccidx; + // process cell list + while (curci <> -1) do + begin + cc := @mCells[curci]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + px := @mProxies[cc.bodies[f]]; + ptag := px.mTag; + if ((ptag and TagDisabled) = 0) and ((ptag and tagmask) <> 0) and (px.mQueryMark <> lq) then + begin + px.mQueryMark := lq; // mark as processed + if assigned(cb) then + begin + if cb(px.mObj, ptag) then begin result := px.mObj; mInQuery := false; exit; end; + end + else + begin + result := px.mObj; + mInQuery := false; + exit; + end; + end; + end; + // next cell + curci := cc.next; + end; + // nothing more interesting in this cell + ccidx := -1; + end; + // move to cell edge, as we have nothing to trace here anymore + if (stx < 0) then xdist := xd and (not (mTileSize-1)) else xdist := xd or (mTileSize-1); + if (sty < 0) then ydist := yd and (not (mTileSize-1)) else ydist := yd or (mTileSize-1); + //e_LogWritefln('0: swapped=%d; xd=%d; yd=%d; stx=%d; sty=%d; e=%d; dx2=%d; dy2=%d; term=%d; xdist=%d; ydist=%d', [swapped, xd, yd, stx, sty, e, dx2, dy2, term, xdist, ydist]); + while (xd <> xdist) and (yd <> ydist) do + begin + // step + xd += stx; + if (e >= 0) then begin yd += sty; e -= dx2; end else e += dy2; + //e_LogWritefln(' xd=%d; yd=%d', [xd, yd]); + if (xd = term) then break; + end; + //e_LogWritefln('1: swapped=%d; xd=%d; yd=%d; stx=%d; sty=%d; e=%d; dx2=%d; dy2=%d; term=%d; xdist=%d; ydist=%d', [swapped, xd, yd, stx, sty, e, dx2, dy2, term, xdist, ydist]); + if (xd = term) then break; + //putPixel(xptr^, yptr^); + // move coords + if (e >= 0) then begin yd += sty; e -= dx2; end else e += dy2; + xd += stx; + end; + + mInQuery := false; +end; + + +{.$DEFINE D2F_DEBUG_OTR} +function TBodyGridBase.traceOrthoRayWhileIn (out ex, ey: Integer; ax0, ay0, ax1, ay1: Integer; tagmask: Integer=-1): Boolean; +var + ccidx: Integer; + cc: PGridCell; + px: PBodyProxyRec; + ptag: Integer; + minx, miny: Integer; + f, c0, c1: Integer; + x0, y0, x1, y1: Integer; + celly0, celly1: Integer; + dy: Integer; + filled: array[0..mTileSize-1] of Byte; + {$IF DEFINED(D2F_DEBUG_OTR)} + s: AnsiString = ''; + {$ENDIF} +begin + result := false; + ex := ax1; + ey := ay1; + if not ((ax0 = ax1) or (ay0 = ay1)) then raise Exception.Create('orthoray is not orthogonal'); + + tagmask := tagmask and TagFullMask; + if (tagmask = 0) then exit; + + if (forEachAtPoint(ax0, ay0, nil, tagmask) = nil) then exit; + + minx := mMinX; + miny := mMinY; + + // offset query coords to (0,0)-based + x0 := ax0-minx; + y0 := ay0-miny; + x1 := ax1-minx; + y1 := ay1-miny; + + if (x0 = x1) then + begin + if (x0 < 0) or (x0 >= mWidth*mTileSize) then exit; // oops + // vertical + if (y0 < y1) then + begin + // down + if (y1 < 0) or (y0 >= mHeight*mTileSize) then exit; + //if (ay0 < 0) then ay0 := 0; + if (y0 < 0) then exit; + if (y1 >= mHeight*mTileSize) then y1 := mHeight*mTileSize-1; + dy := 1; + end + else + begin + // up + if (y0 < 0) or (y1 >= mHeight*mTileSize) then exit; + //if (ay1 < 0) then ay1 := 0; + if (y1 < 0) then exit; + if (y0 >= mHeight*mTileSize) then y0 := mHeight*mTileSize-1; + dy := -1; + end; + // check tile + while true do + begin + ccidx := mGrid[(y0 div mTileSize)*mWidth+(x0 div mTileSize)]; + FillChar(filled, sizeof(filled), 0); + celly0 := y0 and (not (mTileSize-1)); + celly1 := celly0+mTileSize-1; + while (ccidx <> -1) do + begin + cc := @mCells[ccidx]; + for f := 0 to GridCellBucketSize-1 do + begin + if (cc.bodies[f] = -1) then break; + px := @mProxies[cc.bodies[f]]; + ptag := px.mTag; + if ((ptag and TagDisabled) = 0) and ((ptag and tagmask) <> 0) and + (ax0 >= px.x0) and (ax0 <= px.x1) then + begin + // bound c0 and c1 to cell + c0 := nclamp(px.y0-miny, celly0, celly1); + c1 := nclamp(px.y1-miny, celly0, celly1); + // fill the thing + {$IF DEFINED(D2F_DEBUG_OTR)} + e_LogWritefln('**px.y0=%s; px.y1=%s; c0=%s; c1=%s; celly0=%s; celly1=%s; [%s..%s]', [px.y0-miny, px.y1-miny, c0, c1, celly0, celly1, c0-celly0, (c0-celly0)+(c1-c0)]); + {$ENDIF} + //assert(c0 <= c1); + FillChar(filled[c0-celly0], c1-c0+1, 1); + end; + end; + // next cell + ccidx := cc.next; + end; + {$IF DEFINED(D2F_DEBUG_OTR)} + s := formatstrf(' x=%s; ay0=%s; ay1=%s; y0=%s; celly0=%s; celly1=%s; dy=%s; [', [ax0, ay0, ay1, y0, celly0, celly1, dy]); + for f := 0 to High(filled) do if (filled[f] <> 0) then s += '1' else s += '0'; + s += ']'; + e_LogWriteln(s); + {$ENDIF} + // now go till we hit cell boundary or empty space + if (dy < 0) then + begin + // up + while (y0 >= celly0) and (filled[y0-celly0] <> 0) do + begin + {$IF DEFINED(D2F_DEBUG_OTR)} + e_LogWritefln(' filled: cdy=%s; y0=%s; celly0=%s; ay0=%s; ay1=%s', [y0-celly0, y0, celly0, ay0, ay1]); + {$ENDIF} + Dec(y0); + Dec(ay0); + end; + {$IF DEFINED(D2F_DEBUG_OTR)} + e_LogWritefln(' span done: cdy=%s; y0=%s; celly0=%s; ay0=%s; ay1=%s', [y0-celly0, y0, celly0, ay0, ay1]); + {$ENDIF} + if (ay0 <= ay1) then begin ey := ay1; result := false; exit; end; + if (y0 >= celly0) then begin ey := ay0+1; {assert(forEachAtPoint(ex, ey, nil, tagmask) <> nil);} result := true; exit; end; + end + else + begin + // down + while (y0 <= celly1) and (filled[y1-celly0] <> 0) do begin Inc(y0); Inc(ay0); end; + if (ay0 >= ay1) then begin ey := ay1; result := false; exit; end; + if (y0 <= celly1) then begin ey := ay0-1; result := true; exit; end; + end; + end; + end + else + begin + // horizontal + assert(false); + end; end;