1

/* $Id$ */

     2

     3

/** @file map.cpp */

     4

     5

#include "stdafx.h"

     6

#include "openttd.h"

     7

#include "debug.h"

     8

#include "functions.h"

     9

#include "macros.h"

    10

#include "map.h"

    11

#include "direction.h"

    12

#include "helpers.hpp"

    13

    14

#if defined(_MSC_VER) && _MSC_VER >= 1400 /* VStudio 2005 is stupid! */

    15

/* Why the hell is that not in all MSVC headers?? */

    16

extern "C" _CRTIMP void __cdecl _assert(void *, void *, unsigned);

    17

#endif

    18

    19

uint _map_log_x;     ///< 2^_map_log_x == _map_size_x

    20

uint _map_size_x;    ///< Size of the map along the X

    21

uint _map_size_y;    ///< Size of the map along the Y

    22

uint _map_size;      ///< The number of tiles on the map

    23

uint _map_tile_mask; ///< _map_size - 1 (to mask the mapsize)

    24

    25

Tile *_m = NULL;          ///< Tiles of the map

    26

TileExtended *_me = NULL; ///< Extended Tiles of the map

    27

    28

    29

/**

    30

 * (Re)allocates a map with the given dimension

    31

 * @param size_x the width of the map along the NE/SW edge

    32

 * @param size_y the 'height' of the map along the SE/NW edge

    33

 */

    34

void AllocateMap(uint size_x, uint size_y)

    35

{

    36

	/* Make sure that the map size is within the limits and that

    37

	 * the x axis size is a power of 2. */

    38

	if (size_x < 64 || size_x > 2048 ||

    39

			size_y < 64 || size_y > 2048 ||

    40

			(size_x & (size_x - 1)) != 0 ||

    41

			(size_y & (size_y - 1)) != 0)

    42

		error("Invalid map size");

    43

    44

	DEBUG(map, 1, "Allocating map of size %dx%d", size_x, size_y);

    45

    46

	_map_log_x = FindFirstBit(size_x);

    47

	_map_size_x = size_x;

    48

	_map_size_y = size_y;

    49

	_map_size = size_x * size_y;

    50

	_map_tile_mask = _map_size - 1;

    51

    52

	free(_m);

    53

	free(_me);

    54

    55

	_m = CallocT<Tile>(_map_size);

    56

	_me = CallocT<TileExtended>(_map_size);

    57

    58

	/* XXX @todo handle memory shortage more gracefully

    59

	 * Maybe some attemps could be made to try with smaller maps down to 64x64

    60

	 * Maybe check for available memory before doing the calls, after all, we know how big

    61

	 * the map is */

    62

	if ((_m == NULL) || (_me == NULL)) error("Failed to allocate memory for the map");

    63

}

    64

    65

    66

#ifdef _DEBUG

    67

TileIndex TileAdd(TileIndex tile, TileIndexDiff add,

    68

	const char *exp, const char *file, int line)

    69

{

    70

	int dx;

    71

	int dy;

    72

	uint x;

    73

	uint y;

    74

    75

	dx = add & MapMaxX();

    76

	if (dx >= (int)MapSizeX() / 2) dx -= MapSizeX();

    77

	dy = (add - dx) / (int)MapSizeX();

    78

    79

	x = TileX(tile) + dx;

    80

	y = TileY(tile) + dy;

    81

    82

	if (x >= MapSizeX() || y >= MapSizeY()) {

    83

		char buf[512];

    84

    85

		snprintf(buf, lengthof(buf), "TILE_ADD(%s) when adding 0x%.4X and 0x%.4X failed",

    86

			exp, tile, add);

    87

#if !defined(_MSC_VER)

    88

		fprintf(stderr, "%s:%d %s\n", file, line, buf);

    89

#else

    90

		_assert(buf, (char*)file, line);

    91

#endif

    92

	}

    93

    94

	assert(TileXY(x, y) == TILE_MASK(tile + add));

    95

    96

	return TileXY(x, y);

    97

}

    98

#endif

    99

   100

/**

   101

 * Scales the given value by the map size, where the given value is

   102

 * for a 256 by 256 map

   103

 * @param n the value to scale

   104

 * @return the scaled size

   105

 */

   106

uint ScaleByMapSize(uint n)

   107

{

   108

	/* First shift by 12 to prevent integer overflow for large values of n.

   109

	 * >>12 is safe since the min mapsize is 64x64

   110

	 * Add (1<<4)-1 to round upwards. */

   111

	return (n * (MapSize() >> 12) + (1 << 4) - 1) >> 4;

   112

}

   113

   114

   115

