lattices v1.0.16

lattices

lattices.js helps making 2-d symmetric square and hexagonal lattice objects with optional periodic boundary conditions, size and scale options and easy access to nodes, node neighbors and cell geometries.

This can be useful for programming visualizations and agent based simulations in d3js.

Installation

Install the package as a node module:

npm install lattices

or clone this repository and install:

git clone https://github.com/dirkbrockmann/lattices.git
cd lattices
npm install
npm run build

View examples:

npm run examples

or start

http-server ./dist/

and open http://localhost:8080/examples/ in the browser.

Usage

Either load the package as a remote resource like so:

<script src="https://unpkg.com/lattices"></script>

<script>

const a = lattices.square(123)

</script>

Alternatively, use a local copy (dist/lattices.js) and include it in your html-file like so:

<script src="lattices.js"></script>

<script>
	
	const sq = lattices.square(5).boundary("dirichlet")
	
</script>

If you want to use it as part of your own project as a module import @dirkbrockmann/lattices like so:

import * as lattices from "lattices"

Setting up a lattice

You can generate square and hexagonal lattices, both either with periodic boundary conditions or not (Dirichlet boundaries).

Square Lattice

const N = 10
const sq = lattices.square(N)

creates a simple square lattice and stores it in sq. $N$ is the linear half-size of the lattice. Along each dimension the lattice has $2N+1$ nodes, the total number of nodes is therefore $(2N+1)^2$.

By default the lattice has periodic boundaries.

By default the lattice has a spatial scale $L=1$, so it fits into a square of spatial dimensions $L\times L$. So by default size of a node patch is ($dx\times dy$) with $dx=dy=(L/(2N+1))$.

The (x,y)-coordinates of the lattice range from -L/2 to L/2 in both dimensions:

$$x=-L/2+dx/2+n\, dx$$

$$y=-L/2+dy/2+m\, dy$$

with $n,m=0,...,2N$.

Hexagonal Lattices

const N = 10
const hx = lattices.hex(N) 

creates a hexagonal lattice stores it in hx. $2*N+1$ is the number of nodes along the horizontal axis. The total number of nodes in the network is $1+3N(N+1)$.

By default the lattice has periodic boundaries. Periodic boundary conditions are a bit tricky in hex lattices, but it can be worked out, see https://www.redblobgames.com/grids/hexagons/.

By default the lattice has a spatial scale $L=1$, so it fits into a rectangle of spatial dimensions $L\times \sqrt{3}L/2$. Each patch is a hexagon with a corner to corner distance $dx=(L/(2N+1))$.

Lattice properties

Both, square and hexagonal lattice have the following property fields. Say you have defined a lattice G = lattices.hex(10) or G = lattices.square(10) then

• G.N : is the linear node range $N$ (in this case 10)
• G.size : is the number of nodes in the lattice
• G.L: the linear physical extent
• G.type : returns its type, i.e. "hexagonal" / "square"

Nodes

Nodes can be accessed by G.nodes. Each node in the array has x and y coordinates and a list of neighbors array. e.g.

• G.nodes[i].x : returns node i's x-coordinate
• G.nodes[i].y : returns node i's y-coordinate
• G.nodes[i].neighbors : returns an array of all of i's neighbors
• G.nodes[i].cell(): returns a list of coordinates to draw the nodes cell

Scale

G.scale(s)

sets the spatial scale of the lattice to s. The default value is 1. This effects the x and y coordinates of the nodes and their boundaries.

Without argument, the current scale of G is returned.

Boundary Condition

G.boundary(["periodic"|"dirichlet"])

sets the type of boundary condition. Default is "periodic". For dirichlet the nodes on the boundary have a different number of neighbors.

Without an argument, returns the lattice's boundary type.

Square lattice n8 and n4 neighborhoods

For a square lattice you can chose a node's neighborhood to have a neighborhood of the 8 surrounding lattice nodes or only the 4 nodes corresponding to above, below, left and right. Like so

const sq = lattices.square(10).hood("n4")

or

const sq = lattices.square(10).hood("n8")

Without arguments hood() returns the neighborhood type.