(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     48216,       1272]
NotebookOptionsPosition[     43407,       1116]
NotebookOutlinePosition[     43958,       1140]
CellTagsIndexPosition[     43889,       1135]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["FusionAtlas package", "Subtitle",
 CellChangeTimes->{{3.4202909924834576`*^9, 3.4202909948668847`*^9}}],

Cell["\<\
Version 0, May 20 2008, Scott Morrison, Emily Peters, Noah Snyder\
\>", "Text",
 CellChangeTimes->{{3.420290997931291*^9, 3.420291017779832*^9}}],

Cell["Introduction", "Section"],

Cell[CellGroupData[{

Cell["Implementation", "Section"],

Cell[CellGroupData[{

Cell["Declaring subpackages", "Subsubsection",
 CellChangeTimes->{{3.420291037988891*^9, 3.420291045960354*^9}, {
  3.4202953396844254`*^9, 3.4202953402652607`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`Bigraphs`\>\"", ",", 
    RowBox[{"{", 
     RowBox[{
     "\"\<GradedGraph\>\"", ",", "\"\<GradedBigraph\>\"", ",", 
      "\"\<DisplayBigraph\>\"", ",", "\"\<GraphRank\>\"", ",", 
      "\"\<GraphEvenRank\>\"", ",", "\"\<GraphOddRank\>\"", ",", 
      "\"\<GraphDepth\>\"", ",", "\"\<RankAtDepth\>\"", ",", 
      "\"\<GraphAdjacencyMatrix\>\"", ",", "\"\<EvenPart\>\"", ",", 
      "\"\<GraphIndex\>\"", ",", "\"\<DimensionOfGenerator\>\"", ",", 
      "\"\<ReducedDimensionOfGenerator\>\"", ",", 
      "\"\<NumericDimensionOfGenerator\>\"", ",", "\"\<DimensionsByDepth\>\"",
       ",", "\"\<NumericDimensionsByDepth\>\"", ",", 
      "\"\<ReducedDimensionsByDepth\>\"", ",", "\"\<GlobalEvenDimension\>\"", 
      ",", "\"\<DimensionAtMostQ\>\"", ",", 
      "\"\<DimensionOfLowWeightSpace\>\"", ",", 
      "\"\<StableLowWeightSpaceDimensions\>\"", ",", 
      "\"\<AnnularTanglesSubgraphTest\>\"", ",", "\"\<AnnularTanglesTest\>\"",
       ",", "\"\<FindBigraphExtensions\>\"", ",", 
      "\"\<FindBigraphExtensionsUpToRank\>\"", ",", 
      "\"\<FindBigraphExtensionsUpToRankAndDepth\>\"", ",", 
      "\"\<CanonicalizeBigraph\>\"", ",", "\"\<AnBigraph\>\"", ",", 
      "\"\<DnBigraph\>\"", ",", "\"\<EnBigraph\>\"", ",", 
      "\"\<trivalentBigraph\>\"", ",", "\"\<haagerupFamilyBigraph\>\"", ",", 
      "\"\<dualHaagerupFamilyBigraph\>\"", ",", "\"\<HaagerupBigraph\>\"", 
      ",", "\"\<DualHaagerupBigraph\>\"", ",", "\"\<IzumiStarBigraph\>\"", 
      ",", "\"\<ExtendedDnBigraph\>\"", ",", "\"\<HexagonBigraph\>\"", ",", 
      "\"\<DualHexagonBigraph\>\"", ",", "\"\<HaagerupAsaedaBigraph\>\"", 
      ",", "\"\<DualHaagerupAsaedaBigraph\>\"", ",", 
      "\"\<ExtendedEnBigraph\>\"", ",", "\"\<S4S5Bigraph\>\"", ",", " ", 
      "\"\<A4A5Bigraph\>\"", ",", "\"\<A5A6Bigraph\>\"", ",", " ", 
      "\"\<GraphToString\>\"", ",", "\"\<GraphFromString\>\"", ",", 
      "\"\<ClearCachedGraphFromString\>\"", ",", 
      "\"\<CachedGraphFromString\>\"", ",", "\"\<DepthOfBranchPoint\>\"", 
      ",", "\"\<HasSimpleBranchQ\>\"", ",", "\"\<Truncate\>\"", ",", 
      "\"\<BranchFactor\>\""}], "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.420291117603371*^9, 3.420291118895229*^9}, {
   3.4202913849377794`*^9, 3.420291388132373*^9}, {3.4202920292742896`*^9, 
   3.4202920418123183`*^9}, {3.4202927560893984`*^9, 
   3.4202927623383837`*^9}, {3.42045421068404*^9, 3.4204542737647457`*^9}, {
   3.422309142772846*^9, 3.422309143143379*^9}, {3.4223200548937054`*^9, 
   3.4223200573272047`*^9}, {3.422375248668456*^9, 3.422375252814418*^9}, {
   3.4223758021242867`*^9, 3.422375814021394*^9}, {3.4325166524374104`*^9, 
   3.432516656923388*^9}, {3.432761815099677*^9, 3.4327618166418943`*^9}, {
   3.432764768646669*^9, 3.43276477939212*^9}, {3.4330380029733086`*^9, 
   3.4330380035441294`*^9}, {3.433967209088317*^9, 3.4339672122328386`*^9}, {
   3.434302768103797*^9, 3.4343027688749056`*^9}, {3.4343096082895007`*^9, 
   3.4343096097015314`*^9}, {3.4343930690519533`*^9, 3.434393073788765*^9}, {
   3.4345639820139456`*^9, 3.4345639841069555`*^9}, {3.434744344966757*^9, 
   3.4347443466792192`*^9}, {3.4354440045850115`*^9, 
   3.4354440057466817`*^9}, {3.4457949428591948`*^9, 3.445794945712213*^9}, {
   3.4459128520720263`*^9, 3.445912852970158*^9}, 3.446846035638986*^9, {
   3.4475529354378366`*^9, 3.447552946262328*^9}, {3.452449185709619*^9, 
   3.452449190202622*^9}, {3.46284809518544*^9, 3.4628481036430807`*^9}, {
   3.463943116119879*^9, 3.463943117779564*^9}, {3.466796765560052*^9, 
   3.466796768949822*^9}, {3.469391560478197*^9, 3.469391563293811*^9}, {
   3.469391731593257*^9, 3.469391736056587*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`Java`\>\"", ",", 
    RowBox[{"{", 
     RowBox[{
     "\"\<FusionAtlasJavaDirectory\>\"", ",", 
      "\"\<FusionAtlasScalaDirectory\>\"", ",", "\"\<RestartJava\>\"", ",", 
      "\"\<GraphToScalaObject\>\"", ",", "\"\<AsScalaSet\>\"", ",", 
      "\"\<AsScalaList\>\"", ",", "\"\<UnpackJavaLists\>\"", ",", 
      "\"\<UnpackJavaObjects;\>\""}], "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.494961336398065*^9, 3.49496135265196*^9}, {
  3.494981177301717*^9, 3.494981179238475*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`JavaOdometer`\>\"", ",", 
    RowBox[{"{", 
     RowBox[{
     "\"\<FindBigraphExtensionsUpToRankJava\>\"", ",", 
      "\"\<FindBigraphExtensionsUpToRankAndDepthJava\>\"", ",", 
      "\"\<FindBigraphExtensionsUpToRankAndDepthJavaIterator\>\"", ",", 
      "\"\<FindBigraphExtensionsMatching\>\""}], "}"}]}], "]"}], 
  ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4327784707593393`*^9, 3.432778495064288*^9}, {
  3.434392594850085*^9, 3.4343925958415103`*^9}, {3.4347609467483344`*^9, 
  3.434760949822755*^9}, {3.4420723723011217`*^9, 3.4420723735068293`*^9}}],

Cell[BoxData[
 RowBox[{"DeclarePackage", "[", 
  RowBox[{"\"\<FusionAtlas`FindAllBigraphs`\>\"", ",", 
   RowBox[{"{", 
    RowBox[{
    "\"\<AlgebraicIntegerTest\>\"", ",", "\"\<PositivityTest\>\"", ",", 
     "\"\<GHJTest\>\"", ",", "\"\<BigraphConditions\>\"", ",", 
     "\"\<AllBigraphs\>\"", ",", "\"\<FindAllBigraphs\>\"", ",", 
     "\"\<RequestAllBigraphs\>\"", ",", "\"\<KnownObstructions\>\"", ",", 
     "\"\<FindAllObstructions\>\"", ",", "\"\<RequestMoreBigraphs\>\""}], 
    "}"}]}], "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.43242799348228*^9, 3.432428026930376*^9}, {
   3.4329521780855713`*^9, 3.432952183433261*^9}, {3.4337973756456947`*^9, 
   3.4337973767372637`*^9}, {3.434046998429755*^9, 3.434047003957704*^9}, {
   3.434047069071333*^9, 3.4340470698724847`*^9}, {3.4340485038744783`*^9, 
   3.434048504685645*^9}, {3.4340489688130274`*^9, 3.4340489709460945`*^9}, {
   3.4347614220818305`*^9, 3.4347614231633854`*^9}, 3.435618442494402*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`PackageData`\>\"", ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<ValuesAsString\>\"", ",", "\"\<PackageData\>\"", ",", 
      "\"\<MatchingValues\>\"", ",", "\"\<DefiniteValues\>\"", ",", 
      "\"\<PackageEverything\>\"", ",", "\"\<PackageCachedRootReduce\>\"", 
      ",", "\"\<PackageDimensions\>\"", ",", 
      "\"\<PackageBigraphExtensions\>\"", ",", 
      "\"\<PackageBifusionAlgebras\>\"", ",", 
      "\"\<PackageLowWeightSpaces\>\"", ",", 
      "\"\<PackageAlgebraicIntegerTest\>\"", ",", 
      "\"\<PackageObstructions\>\"", ",", "\"\<PackageKnapsackTest\>\"", ",", 
      "\"\<PackageFindGraphPartners\>\"", ",", 
      "\"\<PackageLowestWeightEigenspaces\>\"", ",", 
      "\"\<PackageIdempotents\>\"", ",", "\"\<PackagePairExtensions\>\""}], 
     "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4211987025850735`*^9, 3.421198732668331*^9}, {
   3.4217687607767344`*^9, 3.4217687621286783`*^9}, {3.422313404661141*^9, 
   3.422313406173315*^9}, {3.4223144034873824`*^9, 3.4223144048393264`*^9}, {
   3.4223172035937376`*^9, 3.4223172095723343`*^9}, {3.4316570849647026`*^9, 
   3.431657093346755*^9}, {3.432252673775091*^9, 3.432252674696416*^9}, {
   3.432252868835574*^9, 3.432252869626712*^9}, {3.4324314057288494`*^9, 
   3.432431406670203*^9}, 3.432433620623712*^9, {3.432532945756606*^9, 
   3.432532946517701*^9}, {3.4329526458581953`*^9, 3.432952660569349*^9}, {
   3.4468455105303917`*^9, 3.4468455211718063`*^9}, 3.446846676394143*^9, {
   3.44700322268544*^9, 3.447003223523093*^9}, {3.447532126164963*^9, 
   3.4475321368589067`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`TensorSolver`\>\"", ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<FindFusionAlgebras\>\"", ",", "\"\<FindBifusionAlgebras\>\"", ",", 
      "\"\<FusionAlgebra\>\"", ",", "\"\<BifusionAlgebra\>\"", ",", 
      "\"\<TimeConstrainedFindBifusionAlgebras\>\"", ",", 
      "\"\<TimeConstrainedFindFusionAlgebras\>\"", ",", 
      "\"\<LazyFindFusionAlgebras\>\"", ",", 
      "\"\<FindBifusionAlgebras2\>\""}], "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.421268743427096*^9, 3.4212687695746946`*^9}, {
   3.422317812329056*^9, 3.4223178208813534`*^9}, {3.432863444351768*^9, 
   3.4328634617668095`*^9}, {3.4338737890491314`*^9, 3.4338737898803263`*^9}, 
   3.433874479241579*^9, {3.4339653953502893`*^9, 3.4339653965520177`*^9}, {
   3.4339662905074625`*^9, 3.4339662914488163`*^9}, {3.4342173889185743`*^9, 
   3.434217390651066*^9}, {3.4342196931919565`*^9, 3.434219696867242*^9}, {
   3.434404383180872*^9, 3.4344043839519806`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{
   "\"\<FusionAtlas`FindGraphPartners`\>\"", ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<FindGraphPartners\>\"", ",", 
      "\"\<FindUnobstructedGraphPartners\>\"", ",", 
      "\"\<LazyFindGraphPartners\>\"", ",", "\"\<FindAllGraphPartners\>\"", 
      ",", "\"\<FindGraphPartnersTimeConstraint\>\""}], "}"}]}], "]"}], 
  ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.421268743427096*^9, 3.4212687695746946`*^9}, {
  3.422317812329056*^9, 3.4223178208813534`*^9}, {3.432863444351768*^9, 
  3.4328634617668095`*^9}, {3.4328641733099585`*^9, 3.4328641818221984`*^9}, {
  3.4339726736860223`*^9, 3.4339726744371023`*^9}, {3.4340472883065777`*^9, 
  3.4340472891377726`*^9}, {3.4340473475918255`*^9, 3.4340473490038557`*^9}, {
  3.4354298293319983`*^9, 3.4354298302633376`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{
   "\"\<FusionAtlas`CyclotomicTest`\>\"", ",", "\[IndentingNewLine]", 
    RowBox[{"{", "\"\<CyclotomicTest\>\"", "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.421268743427096*^9, 3.4212687695746946`*^9}, {
  3.422317812329056*^9, 3.4223178208813534`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`KnapsackTest`\>\"", ",", "\[IndentingNewLine]", 
    RowBox[{"{", "\"\<KnapsackTest\>\"", "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4326058911572943`*^9, 3.4326058967353153`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{
   "\"\<FusionAtlas`QuadraticTanglesTest`\>\"", ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"\"\<ParityTest\>\"", ",", "\"\<HaagerupTypeBranchTest\>\""}], 
     "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.421268743427096*^9, 3.4212687695746946`*^9}, {
   3.422317812329056*^9, 3.4223178208813534`*^9}, 3.4325190046011143`*^9, {
   3.442870008289198*^9, 3.442870013943906*^9}, 3.44684603131995*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`dTest`\>\"", ",", 
    RowBox[{"{", "\"\<DTest\>\"", "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4336943604918747`*^9, 3.433694387078576*^9}, {
  3.433695220719633*^9, 3.43369522089264*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`Debugging`\>\"", ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<DebugEcho\>\"", ",", "\"\<DebugPrint\>\"", ",", 
      "\"\<DebugEvaluate\>\"", ",", "\"\<SetDebugMode\>\"", ",", 
      "\"\<SetDebugQueueInterval\>\"", ",", "\"\<SetDebugFile\>\"", ",", 
      "\"\<DebugPrintHeld\>\"", ",", "\"\<ToStringHeld\>\""}], "}"}]}], "]"}],
   ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.421268743427096*^9, 3.4212687695746946`*^9}, {
   3.433201472620451*^9, 3.433201475434498*^9}, 3.433201618620389*^9, {
   3.4343084187289968`*^9, 3.4343084198706384`*^9}, {3.4458128612783012`*^9, 
   3.4458128623303967`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`RemoteInterface`\>\"", ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<RecordObstructions\>\"", ",", "\"\<RequestObstructions\>\"", ",", 
      "\"\<LookupObstructions\>\"", ",", "\"\<LookupCanonicalForms\>\"", ",", 
      "\"\<ProcessObstructionsQueue\>\"", ",", 
      "\"\<WorkOnObstructionsQueue\>\"", ",", 
      "\"\<ClearObstructionsQueue\>\"", ",", 
      "\"\<ObstructionsQueueLength\>\"", ",", 
      "\"\<EstimateQueueCompletionTime\>\"", ",", "\"\<RetrieveGraphs\>\"", 
      ",", "\"\<LookupGraphs\>\"", ",", "\"\<TakeGraphs\>\"", ",", 
      "\"\<RecordCanonicalForms\>\"", ",", "\"\<RecordFindGraphPartners\>\"", 
      ",", "\"\<GraphsNeedingPartners\>\"", ",", 
      "\"\<ReapplyTestToDatabase\>\"", ",", 
      "\"\<RecordRequestAllBigraphs\>\"", ",", "\"\<BigraphRequests\>\"", 
      ",", "\"\<MaximalBigraphRequests\>\"", ",", 
      "\"\<MaximalIncompleteBigraphRequests\>\"", ",", 
      "\"\<DisplayBigraphRequests\>\"", ",", "\"\<QueueExpression\>\""}], 
     "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.434310388541448*^9, 3.4343104010995054`*^9}, {
   3.4343113325288353`*^9, 3.4343113388879795`*^9}, {3.434315474344474*^9, 
   3.434315477458952*^9}, {3.4343993407802563`*^9, 3.4343993415012927`*^9}, {
   3.43440255525244*^9, 3.4344025565442977`*^9}, {3.4345616588233614`*^9, 
   3.434561659734672*^9}, {3.434562994834451*^9, 3.4345629958258767`*^9}, {
   3.434564956455123*^9, 3.4345649584079313`*^9}, {3.4345718103405266`*^9, 
   3.434571811682456*^9}, {3.4345912761610017`*^9, 3.4345912778133774`*^9}, {
   3.4346550182436466`*^9, 3.434655025734418*^9}, {3.434655213664648*^9, 
   3.4346552144658003`*^9}, 3.4346821631975794`*^9, {3.435433552625832*^9, 
   3.435433553637286*^9}, {3.4354393239646096`*^9, 3.435439324575488*^9}, {
   3.435508322319064*^9, 3.4355083250129375`*^9}, {3.4356134448281107`*^9, 
   3.4356134480527472`*^9}, {3.435613768092942*^9, 3.4356137805708847`*^9}, {
   3.43561567221093*^9, 3.4356156767474527`*^9}, {3.4356165723753023`*^9, 
   3.435616573797347*^9}, {3.4358825131991296`*^9, 3.4358825140002813`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`GraphPairs`\>\"", ",", 
    RowBox[{"{", 
     RowBox[{
     "\"\<BigraphWithDuals\>\"", ",", "\"\<DualData\>\"", ",", 
      "\"\<OcneanuDisplay\>\"", ",", "\"\<HaagerupWithDuals\>\"", ",", 
      "\"\<DualHaagerupWithDuals\>\"", ",", "\"\<HaagerupAsaedaWithDuals\>\"",
       ",", "\"\<DualHaagerupAsaedaWithDuals\>\"", ",", 
      "\"\<TwistBigraph\>\"", ",", "\"\<FullOcneanuTest\>\"", ",", 
      "\"\<PartialOcneanuTest\>\"", ",", "\"\<DualDataList\>\"", ",", 
      "\"\<CompatibleDualDataList\>\"", ",", "\"\<CheckDualDimensions\>\"", 
      ",", "\"\<GraphPairsIsomorphicQ\>\"", ",", 
      "\"\<ClearGraphPairsCachedValues\>\"", ",", 
      "\"\<PermuteToReduceCrossings\>\""}], "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4434668365282307`*^9, 3.443466892940999*^9}, {
  3.443467181474247*^9, 3.443467197704023*^9}, {3.4434675965538397`*^9, 
  3.443467607917555*^9}, {3.447606347932467*^9, 3.447606352737217*^9}, {
  3.497328027375166*^9, 3.4973280281671867`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`PairOdometer`\>\"", ",", 
    RowBox[{"{", 
     RowBox[{
     "\"\<FindBigraphPairExtensionsUpToDepth\>\"", ",", 
      "\"\<FamiliesOfBigraphPairs\>\"", ",", 
      "\"\<DisplayFamiliesOfBigraphPairs\>\"", ",", 
      "\"\<ChooseIsomorphismRepresentatives\>\""}], "}"}]}], "]"}], 
  ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4327784707593393`*^9, 3.432778495064288*^9}, {
  3.434392594850085*^9, 3.4343925958415103`*^9}, {3.4347609467483344`*^9, 
  3.434760949822755*^9}, {3.4420723723011217`*^9, 3.4420723735068293`*^9}, {
  3.447362492494246*^9, 3.447362516146487*^9}, {3.48447538940093*^9, 
  3.4844753906293287`*^9}}],

Cell[BoxData[
 RowBox[{"DeclarePackage", "[", 
  RowBox[{"\"\<FusionAtlas`GraphPlanarAlgebra`\>\"", ",", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{
    "\"\<Vertex\>\"", ",", "\[IndentingNewLine]", "\"\<GPAElement\>\"", ",", 
     "\"\<GPAMatrix\>\"", ",", "\n", "\"\<NumberFieldGenerator\>\"", ",", 
     "\"\<EnlargeNumberField\>\"", ",", "\n", "\"\<LopsidedDimension\>\"", 
     ",", "\"\<EmptyGPAElement\>\"", ",", "\"\<GPACoefficients\>\"", ",", 
     "\"\<GPAMultiply\>\"", ",", "\n", "\"\<RotateOneClick\>\"", ",", 
     "\"\<RotateTwoClicks\>\"", ",", "\"\<CapTopLeft\>\"", ",", 
     "\"\<CapAt\>\"", ",", "\"\<AddStrandOnRight\>\"", ",", " ", 
     "\"\<AddCupOnTop\>\"", ",", "\"\<JonesWenzlIdempotent\>\"", ",", "\n", 
     "\"\<LowestWeightEigenvectorConditions\>\"", ",", "\n", 
     "\"\<LowestWeightEigenspace\>\""}], "\[IndentingNewLine]", "}"}]}], 
  "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.446845641507465*^9, 3.4468456434488688`*^9}, {
  3.446845803228867*^9, 3.446845833163484*^9}, {3.4468463123606853`*^9, 
  3.4468463178259287`*^9}, {3.446849955796545*^9, 3.446849961808324*^9}, {
  3.446850089214408*^9, 3.4468500939946737`*^9}, {3.4468505694305267`*^9, 
  3.446850618370099*^9}, {3.446856990391223*^9, 3.446856991403557*^9}, {
  3.446912818498979*^9, 3.44691282042977*^9}, {3.483034922782889*^9, 
  3.4830349355462112`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`RelativeDimensions`\>\"", ",", 
    RowBox[{"{", 
     RowBox[{
     "\"\<RelativeDimensionsByDepthEquations\>\"", ",", 
      "\"\<SolveRelativeDimensionsByDepthEquations\>\"", ",", 
      "\"\<RelativeBranchFactor\>\"", ",", 
      "\"\<RelativeDimensionsByDepth\>\"", ",", 
      "\"\<RelativeDimensionsConditions\>\"", ",", 
      "\"\<PositiveOmegaInequality\>\""}], "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.469391597841446*^9, 3.469391698539899*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`TriplePointTest`\>\"", ",", 
    RowBox[{"{", 
     RowBox[{
     "\"\<TriplePointObstruction\>\"", ",", 
      "\"\<EvenTriplePointObstruction\>\"", ",", 
      "\"\<OddTriplePointObstruction\>\""}], "}"}]}], "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.469391597841446*^9, 3.469391698539899*^9}, {
  3.4783066414764223`*^9, 3.4783066671764927`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"DeclarePackage", "[", 
   RowBox[{"\"\<FusionAtlas`Utilities`IntersectSubspaces`\>\"", ",", 
    RowBox[{"{", "\"\<IntersectSubspaces\>\"", "}"}]}], "]"}], ";"}]], "Input",\

 InitializationCell->True,
 CellChangeTimes->{{3.469391597841446*^9, 3.469391698539899*^9}, {
  3.4783066414764223`*^9, 3.4783066671764927`*^9}, {3.483185632800359*^9, 
  3.483185650453915*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Print loading statement", "Subsubsection",
 CellChangeTimes->{{3.4202910743011055`*^9, 3.4202910795386367`*^9}}],

Cell[BoxData[
 RowBox[{"Print", "[", "\[IndentingNewLine]", 
  RowBox[{
  "\"\<Loading FusionAtlas` version 0\\n\>\"", ",", "\[IndentingNewLine]", 
   "\"\<Read more at http://tqft.net/wiki/Atlas_of_subfactors\>\""}], 
  "\[IndentingNewLine]", "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.420291094860669*^9, 3.4202911096118803`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["The core package", "Subsubsection",
 CellChangeTimes->{{3.420291084395621*^9, 3.4202910873398542`*^9}, {
  3.420291344749992*^9, 3.420291345451*^9}, {3.420291407660453*^9, 
  3.420291409573203*^9}}],

Cell["\<\
This doesn' t actually contain much; just some common functions. Everything \
interesting is in a subpackage.\
\>", "Text",
 CellChangeTimes->{{3.420291411676227*^9, 3.420291435810931*^9}}],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"$VersionNumber", "<", "6"}], ",", "\[IndentingNewLine]", 
   RowBox[{"BeginPackage", "[", 
    RowBox[{"\"\<FusionAtlas`\>\"", ",", 
     RowBox[{"{", "\"\<LinearAlgebra`MatrixManipulation`\>\"", "}"}]}], "]"}],
    ",", "\[IndentingNewLine]", 
   RowBox[{"BeginPackage", "[", "\"\<FusionAtlas`\>\"", "]"}]}], 
  "\[IndentingNewLine]", "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4202911560186095`*^9, 3.4202911587124834`*^9}, {
  3.420295014256483*^9, 3.4202950158788157`*^9}, {3.4202954024747133`*^9, 
  3.42029543357944*^9}, {3.4202954793853054`*^9, 3.4202955565763006`*^9}, {
  3.4202956106340322`*^9, 3.420295640847477*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"FusionAtlasDirectory", "::", "usage"}], "=", 
   "\"\<FusionAtlasDirectory[] should hopefully return the location the \
FusionAtlas` package was loaded from.\>\""}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"FusionAtlasDataDirectory", "::", "usage"}], "=", 
   "\"\<FusionAtlasDataDirectory[] specifies were the FusionAtlas` package \
should look for, and save, precomputed data.\>\""}], ";"}]}], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4211987935558834`*^9, 3.4211988124230127`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"FusionAtlas", "::", "loading"}], "=", 
   "\"\<Loading precomputed data in `1`.\>\""}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.421200368931162*^9, 3.4212003847038417`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"cachedRootReduce", "::", "usage"}], "=", "\"\<\>\""}], 
  ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.420454668872883*^9, 3.4204546717370014`*^9}, {
  3.4223130452443247`*^9, 3.4223130471670895`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"UnsortedUnion", "::", "usage"}], "=", 
   "\"\<UnsortedUnion[list] a list of all the unique elements in list, in the \
order that they first appear.\>\""}], ";"}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"FindFirst", "::", "usage"}], "=", "\"\<\>\""}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4202930097040787`*^9, 3.4202930126683407`*^9}, {
  3.4202947327316704`*^9, 3.420294733723096*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"RealMax", "::", "usage"}], "=", "\"\<\>\""}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"RealMaxPosition", "::", "usage"}], "=", "\"\<\>\""}], 
  ";"}]}], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4202947266929874`*^9, 3.420294741023594*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"NaturalQ", "::", "usage"}], "=", 
   "\"\<NaturalNumberQ[n] tests if n is a non-negative integer.\>\""}], 
  ";"}]], "Input",
 InitializationCell->True],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"qInteger", "::", "usage"}], "=", 
   "\"\<qInteger[n][q] computes the quantum integer n with the variable q; in \
our conventions qInteger[2][q]==q+\!\(\*SuperscriptBox[\(q\), \
\(-1\)]\)\>\""}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"{", 
   RowBox[{"qFactorial", ",", "qBinomial"}], "}"}], ";"}]}], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4202912687707396`*^9, 3.420291287768056*^9}}],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"$VersionNumber", "<", "7"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"SplitBy", "[", 
     RowBox[{"x_", ",", "f_"}], "]"}], ":=", 
    RowBox[{"Split", "[", 
     RowBox[{"x", ",", 
      RowBox[{
       RowBox[{
        RowBox[{"f", "[", "#1", "]"}], "===", 
        RowBox[{"f", "[", "#2", "]"}]}], "&"}]}], "]"}]}]}], 
  "\[IndentingNewLine]", "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.443479169922008*^9, 3.4434791770836267`*^9}}],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"$VersionNumber", "<", "6"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"UnitVector", "::", "usage"}], "=", 
     "\"\<UnitVector[n,i] returns the i-th n-dimensional unit vector, if i is \
an integer between 1 and n, and the n-dimensional zero vector \
otherwise.\>\""}], ";"}]}], "\[IndentingNewLine]", "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4094369769829683`*^9, 3.4094369851747475`*^9}}],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"$VersionNumber", "\[GreaterEqual]", "6"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"BlockMatrix", "=", "ArrayFlatten"}], ";", "\[IndentingNewLine]", 
    
    RowBox[{
     RowBox[{"ZeroMatrix", "[", "n_", "]"}], ":=", 
     RowBox[{"ZeroMatrix", "[", 
      RowBox[{"n", ",", "n"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ZeroMatrix", "[", 
      RowBox[{"n_", ",", "m_"}], "]"}], ":=", 
     RowBox[{"ConstantArray", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", 
        RowBox[{"n", ",", "m"}], "}"}]}], "]"}]}]}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4094369769829683`*^9, 3.4094369851747475`*^9}, {
  3.42029346825344*^9, 3.4202934790289345`*^9}, {3.4202935339979763`*^9, 
  3.4202935487091293`*^9}, {3.4202943339883056`*^9, 3.420294351593621*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ZeroVector", "::", "usage"}], "=", 
   "\"\<ZeroVector[n] returns the n-dimensional 0 vector.\>\""}], 
  ";"}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"UnitVectorQ", "::", "usage"}], "=", 
   "\"\<UnitVectorQ[v] tests if v is a unit coordinate vector.\>\""}], 
  ";"}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ZeroVectorQ", "::", "usage"}], "=", 
   "\"\<ZeroVectorQ[v] tests if v is the zero vector.\>\""}], ";"}]], "Input",\

 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"ZeroMatrixQ", "::", "usage"}], "=", "\"\<\>\""}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4202941778337665`*^9, 3.4202941799367905`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"qFromd", "::", "Usage"}], "=", 
   "\"\<Computes q given d=q+q^(-1)\>\""}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.466791516339015*^9, 3.4667915397614098`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Begin", "[", "\"\<`Private`\>\"", "]"}], ";"}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"FusionAtlasDirectory", "[", "]"}], ":=", 
  RowBox[{
   RowBox[{"FusionAtlasDirectory", "[", "]"}], "=", 
   RowBox[{"StringDrop", "[", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"File", "/.", 
       RowBox[{"Flatten", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"FileInformation", "[", 
           RowBox[{"ToFileName", "[", 
            RowBox[{"#", ",", "\"\<FusionAtlas\>\""}], "]"}], "]"}], "&"}], "/@", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Reverse", "[", "$Path", "]"}], "/.", 
           RowBox[{"\"\<.\>\"", "\[Rule]", 
            RowBox[{"Directory", "[", "]"}]}]}], ")"}]}], "]"}]}], ")"}], ",", 
     RowBox[{"-", "12"}]}], "]"}]}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.421198815136915*^9, 3.4211988188322287`*^9}, {
   3.4211993678116226`*^9, 3.421199367971853*^9}, 3.494103885026681*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"might", " ", "be", " ", "dangerous", " ", "if", " ", 
    RowBox[{"FusionAtlasDirectory", "[", "]"}], " ", "is", " ", "somehow", 
    " ", 
    RowBox[{"incorrect", "!"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{"!", 
     RowBox[{"MemberQ", "[", 
      RowBox[{"$Path", ",", 
       RowBox[{"FusionAtlasDirectory", "[", "]"}]}], "]"}]}], ",", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"$Path", ",", 
      RowBox[{"FusionAtlasDirectory", "[", "]"}]}], "]"}]}], "]"}]}]], "Input",\

 InitializationCell->True,
 CellChangeTimes->{{3.4211988253315744`*^9, 3.421198831550517*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"StringTake", "[", 
      RowBox[{
       RowBox[{"FusionAtlasDirectory", "[", "]"}], ",", 
       RowBox[{"-", "7"}]}], "]"}], "==", "\"\<package\>\""}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"FusionAtlasDataDirectory", "[", "]"}], ":=", 
      RowBox[{
       RowBox[{"StringDrop", "[", 
        RowBox[{
         RowBox[{"FusionAtlasDirectory", "[", "]"}], ",", 
         RowBox[{"-", "7"}]}], "]"}], "<>", "\"\<data\>\""}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{"\"\<Found precomputed data in \>\"", ",", 
       RowBox[{"FusionAtlasDataDirectory", "[", "]"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"!", 
        RowBox[{"MemberQ", "[", 
         RowBox[{"$Path", ",", 
          RowBox[{"FusionAtlasDataDirectory", "[", "]"}]}], "]"}]}], ",", 
       RowBox[{"AppendTo", "[", 
        RowBox[{"$Path", ",", 
         RowBox[{"FusionAtlasDataDirectory", "[", "]"}]}], "]"}]}], "]"}]}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{
    "Print", "[", 
     "\"\<Remember to set FusionAtlasDataDirectory[] to the appropriate path, \
if you've downloaded precomputed data.\>\"", "]"}]}], "\[IndentingNewLine]", 
   "]"}], ";"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4211988348152113`*^9, 3.4211988575078416`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"UnsortedUnion", "[", "x_", "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", "f", "}"}], ",", 
    RowBox[{
     RowBox[{
      RowBox[{"f", "[", "y_", "]"}], ":=", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"f", "[", "y", "]"}], "=", 
         RowBox[{"Sequence", "[", "]"}]}], ";", "y"}], ")"}]}], ";", 
     RowBox[{"f", "/@", "x"}]}]}], "]"}]}]], "Input",
 InitializationCell->True,
 CellTags->"Union"],