/**

   116

 * Scales the given value by the maps circumference, where the given

   117

 * value is for a 256 by 256 map

   118

 * @param n the value to scale

   119

 * @return the scaled size

   120

 */

   121

uint ScaleByMapSize1D(uint n)

   122

{

   123

	/* Normal circumference for the X+Y is 256+256 = 1<<9

   124

	 * Note, not actually taking the full circumference into account,

   125

	 * just half of it.

   126

	 * (1<<9) - 1 is there to scale upwards. */

   127

	return (n * (MapSizeX() + MapSizeY()) + (1 << 9) - 1) >> 9;

   128

}

   129

   130

   131

/**

   132

 * This function checks if we add addx/addy to tile, if we

   133

 *  do wrap around the edges. For example, tile = (10,2) and

   134

 *  addx = +3 and addy = -4. This function will now return

   135

 *  INVALID_TILE, because the y is wrapped. This is needed in

   136

 *  for example, farmland. When the tile is not wrapped,

   137

 *  the result will be tile + TileDiffXY(addx, addy)

   138

 * @param tile the 'starting' point of the adding

   139

 * @param addx the amount of tiles in the X direction to add

   140

 * @param addy the amount of tiles in the Y direction to add

   141

 * @return translated tile, or INVALID_TILE when it would've wrapped.

   142

 */

   143

uint TileAddWrap(TileIndex tile, int addx, int addy)

   144

{

   145

	uint x = TileX(tile) + addx;

   146

	uint y = TileY(tile) + addy;

   147

   148

	/* Are we about to wrap? */

   149

	if (x < MapMaxX() && y < MapMaxY())

   150

		return tile + TileDiffXY(addx, addy);

   151

   152

	return INVALID_TILE;

   153

}

   154

   155

/** 'Lookup table' for tile offsets given a DiagDirection */

   156

extern const TileIndexDiffC _tileoffs_by_diagdir[] = {

   157

	{-1,  0}, ///< DIAGDIR_NE

   158

	{ 0,  1}, ///< DIAGDIR_SE

   159

	{ 1,  0}, ///< DIAGDIR_SW

   160

	{ 0, -1}  ///< DIAGDIR_NW

   161

};

   162

   163

/** 'Lookup table' for tile offsets given a Direction */

   164

extern const TileIndexDiffC _tileoffs_by_dir[] = {

   165

	{-1, -1}, ///< DIR_N

   166

	{-1,  0}, ///< DIR_NE

   167

	{-1,  1}, ///< DIR_E

   168

	{ 0,  1}, ///< DIR_SE

   169

	{ 1,  1}, ///< DIR_S

   170

	{ 1,  0}, ///< DIR_SW

   171

	{ 1, -1}, ///< DIR_W

   172

	{ 0, -1}  ///< DIR_NW

   173

};

   174

   175

/**

   176

 * Gets the Manhattan distance between the two given tiles.

   177

 * The Manhattan distance is the sum of the delta of both the

   178

 * X and Y component.

   179

 * Also known as L1-Norm

   180

 * @param t0 the start tile

   181

 * @param t1 the end tile

   182

 * @return the distance

   183

 */

   184

uint DistanceManhattan(TileIndex t0, TileIndex t1)

   185

{

   186

	const uint dx = delta(TileX(t0), TileX(t1));

   187

	const uint dy = delta(TileY(t0), TileY(t1));

   188

	return dx + dy;

   189

}

   190

   191

   192

/**

   193

 * Gets the 'Square' distance between the two given tiles.

   194

 * The 'Square' distance is the square of the shortest (straight line)

   195

 * distance between the two tiles.

   196

 * Also known as euclidian- or L2-Norm squared.

   197

 * @param t0 the start tile

   198

 * @param t1 the end tile

   199

 * @return the distance

   200

 */

   201

uint DistanceSquare(TileIndex t0, TileIndex t1)

   202

{

   203

	const int dx = TileX(t0) - TileX(t1);

   204

	const int dy = TileY(t0) - TileY(t1);

   205

	return dx * dx + dy * dy;

   206

}

   207

   208

   209

/**

   210

 * Gets the biggest distance component (x or y) between the two given tiles.

   211

 * Also known as L-Infinity-Norm.

   212

 * @param t0 the start tile

   213

 * @param t1 the end tile

   214

 * @return the distance

   215

 */

   216

