| 1 | #ifndef __GEOWORLD_H__
|
|---|
| 2 | #define __GEOWORLD_H__
|
|---|
| 3 |
|
|---|
| 4 | #include <mars_grids.h>
|
|---|
| 5 | #include <mars_lookup.h>
|
|---|
| 6 | #include <mars_linalg.h>
|
|---|
| 7 | #include <mars_log.h>
|
|---|
| 8 | #include <mars_events.h>
|
|---|
| 9 | #include <mars_ptr.h>
|
|---|
| 10 | #include <mars_util.h>
|
|---|
| 11 | #include <mars_calc.h>
|
|---|
| 12 |
|
|---|
| 13 | #include <fstream>
|
|---|
| 14 | #include <boost/thread.hpp>
|
|---|
| 15 | #include <list>
|
|---|
| 16 | #include <string>
|
|---|
| 17 | #include <utility>
|
|---|
| 18 | #include <set>
|
|---|
| 19 |
|
|---|
| 20 | #include "sphpack.h"
|
|---|
| 21 | #include "worldfile.h"
|
|---|
| 22 | #include "gwutil.h"
|
|---|
| 23 | #include "PNGDEM.h"
|
|---|
| 24 | #include "geomorph.h"
|
|---|
| 25 | #include "oregen.h"
|
|---|
| 26 |
|
|---|
| 27 | namespace geoworld
|
|---|
| 28 | {
|
|---|
| 29 | class UnlockEx : public std::exception
|
|---|
| 30 | {
|
|---|
| 31 | public:
|
|---|
| 32 | const mars::BBox< long > bbox;
|
|---|
| 33 |
|
|---|
| 34 | UnlockEx (const mars::BBox< long > & bbox)
|
|---|
| 35 | : bbox(bbox) {}
|
|---|
| 36 | };
|
|---|
| 37 |
|
|---|
| 38 | class PGWFileEx : public std::exception
|
|---|
| 39 | {
|
|---|
| 40 | public:
|
|---|
| 41 | const std::string message;
|
|---|
| 42 |
|
|---|
| 43 | PGWFileEx(const char * szMsg) : message(szMsg) {}
|
|---|
| 44 | };
|
|---|
| 45 |
|
|---|
| 46 | class PGWIllegal : public std::exception
|
|---|
| 47 | {
|
|---|
| 48 | public:
|
|---|
| 49 | const std::string message;
|
|---|
| 50 |
|
|---|
| 51 | PGWIllegal(const char * szMsg) : message(szMsg) {}
|
|---|
| 52 | };
|
|---|
| 53 |
|
|---|
| 54 | class LockedGWSection
|
|---|
| 55 | {
|
|---|
| 56 | public:
|
|---|
| 57 | static const unsigned short MAX_LOCK_DIM = 5000;
|
|---|
| 58 | enum StratumSyncOp
|
|---|
| 59 | {
|
|---|
| 60 | SSync_None = 0,
|
|---|
| 61 | SSync_Perfect = 1
|
|---|
| 62 | };
|
|---|
| 63 |
|
|---|
| 64 | class View
|
|---|
| 65 | {
|
|---|
| 66 | private:
|
|---|
| 67 | struct
|
|---|
| 68 | {
|
|---|
| 69 | const GeoHeightMap::View * heightmap;
|
|---|
| 70 | const Stratum::View * stratum;
|
|---|
| 71 | } const_view;
|
|---|
| 72 |
|
|---|
| 73 | struct
|
|---|
| 74 | {
|
|---|
| 75 | GeoHeightMap::View * heightmap;
|
|---|
| 76 | Stratum::View * stratum;
|
|---|
| 77 | } mutable_view;
|
|---|
| 78 |
|
|---|
| 79 | // What? do I look like a Staples store to you?
|
|---|
| 80 | View (const View & copy);
|
|---|
| 81 |
|
|---|
| 82 | public:
|
|---|
| 83 | const unsigned short
|
|---|
| 84 | width, height;
|
|---|
| 85 |
|
|---|
| 86 | View (GeoHeightMap::View * heightmap, Stratum::View * stratum)
|
|---|
| 87 | : width(heightmap->width), height(heightmap->height)
|
|---|
| 88 | {
|
|---|
| 89 | mutable_view.heightmap = heightmap;
|
|---|
| 90 | mutable_view.stratum = stratum;
|
|---|
| 91 | const_view.heightmap = heightmap;
|
|---|
| 92 | const_view.stratum = stratum;
|
|---|
| 93 | }
|
|---|
| 94 | View (const GeoHeightMap::View * heightmap, const Stratum::View * stratum)
|
|---|
| 95 | : width(heightmap->width), height(heightmap->height)
|
|---|
| 96 | {
|
|---|
| 97 | mutable_view.heightmap = NULL;
|
|---|
| 98 | mutable_view.stratum = NULL;
|
|---|
| 99 | const_view.heightmap = heightmap;
|
|---|
| 100 | const_view.stratum = stratum;
|
|---|
| 101 | }
|
|---|
| 102 |
|
|---|
| 103 | inline GeoHeightMap::View * getHeightMap ()
|
|---|
| 104 | {
|
|---|
| 105 | assert(mutable_view.heightmap != NULL);
|
|---|
| 106 | return mutable_view.heightmap;
|
|---|
| 107 | }
|
|---|
| 108 | inline const GeoHeightMap::View * getHeightMap () const { return const_view.heightmap; }
|
|---|
| 109 | inline Stratum::View * getStratum ()
|
|---|
| 110 | {
|
|---|
| 111 | assert(mutable_view.stratum != NULL);
|
|---|
| 112 | return mutable_view.stratum;
|
|---|
| 113 | }
|
|---|
| 114 | inline const Stratum::View * getStratum () const { return const_view.stratum; }
|
|---|
| 115 | };
|
|---|
| 116 |
|
|---|
| 117 | const mars::BBox< long >
|
|---|
| 118 | bbox;
|
|---|
| 119 |
|
|---|
| 120 | LockedGWSection( const mars::BBox< long > & bbox, const long nDepth, const unsigned short nDEMLOD, const size_t nLayers, const unsigned short nStratumLOD );
|
|---|
| 121 | LockedGWSection( const GeoHeightMap & hm, const unsigned short nDEMLOD, const Stratum & stratum );
|
|---|
| 122 |
|
|---|
| 123 | static bool isValid (const mars::BBox< long > & bbox, const unsigned short nTileDim, const bool bPad = true );
|
|---|
| 124 |
|
|---|
| 125 | static mars::BBox< long > computeLockedSection (const mars::BBox< long > & bbox, const unsigned short nTileDim, const bool bPad = true );
|
|---|
| 126 | inline static unsigned short getMaxViewportDim (const unsigned short nTileDim, const bool bPad = true )
|
|---|
| 127 | { return getActualViewportDim(MAX_LOCK_DIM, nTileDim, bPad); }
|
|---|
| 128 | inline static unsigned short getActualViewportDim (const unsigned short nDimTest, const unsigned short nTileDim, const bool bPad = true)
|
|---|
| 129 | { return nDimTest - nTileDim * (bPad ? 2 : 1); } // DEPS: LGWSBBA
|
|---|
| 130 |
|
|---|
| 131 | inline unsigned short getActualWidth () const { return _hm.width; }
|
|---|
| 132 | inline unsigned short getActualHeight () const { return _hm.height; }
|
|---|
| 133 |
|
|---|
| 134 | inline mars::Range< GeoHeightMap::Precision > getElevationRange () const { return _hm.range(); }
|
|---|
| 135 |
|
|---|
| 136 | inline void setStratumSyncOp (const StratumSyncOp & ensso)
|
|---|
| 137 | { _ssop = ensso; }
|
|---|
| 138 | inline const StratumSyncOp & getStratumSyncOp () const
|
|---|
| 139 | { return _ssop; }
|
|---|
| 140 |
|
|---|
| 141 | inline mars::ScalarFieldView< GeoHeightMap::Precision >::Offset at (const long x, const long y)
|
|---|
| 142 | { return GeoHeightMap::Offset(&_hm, static_cast< signed int > (x - bbox.left), static_cast< signed int > (y - bbox.top)); }
|
|---|
| 143 |
|
|---|
| 144 | void operator << (const View & view);
|
|---|
| 145 | void operator >> (View && view) const;
|
|---|
| 146 |
|
|---|
| 147 | inline void operator << (const Stratum::Offset offs)
|
|---|
| 148 | { _stratum << offs; }
|
|---|
| 149 |
|
|---|
| 150 | void blitFrom (const FileSource::MainDEM::Block & block, const unsigned short nTileX, const unsigned short nTileY);
|
|---|
| 151 | void blitTo (FileSource::MainDEM::Block & block, const unsigned short nTileX, const unsigned short nTileY) const;
|
|---|
| 152 | inline void processSyncOp ()
|
|---|
| 153 | {
|
|---|
| 154 | switch (_ssop)
|
|---|
| 155 | {
|
|---|
| 156 | case SSync_Perfect:
|
|---|
| 157 | _stratum.applyDEM(_hm, _nDEMLOD);
|
|---|
| 158 | break;
|
|---|
| 159 | case SSync_None:
|
|---|
| 160 | break;
|
|---|
| 161 | }
|
|---|
| 162 | }
|
|---|
| 163 |
|
|---|
| 164 | void operator >> (PNGDEM & pngdem) const;
|
|---|
| 165 | void operator << (const PNGDEM & pngdem);
|
|---|
| 166 |
|
|---|
| 167 | View * createView (const long nLeft, const long nTop, const long nRight, const long nBottom);
|
|---|
| 168 | const View * createView (const long nLeft, const long nTop, const long nRight, const long nBottom) const;
|
|---|
| 169 |
|
|---|
| 170 | inline void releaseView (const View * pView) const
|
|---|
| 171 | {
|
|---|
| 172 | LOG(mars::Log::Debug) << pView->getHeightMap();
|
|---|
| 173 | _hm.releaseView(pView->getHeightMap());
|
|---|
| 174 | _stratum.releaseView(pView->getStratum());
|
|---|
| 175 | delete pView;
|
|---|
| 176 | }
|
|---|
| 177 |
|
|---|
| 178 | template< typename J, typename JST >
|
|---|
| 179 | inline LockedGWSection & operator += (const mars::ScalarFieldView< J, JST > & other)
|
|---|
| 180 | {
|
|---|
| 181 | _hm += other;
|
|---|
| 182 | return *this;
|
|---|
| 183 | }
|
|---|
| 184 | inline LockedGWSection & operator += (const GeoHeightMap::View & other)
|
|---|
| 185 | {
|
|---|
| 186 | _hm += other;
|
|---|
| 187 | return *this;
|
|---|
| 188 | }
|
|---|
| 189 | template< typename J, typename JST >
|
|---|
| 190 | inline LockedGWSection & operator -= (const mars::ScalarFieldView< J, JST > & other)
|
|---|
| 191 | {
|
|---|
| 192 | _hm -= other;
|
|---|
| 193 | return *this;
|
|---|
| 194 | }
|
|---|
| 195 | inline LockedGWSection & operator -= (const GeoHeightMap::View & other)
|
|---|
| 196 | {
|
|---|
| 197 | _hm -= other;
|
|---|
| 198 | return *this;
|
|---|
| 199 | }
|
|---|
| 200 |
|
|---|
| 201 | inline const GeoHeightMap::Precision & heightmap (const long x, const long y) const
|
|---|
| 202 | {
|
|---|
| 203 | assert (bbox.inside(x, y));
|
|---|
| 204 | return _hm.getw(x - bbox.left, y - bbox.top);
|
|---|
| 205 | }
|
|---|
| 206 | inline GeoHeightMap::Precision & heightmap (const long x, const long y)
|
|---|
| 207 | {
|
|---|
| 208 | assert (bbox.inside(x, y));
|
|---|
| 209 | return _hm.getw(x - bbox.left, y - bbox.top);
|
|---|
| 210 | }
|
|---|
| 211 |
|
|---|
| 212 | inline const GeoHeightMap::Precision getSurfaceMean (const long x, const long y) const
|
|---|
| 213 | {
|
|---|
| 214 | assert(bbox.inside(x - 1, y - 1) && bbox.inside(x + 1, y + 1));
|
|---|
| 215 | return _hm.mean (x - bbox.left, y - bbox.top);
|
|---|
| 216 | }
|
|---|
| 217 | inline const GeoHeightMap::Precision getSurfaceMeanVN( const long x, const long y ) const
|
|---|
| 218 | {
|
|---|
| 219 | assert(bbox.inside(x - 1, y - 1) && bbox.inside(x + 1, y + 1));
|
|---|
| 220 | return _hm.meanVN (x - bbox.left, y - bbox.top);
|
|---|
| 221 | }
|
|---|
| 222 |
|
|---|
| 223 | inline GWCoords createCoordsRVN(const long x, const long y, const unsigned short nScale) const
|
|---|
| 224 | {
|
|---|
| 225 | typedef mars::Numerics< GeoHeightMap::Precision >::UP R;
|
|---|
| 226 |
|
|---|
| 227 | return GWCoords (x, y,
|
|---|
| 228 | static_cast< GWCoords::Precision > (
|
|---|
| 229 | (
|
|---|
| 230 | static_cast< R > (_hm.meanVN(x - nScale - bbox.left, y - nScale - bbox.top)) +
|
|---|
| 231 | static_cast< R > (_hm.meanVN(x + nScale - bbox.left, y - nScale - bbox.top)) +
|
|---|
| 232 | static_cast< R > (_hm.meanVN(x - nScale - bbox.left, y + nScale - bbox.top)) +
|
|---|
| 233 | static_cast< R > (_hm.meanVN(x + nScale - bbox.left, y + nScale - bbox.top)) +
|
|---|
| 234 | static_cast< R > (_hm.meanVN(x - bbox.left, y - bbox.top))
|
|---|
| 235 | ) / 5
|
|---|
| 236 | )
|
|---|
| 237 | );
|
|---|
| 238 | }
|
|---|
| 239 |
|
|---|
| 240 | inline GWCoords createCoords(const long x, const long y) const
|
|---|
| 241 | { return GWCoords (x, y, _hm.meanVN(x - bbox.left, y - bbox.top)); }
|
|---|
| 242 |
|
|---|
| 243 | inline GWCoords createCoords(const GWSurfacePos & pos) const
|
|---|
| 244 | { return createCoords(pos.x, pos.y); }
|
|---|
| 245 |
|
|---|
| 246 | inline GWCoords createCoordsRVN(const GWSurfacePos & pos, const unsigned short nScale) const
|
|---|
| 247 | { return createCoordsRVN(pos.x, pos.y, nScale); }
|
|---|
| 248 |
|
|---|
| 249 | private:
|
|---|
| 250 | GeoHeightMap _hm;
|
|---|
| 251 | unsigned short _nDEMLOD;
|
|---|
| 252 | Stratum _stratum;
|
|---|
| 253 | StratumSyncOp _ssop;
|
|---|
| 254 |
|
|---|
| 255 | /* LGWSBBA: Section bbox computation algorithm:
|
|---|
| 256 | * - Aligned along tile-boundary
|
|---|
| 257 | * - Iff padding:
|
|---|
| 258 | * - Minimum extents minus one-half tile dimension
|
|---|
| 259 | * - Maximum extents plus one-half tile dimension
|
|---|
| 260 | * - Maximum extents plus one tile dimension
|
|---|
| 261 | */
|
|---|
| 262 | };
|
|---|
| 263 |
|
|---|
| 264 | class RegionContainer : private std::list< mars::BBox< long > >
|
|---|
| 265 | {
|
|---|
| 266 | private:
|
|---|
| 267 | const unsigned long
|
|---|
| 268 | _nWWidth, _nWHeight;
|
|---|
| 269 |
|
|---|
| 270 | public:
|
|---|
| 271 | RegionContainer(const unsigned long nWWidth, const unsigned long nWHeight)
|
|---|
| 272 | : _nWWidth(nWWidth), _nWHeight(nWHeight) {}
|
|---|
| 273 |
|
|---|
| 274 | inline void addRegion (const mars::BBox< long > & bbox)
|
|---|
| 275 | {
|
|---|
| 276 | assert(!intersects(bbox));
|
|---|
| 277 | push_back(bbox);
|
|---|
| 278 | }
|
|---|
| 279 | bool removeRegion (const mars::BBox< long > & bbox);
|
|---|
| 280 | bool intersects (const mars::BBox< long > & bbox) const;
|
|---|
| 281 | inline bool empty () const { return std::list< mars::BBox< long > >::empty(); }
|
|---|
| 282 | };
|
|---|
| 283 |
|
|---|
| 284 | class QTMorphs : private mars::QuadTree_BBoxIndexed< GeoMorph, long >
|
|---|
| 285 | {
|
|---|
| 286 | private:
|
|---|
| 287 | unsigned long
|
|---|
| 288 | _nWWidth, _nWHeight;
|
|---|
| 289 |
|
|---|
| 290 | public:
|
|---|
| 291 | using MorphIterator = mars::QuadTree_BBoxIndexed< GeoMorph, long >::EntityRectIterator;
|
|---|
| 292 | using Region = mars::QuadTree_BBoxIndexed< GeoMorph, long >::Region;
|
|---|
| 293 | using VECTOR = mars::QuadTree_BBoxIndexed< GeoMorph, long >::VECTOR;
|
|---|
| 294 |
|
|---|
| 295 | QTMorphs (const unsigned long nWWidth, const unsigned long nWHeight)
|
|---|
| 296 | : QuadTree_BBoxIndexed(std::max(nWWidth, nWHeight)), _nWWidth(nWWidth), _nWHeight(nWHeight) {}
|
|---|
| 297 |
|
|---|
| 298 | bool addMorph (mars::ptr <GeoMorph> gm);
|
|---|
| 299 | inline MorphIterator morphs (const mars::BBox< long > & bbox)
|
|---|
| 300 | { return morphs(Region(bbox.getMinimum(), bbox.getMaximum())); }
|
|---|
| 301 | inline MorphIterator morphs (const long nLeft, const long nTop, const long nRight, const long nBottom)
|
|---|
| 302 | { return morphs(Region(VECTOR(nLeft, nTop), VECTOR(nRight, nBottom))); }
|
|---|
| 303 | inline MorphIterator morphs (const Region & rg)
|
|---|
| 304 | { return contents(rg); }
|
|---|
| 305 | };
|
|---|
| 306 |
|
|---|
| 307 | class AgentMap
|
|---|
| 308 | {
|
|---|
| 309 | private:
|
|---|
| 310 | class FieldTypeWrapper
|
|---|
| 311 | {
|
|---|
| 312 | public:
|
|---|
| 313 | FieldType type;
|
|---|
| 314 |
|
|---|
| 315 | FieldTypeWrapper(const FieldType enftType)
|
|---|
| 316 | : type(enftType) {}
|
|---|
| 317 |
|
|---|
| 318 | inline operator FieldType () const { return type; }
|
|---|
| 319 | };
|
|---|
| 320 |
|
|---|
| 321 | public:
|
|---|
| 322 | FieldType type;
|
|---|
| 323 | const unsigned long width, height;
|
|---|
| 324 | const float base;
|
|---|
| 325 |
|
|---|
| 326 | AgentMap (const FieldType enftType, const unsigned long nWidth, const unsigned long nHeight, const float fBase = 0.0f)
|
|---|
| 327 | : type(enftType), width(nWidth), height(nHeight), _qtField(std::max(nWidth, nHeight)), base(fBase) {}
|
|---|
| 328 |
|
|---|
| 329 | inline void addSpherePack (const mars::ptr< GeoSpherePack > & gsp)
|
|---|
| 330 | { gsp->store(_qtField, mars::ptr< FieldTypeWrapper > (new FieldTypeWrapper(type))); }
|
|---|
| 331 |
|
|---|
| 332 | inline float sample (const long x, const long y, const long z) const
|
|---|
| 333 | { return sample (GWCoords(x, y, z)); }
|
|---|
| 334 | inline mars::vector3Df gradient (const long x, const long y, const long z) const
|
|---|
| 335 | { return gradient (GWCoords(x, y, z)); }
|
|---|
| 336 |
|
|---|
| 337 | float sample (const GWCoords & coords) const;
|
|---|
| 338 | mars::vector3Df gradient (const GWCoords & coords) const;
|
|---|
| 339 |
|
|---|
| 340 | private:
|
|---|
| 341 | typedef mars::QuadTree_SphereIndexed< FieldTypeWrapper, GeoSpherePack::SCALAR > QuadTreeSpherePacks; // TODO: Name collision in worldfile.h
|
|---|
| 342 | typedef std::vector< std::pair< std::string, mars::ptr< GeoSpherePack > > > SphPackList;
|
|---|
| 343 |
|
|---|
| 344 | SphPackList _sphpacks;
|
|---|
| 345 | QuadTreeSpherePacks _qtField;
|
|---|
| 346 | };
|
|---|
| 347 |
|
|---|
| 348 | class AllAgentMaps
|
|---|
| 349 | {
|
|---|
| 350 | public:
|
|---|
| 351 | AgentMap density;
|
|---|
| 352 |
|
|---|
| 353 | AllAgentMaps (const unsigned long nWidth, const unsigned long nHeight);
|
|---|
| 354 |
|
|---|
| 355 | AgentMap & operator [] ( const FieldType enftType );
|
|---|
| 356 | const AgentMap & operator [] ( const FieldType enftType ) const;
|
|---|
| 357 | };
|
|---|
| 358 |
|
|---|
| 359 | // TODO: Move into DepositSynthesizer
|
|---|
| 360 | class MineralsContainer
|
|---|
| 361 | {
|
|---|
| 362 | public:
|
|---|
| 363 | inline void addSpherePack (const std::string & sName, const mars::ptr< GeoSpherePack > & rsp)
|
|---|
| 364 | { _sphpacks.push_back(std::pair< std::string, mars::ptr< GeoSpherePack > > (sName, rsp)); }
|
|---|
| 365 |
|
|---|
| 366 | void flushTo (FileSource::MineralMap & minmap);
|
|---|
| 367 |
|
|---|
| 368 | private:
|
|---|
| 369 | typedef std::vector< std::pair< std::string, mars::ptr< GeoSpherePack > > > SphPackList;
|
|---|
| 370 |
|
|---|
| 371 | SphPackList _sphpacks;
|
|---|
| 372 | };
|
|---|
| 373 |
|
|---|
| 374 | class IGeoWorldListener
|
|---|
| 375 | {
|
|---|
| 376 | public:
|
|---|
| 377 | virtual void runningMorph (const mars::ptr< GeoMorph > & pGM, const size_t i, const size_t nTotal) {};
|
|---|
| 378 | virtual void runningMorphs (MorphList::const_iterator i0, MorphList::const_iterator iN, const size_t i, const size_t nCount, const size_t nTotal) {};
|
|---|
| 379 | virtual void rejectedMorph (const mars::ptr< GeoMorph > & pGM) {};
|
|---|
| 380 | };
|
|---|
| 381 |
|
|---|
| 382 | class PagedGeoWorld : public mars::Observable< IGeoWorldListener >
|
|---|
| 383 | {
|
|---|
| 384 | private:
|
|---|
| 385 | unsigned short
|
|---|
| 386 | _nNumProcs;
|
|---|
| 387 |
|
|---|
| 388 | GW_MUTEX(mutex);
|
|---|
| 389 |
|
|---|
| 390 | mutable FileSource * _pFile;
|
|---|
| 391 | mutable RegionContainer * _plsbboxLockedMutableRegions;
|
|---|
| 392 |
|
|---|
| 393 | DepositSynthesizer
|
|---|
| 394 | * _pDepSynth;
|
|---|
| 395 |
|
|---|
| 396 | struct Pass
|
|---|
| 397 | {
|
|---|
| 398 | QTMorphs quadtree;
|
|---|
| 399 | MorphList list;
|
|---|
| 400 |
|
|---|
| 401 | Pass(const unsigned long nWWidth, const unsigned long nWHeight)
|
|---|
| 402 | : quadtree(nWWidth, nWHeight) {}
|
|---|
| 403 | };
|
|---|
| 404 | typedef std::vector< Pass * > PassList;
|
|---|
| 405 |
|
|---|
| 406 | PassList _passes;
|
|---|
| 407 |
|
|---|
| 408 | AllAgentMaps * _pAgentMaps;
|
|---|
| 409 |
|
|---|
| 410 | protected:
|
|---|
| 411 | typedef std::set< mars::ptr< GeoMorph > > MorphSet;
|
|---|
| 412 |
|
|---|
| 413 | enum Mode
|
|---|
| 414 | {
|
|---|
| 415 | Read,
|
|---|
| 416 | Write
|
|---|
| 417 | };
|
|---|
| 418 |
|
|---|
| 419 | class Picks : public MorphSet
|
|---|
| 420 | {
|
|---|
| 421 | private:
|
|---|
| 422 | mars::BBox< long > _bbox;
|
|---|
| 423 |
|
|---|
| 424 | public:
|
|---|
| 425 | inline const mars::BBox< long > & getBBox() const { return _bbox; }
|
|---|
| 426 | inline void setBBox (const mars::BBox< long > & bbox) { _bbox = bbox; }
|
|---|
| 427 | void computeBBox ();
|
|---|
| 428 | void consume( const IMineralQuery * pQMin, LockedGWSection * pSection, MineralsContainer & ctrMins );
|
|---|
| 429 | };
|
|---|
| 430 |
|
|---|
| 431 | struct WorkItem
|
|---|
| 432 | {
|
|---|
| 433 | mars::ptr< Picks > picks;
|
|---|
| 434 | LockedGWSection * section;
|
|---|
| 435 | MineralsContainer minerals;
|
|---|
| 436 |
|
|---|
| 437 | WorkItem ()
|
|---|
| 438 | : picks(NULL), section(NULL) {}
|
|---|
| 439 |
|
|---|
| 440 | WorkItem(const mars::ptr< Picks > & pPicks)
|
|---|
| 441 | : picks(pPicks), section(NULL) {}
|
|---|
| 442 |
|
|---|
| 443 | WorkItem(const mars::ptr< GeoMorph > & pGM, const mars::BBox< long > & bbox = mars::BBox< long > ())
|
|---|
| 444 | : picks(new Picks())
|
|---|
| 445 | {
|
|---|
| 446 | picks->insert(pGM);
|
|---|
| 447 | if (bbox.empty())
|
|---|
| 448 | picks->computeBBox();
|
|---|
| 449 | else
|
|---|
| 450 | picks->setBBox(bbox);
|
|---|
| 451 | }
|
|---|
| 452 | };
|
|---|
| 453 |
|
|---|
| 454 | struct WorkQueue
|
|---|
| 455 | {
|
|---|
| 456 | typedef std::deque< WorkItem > ItemQueue;
|
|---|
| 457 |
|
|---|
| 458 | ItemQueue inbox, outbox, ready;
|
|---|
| 459 | };
|
|---|
| 460 |
|
|---|
| 461 | class IWorkQueueListener
|
|---|
| 462 | {
|
|---|
| 463 | public:
|
|---|
| 464 | virtual void readyQueue (const WorkQueue::ItemQueue & readyq) {}
|
|---|
| 465 | virtual void collisions (const WorkQueue::ItemQueue & colq) {}
|
|---|
| 466 | virtual void unlockQueue (const WorkQueue::ItemQueue & unlq) {}
|
|---|
| 467 | };
|
|---|
| 468 |
|
|---|
| 469 | class MasterSet : private MorphSet
|
|---|
| 470 | {
|
|---|
| 471 | private:
|
|---|
| 472 | QTMorphs * _pQTMorphs;
|
|---|
| 473 |
|
|---|
| 474 | GW_MUTEX(mutex);
|
|---|
| 475 |
|
|---|
| 476 | public:
|
|---|
| 477 | MasterSet (QTMorphs * pQTMorphs, MorphList::const_iterator iMorphs0, MorphList::const_iterator iMorphsN);
|
|---|
| 478 |
|
|---|
| 479 | bool empty ();
|
|---|
| 480 | mars::ptr< Picks > collect( const QTMorphs::Region & rgBBox );
|
|---|
| 481 | void push (const mars::ptr< Picks > & pPicks);
|
|---|
| 482 |
|
|---|
| 483 | MorphSet::size_type size() const { return MorphSet::size(); }
|
|---|
| 484 |
|
|---|
| 485 | const_iterator begin () const { return MorphSet::begin(); }
|
|---|
| 486 | const_iterator end () const { return MorphSet::end(); }
|
|---|
| 487 | };
|
|---|
| 488 |
|
|---|
| 489 | LockedGWSection * doLock (const mars::BBox< long > & bbox, const bool bPad = true, const bool bReadOnly = false) const;
|
|---|
| 490 |
|
|---|
| 491 | template< Mode MODE, typename TYPE_LockedGWSection >
|
|---|
| 492 | void performIOOperation (TYPE_LockedGWSection * pSection) const;
|
|---|
| 493 |
|
|---|
| 494 | template< Mode MODE >
|
|---|
| 495 | inline void blockOp ( LockedGWSection * pSection, FileSource::MainDEM::Block & block, const unsigned short i, const unsigned short j ) const
|
|---|
| 496 | {
|
|---|
| 497 | switch (MODE)
|
|---|
| 498 | {
|
|---|
| 499 | case Read: readBlock(pSection, block, i, j); break;
|
|---|
| 500 | case Write: writeBlock(pSection, block, i, j); break;
|
|---|
| 501 | }
|
|---|
| 502 | }
|
|---|
| 503 | template< Mode MODE >
|
|---|
| 504 | inline void seekTileOp (const int nTX, const int nTY, const std::ios_base::seekdir dir = std::ios_base::cur) const
|
|---|
| 505 | {
|
|---|
| 506 | switch (MODE)
|
|---|
| 507 | {
|
|---|
| 508 | case Read: _pFile->tilemap.seekg(nTX, nTY, dir); break;
|
|---|
| 509 | case Write: _pFile->tilemap.seekp(nTX, nTY, dir); break;
|
|---|
| 510 | }
|
|---|
| 511 | }
|
|---|
| 512 |
|
|---|
| 513 | void readBlock( LockedGWSection * pSection, FileSource::MainDEM::Block & block, const unsigned short i, const unsigned short j ) const;
|
|---|
| 514 | void writeBlock( const LockedGWSection * pSection, FileSource::MainDEM::Block & block, const unsigned short i, const unsigned short j ) const;
|
|---|
| 515 |
|
|---|
| 516 | void processWorkQueue (PagedGeoWorld::WorkQueue & q, IWorkQueueListener * pListener = NULL);
|
|---|
| 517 |
|
|---|
| 518 | void flushMinerals (MineralsContainer & ctrMins);
|
|---|
| 519 | void flushOutbox( WorkQueue::ItemQueue &outbox );
|
|---|
| 520 | bool isPickValid (const Picks & p) const;
|
|---|
| 521 | bool isValidSectionBBox (const mars::BBox< long > & bbox) const;
|
|---|
| 522 |
|
|---|
| 523 | bool checkGeoMorph( const mars::ptr<GeoMorph> &gm );
|
|---|
| 524 |
|
|---|
| 525 | public:
|
|---|
| 526 | static const unsigned short DEFAULT_NUM_PROCS = 16;
|
|---|
| 527 |
|
|---|
| 528 | PagedGeoWorld (DepositSynthesizer * pDepSynth, FileSource * pFile, const unsigned short nPasses);
|
|---|
| 529 | ~PagedGeoWorld();
|
|---|
| 530 |
|
|---|
| 531 | inline void setNumProcs(const unsigned short nNumProcs) { _nNumProcs = nNumProcs; }
|
|---|
| 532 |
|
|---|
| 533 | inline const unsigned short getLOD () const { return 0; } // TODO: Get rid of LOD for PagedGeoWorld
|
|---|
| 534 | unsigned short getMaxGeoMorphWidth () const;
|
|---|
| 535 | unsigned short getMaxGeoMorphHeight () const;
|
|---|
| 536 |
|
|---|
| 537 | void addMorphs (MorphList & morphs, MorphList & rejects, const unsigned short nPass);
|
|---|
| 538 | bool addMorph (mars::ptr <GeoMorph> & gm, const unsigned short nPass);
|
|---|
| 539 |
|
|---|
| 540 | inline MorphList::const_iterator beginPass (const unsigned short nPass) const { return _passes[nPass]->list.begin(); }
|
|---|
| 541 | inline MorphList::iterator beginPass (const unsigned short nPass) { return _passes[nPass]->list.begin(); }
|
|---|
| 542 | inline MorphList::const_iterator endPass (const unsigned short nPass) const { return _passes[nPass]->list.end(); }
|
|---|
| 543 | inline MorphList::iterator endPass (const unsigned short nPass) { return _passes[nPass]->list.end(); }
|
|---|
| 544 | inline unsigned short getPassCount () const { return _passes.size(); }
|
|---|
| 545 | void runMorphs (const unsigned short nPass, MorphList & rejects);
|
|---|
| 546 |
|
|---|
| 547 | inline const AllAgentMaps & getAgentMaps () const
|
|---|
| 548 | {
|
|---|
| 549 | assert(_pAgentMaps != NULL);
|
|---|
| 550 | return *_pAgentMaps;
|
|---|
| 551 | }
|
|---|
| 552 | inline AllAgentMaps & getAgentMaps ()
|
|---|
| 553 | {
|
|---|
| 554 | assert(_pAgentMaps != NULL);
|
|---|
| 555 | return *_pAgentMaps;
|
|---|
| 556 | }
|
|---|
| 557 |
|
|---|
| 558 | inline LockedGWSection * lock (const mars::BBox< long > & bbox, const bool bPad = true)
|
|---|
| 559 | { return doLock(bbox, bPad, false); }
|
|---|
| 560 | inline const LockedGWSection * lockro (const mars::BBox< long > & bbox, const bool bPad = true) const
|
|---|
| 561 | { return doLock(bbox, bPad, true); }
|
|---|
| 562 |
|
|---|
| 563 | inline LockedGWSection * lock (const long x, const long y, const unsigned short nWidth, const unsigned short nHeight, const bool bPad = true)
|
|---|
| 564 | { return doLock (mars::BBox< long > (x, y, x + nWidth - 1, y + nHeight - 1), bPad, false); }
|
|---|
| 565 | inline const LockedGWSection * lockro (const long x, const long y, const unsigned short nWidth, const unsigned short nHeight, const bool bPad = true) const
|
|---|
| 566 | { return doLock (mars::BBox< long > (x, y, x + nWidth - 1, y + nHeight - 1), bPad, true); }
|
|---|
| 567 |
|
|---|
| 568 | void unlock (LockedGWSection * pSection);
|
|---|
| 569 | void unlock (const LockedGWSection * pSection) const;
|
|---|
| 570 |
|
|---|
| 571 | inline void operator >> (mars::ScalarFieldView< GeoHeightMap::Precision > & view) const
|
|---|
| 572 | {
|
|---|
| 573 | assert(_pFile != NULL);
|
|---|
| 574 | _pFile->tilemap >> view;
|
|---|
| 575 | }
|
|---|
| 576 |
|
|---|
| 577 | inline const unsigned long getWorldWidth () const { return static_cast< unsigned long > (_pFile->tilemap.getTilesWide()) * static_cast< unsigned long > (_pFile->tilemap.getVirtualTileDim()); }
|
|---|
| 578 | inline const unsigned long getWorldHeight () const { return static_cast< unsigned long > (_pFile->tilemap.getTilesHigh()) * static_cast< unsigned long > (_pFile->tilemap.getVirtualTileDim()); }
|
|---|
| 579 | inline const unsigned short getTileDimension () const { return _pFile->tilemap.getVirtualTileDim(); }
|
|---|
| 580 | };
|
|---|
| 581 |
|
|---|
| 582 | class GeoWorld : public mars::Observable< IGeoWorldListener >
|
|---|
| 583 | {
|
|---|
| 584 | private:
|
|---|
| 585 | unsigned short
|
|---|
| 586 | _nLODLowest, _nLODMid,
|
|---|
| 587 | _nTileDim;
|
|---|
| 588 |
|
|---|
| 589 | GeoHeightMap
|
|---|
| 590 | * _phmLowest;
|
|---|
| 591 | DepositSynthesizer
|
|---|
| 592 | * _pDepSynth;
|
|---|
| 593 |
|
|---|
| 594 | MineralsContainer
|
|---|
| 595 | _ctrMins;
|
|---|
| 596 | AllAgentMaps
|
|---|
| 597 | * _pAllAgentMaps;
|
|---|
| 598 |
|
|---|
| 599 | unsigned long
|
|---|
| 600 | _nLowWidth, _nLowHeight;
|
|---|
| 601 |
|
|---|
| 602 | FileSource * _pFile;
|
|---|
| 603 |
|
|---|
| 604 | MorphList _morphs;
|
|---|
| 605 | mutable RegionContainer * _plsbboxLockedMutableRegions;
|
|---|
| 606 |
|
|---|
| 607 | GW_MUTEX(mutex);
|
|---|
| 608 |
|
|---|
| 609 | protected:
|
|---|
| 610 | LockedGWSection * lock (const mars::BBox< long > & bbox, const bool bReadOnly) const;
|
|---|
| 611 | LockedGWSection * lock (const bool bReadOnly) const;
|
|---|
| 612 |
|
|---|
| 613 | void runMorph( const MineralsDB & mindb, mars::ptr< GeoMorph > pGM, LockedGWSection * pLockedSection );
|
|---|
| 614 | bool checkGeoMorph( const mars::ptr<GeoMorph> &gm );
|
|---|
| 615 |
|
|---|
| 616 | public:
|
|---|
| 617 | static const unsigned short MIN_LOWEST_LOD = 2;
|
|---|
| 618 |
|
|---|
| 619 | GeoWorld (
|
|---|
| 620 | DepositSynthesizer * pDepSynth,
|
|---|
| 621 | const unsigned short nLODLowest,
|
|---|
| 622 | const unsigned short nLODMid,
|
|---|
| 623 | FileSource * pFile
|
|---|
| 624 | );
|
|---|
| 625 | ~GeoWorld();
|
|---|
| 626 |
|
|---|
| 627 | inline AllAgentMaps & getAgentMaps () { return *_pAllAgentMaps; }
|
|---|
| 628 | inline const AllAgentMaps & getAgentMaps () const { return *_pAllAgentMaps; }
|
|---|
| 629 |
|
|---|
| 630 | inline unsigned long getWidth () const { return _nLowWidth; }
|
|---|
| 631 | inline unsigned long getHeight () const { return _nLowHeight; }
|
|---|
| 632 | inline unsigned short getLOD () const { return _nLODLowest; }
|
|---|
| 633 | inline unsigned short getDeltaLOD () const { return _nLODLowest - _nLODMid; }
|
|---|
| 634 | unsigned short getMaxGeoMorphWidth () const;
|
|---|
| 635 | unsigned short getMaxGeoMorphHeight () const;
|
|---|
| 636 |
|
|---|
| 637 | void addMorphs( MorphList & morphs, MorphList & rejects );
|
|---|
| 638 | bool addMorph (mars::ptr <GeoMorph> & gm);
|
|---|
| 639 |
|
|---|
| 640 | inline MorphList::const_iterator beginMorphs () const { return _morphs.begin(); }
|
|---|
| 641 | inline MorphList::const_iterator endMorphs () const { return _morphs.end(); }
|
|---|
| 642 | inline MorphList::iterator beginMorphs () { return _morphs.begin(); }
|
|---|
| 643 | inline MorphList::iterator endMorphs () { return _morphs.end(); }
|
|---|
| 644 |
|
|---|
| 645 | inline LockedGWSection * lock ()
|
|---|
| 646 | { return lock(false); }
|
|---|
| 647 | inline const LockedGWSection * lock () const
|
|---|
| 648 | { return lock(true); }
|
|---|
| 649 | inline LockedGWSection * lock (const mars::BBox< long > & bbox)
|
|---|
| 650 | { return lock(bbox, false); }
|
|---|
| 651 | inline const LockedGWSection * lock (const mars::BBox< long > & bbox) const
|
|---|
| 652 | { return lock(bbox, true); }
|
|---|
| 653 |
|
|---|
| 654 | inline LockedGWSection * lock (const long x, const long y, const unsigned short nWidth, const unsigned short nHeight)
|
|---|
| 655 | { return lock (mars::BBox< long > (x, y, x + nWidth - 1, y + nHeight - 1)); }
|
|---|
| 656 | inline const LockedGWSection * lock (const long x, const long y, const unsigned short nWidth, const unsigned short nHeight) const
|
|---|
| 657 | { return lock (mars::BBox< long > (x, y, x + nWidth - 1, y + nHeight - 1)); }
|
|---|
| 658 |
|
|---|
| 659 | void unlock (LockedGWSection * pSection);
|
|---|
| 660 | void unlock (const LockedGWSection * pSection) const;
|
|---|
| 661 |
|
|---|
| 662 | void runMorphs (const MineralsDB & mindb, MorphList & rejects);
|
|---|
| 663 |
|
|---|
| 664 | void write(const float fCoarseness);
|
|---|
| 665 |
|
|---|
| 666 | inline const GeoHeightMap & getLowHM () const { return *_phmLowest; }
|
|---|
| 667 |
|
|---|
| 668 | inline void operator << (const mars::ScalarFieldView< GeoHeightMap::Precision > & view)
|
|---|
| 669 | { *_phmLowest << view; }
|
|---|
| 670 | inline void operator >> (mars::ScalarFieldView< GeoHeightMap::Precision > & view) const
|
|---|
| 671 | { *_phmLowest >> view; }
|
|---|
| 672 | };
|
|---|
| 673 | }
|
|---|
| 674 |
|
|---|
| 675 | #endif
|
|---|