docs/LLLplusTutorial.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# LLLplus Tutorial"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Lattice reduction and related lattice tools are used in cryptography, digital communication, and integer programming. The [LLLplus.jl](https://github.com/christianpeel/LLLplus.jl) package includes several lattice reduction tools, a CVP problem solver, and more. This notebook gives an introductory tutorial on lattices and lattice reduction and then shows the application of LLLplus tools in wireless communications, cryptography, and integer programming.\n",
    "\n",
    "### Intro to lattices\n",
    "\n",
    "What is a lattice? For an $m$ by $n$ matrix $H$ with $m\\ge n$, then the set of points $\\mathbf{y}=H\\mathbf{x}$ for any integer valued vector $\\mathbf{x}$ is a lattice.\n",
    "\n",
    "Let's illustrate this with an $m=n=2$ case. Rather than looking all possible vectors $\\mathbf{x}$, let's look at vectors with values from $-2,\\ldots,2$, represented by the columns of the following matrix $X$:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "X = [ -2  -2  -2  -2  -2  -1  -1  -1  -1  -1   0   0  0  0  0   1   1  1  1  1   2   2  2  2  2\n",
    "      -2  -1   0   1   2  -2  -1   0   1   2  -2  -1  0  1  2  -2  -1  0  1  2  -2  -1  0  1  2];"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's see what these points look like:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\"?>\n",
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600\" height=\"400\" viewBox=\"0 0 2400 1600\">\n",
       "<defs>\n",
       "  <clipPath id=\"clip0200\">\n",
       "    <rect x=\"0\" y=\"0\" width=\"2400\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polygon clip-path=\"url(#clip0200)\" points=\"\n",
       "0,1600 2400,1600 2400,0 0,0 \n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip0201\">\n",
       "    <rect x=\"480\" y=\"0\" width=\"1681\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polygon clip-path=\"url(#clip0200)\" points=\"\n",
       "149.361,1503.47 2321.26,1503.47 2321.26,47.2441 149.361,47.2441 \n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip0202\">\n",
       "    <rect x=\"149\" y=\"47\" width=\"2173\" height=\"1457\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  210.829,1503.47 210.829,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  723.07,1503.47 723.07,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1235.31,1503.47 1235.31,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1747.55,1503.47 1747.55,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  2259.79,1503.47 2259.79,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,1462.26 2321.26,1462.26 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,1118.81 2321.26,1118.81 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,775.359 2321.26,775.359 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,431.909 2321.26,431.909 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,88.4582 2321.26,88.4582 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1503.47 2321.26,1503.47 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1503.47 149.361,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  210.829,1503.47 210.829,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  723.07,1503.47 723.07,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1235.31,1503.47 1235.31,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1747.55,1503.47 1747.55,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  2259.79,1503.47 2259.79,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1462.26 181.939,1462.26 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1118.81 181.939,1118.81 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,775.359 181.939,775.359 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,431.909 181.939,431.909 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0200)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,88.4582 181.939,88.4582 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 210.829, 1557.47)\" x=\"210.829\" y=\"1557.47\">-2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 723.07, 1557.47)\" x=\"723.07\" y=\"1557.47\">-1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1235.31, 1557.47)\" x=\"1235.31\" y=\"1557.47\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1747.55, 1557.47)\" x=\"1747.55\" y=\"1557.47\">1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 2259.79, 1557.47)\" x=\"2259.79\" y=\"1557.47\">2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1479.76)\" x=\"125.361\" y=\"1479.76\">-2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1136.31)\" x=\"125.361\" y=\"1136.31\">-1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 792.859)\" x=\"125.361\" y=\"792.859\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 449.409)\" x=\"125.361\" y=\"449.409\">1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0200)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 105.958)\" x=\"125.361\" y=\"105.958\">2</text>\n",
       "</g>\n",
       "<polyline clip-path=\"url(#clip0202)\" style=\"stroke:#009af9; stroke-width:4; stroke-opacity:0; fill:none\" points=\"\n",
       "  210.829,1462.26 210.829,1118.81 210.829,775.359 210.829,431.909 210.829,88.4582 723.07,1462.26 723.07,1118.81 723.07,775.359 723.07,431.909 723.07,88.4582 \n",
       "  1235.31,1462.26 1235.31,1118.81 1235.31,775.359 1235.31,431.909 1235.31,88.4582 1747.55,1462.26 1747.55,1118.81 1747.55,775.359 1747.55,431.909 1747.55,88.4582 \n",
       "  2259.79,1462.26 2259.79,1118.81 2259.79,775.359 2259.79,431.909 2259.79,88.4582 \n",
       "  \"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"1462.26\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"1462.26\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"1118.81\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"1118.81\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"431.909\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"431.909\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"88.4582\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"88.4582\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"1462.26\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"1462.26\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"1118.81\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"1118.81\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"431.909\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"431.909\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"88.4582\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"88.4582\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1462.26\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1462.26\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1118.81\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1118.81\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"431.909\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"431.909\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"88.4582\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"88.4582\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"1462.26\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"1462.26\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"1118.81\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"1118.81\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"431.909\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"431.909\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"88.4582\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"88.4582\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"1462.26\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"1462.26\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"1118.81\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"1118.81\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"431.909\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"431.909\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"88.4582\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0202)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"88.4582\" r=\"14\"/>\n",
       "</svg>\n"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "using Plots\n",
    "plot(X[1,:],X[2,:], linealpha = 0,markershape=:circle,legend=false)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Hopefully this plot looks like a subset of $\\mathbb{Z}^2$ to you.\n",
    "\n",
    "So what about a lattice?  Let's look at a specially-chosen lattice, illustrated by the matrix $H$ below, and plot it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2×25 Array{Float64,2}:\n",
       " -2.0  -3.5       -5.0      -6.5      …  6.5      5.0      3.5       2.0\n",
       "  0.0  -0.866025  -1.73205  -2.59808     2.59808  1.73205  0.866025  0.0"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "H = [ 2.5       -1.5     \n",
    " 0.866025  -0.866025]\n",
    "Y = H*X"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\"?>\n",
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600\" height=\"400\" viewBox=\"0 0 2400 1600\">\n",
       "<defs>\n",
       "  <clipPath id=\"clip0600\">\n",
       "    <rect x=\"0\" y=\"0\" width=\"2400\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polygon clip-path=\"url(#clip0600)\" points=\"\n",
       "0,1600 2400,1600 2400,0 0,0 \n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip0601\">\n",
       "    <rect x=\"480\" y=\"0\" width=\"1681\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polygon clip-path=\"url(#clip0600)\" points=\"\n",
       "149.361,1503.47 2321.26,1503.47 2321.26,47.2441 149.361,47.2441 \n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip0602\">\n",
       "    <rect x=\"149\" y=\"47\" width=\"2173\" height=\"1457\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  466.95,1503.47 466.95,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  851.13,1503.47 851.13,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1235.31,1503.47 1235.31,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1619.49,1503.47 1619.49,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  2003.67,1503.47 2003.67,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,1370.23 2321.26,1370.23 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,1171.94 2321.26,1171.94 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,973.651 2321.26,973.651 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,775.359 2321.26,775.359 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,577.068 2321.26,577.068 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,378.777 2321.26,378.777 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,180.485 2321.26,180.485 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1503.47 2321.26,1503.47 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1503.47 149.361,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  466.95,1503.47 466.95,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  851.13,1503.47 851.13,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1235.31,1503.47 1235.31,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1619.49,1503.47 1619.49,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  2003.67,1503.47 2003.67,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1370.23 181.939,1370.23 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1171.94 181.939,1171.94 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,973.651 181.939,973.651 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,775.359 181.939,775.359 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,577.068 181.939,577.068 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,378.777 181.939,378.777 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0600)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,180.485 181.939,180.485 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 466.95, 1557.47)\" x=\"466.95\" y=\"1557.47\">-6</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 851.13, 1557.47)\" x=\"851.13\" y=\"1557.47\">-3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1235.31, 1557.47)\" x=\"1235.31\" y=\"1557.47\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1619.49, 1557.47)\" x=\"1619.49\" y=\"1557.47\">3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 2003.67, 1557.47)\" x=\"2003.67\" y=\"1557.47\">6</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1387.73)\" x=\"125.361\" y=\"1387.73\">-3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1189.44)\" x=\"125.361\" y=\"1189.44\">-2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 991.151)\" x=\"125.361\" y=\"991.151\">-1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 792.859)\" x=\"125.361\" y=\"792.859\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 594.568)\" x=\"125.361\" y=\"594.568\">1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 396.277)\" x=\"125.361\" y=\"396.277\">2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0600)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 197.985)\" x=\"125.361\" y=\"197.985\">3</text>\n",
       "</g>\n",
       "<polyline clip-path=\"url(#clip0602)\" style=\"stroke:#009af9; stroke-width:4; stroke-opacity:0; fill:none\" points=\"\n",
       "  979.19,775.359 787.1,947.085 595.01,1118.81 402.92,1290.54 210.829,1462.26 1299.34,603.634 1107.25,775.359 915.16,947.085 723.07,1118.81 530.98,1290.54 \n",
       "  1619.49,431.909 1427.4,603.634 1235.31,775.359 1043.22,947.085 851.13,1118.81 1939.64,260.183 1747.55,431.909 1555.46,603.634 1363.37,775.359 1171.28,947.085 \n",
       "  2259.79,88.4582 2067.7,260.183 1875.61,431.909 1683.52,603.634 1491.43,775.359 \n",
       "  \"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"979.19\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"979.19\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"787.1\" cy=\"947.085\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"787.1\" cy=\"947.085\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"595.01\" cy=\"1118.81\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"595.01\" cy=\"1118.81\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"402.92\" cy=\"1290.54\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"402.92\" cy=\"1290.54\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"1462.26\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"210.829\" cy=\"1462.26\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1299.34\" cy=\"603.634\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1299.34\" cy=\"603.634\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1107.25\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1107.25\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"915.16\" cy=\"947.085\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"915.16\" cy=\"947.085\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"1118.81\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"723.07\" cy=\"1118.81\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"530.98\" cy=\"1290.54\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"530.98\" cy=\"1290.54\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1619.49\" cy=\"431.909\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1619.49\" cy=\"431.909\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1427.4\" cy=\"603.634\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1427.4\" cy=\"603.634\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1043.22\" cy=\"947.085\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1043.22\" cy=\"947.085\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"851.13\" cy=\"1118.81\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"851.13\" cy=\"1118.81\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1939.64\" cy=\"260.183\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1939.64\" cy=\"260.183\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"431.909\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1747.55\" cy=\"431.909\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1555.46\" cy=\"603.634\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1555.46\" cy=\"603.634\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1363.37\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1363.37\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1171.28\" cy=\"947.085\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1171.28\" cy=\"947.085\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"88.4582\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2259.79\" cy=\"88.4582\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2067.7\" cy=\"260.183\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2067.7\" cy=\"260.183\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1875.61\" cy=\"431.909\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1875.61\" cy=\"431.909\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1683.52\" cy=\"603.634\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1683.52\" cy=\"603.634\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1491.43\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip0602)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1491.43\" cy=\"775.359\" r=\"14\"/>\n",
       "</svg>\n"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plot(Y[1,:],Y[2,:], linealpha = 0,markershape=:circle,legend=false)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Despite the different shapes, the first two figures above contain the *same* number of points.\n",
    "\n",
    "A *lattice* really means that we are spanning the entire $\\mathbb{Z}^2$ space; in the figure above we're only looking a subset of those points.  To get a better feel for what's going on, we'll define the function 'count' below. The only thing you need to know is that it allows us to look at a larger set of 'lattice' points without writting the explicit `X` vector down."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2×169 Array{Float64,2}:\n",
       " -6.0  -7.5       -9.0      -10.5      …  10.5      9.0      7.5       6.0\n",
       "  0.0  -0.866025  -1.73205   -2.59808      2.59808  1.73205  0.866025  0.0"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "function count(M,N)\n",
    "    s = zeros(M,N^M)\n",
    "    z = (1:N^M) .-1\n",
    "    for n = 0:M-1\n",
    "        m = floor.(z/N)\n",
    "        s[M-n,:] = z .- N.*m\n",
    "        z = m\n",
    "    end\n",
    "    return s\n",
    "end\n",
    "X = Int.(count(2,13) .- 6)\n",
    "Y = H*X"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We'll plot the same lattice, with the same plot window (same x and y domain), just with the larger number of points in X."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\"?>\n",
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600\" height=\"400\" viewBox=\"0 0 2400 1600\">\n",
       "<defs>\n",
       "  <clipPath id=\"clip1000\">\n",
       "    <rect x=\"0\" y=\"0\" width=\"2400\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polygon clip-path=\"url(#clip1000)\" points=\"\n",
       "0,1600 2400,1600 2400,0 0,0 \n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip1001\">\n",
       "    <rect x=\"480\" y=\"0\" width=\"1681\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polygon clip-path=\"url(#clip1000)\" points=\"\n",
       "149.361,1503.47 2321.26,1503.47 2321.26,47.2441 149.361,47.2441 \n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip1002\">\n",
       "    <rect x=\"149\" y=\"47\" width=\"2173\" height=\"1457\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  420.848,1503.47 420.848,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  828.079,1503.47 828.079,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1235.31,1503.47 1235.31,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1642.54,1503.47 1642.54,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  2049.77,1503.47 2049.77,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,1399.46 2321.26,1399.46 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,1191.43 2321.26,1191.43 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,983.392 2321.26,983.392 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,775.359 2321.26,775.359 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,567.326 2321.26,567.326 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,359.294 2321.26,359.294 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,151.261 2321.26,151.261 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1503.47 2321.26,1503.47 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1503.47 149.361,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  420.848,1503.47 420.848,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  828.079,1503.47 828.079,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1235.31,1503.47 1235.31,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1642.54,1503.47 1642.54,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  2049.77,1503.47 2049.77,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1399.46 181.939,1399.46 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1191.43 181.939,1191.43 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,983.392 181.939,983.392 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,775.359 181.939,775.359 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,567.326 181.939,567.326 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,359.294 181.939,359.294 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1000)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,151.261 181.939,151.261 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 420.848, 1557.47)\" x=\"420.848\" y=\"1557.47\">-6</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 828.079, 1557.47)\" x=\"828.079\" y=\"1557.47\">-3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1235.31, 1557.47)\" x=\"1235.31\" y=\"1557.47\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1642.54, 1557.47)\" x=\"1642.54\" y=\"1557.47\">3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 2049.77, 1557.47)\" x=\"2049.77\" y=\"1557.47\">6</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1416.96)\" x=\"125.361\" y=\"1416.96\">-3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1208.93)\" x=\"125.361\" y=\"1208.93\">-2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1000.89)\" x=\"125.361\" y=\"1000.89\">-1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 792.859)\" x=\"125.361\" y=\"792.859\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 584.826)\" x=\"125.361\" y=\"584.826\">1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 376.794)\" x=\"125.361\" y=\"376.794\">2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1000)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 168.761)\" x=\"125.361\" y=\"168.761\">3</text>\n",
       "</g>\n",
       "<polyline clip-path=\"url(#clip1002)\" style=\"stroke:#009af9; stroke-width:4; stroke-opacity:0; fill:none\" points=\"\n",
       "  420.848,775.359 217.232,955.521 13.6169,1135.68 -189.999,1315.84 -393.614,1496.01 -597.23,1676.17 -800.845,1856.33 -1004.46,2036.49 -1208.08,2216.65 -1411.69,2396.82 \n",
       "  -1615.31,2576.98 -1818.92,2757.14 -2022.54,2937.3 760.207,595.198 556.592,775.359 352.976,955.521 149.361,1135.68 -54.2549,1315.84 -257.87,1496.01 -461.486,1676.17 \n",
       "  -665.102,1856.33 -868.717,2036.49 -1072.33,2216.65 -1275.95,2396.82 -1479.56,2576.98 -1683.18,2757.14 1099.57,415.036 895.951,595.198 692.335,775.359 488.72,955.521 \n",
       "  285.104,1135.68 81.4888,1315.84 -122.127,1496.01 -325.742,1676.17 -529.358,1856.33 -732.973,2036.49 -936.589,2216.65 -1140.2,2396.82 -1343.82,2576.98 1438.93,234.874 \n",
       "  1235.31,415.036 1031.69,595.198 828.079,775.359 624.464,955.521 420.848,1135.68 217.232,1315.84 13.6169,1496.01 -189.999,1676.17 -393.614,1856.33 -597.23,2036.49 \n",
       "  -800.845,2216.65 -1004.46,2396.82 1778.29,54.7125 1574.67,234.874 1371.05,415.036 1167.44,595.198 963.823,775.359 760.207,955.521 556.592,1135.68 352.976,1315.84 \n",
       "  149.361,1496.01 -54.2549,1676.17 -257.87,1856.33 -461.486,2036.49 -665.102,2216.65 2117.64,-125.449 1914.03,54.7125 1710.41,234.874 1506.8,415.036 1303.18,595.198 \n",
       "  1099.57,775.359 895.951,955.521 692.335,1135.68 488.72,1315.84 285.104,1496.01 81.4888,1676.17 -122.127,1856.33 -325.742,2036.49 2457,-305.611 2253.39,-125.449 \n",
       "  2049.77,54.7125 1846.16,234.874 1642.54,415.036 1438.93,595.198 1235.31,775.359 1031.69,955.521 828.079,1135.68 624.464,1315.84 420.848,1496.01 217.232,1676.17 \n",
       "  13.6169,1856.33 2796.36,-485.773 2592.75,-305.611 2389.13,-125.449 2185.52,54.7125 1981.9,234.874 1778.29,415.036 1574.67,595.198 1371.05,775.359 1167.44,955.521 \n",
       "  963.823,1135.68 760.207,1315.84 556.592,1496.01 352.976,1676.17 3135.72,-665.934 2932.11,-485.773 2728.49,-305.611 2524.88,-125.449 2321.26,54.7125 2117.64,234.874 \n",
       "  1914.03,415.036 1710.41,595.198 1506.8,775.359 1303.18,955.521 1099.57,1135.68 895.951,1315.84 692.335,1496.01 3475.08,-846.096 3271.47,-665.934 3067.85,-485.773 \n",
       "  2864.23,-305.611 2660.62,-125.449 2457,54.7125 2253.39,234.874 2049.77,415.036 1846.16,595.198 1642.54,775.359 1438.93,955.521 1235.31,1135.68 1031.69,1315.84 \n",
       "  3814.44,-1026.26 3610.83,-846.096 3407.21,-665.934 3203.59,-485.773 2999.98,-305.611 2796.36,-125.449 2592.75,54.7125 2389.13,234.874 2185.52,415.036 1981.9,595.198 \n",
       "  1778.29,775.359 1574.67,955.521 1371.05,1135.68 4153.8,-1206.42 3950.18,-1026.26 3746.57,-846.096 3542.95,-665.934 3339.34,-485.773 3135.72,-305.611 2932.11,-125.449 \n",
       "  2728.49,54.7125 2524.88,234.874 2321.26,415.036 2117.64,595.198 1914.03,775.359 1710.41,955.521 4493.16,-1386.58 4289.54,-1206.42 4085.93,-1026.26 3882.31,-846.096 \n",
       "  3678.7,-665.934 3475.08,-485.773 3271.47,-305.611 3067.85,-125.449 2864.23,54.7125 2660.62,234.874 2457,415.036 2253.39,595.198 2049.77,775.359 \n",
       "  \"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"217.232\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"217.232\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"352.976\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"352.976\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"285.104\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"285.104\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"217.232\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"217.232\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"352.976\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"352.976\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"285.104\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"285.104\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2185.52\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2185.52\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2321.26\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2321.26\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2117.64\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2117.64\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2253.39\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2253.39\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2185.52\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2185.52\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2321.26\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2321.26\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2117.64\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2117.64\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2253.39\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2253.39\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1002)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"775.359\" r=\"14\"/>\n",
       "</svg>\n"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plot(Y[1,:],Y[2,:], linealpha = 0,markershape=:circle,legend=false,xlims=(-8,8),ylims=(-3.5,3.5))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Lattice Reduction\n",
    "Since the lattice is over the entire $\\mathbb{Z}^2$ space, it's natural to ask if there's a different matrix which will generate the same points. Hopefully, this new lattice will be less skew.\n",
    "\n",
    "This is what [Lenstra, Lenstra, Lovász](https://en.wikipedia.org/wiki/Lenstra%E2%80%93Lenstra%E2%80%93Lov%C3%A1sz_lattice_basis_reduction_algorithm) (LLL) lattice reduction does.  It decomposes a matrix `H` into a matrix `B` that is still a basis for the lattice that `H` makes, yet has shorter vectors and is closer to orthogonal. It also returns a unimodular (determinant= +/-1) integer valued matrix `T` such that `H T= B`.   \n",
    "\n",
    "Let's take the LLL decomposition of our `H` from above and plot the resulting lattice."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2×2 Array{Float64,2}:\n",
       " 1.0  -0.5     \n",
       " 0.0   0.866025"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "using LLLplus\n",
    "B,T = lll(H); B"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\"?>\n",
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600\" height=\"400\" viewBox=\"0 0 2400 1600\">\n",
       "<defs>\n",
       "  <clipPath id=\"clip1400\">\n",
       "    <rect x=\"0\" y=\"0\" width=\"2400\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polygon clip-path=\"url(#clip1400)\" points=\"\n",
       "0,1600 2400,1600 2400,0 0,0 \n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip1401\">\n",
       "    <rect x=\"480\" y=\"0\" width=\"1681\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polygon clip-path=\"url(#clip1400)\" points=\"\n",
       "149.361,1503.47 2321.26,1503.47 2321.26,47.2441 149.361,47.2441 \n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip1402\">\n",
       "    <rect x=\"149\" y=\"47\" width=\"2173\" height=\"1457\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  420.848,1503.47 420.848,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  828.079,1503.47 828.079,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1235.31,1503.47 1235.31,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1642.54,1503.47 1642.54,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  2049.77,1503.47 2049.77,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,1399.46 2321.26,1399.46 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,1191.43 2321.26,1191.43 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,983.392 2321.26,983.392 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,775.359 2321.26,775.359 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,567.326 2321.26,567.326 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,359.294 2321.26,359.294 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  149.361,151.261 2321.26,151.261 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1503.47 2321.26,1503.47 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1503.47 149.361,47.2441 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  420.848,1503.47 420.848,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  828.079,1503.47 828.079,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1235.31,1503.47 1235.31,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1642.54,1503.47 1642.54,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  2049.77,1503.47 2049.77,1481.63 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1399.46 181.939,1399.46 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,1191.43 181.939,1191.43 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,983.392 181.939,983.392 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,775.359 181.939,775.359 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,567.326 181.939,567.326 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,359.294 181.939,359.294 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip1400)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  149.361,151.261 181.939,151.261 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 420.848, 1557.47)\" x=\"420.848\" y=\"1557.47\">-6</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 828.079, 1557.47)\" x=\"828.079\" y=\"1557.47\">-3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1235.31, 1557.47)\" x=\"1235.31\" y=\"1557.47\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1642.54, 1557.47)\" x=\"1642.54\" y=\"1557.47\">3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 2049.77, 1557.47)\" x=\"2049.77\" y=\"1557.47\">6</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1416.96)\" x=\"125.361\" y=\"1416.96\">-3</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1208.93)\" x=\"125.361\" y=\"1208.93\">-2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 1000.89)\" x=\"125.361\" y=\"1000.89\">-1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 792.859)\" x=\"125.361\" y=\"792.859\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 584.826)\" x=\"125.361\" y=\"584.826\">1</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 376.794)\" x=\"125.361\" y=\"376.794\">2</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip1400)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 125.361, 168.761)\" x=\"125.361\" y=\"168.761\">3</text>\n",
       "</g>\n",
       "<polyline clip-path=\"url(#clip1402)\" style=\"stroke:#009af9; stroke-width:4; stroke-opacity:0; fill:none\" points=\"\n",
       "  828.079,1856.33 760.207,1676.17 692.335,1496.01 624.464,1315.84 556.592,1135.68 488.72,955.521 420.848,775.359 352.976,595.198 285.104,415.036 217.232,234.874 \n",
       "  149.361,54.7125 81.4888,-125.449 13.6169,-305.611 963.823,1856.33 895.951,1676.17 828.079,1496.01 760.207,1315.84 692.335,1135.68 624.464,955.521 556.592,775.359 \n",
       "  488.72,595.198 420.848,415.036 352.976,234.874 285.104,54.7125 217.232,-125.449 149.361,-305.611 1099.57,1856.33 1031.69,1676.17 963.823,1496.01 895.951,1315.84 \n",
       "  828.079,1135.68 760.207,955.521 692.335,775.359 624.464,595.198 556.592,415.036 488.72,234.874 420.848,54.7125 352.976,-125.449 285.104,-305.611 1235.31,1856.33 \n",
       "  1167.44,1676.17 1099.57,1496.01 1031.69,1315.84 963.823,1135.68 895.951,955.521 828.079,775.359 760.207,595.198 692.335,415.036 624.464,234.874 556.592,54.7125 \n",
       "  488.72,-125.449 420.848,-305.611 1371.05,1856.33 1303.18,1676.17 1235.31,1496.01 1167.44,1315.84 1099.57,1135.68 1031.69,955.521 963.823,775.359 895.951,595.198 \n",
       "  828.079,415.036 760.207,234.874 692.335,54.7125 624.464,-125.449 556.592,-305.611 1506.8,1856.33 1438.93,1676.17 1371.05,1496.01 1303.18,1315.84 1235.31,1135.68 \n",
       "  1167.44,955.521 1099.57,775.359 1031.69,595.198 963.823,415.036 895.951,234.874 828.079,54.7125 760.207,-125.449 692.335,-305.611 1642.54,1856.33 1574.67,1676.17 \n",
       "  1506.8,1496.01 1438.93,1315.84 1371.05,1135.68 1303.18,955.521 1235.31,775.359 1167.44,595.198 1099.57,415.036 1031.69,234.874 963.823,54.7125 895.951,-125.449 \n",
       "  828.079,-305.611 1778.29,1856.33 1710.41,1676.17 1642.54,1496.01 1574.67,1315.84 1506.8,1135.68 1438.93,955.521 1371.05,775.359 1303.18,595.198 1235.31,415.036 \n",
       "  1167.44,234.874 1099.57,54.7125 1031.69,-125.449 963.823,-305.611 1914.03,1856.33 1846.16,1676.17 1778.29,1496.01 1710.41,1315.84 1642.54,1135.68 1574.67,955.521 \n",
       "  1506.8,775.359 1438.93,595.198 1371.05,415.036 1303.18,234.874 1235.31,54.7125 1167.44,-125.449 1099.57,-305.611 2049.77,1856.33 1981.9,1676.17 1914.03,1496.01 \n",
       "  1846.16,1315.84 1778.29,1135.68 1710.41,955.521 1642.54,775.359 1574.67,595.198 1506.8,415.036 1438.93,234.874 1371.05,54.7125 1303.18,-125.449 1235.31,-305.611 \n",
       "  2185.52,1856.33 2117.64,1676.17 2049.77,1496.01 1981.9,1315.84 1914.03,1135.68 1846.16,955.521 1778.29,775.359 1710.41,595.198 1642.54,415.036 1574.67,234.874 \n",
       "  1506.8,54.7125 1438.93,-125.449 1371.05,-305.611 2321.26,1856.33 2253.39,1676.17 2185.52,1496.01 2117.64,1315.84 2049.77,1135.68 1981.9,955.521 1914.03,775.359 \n",
       "  1846.16,595.198 1778.29,415.036 1710.41,234.874 1642.54,54.7125 1574.67,-125.449 1506.8,-305.611 2457,1856.33 2389.13,1676.17 2321.26,1496.01 2253.39,1315.84 \n",
       "  2185.52,1135.68 2117.64,955.521 2049.77,775.359 1981.9,595.198 1914.03,415.036 1846.16,234.874 1778.29,54.7125 1710.41,-125.449 1642.54,-305.611 \n",
       "  \"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"352.976\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"352.976\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"285.104\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"285.104\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"217.232\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"217.232\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"352.976\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"352.976\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"285.104\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"285.104\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"488.72\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"420.848\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"624.464\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"556.592\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"760.207\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"692.335\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"895.951\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"828.079\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1031.69\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"963.823\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1167.44\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1099.57\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1303.18\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1235.31\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1438.93\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1371.05\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1574.67\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1506.8\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2185.52\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2185.52\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2117.64\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2117.64\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1710.41\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1642.54\" cy=\"54.7125\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2321.26\" cy=\"1496.01\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2321.26\" cy=\"1496.01\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2253.39\" cy=\"1315.84\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2253.39\" cy=\"1315.84\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2185.52\" cy=\"1135.68\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2185.52\" cy=\"1135.68\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2117.64\" cy=\"955.521\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2117.64\" cy=\"955.521\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"775.359\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"2049.77\" cy=\"775.359\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"595.198\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1981.9\" cy=\"595.198\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"415.036\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1914.03\" cy=\"415.036\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"234.874\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1846.16\" cy=\"234.874\" r=\"14\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#000000; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"54.7125\" r=\"18\"/>\n",
       "<circle clip-path=\"url(#clip1402)\" style=\"fill:#009af9; stroke:none; fill-opacity:1\" cx=\"1778.29\" cy=\"54.7125\" r=\"14\"/>\n",
       "</svg>\n"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Y = B*X\n",
    "plot(Y[1,:],Y[2,:], linealpha = 0,markershape=:circle,legend=false,xlims=(-8,8),ylims=(-3.5,3.5))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If you compare the two figures above, you'll see that they cover the same lattice points.  The lattice basis is also less skew. \n",
    "\n",
    "What does that matrix `T` contain?  Does `T` really have determinant of +/- 1?  The lines below show that this is true.  You can check yourself that the inverse of `T` is also integer-valued, and the determinant of the inverse is also +/- 1."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2×2 Array{Int64,2}:\n",
       " 1  -2\n",
       " 1  -3"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "T"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "-1.0"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "det(T)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### CVP\n",
    "\n",
    "The \"closest vector problem\" is to find the closest point in a lattice to a point in the same underlying space.  Specifically, given the lattice induced by a  matrix `H`, and a vector `y`, find the integer-valued vector `x` that such that `Hx` is closest to `y`.  \n",
    "\n",
    "$$\\hat{\\mathbf{x}} = \\arg\\min_{x \\in \\mathbf{Z}^n} ||y-Hx||_2$$\n",
    "\n",
    "Let's look at an example:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2-element Array{Float64,1}:\n",
       " 0.029940685725328726\n",
       " 2.1336665272450257  "
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "H=[1 2; 3 4]\n",
    "x=[2,-1]\n",
    "y=H*x + randn(2)/10"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Without noise `Hx` would be `[0,2]`; the addition of noise moves it away from that baseline.\n",
    "\n",
    "The `cvp` function in LLLplus requires that we use a QR decomposition to prepare the data. Runing the code, we find we can recover the initial `x` vector."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2-element Array{Float64,1}:\n",
       "  2.0\n",
       " -1.0"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Q,R=qr(H)\n",
    "x̂ = cvp(Q'*y,R)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Application to Digital communication\n",
    "\n",
    "A wireless system with multiple antennas at both the transmitter and receiver is often referred to as a multiple-input, multiple-output (MIMO) system. Some engineers and researchers have proposed \"Massive-MIMO\" systems with hundreds of antennas. Let's look at decoding a binary signal sent on each antenna of such a system.  *Some details that everyone but wireless gurus can ignore: we'll assume a synchronized, narrow-band system, with no correlation between antennas and no impairments other than thermal noise*. Let's assume `n` antennas at each transmitter and receiver, with the propogation channel between each transmit and each receive antenna represented by one entry in a `n` by `n` matrix `H`.  The propgation matrix `H` reflects on how the transmitted signal bounced off of buildings and other metal surfaces between the transmitter and receiver. We will also add gaussian thermal noise `w`. The overal model is \n",
    "    $$\\mathbf{y} = H\\mathbf{x} + \\mathbf{w}$$\n",
    "Let's simulate this system:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "n=100                   # Let's use 100 antennas\n",
    "x=Int.(rand(0:1,n))     # A binary signal is transmitted from each tx anteanna\n",
    "H=randn(n,n)            # channel coefficients btwn each tx and rx antenna\n",
    "\n",
    "# send the signal through the channel\n",
    "w = rand(n)/10          # thermal noise\n",
    "y=H*x+w                 # received data\n",
    "\n",
    "# decode the signal\n",
    "Q,R=qr(H)               # The receiver is assumed to know H via training (not modeled)\n",
    "x̂=cvp(Q'*y,R)\n",
    "\n",
    "# check if we decoded the signal correctly\n",
    "sum(abs.(x-x̂))          # Success!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The signal is decoded correctly.\n",
    "\n",
    "Use of lattice tools such as a cvp solver, or using lattice-reduction typically gives much better performance to a broadband wireless system than simple linear receivers (such as `round(inv(H)*y)` ). These principles apply to LTE, WiFi, and many other digital communications problems. It's possible that the difference between the figures in this [article](https://9to5mac.com/2016/10/20/iphone-7-qualcomm-intel-modem-cellular-performance/) is that Qualcomm is using a nonlinear (CVP, LLL, etc...) receiver, while Intel is using a simpler, but poorer-performance linear receiver."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Subset-sum problem\n",
    "\n",
    "In the subset-sum problem, one is given a set of numbers `a`, a sum `s`, with the goal of finding whether some subset of `a` sums to `s`.  Specifically, if `a` is a vector, then is there a binary vector `x` such that `x'*a==s`?  The `subsetsum` function implements the Lagarias and Odlyzko (1983) technique which uses the LLL algorithm to find a polynomial-time solution to the problem. The LLL has been used to crack the [Merkle-Hellman](https://en.wikipedia.org/wiki/Merkle%E2%80%93Hellman_knapsack_cryptosystem) cryptosystem, and has been part of many other cryptoanalytic attacks.\n",
    "\n",
    "We do not dive into the details of the `subsetsum` function, rather we give two examples, starting with an example that can be visually inspected to be correct."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "6-element Array{Float64,1}:\n",
       " 1.0\n",
       " 0.0\n",
       " 0.0\n",
       " 1.0\n",
       " 1.0\n",
       " 0.0"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a=[1,2,4,9,20,38]\n",
    "s=30\n",
    "x=subsetsum(a,s)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "true"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# check the result\n",
    "s==x'*a"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The second example sets the vector `a` to be integers sampled uniformly from the range `[1:2^256]`.  Despite the large bit-depth of the integers, the LLL-based technique solves the problem quickly."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "40-element Array{BigInt,1}:\n",
       "  52966019994847903947942438279253224710179175905054670645485359246712459937580\n",
       "  60284285807960226905685998406630132778973408018084486443877102388277646793029\n",
       "  88088434549054603795035220275158003859278583319576200298335243431186049890688\n",
       " 102216526612158525608427694592741849895831500557382123931052317507192769873414\n",
       " 111366934726704812336357159967275227136813435253283639477131741118741497414209\n",
       "  68312097531131194452312516287134976387408608625158155795157983708319827992525\n",
       "  55125975044263260301331673399718883155408719786422595646401499451424633490680\n",
       "  47238834081716529338009400277108320358882359342047456245665496785110661486235\n",
       "  72736875888792828199967782662535352416765087576704047857802048708869240140699\n",
       "  16782831561841572851701628030674265972955034348219805218945406149166119288657\n",
       " 111518864133851420097778624452034324362035413739427444854597266606003627923669\n",
       "  57831768993344032657816302725040745056254556148878139499279348917001800710818\n",
       "  97098252782077376706106816615328949582935885760188519504856906240934843729036\n",
       "                                                                              ⋮\n",
       "  43617756509028935268065153506657982901346484655194542378962659172263286874567\n",
       "  24514787820482777097850819996878336822470282003929493899369788401994115725049\n",
       "    659019522081938406151199661909560306826152618189253797350605367892643425394\n",
       "  18524219368477146697571621538975472659957913040167655493576176702210093300777\n",
       "  19991360011325671962255145499577849636232056395638210789274105428863669504665\n",
       "  88153651723799965234202277829535942140168079285597169892962156845840816046269\n",
       "  28353066670915811088074600658821193560234484933732011087030237533483774640352\n",
       "  82601474594714678935549600454921758516384400812400780254743103044182810355919\n",
       "  63154568086128685299021557884272278280055552997286641122525772128245207399865\n",
       "   6834567513077213689931962590272659344149774329240293527553631454605963273832\n",
       "  69246004972092376302010948528430353858523542479157092028493846974578313098396\n",
       "  94103602193652048517365441015548431175193581166147658119497145544211926905650"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "N=40\n",
    "a=rand(1:2^BigInt(256),N)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "true"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "xtrue=rand(Bool,N)\n",
    "s=a'*xtrue\n",
    "setprecision(BigFloat,300)\n",
    "x=subsetsum(a,s)\n",
    "s==x'*a"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Integer Programming Feasibility Problem\n",
    "\n",
    "The integer programming feasibility problem is: given an `m` by `n` matrix `A`, and an `m`-dimensional vector `d`, determine if there is an `x` such that `Ax≤d`.  For a fixed `n`, H.W. Lenstra (1983) gave a polynomial time solution, which we illustrate in the `integerfeasibility` function. We do not go into the details of the technique, rather illustrate it with a simple example."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "true"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "A=[10 1 -9; 1 8 8]\n",
    "xtrue=[0; 2; 9]\n",
    "d=A*xtrue;\n",
    "x̂= integerfeasibility(A,d)\n",
    "xtrue==x̂"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In a more complex example, we set the values of the `20x30` matrix `A` to be uniform from `-10:10` and the `30`-dimensional vector `x` to have values taken uniformly from `0:10`. The `integerfeasibility` function finds the solution quickly."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "20×30 Array{Int64,2}:\n",
       "  -1    9   -2    4  -9   4   -4  -8  …    9    1   3  -3  -5    3   8    6\n",
       "   4   -5    5    8   6  -3   -3  -7      -9    2   1  -2  -6    7   3   -4\n",
       "  10   -2    6    9  -5  -5  -10   7      -1  -10   5  -4  -3   -1   1  -10\n",
       "  -7  -10   -5    9  -7  10    9  -3      10   -6  -5  -4   1    2   0    5\n",
       "   8   -7    7    7   0  -7    2  -2      -2    3   4   1  -2    4  -7   -1\n",
       "  -3   -4  -10    0   9   5    2  -7  …    6    7   1  -3   3    3  10   -5\n",
       "   4   -4    1  -10   6   2   -9  -3      -2    2  -8   7  -3   -5  -7    9\n",
       "  -4   -6   -2   -6  -1   3   10  -4      -9    4   4  -7  -6    9  -5    8\n",
       "   3   -6   -2    8   5   1   -4  -2       3    8   9  -1  -7   10   0   -8\n",
       "  -5    2    0    7   0   9    0   0      -3   -8  -2   7  -7   -7   6    4\n",
       "   5   -4   -3    8   4  -2    7  -7  …   -5    6  -7  -9   5   -5  -6   -6\n",
       " -10    3    6    9   0  -6    8   6      -8    8  -7   5  -7   -1  -3   -6\n",
       "  -3  -10    5    8   6   9   10   4      -5   -8   9   1  -8    8  10   -5\n",
       "   9    0   -4   -5   1  -3    8  10     -10   -4  -9  -3  -4    6   8   -3\n",
       "   2    5    6   -7  -6  -3    3   2      -3    3   6   8   5   -1  -9   -4\n",
       "  10    3  -10   -1  -9  -1   -8   5  …   -8    6   8   8  10   -6   9   -5\n",
       "   4  -10   -7   -7  -4  -9    4   1      -5   -4   6   1  -7    2   9   -3\n",
       "   5    8    6   -1   5   0    5   4      -3    0  -9   0  10    3   1  -10\n",
       "   6   -5    7   10  -7  -3    9   4      -3    1  -6   4  -2  -10  -7   -6\n",
       "   9  -10   -5    0   6   9   10   3     -10  -10   9  10  -9    9  -7    9"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "n=20\n",
    "m=30\n",
    "xtrue = rand(0:10,m)\n",
    "A = rand(-10:10,n,m)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "true"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "d=A*xtrue;\n",
    "x̂ = integerfeasibility(A,d)\n",
    "xtrue == x̂"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Why Julia?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Why write LLLplus in Julia? The C++ library [fplll](https://github.com/fplll/fplll) is a widely-used, fast, and robust lattice toolbox and is a reasonable choice for many. Even so, there are several advantages to native Julia code. For example, the same Julia code can accept matrices of integers, floats, [DoubleFloats](https://github.com/JuliaMath/DoubleFloats.jl), Float128s from [Quadmath](https://github.com/JuliaMath/Quadmath.jl), BigInts, BigFloats, and likely others. Hopefully LLLplus can also be used in other ways as part of a composable Julia  program.\n",
    "\n",
    "Want to learn more about lattice techniques? If you want a broad survey, try [\"The LLL Algorithm\"](https://www.springer.com/us/book/9783642022944).  The algorithm pseudocode in the [monograph by Bremner](https://www.amazon.com/Lattice-Basis-Reduction-Introduction-Applications/dp/1439807027) and the [survey paper](http://www.ant.uni-bremen.de/sixcms/media.php/102/10740/SPM_2011_Wuebben.pdf) were helpful to me. "
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Julia 1.1.0",
   "language": "julia",
   "name": "julia-1.1"
  },
  "language_info": {
   "file_extension": ".jl",
   "mimetype": "application/julia",
   "name": "julia",
   "version": "1.1.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}