uint DistanceMax(TileIndex t0, TileIndex t1)

   217

{

   218

	const uint dx = delta(TileX(t0), TileX(t1));

   219

	const uint dy = delta(TileY(t0), TileY(t1));

   220

	return dx > dy ? dx : dy;

   221

}

   222

   223

   224

/**

   225

 * Gets the biggest distance component (x or y) between the two given tiles

   226

 * plus the Manhattan distance, i.e. two times the biggest distance component

   227

 * and once the smallest component.

   228

 * @param t0 the start tile

   229

 * @param t1 the end tile

   230

 * @return the distance

   231

 */

   232

uint DistanceMaxPlusManhattan(TileIndex t0, TileIndex t1)

   233

{

   234

	const uint dx = delta(TileX(t0), TileX(t1));

   235

	const uint dy = delta(TileY(t0), TileY(t1));

   236

	return dx > dy ? 2 * dx + dy : 2 * dy + dx;

   237

}

   238

   239

/**

   240

 * Param the minimum distance to an edge

   241

 * @param tile the tile to get the distance from

   242

 * @return the distance from the edge in tiles

   243

 */

   244

uint DistanceFromEdge(TileIndex tile)

   245

{

   246

	const uint xl = TileX(tile);

   247

	const uint yl = TileY(tile);

   248

	const uint xh = MapSizeX() - 1 - xl;

   249

	const uint yh = MapSizeY() - 1 - yl;

   250

	const uint minl = xl < yl ? xl : yl;

   251

	const uint minh = xh < yh ? xh : yh;

   252

	return minl < minh ? minl : minh;

   253

}

   254

   255

/**

   256

 * Function performing a search around a center tile and going outward, thus in circle.

   257

 * Although it really is a square search...

   258

 * Every tile will be tested by means of the callback function proc,

   259

 * which will determine if yes or no the given tile meets criteria of search.

   260

 * @param tile to start the search from

   261

 * @param size: number of tiles per side of the desired search area

   262

 * @param proc: callback testing function pointer.

   263

 * @param data to be passed to the callback function. Depends on the implementation

   264

 * @return result of the search

   265

 * @pre proc != NULL

   266

 * @pre size > 0

   267

 */

   268

bool CircularTileSearch(TileIndex tile, uint size, TestTileOnSearchProc proc, uint32 data)

   269

{

   270

	uint n, x, y;

   271

	DiagDirection dir;

   272

   273

	assert(proc != NULL);

   274

	assert(size > 0);

   275

   276

	x = TileX(tile);

   277

	y = TileY(tile);

   278

   279

	if (size % 2 == 1) {

   280

		/* If the length of the side is uneven, the center has to be checked

   281

		 * separately, as the pattern of uneven sides requires to go around the center */

   282

		n = 2;

   283

		if (proc(TileXY(x, y), data)) return true;

   284

   285

		/* If tile test is not successfull, get one tile down and left,

   286

		 * ready for a test in first circle around center tile */

   287

		x += _tileoffs_by_dir[DIR_W].x;

   288

		y += _tileoffs_by_dir[DIR_W].y;

   289

	} else {

   290

		n = 1;

   291

		/* To use _tileoffs_by_diagdir's order, we must relocate to

   292

		 * another tile, as we now first go 'up', 'right', 'down', 'left'

   293

		 * instead of 'right', 'down', 'left', 'up', which the calling

   294

		 * function assume. */

   295

		x++;

   296

	}

   297

   298

	for (; n < size; n += 2) {

   299

		for (dir = DIAGDIR_NE; dir < DIAGDIR_END; dir++) {

   300

			uint j;

   301

			for (j = n; j != 0; j--) {

   302

				if (x <= MapMaxX() && y <= MapMaxY() && ///< Is the tile within the map?

   303

						proc(TileXY(x, y), data)) {     ///< Is the callback successfulll?

   304

					return true;                        ///< then stop the search

   305

				}

   306

   307

				/* Step to the next 'neighbour' in the circular line */

   308

				x += _tileoffs_by_diagdir[dir].x;

   309

				y += _tileoffs_by_diagdir[dir].y;

   310

			}

   311

		}

   312

		/* Jump to next circle to test */

   313

		x += _tileoffs_by_dir[DIR_W].x;

   314

		y += _tileoffs_by_dir[DIR_W].y;

   315

	}

   316

	return false;

   317

}