Reinforcement Learning and
Artificial
Intelligence (RLAI) 

Tile
Coding Software  Reference Manual, Version 2.1 (see version 3) 
This version 2.1 of the tiling software provides a optional variation of the main routine that prevents hashing collisions, so as to exactly simulate infinite tilings (until it runs out of memory). This variation can also be used to gather statistics on hashing collisions even if they are be ignored (as in the standard version). The variation requires that an additional hashing table be passed in, and there is greater computational expense.
Next we describe the use of the main tile coding routine.
At the core of the tile coding software is a single, memoryless
routine, "tiles
":
void tiles(
int tiles[], // provided array will contain the returned tile indices
int num_tilings, // number of tile indices to be returned
int memory_size, // total number of possible tiles
float floats[], // array of floatingpoint vars to be converted to tiles
int num_floats, // number of floatingpoint variables
int ints[], // array of integer variables to be converted to tiles
int num_ints) // number of integer variables
The first argument (tiles
) will contain the output of
the routine. This array is provided by the caller and
filled up with tile indices (integers in [0,memory_size
)).
The second argument is the number of tile indices returned, so the
provided
array had better be of at least this length. It is recommended by
the folks at UNH (who designed the key algorithm here) that num_tilings
be a power of two.
The third argument, memory_size
, specifies the largest
permitted
tile
index.
The next two arguments are the primary
inputs to tiles
. They specify the floatingpoint
variables which are to be tiled, in an array, with an accompanying
integer argument to indicate the length of the array.
The routine will overlay gridtilings on these variables, generalizing
across values that differ by less than 1.0. This argument should also
be a power of two.
The last two arguments are optional. For each different
value of the num_ints
integers in ints
, the
tiles will be hashed totally differently, such that there
is no generalization to any other value of the integers. These
inputs are often used to ensure different hashings for
different groups of tilings, e.g., for different actions in
reinforcement learning.
If there are three or fewer such integer inputs, you can simply provide
them instead of the integer array, as simple integer arguments, and the
magic of C++ will figure it out.
Next we describe examples, gradually increasing in complexity, of the use of the tiles routine.
tiles
assumes unit generalization in its inputs, so we arrange for .25
generalization by scaling x and y and stuffing the
results into a float array float_array
to be input to tiles:
float float_array[2];
float_array[0] = x / 0.25;
float_array[1] = y / 0.25;
Despite this fairly
broad generalization, we would like the ability to make fine
distinctions if required.
For this we need many tilings, say 32 (a power of two is good here).
This will give us an ultimate resolution of .25/32=.0078, or better
than 1%.
To arrange for this we pass the tiles routine an integer array, tiles_array
,
of at least 32 elements, to hold the returned
results, and pass 32 as the num_tilings
argument.
int tiles_array[32];
Finally, we have to select memory_size
, the range
of tile indices that will be hashed down into. First we consider the
"natural" memory size, as if we were doing no hashing. This is simply
the total number of tiles in all the tilings,
or 4*4*32=512 in our case. This is not a large number and might be used
directly as the memory_size
.
Note, however, that this does not mean there will be no hashing or hash
collisions.
This routine always hashes and so collisions are possible. You could
make them very unlikely
by using a large memory size, say 5120, or even 512000 if you have
enough memory for
how the tiles are to be used. More often the natural memory size is too
large, say 10^10, and we must hash into a smaller memory size.
Thus, our final call to the tiles
routine might be
tiles(tiles_array,32,512,float_array,2);
The returned tiles can be used in many ways, as is discussed elsewhere. But let
us complete this simple example by showing one use. To do simple linear
function approximation
one uses an array of weights of number equal to the total number of
possible tiles
(memory_size
).
The approximate function value for the input vector of real values is
obtained
simply by summing the weights of the corresponding returned tiles:
float weights[512];
float result = 0.0;
for (i=0;i<32;i++) \
result += weights[tiles_array[i]];
target
that
we wish to adjust the approximation result towards:
float target;
float alpha = 0.1 / 32; // stepsize parameter
for (i=0;i<32;i++) \\
weight[tiles_array[i]] += alpha * (target  result);
For simple uses of tiles
, packing into arrays
is
tiresome,
so we provide routines to do the packing for you for the cases
with 1 or 2 floatingpoint state variables, called tiles1
and tiles2
.
For example, the 4x4 grid case described above could be done without
explicitly
packing into the float_array
by calling
tiles2(tiles_array,32,512,x/0.25,y/0.25);
There are a number of useful simple variations on the example
described
above.
Obviously, one could scale the two dimensions differently to get
rectangular
rather than square tiles. One could also stretch and distort the
dimensions, for example,
making one logarithmic. Or one could obtain a rotated grid by linear
transformations
on the original x and y variables. For example, to get
a grid oriented along the diagonals, we construct the float_array
like this:
float_array[0] = (x+y) / 0.25;
float_array[1] = (xy) / 0.25;
Note that the tiling produced in all these cases is not
restricted to
any region of
the space. It extends in all directions to infinity; it is a tiling of
the infinite
plane. These infinite number of tiles are reduced to memory_size
in number
by the hashing.
If you want to restrict the input vector to a limited region,
then
you must do so
yourself before calling tiles
. If you want to achieve
a wraparound
effect, then you should cause wraparound when you construct the float_array
,
e.g., by dividing by the width of the region and filling float_array
with the remainder.
Sometimes one can get by with a single kind of tiling, as
produced by a
single call
to tiles
(note that a single call creates multiple
tillings, but these are
all simple offset copies of each other). But for many other cases one
needs to be a little more sophisticated.
A great deal can be done to control the generalization and to limit the
memory requirements
by using multiple tilings through multiple calls to tiles
.
We consider some of
these below.
Suppose as previously you have input vectors in the unit square,
but
you want to generalize
between input vectors that are similar along the x dimension OR
similar along the y
dimension. To do this you need two kinds of tilings, one
producing stripes horizontally and one vertically. Suppose you want to
make 32 tilings of each type, and combine the results
in one 64element tiles_array
:
int tiles_array[64];
tiles1(tiles_array,32,4*32*2,x/0.25,0);
tiles1(&tiles_array[32],32,4*32*2,y/0.25,1);
tiles
passes in, in effect, the
second half of the full tiles_array
. The tiles from the
second group of tilings will go there. Notice that the memory_size
input
is the same for both calls and set roughly equal to the memory
requirements for both
sets of tiles: 4 (stripes over the unit interval) * 32 (number of
tilings within a stripe) * 2
(because there are two sets of tilings over this memory space). As
always, the memory_size
can be chosen within wide limits,
but it makes sense to chose it the same for all the groups of
tilings. Finally, notice that we change the final
integer input from 0 to 1 in the second call. Whenever using multiple
calls like this, it is a good idea to give them each a different
setting for the integer arguments.
There are a lot of possibilities here. One could use two groups of tilings, one for each single dimension, as illustrated above, PLUS an additional group of tilings that depends on both dimensions. This would create a system that generalizes partly along each dimension separately, but also on the two dimensions conjunctively, so that it could if needed learn to give a special response to a particular pair of x,y values. Another possibility would be to choose some groups of tilings with coarse generalization and some with fine.
One can control the tendency to generalize in the different ways given by each group of tilings by controlling the number of tilings in each group. For example, if you wanted discrimination to be stronger along x rather than y, then you would ask for more x tillings (more than 32) in the above example, or fewer y tilings.
In many applications one wants to make sure hashing, a random process, is done in a consistent way, e.g., across different runs of a program. The hashing done by the routines described here is determined by some tables filled with random numbers the first time the routine is called. To get consistent hashing you need only ensure that the random number generator is in a consistent state for the first call. A convenient way to do this is to make a dummy first call during an initialization portion of your program. For example, an early portion of your C++ code could include:
int dummy_tiles[1];
float dummy_vars[1];
srand(0);
tiles(dummy_tiles,1,1,dummy_vars,0); // initializes tiling code
*randomstate*
in
a reproducable manner, then call (tiles nil 1 1)
.random.seed(number)
, and then import the
tiles
module.The LISP function "tiles
" returns a list of
tile
indices. It is called as follows:
(tiles
numtilings ; the number of tilings desired
memorysize ; the number of possible tile indices
floats ; a list of real values making up the input vector
ints) ; optional list of integers to get different hashings
All the inputs have the same meanings and uses as described above.
Multiple calls
are combined by appending (or nconcing) the resultant lists into one
list.
The function "loadtiles
" is also available as a
maximally efficient interface.
It loads a portion of a provided array with the returned tile indices.
It is called as follows:
(loadtiles
tiles ; a provided array for the tile indices to go into
startingelement ; first element of "tiles" to be changed (typically 0)
numtilings ; the number of tilings desired
memorysize ; the number of possible tile indices
floats ; a list of real values making up the input vector
ints) ; list of optional inputs to get different hashings
The Python function "tiles
" returns a tuple of
tile
indices. It is called as follows:
tiles(numtilings, memsize, floats, ints)
All the inputs have the same meanings and uses as described above.
Multiple calls
are combined by concatenating the resultant tuples into one.
The function "loadtiles
" is also available as a
maximally efficient interface.
It loads a portion of a provided array with the returned tile indices.
It is called as follows:
loadtiles(tiles, startelement, numtilings, memsize, floats, ints)
with meanings as in the lisp calls.
In order to improve the efficiency of the tiling functions, but
still have the flexibility of
Python, there is also a version available for Python which calls the C
version of tiles. The calling
sequences are exactly the same as for the Python versions.
With the code described so far, there is a small probability that unrelated inputs will hash into some of the same tiles. In a group of tilings, usually there will be no more than one such "collision", so that it is not a big problem; the learning process will sort it out. There will not be a big effect on performance unless the memory is too small or the hash functions are poorly designed. Nevertheless, one sometimes wants to rule this out. (When one's program doesn't work, there is a tendency, deserved or not, to suspect a failure of the hashing function.)
To alleviate this kind of concern, we provide another "safe"
version
of the tiles routine that
does the additional work to prevent hashing collisions  or rather to
detect them and
take the necessary corrective steps. The result is a guaranteed exact
emulation
of the infinite tilings. Every different tile used will have a
different index. If more tiles are used than there is space within the memory_size
,
then an error
is signaled. Of course there is also a cost in terms of additional
memory and computation.
The additional memory is contained in an additional data structure, called a "collision table". To create and use a collision table in the original 4x4 example:
collision_table ct(512,1);
float float_array[2];
int tiles_array[32];
float_array[0] = x / 0.25;
float_array[1] = y / 0.25;
tiles(tiles_array,32,&ct,float_array,2);
In the declaration of the collision table, the 512 is the size of the
table. The size must be a power of two. The 1 in the declaration
indicates that this is to be a "safe" table; provide a 0 here if you
just
want to collect statistics on a regular table. Notice that the
collision
table is then used in place of the memory_size
in calls
to tiles
.
This is possible because of course the collision table has the memory_size
stored within it.
You can call ct.usage()
to determine how many
entries
in the table
are currently in use. You can also examine ct.calls
, ct.clearhits
,
and ct.collisions
to determine the number of times this
collision
table has been used (calls), has returned a result without a collision
(clearhits),
and has had collisions (collisions, which may be more than one per
call).
ct.reset
clears the collision table, making room
for a
new set of entries.
In Lisp, you construct a collisiontable by calling makect
:
(makect
memorysize ; the number of possible tile
indices
safetylevel) ; one of :safe, :unsafe, or
:supersafe
which returns the newly initialized collisiontable.
The memorysize
here must be a power of two, e.g., 1024;
otherwise it is
the same memorysize
as before. If the safetylevel
is :unsafe
, then
the collision table is used only to collect statistics; collisions are
still allowed
and not corrected for. Normally, a safetylevel
of :safe
is adequate, but :supersafe
is available to be
absolutely sure.
Once a collision table has been constructed, it is used in any
of
the routines
previously described simply by passing it in in place of the memorysize
input to the routine. For example, a lisp version of the 32 4x4 tilings
example, using the
safe version, would look
like
this:
(setq ctable (makect 512 :safe)) ; this is done once
(tiles 32 ctable (list (/ x .25) (/ y .25))) ; done for each x,y
In Python, you obtain a CollisionTable by calling CollisionTable(memsize,
safetyval)
, with meanings as above except safetyval
is one of 'safe', 'unsafe', or 'super safe'.
Once a collision table has been constructed, it is used in any of the
routines
previously described simply by passing it in in place of the memsize
input to the routine. For example, a python version of the 32 4x4
tilings example,
using the safe version, would
look like this:
ctable = CollisionTable(512, 'safe') # this is done once
tiles(32, ctable, (x/0.25, y/0.25)) # this is done for each x,y
The version where Python calls C does not have the 'super safe'
collision table options, since it has not been implemented in the C
version (the required comparisons would slow things down considerably).The tilings we have discussed so far stretch out to infinity with no need to specify a range for them. This is cool, but sometimes not what you want. Sometimes you want the variables to wraparound over some range. For example, you may have an angle variable that goes from 0 to 2pi and then should wrap around, that is, you would like generalization to occur between angles just greater than 0 and angles that are nearly 2pi. To accomodate this, we provide some alternative, wraparound versions of the tiling routines.
These versions take an additional input, wrapwidths
,
which parallels the float structure (array or list), and which
specifies for each float the width of the range over which it wraps. If
you don't want a float to wrap, it's wrap_width should
be zero (in LISP, nil). The wrap_width is in the same units as the
floats, but should be an integer. This can be confusing, so let's do a
simple 1D example. Suppose you have one real input, an angle theta.
Theta is originally in radians, that is, in [0,2pi), which you would
like to wrap around at 2pi. Remember that these functions all have
their tile boundaries at the integers, so if you passed in theta
directly there would be tile boundaries at 0, 1, and 2, i.e., just a
few tiles over the whole range, which is probably not what you want. So
let's say what we want! Suppose we want tilings with 10 tiles over the
[0,2pi) range. Then we have to scale theta so that it goes from 0 to 10
(instead of 0 to 2pi). One would do this by multiplying theta by
10/2pi. With the new scaled theta, the wrapping is over [0,10), for a
wrap_width of 10. Here is the C code for this case, assuming we want 16
tilings over the original [0,2pi) range with a memory size of 512:
int tiles_array[16];
float float_array[1];
int wrap_array[1];
float_array[0] = theta * 10 / (
2 * 3.14159);
wrap_array[0] = 10;
tileswrap(tiles_array,16,512,float_array,1,wrap_array);
Note that the code would be exactly the same if the original range of theta was [pi,+pi]. Specifying the complete range of wrapping (rather than just the width) is not necessary for the same reason as we did not need to give a range at all in the previous routines.
As another example, suppose you wanted to cover the 2pi x 3 rectangular area with 16 tilings, each with a width of generalization onetenth of the space in each dimension, and with wraparound in the first dimension but not the second. In C you would do
int tiles_array[16];
float float_array[2];
int wrap_array[2];
float_array[0] = x / (3 * 0.1);
wrap_array[0] = 0;
float_array[1] = y / (2 * 3.14159 * 0.1);
wrap_array[1] = 10;
tileswrap(tiles_array,16,512,float_array,2,wrap_array);
tileswrap
has the following interface:
void tileswrap(
int tiles[], // provided array will contain the returned tile indices
int num_tilings, // number of tile indices to be returned
int memory_size, // total number of possible tiles
float floats[], // array of floatingpoint vars to be converted to tiles
int num_floats, // number of floatingpoint variables
int wrap_widths[], // array of widths (length and units as in floats)
int ints[], // array of integer variables to be converted to tiles
int num_ints) // number of integer variables
Previous C versions had a bug because of an oddity of
the C
mod function, but this has been fixed.
In LISP, the wrapping routines are:
(tileswrap
numtilings ; the number of tilings desired
memorysize ; the number of possible tile indices
floats ; a list of real values making up the input vector
wrapwidths ; list of widths for wraparound on floats (ints or nil)
ints) ; optional list of integers to get different hashings
and
(loadtileswrap
tiles ; a provided array for the tile indices to go into
startingelement ; first element of "tiles" to be changed (typically 0)
numtilings ; the number of tilings desired
memorysize ; the number of possible tile indices
floats ; a list of real values making up the input vector
wrapwidths ; list of widths for wraparound on floats (ints or nil)
ints) ; list of optional inputs to get different hashings
In Python (and the Python calling C version), the wrapping routines
are:
tileswrap(numtilings,
memsize, floats, wrapwidths, ints)
loadtileswrap(tiles, startelement, numtilings, memsize, floats,
wrapwidths,
ints)