Cell[BoxData[
 RowBox[{
  RowBox[{"cachedRootReduce", "[", "x_", "]"}], ":=", 
  RowBox[{
   RowBox[{"cachedRootReduce", "[", "x", "]"}], "=", 
   RowBox[{"RootReduce", "[", "x", "]"}]}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4204574886323776`*^9, 3.4204574984364753`*^9}, {
  3.422313055348854*^9, 3.4223130594547586`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"NaturalQ", "[", "n_", "]"}], ":=", 
  RowBox[{
   RowBox[{"NonNegative", "[", "n", "]"}], "&&", 
   RowBox[{"IntegerQ", "[", "n", "]"}]}]}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"FindFirst", "[", 
   RowBox[{"a_", ",", "b_"}], "]"}], ":=", 
  RowBox[{"With", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"p", "=", 
      RowBox[{"Position", "[", 
       RowBox[{"a", ",", "b", ",", 
        RowBox[{"{", "1", "}"}], ",", "1", ",", 
        RowBox[{"Heads", "\[Rule]", "False"}]}], "]"}]}], "}"}], ",", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"Length", "[", "p", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{"-", "1"}], ",", 
      RowBox[{"p", "\[LeftDoubleBracket]", 
       RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}]}], "]"}]}], 
   "]"}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4202928575352707`*^9, 3.420292858706955*^9}, {
   3.432527304477088*^9, 3.432527341153398*^9}, 3.466790210028186*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"RealMax", "[", "a_", "]"}], " ", ":=", " ", 
  RowBox[{"Max", "[", 
   RowBox[{"Cases", "[", 
    RowBox[{"a", ",", 
     RowBox[{"x_", "/;", 
      RowBox[{
       RowBox[{"Im", "[", "x", "]"}], "\[Equal]", "0"}]}]}], "]"}], 
   "]"}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.420292889351019*^9, 3.4202928908431644`*^9}, {
  3.420293030824448*^9, 3.420293043372491*^9}, {3.4202931758529887`*^9, 
  3.4202931829331694`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"RealMaxPosition", "[", "a_", "]"}], " ", ":=", " ", 
  RowBox[{"FindFirst", "[", 
   RowBox[{"a", ",", 
    RowBox[{"RealMax", "[", "a", "]"}]}], "]"}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4202928605896626`*^9, 3.42029286538656*^9}, {
  3.4202947072450223`*^9, 3.420294708216419*^9}}],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"$VersionNumber", "\[GreaterEqual]", "6."}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Unprotect", "[", "IdentityMatrix", "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"IdentityMatrix", "[", "0", "]"}], "=", 
     RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Protect", "[", "IdentityMatrix", "]"}], ";"}]}], 
  "\[IndentingNewLine]", "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4094379343596077`*^9, 3.4094379584642687`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"ZeroVector", "[", "n_", "]"}], ":=", 
  RowBox[{"Table", "[", 
   RowBox[{"0", ",", 
    RowBox[{"{", "n", "}"}]}], "]"}]}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"$VersionNumber", "<", "6"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"UnitVector", "[", 
        RowBox[{"n_", ",", "i_Integer"}], "]"}], "/;", 
       RowBox[{"(", 
        RowBox[{"1", "\[LessEqual]", "i", "\[LessEqual]", "n"}], ")"}]}], ":=", 
      RowBox[{"Module", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"z", "=", 
          RowBox[{"Table", "[", 
           RowBox[{"0", ",", 
            RowBox[{"{", "n", "}"}]}], "]"}]}], "}"}], ",", 
        RowBox[{
         RowBox[{
          RowBox[{"z", "\[LeftDoubleBracket]", "i", "\[RightDoubleBracket]"}],
           "=", "1"}], ";", "z"}]}], "]"}]}], ")"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"UnitVector", "[", 
       RowBox[{"n_", ",", "i_Integer"}], "]"}], ":=", 
      RowBox[{"Table", "[", 
       RowBox[{"0", ",", 
        RowBox[{"{", "n", "}"}]}], "]"}]}], ")"}]}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Unprotect", "[", "UnitVector", "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"UnitVector", "[", 
        RowBox[{"n_", ",", "i_Integer"}], "]"}], "/;", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"i", "<", "1"}], "\[Or]", 
         RowBox[{"i", ">", "n"}]}], ")"}]}], ":=", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"Message", "[", 
         RowBox[{
          RowBox[{"UnitVector", "::", "nokun"}], ",", "n", ",", "i"}], "]"}], 
        ";", 
        RowBox[{"ZeroVector", "[", "n", "]"}]}], ")"}]}], ")"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Protect", "[", "UnitVector", "]"}], ";"}]}], 
  "\[IndentingNewLine]", "]"}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.409437049016547*^9, 3.409437188376938*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"UnitVectorQ", "[", 
   RowBox[{"v_", "?", "VectorQ"}], "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Complement", "[", 
     RowBox[{"v", ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "1"}], "}"}]}], "]"}], "\[Equal]", 
    RowBox[{"{", "}"}]}], "\[And]", 
   RowBox[{
    RowBox[{"Count", "[", 
     RowBox[{"v", ",", "1"}], "]"}], "\[Equal]", "1"}]}]}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"ZeroVectorQ", "[", 
   RowBox[{"v_", "?", "VectorQ"}], "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"Union", "[", "v", "]"}], "===", 
    RowBox[{"{", "0", "}"}]}], "\[Or]", 
   RowBox[{"v", "\[Equal]", 
    RowBox[{"{", "}"}]}]}]}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"ZeroMatrixQ", "[", "m_", "]"}], " ", ":=", " ", 
  RowBox[{"(", 
   RowBox[{"m", "\[Equal]", 
    RowBox[{"ZeroMatrix", "@@", 
     RowBox[{"Dimensions", "[", "m", "]"}]}]}], ")"}]}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"This", " ", "is", " ", "horribly", " ", 
    RowBox[{"broken", ":", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"e", ".", "g", ".", " ", 
        RowBox[{"MinimalPolynomial", "[", 
         RowBox[{
          SqrtBox[
           RowBox[{"Root", "[", 
            RowBox[{
             RowBox[{
              RowBox[{
               RowBox[{"-", "2"}], "+", 
               RowBox[{"12", " ", "#1"}], "-", 
               RowBox[{"7", " ", 
                SuperscriptBox["#1", "2"]}], "+", 
               SuperscriptBox["#1", "3"]}], "&"}], ",", "3"}], "]"}]], ",", 
          "x"}], "]"}]}], " ", "gives"}], " ", "-", 
      RowBox[{"x", "^", 
       RowBox[{
        RowBox[{"3", "!!"}], "!"}]}]}]}]}], " ", "*)"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"Unfortunately", " ", "I", " ", 
    RowBox[{"can", "'"}], "t", " ", "remember", " ", "why", " ", "we", " ", 
    "wrote", " ", 
    RowBox[{"this", ":", " ", 
     RowBox[{
     "it", " ", "must", " ", "have", " ", "been", " ", "to", " ", "speed", 
      " ", "something", " ", 
      RowBox[{"up", ".", " ", "Hopefully"}], " ", 
      RowBox[{"it", "'"}], "s", " ", "just", " ", "not", " ", "a", " ", 
      "problem", " ", "anymore", " ", "to", " ", "use", " ", "the", " ", 
      "default", " ", "implementation", " ", "in", " ", "Mathematica"}]}]}], 
   " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Unprotect", "[", "MinimalPolynomial", "]"}], ";", "\n", 
    RowBox[{
     RowBox[{
      RowBox[{"MinimalPolynomial", "[", 
       RowBox[{"\[Mu]_", ",", "z_"}], "]"}], "/;", 
      RowBox[{
       RowBox[{"Length", "[", 
        RowBox[{"Union", "[", 
         RowBox[{"Cases", "[", 
          RowBox[{"\[Mu]", ",", "_Root", ",", "\[Infinity]"}], "]"}], "]"}], 
        "]"}], "\[Equal]", "1"}]}], ":=", 
     RowBox[{
      RowBox[{"Module", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
         "roots", ",", "z0", ",", "\[Lambda]", ",", "p", ",", "q", ",", "x", 
          ",", "n", ",", "rule", ",", "matrix"}], "}"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"roots", "=", 
          RowBox[{"Cases", "[", 
           RowBox[{"\[Mu]", ",", "_Root", ",", "\[Infinity]", ",", "1"}], 
           "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"p", "=", 
          RowBox[{
           RowBox[{"roots", "\[LeftDoubleBracket]", 
            RowBox[{"1", ",", "1"}], "\[RightDoubleBracket]"}], "[", 
           "\[Lambda]", "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"n", "=", 
          RowBox[{"Exponent", "[", 
           RowBox[{"p", ",", "\[Lambda]"}], "]"}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"rule", "=", 
          RowBox[{
           SuperscriptBox["\[Lambda]", 
            RowBox[{"m_", "/;", 
             RowBox[{"m", "\[GreaterEqual]", "n"}]}]], "\[RuleDelayed]", 
           RowBox[{"Evaluate", "[", 
            RowBox[{
             SuperscriptBox["\[Lambda]", 
              RowBox[{"m", "-", "n"}]], 
             RowBox[{"(", 
              RowBox[{"-", 
               RowBox[{"(", 
                RowBox[{"p", "-", 
                 RowBox[{
                  RowBox[{"Coefficient", "[", 
                   RowBox[{"p", ",", "\[Lambda]", ",", "n"}], "]"}], 
                  SuperscriptBox["\[Lambda]", "n"]}]}], ")"}]}], ")"}]}], 
            "]"}]}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"q", "=", 
          RowBox[{"Table", "[", 
           RowBox[{
            SuperscriptBox["\[Lambda]", "i"], ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", "0", ",", 
              RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"x", "=", 
          RowBox[{"\[Mu]", "/.", 
           RowBox[{
            RowBox[{
            "roots", "\[LeftDoubleBracket]", "1", "\[RightDoubleBracket]"}], 
            "\[Rule]", "\[Lambda]"}]}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"matrix", "=", 
          RowBox[{
           RowBox[{
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{"Coefficient", "[", 
               RowBox[{"#", ",", "\[Lambda]", ",", "k"}], "]"}], ",", 
              RowBox[{"{", 
               RowBox[{"k", ",", "0", ",", 
                RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}], "&"}], "/@", 
           RowBox[{"FixedPoint", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"Expand", "[", 
               RowBox[{"#", "//.", "rule"}], "]"}], "&"}], ",", 
             RowBox[{"q", " ", "x"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"CharacteristicPolynomial", "[", 
           RowBox[{"matrix", ",", "z0"}], "]"}], "/.", 
          RowBox[{"z0", "\[Rule]", "z"}]}]}]}], "\[IndentingNewLine]", "]"}], 
      "\[IndentingNewLine]", 
      RowBox[{"Protect", "[", "MinimalPolynomial", "]"}]}]}], ";"}], 
   "*)"}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.4359540111090126`*^9, 3.435954021884507*^9}, 
   3.4359541158095646`*^9, {3.495995737530705*^9, 3.4959958161145477`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"qFromd", "[", "d_", "]"}], ":=", " ", 
  RowBox[{"RootReduce", "[", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"d", "+", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"d", "^", "2"}], "-", "4"}], "]"}]}], ")"}], "/", "2"}], 
   "]"}]}]], "Input",
 InitializationCell->True,
 CellChangeTimes->{{3.466790488757832*^9, 3.466790500732782*^9}, {
  3.466790782668639*^9, 3.466790810587059*^9}, {3.466790844035144*^9, 
  3.4667908521631937`*^9}, {3.4667971973237543`*^9, 3.466797201954014*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"qInteger", "[", "n_", "]"}], "[", "q_", "]"}], ":=", 
  RowBox[{"RootReduce", "[", 
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{"q", "^", "i"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{
        RowBox[{"-", "n"}], "+", "1"}], ",", 
       RowBox[{"n", "-", "1"}], ",", "2"}], "}"}]}], "]"}], "]"}]}]], "Input",\

 InitializationCell->True,
 CellChangeTimes->{{3.466798284808736*^9, 3.4667983105035467`*^9}, {
  3.466798345426547*^9, 3.4667983748570004`*^9}, {3.46679840595364*^9, 
  3.4667984061296673`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"End", "[", "]"}], ";"}]], "Input",
 InitializationCell->True],

Cell[BoxData[
 RowBox[{
  RowBox[{"EndPackage", "[", "]"}], ";"}]], "Input",
 InitializationCell->True]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
AutoGeneratedPackage->Automatic,
WindowSize->{759, 738},
WindowMargins->{{Automatic, 7}, {Automatic, 13}},
TaggingRules->{"SearchTextTranslated" -> ""},
FrontEndVersion->"7.0 for Mac OS X x86 (32-bit) (February 18, 2009)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{
 "Union"->{
  Cell[30434, 734, 488, 16, 27, "Input",
   InitializationCell->True,
   CellTags->"Union"]}
 }
*)
(*CellTagsIndex
CellTagsIndex->{
 {"Union", 43772, 1128}
 }
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 109, 1, 46, "Subtitle"],
Cell[679, 25, 155, 3, 26, "Text"],
Cell[837, 30, 31, 0, 67, "Section"],
Cell[CellGroupData[{
Cell[893, 34, 33, 0, 37, "Section"],
Cell[CellGroupData[{
Cell[951, 38, 165, 2, 24, "Subsubsection"],
Cell[1119, 42, 3759, 59, 268, "Input",
 InitializationCell->True],
Cell[4881, 103, 607, 13, 73, "Input",
 InitializationCell->True],
Cell[5491, 118, 669, 14, 73, "Input",
 InitializationCell->True],
Cell[6163, 134, 989, 17, 73, "Input",
 InitializationCell->True],
Cell[7155, 153, 1709, 30, 118, "Input",
 InitializationCell->True],
Cell[8867, 185, 1080, 19, 73, "Input",
 InitializationCell->True],
Cell[9950, 206, 893, 18, 58, "Input",
 InitializationCell->True],
Cell[10846, 226, 362, 8, 43, "Input",
 InitializationCell->True],
Cell[11211, 236, 305, 6, 43, "Input",
 InitializationCell->True],
Cell[11519, 244, 521, 11, 43, "Input",
 InitializationCell->True],
Cell[12043, 257, 314, 7, 27, "Input",
 InitializationCell->True],
Cell[12360, 266, 736, 15, 58, "Input",
 InitializationCell->True],
Cell[13099, 283, 2214, 36, 133, "Input",
 InitializationCell->True],
Cell[15316, 321, 1099, 20, 103, "Input",
 InitializationCell->True],
Cell[16418, 343, 733, 16, 58, "Input",
 InitializationCell->True],
Cell[17154, 361, 1386, 23, 163, "Input",
 InitializationCell->True],
Cell[18543, 386, 583, 13, 73, "Input",
 InitializationCell->True],
Cell[19129, 401, 465, 11, 58, "Input",
 InitializationCell->True],
Cell[19597, 414, 404, 9, 27, "Input",
 InitializationCell->True]
}, Open  ]],
Cell[CellGroupData[{
Cell[20038, 428, 118, 1, 24, "Subsubsection"],
Cell[20159, 431, 355, 7, 73, "Input",
 InitializationCell->True]
}, Open  ]],
Cell[CellGroupData[{
Cell[20551, 443, 204, 3, 24, "Subsubsection"],
Cell[20758, 448, 199, 4, 26, "Text"],
Cell[20960, 454, 709, 14, 73, "Input",
 InitializationCell->True],
Cell[21672, 470, 572, 12, 103, "Input",
 InitializationCell->True],
Cell[22247, 484, 246, 6, 27, "Input",
 InitializationCell->True],
Cell[22496, 492, 269, 7, 27, "Input",
 InitializationCell->True],
Cell[22768, 501, 247, 6, 58, "Input",
 InitializationCell->True],
Cell[23018, 509, 259, 6, 27, "Input",
 InitializationCell->True],
Cell[23280, 517, 330, 10, 43, "Input",
 InitializationCell->True],
Cell[23613, 529, 199, 6, 27, "Input",
 InitializationCell->True],
Cell[23815, 537, 465, 11, 77, "Input",
 InitializationCell->True],
Cell[24283, 550, 530, 15, 58, "Input",
 InitializationCell->True],
Cell[24816, 567, 499, 11, 88, "Input",
 InitializationCell->True],
Cell[25318, 580, 917, 23, 88, "Input",
 InitializationCell->True],
Cell[26238, 605, 195, 6, 27, "Input",
 InitializationCell->True],
Cell[26436, 613, 201, 6, 27, "Input",
 InitializationCell->True],
Cell[26640, 621, 191, 6, 27, "Input",
 InitializationCell->True],
Cell[26834, 629, 191, 4, 27, "Input",
 InitializationCell->True],
Cell[27028, 635, 234, 6, 27, "Input",
 InitializationCell->True],
Cell[27265, 643, 119, 3, 27, "Input",
 InitializationCell->True],
Cell[27387, 648, 906, 23, 88, "Input",
 InitializationCell->True],
Cell[28296, 673, 672, 18, 43, "Input",
 InitializationCell->True],
Cell[28971, 693, 1460, 39, 148, "Input",
 InitializationCell->True],
Cell[30434, 734, 488, 16, 27, "Input",
 InitializationCell->True,
 CellTags->"Union"],
Cell[30925, 752, 347, 8, 27, "Input",
 InitializationCell->True],
Cell[31275, 762, 212, 6, 27, "Input",
 InitializationCell->True],
Cell[31490, 770, 818, 22, 43, "Input",
 InitializationCell->True],
Cell[32311, 794, 480, 13, 27, "Input",
 InitializationCell->True],
Cell[32794, 809, 344, 8, 27, "Input",
 InitializationCell->True],
Cell[33141, 819, 560, 14, 88, "Input",
 InitializationCell->True],
Cell[33704, 835, 196, 6, 27, "Input",
 InitializationCell->True],
Cell[33903, 843, 1887, 54, 133, "Input",
 InitializationCell->True],
Cell[35793, 899, 436, 14, 27, "Input",
 InitializationCell->True],
Cell[36232, 915, 309, 10, 27, "Input",
 InitializationCell->True],
Cell[36544, 927, 255, 7, 27, "Input",
 InitializationCell->True],
Cell[36802, 936, 5230, 130, 351, "Input",
 InitializationCell->True],
Cell[42035, 1068, 538, 14, 27, "Input",
 InitializationCell->True],
Cell[42576, 1084, 586, 17, 27, "Input",
 InitializationCell->True],
Cell[43165, 1103, 96, 3, 27, "Input",
 InitializationCell->True],
Cell[43264, 1108, 103, 3, 27, "Input",
 InitializationCell->True]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)