ControlTheory.nb

(* Content-type: application/vnd.wolfram.mathematica *)

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

(* CreatedBy='Mathematica 10.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     40898,       1187]
NotebookOptionsPosition[     36144,       1032]
NotebookOutlinePosition[     37977,       1089]
CellTagsIndexPosition[     37521,       1075]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["ControlTheory Documentation", "Chapter",
 CellChangeTimes->{{3.6333515346579514`*^9, 3.63335154627586*^9}, {
  3.633364320690548*^9, 3.63336432261762*^9}, {3.633781291117386*^9, 
  3.633781292159832*^9}, {3.633781364287074*^9, 3.633781398056273*^9}, {
  3.6337814322542763`*^9, 3.633781434518094*^9}, {3.634304392561145*^9, 
  3.634304394886302*^9}, {3.6343299506479483`*^9, 3.6343299526643972`*^9}}],

Cell[CellGroupData[{

Cell["Preamble", "Subsection",
 CellChangeTimes->{{3.633351558518317*^9, 3.633351560041299*^9}, {
  3.6337812001050663`*^9, 3.633781202247016*^9}, {3.634328305258955*^9, 
  3.634328305340747*^9}, {3.636214360639636*^9, 3.636214362583493*^9}}],

Cell[BoxData[
 RowBox[{"Needs", "[", "\"\<ControlTheory`\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.633351696261045*^9, 3.633351701530016*^9}, {
  3.633781146956828*^9, 3.633781148275402*^9}, {3.633781294331254*^9, 
  3.633781296275483*^9}, {3.633781369011711*^9, 3.633781370275049*^9}, {
  3.633781401233733*^9, 3.633781430288822*^9}, {3.634304397916306*^9, 
  3.6343043997977057`*^9}, {3.6343168373823338`*^9, 3.6343168385706367`*^9}, {
  3.634329956366181*^9, 3.634329957877204*^9}}],

Cell["\<\
The following packages are needed to run some code found in this \
documentation notebook.\
\>", "Text",
 CellChangeTimes->{{3.633781139005548*^9, 3.6337811426309013`*^9}}],

Cell[BoxData[{
 RowBox[{"Needs", "[", "\"\<Tensor`\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{"Needs", "[", "\"\<Visualization`\>\"", "]"}]}], "Input",
 CellChangeTimes->{{3.6343095014776506`*^9, 3.634309506874428*^9}, {
  3.634309840365309*^9, 3.634309841590867*^9}, {3.634309893412298*^9, 
  3.63430989489682*^9}, {3.634310127176152*^9, 3.634310128430011*^9}, {
  3.634310879956975*^9, 3.63431088630688*^9}, {3.634313448232463*^9, 
  3.6343134491585693`*^9}, {3.6343167133891573`*^9, 3.634316719748117*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Source Code", "Subsection",
 CellChangeTimes->{{3.6351824420646763`*^9, 3.635182446927919*^9}}],

Cell[BoxData[
 ButtonBox["\<\"Open Source Code\"\>",
  Appearance->Automatic,
  ButtonFunction:>(Needs["QUDevTools`"]; NotebookOpen[
     FileNameJoin[{QUDevTools`$QUSourcePath, 
       StringJoin[
        FileBaseName[
         NotebookFileName[]], ".m"]}]]),
  Evaluator->Automatic,
  Method->"Preemptive"]], "Output",
 CellChangeTimes->{{3.635182609676621*^9, 3.6351826444407997`*^9}, {
  3.635182842680962*^9, 3.63518284820653*^9}, {3.635182903083259*^9, 
  3.635182906619155*^9}, {3.635182952870352*^9, 3.635182956778173*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Introduction and Overview", "Section",
 CellChangeTimes->{{3.633351563673402*^9, 3.633351567143589*^9}}],

Cell["\<\
This package provides tools commonly used in quantum control theory.\
\>", "Text",
 CellChangeTimes->{{3.6343079942006617`*^9, 3.6343080255678368`*^9}, 
   3.634308056127495*^9, {3.634308094880814*^9, 3.634308095527034*^9}, {
   3.634308126663557*^9, 3.6343082843108683`*^9}, {3.634309233626235*^9, 
   3.6343092358176537`*^9}, {3.634329976232435*^9, 3.63432998670681*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Generating Lie Algebras", "Section",
 CellChangeTimes->{{3.634330005229486*^9, 3.634330012007537*^9}, {
  3.634330526524366*^9, 3.634330530756113*^9}}],

Cell[TextData[{
 StyleBox["GenerateLieAlgebra[", "Input"],
 StyleBox["g1", "Input",
  FontWeight->"Plain",
  FontSlant->"Italic"],
 StyleBox[",", "Input"],
 StyleBox["g2", "Input",
  FontWeight->"Plain",
  FontSlant->"Italic"],
 StyleBox[",", "Input"],
 StyleBox["...", "Input",
  FontWeight->"Plain",
  FontSlant->"Italic"],
 StyleBox["]", "Input"],
 " generates a basis for the minimal matrix Lie algebra containing the given \
generators ",
 StyleBox["g1,g2,...", "Input"],
 ". The output is of the form ",
 StyleBox["{dim, labels, basis}", "Input"],
 " where ",
 StyleBox["dim", "Input"],
 " is the computed Lie dimension, ",
 StyleBox["labels", "Input"],
 " gives a readable representation of what was commuted with what to get the ",
 StyleBox["basis", "Input"],
 ", and ",
 StyleBox["basis", "Input"],
 " is a list of linearly independent matrices which span the matrix Lie \
algebra.\n\n",
 StyleBox["Notes",
  FontWeight->"Bold",
  FontSlant->"Italic",
  FontColor->GrayLevel[0.5]],
 "\n",
 StyleBox["\[Bullet]",
  FontColor->GrayLevel[0]],
 " In general, the elements of ",
 StyleBox["basis", "Input"],
 " will not be equal to operations described in ",
 StyleBox["labels", "Input"],
 "; to remain numerically stable, it is often advantageous to work with an \
orthogonalized and normalized system.\n",
 StyleBox["\[Bullet]",
  FontColor->GrayLevel[0]],
 " In the case of testing for complete control of a quantum system, enter the \
system drift and control Hamiltonians as input generators, and ensure ",
 StyleBox["AssumeHamiltonians", "Input"],
 " is ",
 StyleBox["True", "Input"],
 " (it is by default). \n",
 StyleBox["\[Bullet]",
  FontColor->GrayLevel[0]],
 " To have complete control of a quantum system, it is required that the \
resulting Lie algebra have dimension at least ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    SuperscriptBox["d", "2"], "-", "1"}], TraditionalForm]]],
 ", where ",
 Cell[BoxData[
  FormBox["d", TraditionalForm]]],
 " is the Hilbert space dimension.\n",
 StyleBox["\[Bullet]",
  FontColor->GrayLevel[0]],
 " The free-basis elements are chosen using the Phillip Hall basis algorithm"
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{{3.634330561147408*^9, 3.634330600514605*^9}, {
  3.634330651205419*^9, 3.634330696599441*^9}, {3.634330769771941*^9, 
  3.6343308639556017`*^9}, {3.6343325977198772`*^9, 3.6343326053441553`*^9}, {
  3.634333236971307*^9, 3.634333384797036*^9}, {3.6343334632673597`*^9, 
  3.634333482579919*^9}, {3.634334269518077*^9, 3.634334326309326*^9}, {
  3.63439528485697*^9, 3.634395338649721*^9}, {3.634395407273575*^9, 
  3.634395487239188*^9}, {3.634396258581321*^9, 3.634396298643025*^9}},
 CellTags->"GenerateLieAlgebra::usage"],

Cell[CellGroupData[{

Cell["References", "Subsubsection",
 CellChangeTimes->{{3.63433418924352*^9, 3.634334190219788*^9}}],

Cell["\<\
D\[CloseCurlyQuote]Alessandro, D., 2007. Introduction to Quantum Control and \
Dynamics, 1st ed. Chapman and Hall/CRC.\
\>", "Text",
 CellChangeTimes->{{3.634334206555003*^9, 3.6343342083485613`*^9}}],

Cell["\<\
Phillip Hall algorithm: http://planning.cs.uiuc.edu/node834.html\
\>", "Text",
 CellChangeTimes->{{3.634395346389385*^9, 3.634395358248818*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Options", "Subsection",
 CellChangeTimes->{{3.6343309008417053`*^9, 3.634330907002307*^9}}],

Cell[BoxData[GridBox[{
   {
    StyleBox["Option", "Text",
     FontWeight->Bold], 
    StyleBox[
     RowBox[{"Default", " ", "Value"}], "Text",
     FontWeight->Bold], 
    StyleBox["Description", "Text",
     FontWeight->Bold]},
   {"LieOperation", "Com", 
    StyleBox[
     RowBox[{
     "The", " ", "Lie", " ", "operation", " ", "to", " ", "be", " ", "used", 
      " ", "between", " ", "two", " ", "Lie", " ", "algebra", " ", "elements",
       " ", "in", " ", "their", " ", "matrix", " ", 
      RowBox[{"representation", "."}]}], "Text"]},
   {"Orthogonalize", "True", 
    RowBox[{
     RowBox[{
      StyleBox["If", "Text"], 
      StyleBox[" ", "Text"], 
      StyleBox["True", "Input"]}], 
     StyleBox[",", "Text"], 
     StyleBox[" ", "Text"], 
     RowBox[{
      StyleBox["the", "Text"], 
      StyleBox[" ", "Text"], 
      StyleBox["output", "Text"], 
      StyleBox[" ", "Text"], 
      StyleBox["basis", "Input"], 
      StyleBox[" ", "Text"], 
      StyleBox["is", "Text"], 
      StyleBox[" ", "Text"], 
      StyleBox["orthogonal", "Text"]}], 
     StyleBox[",", "Text"], 
     StyleBox[" ", "Text"], 
     StyleBox[
      RowBox[{
      "and", " ", "the", " ", "algorithm", " ", "uses", " ", "an", " ", 
       "orthogonal", " ", "basis", " ", "in", " ", "its", " ", 
       RowBox[{"calculations", ".", " ", "This"}], " ", "improves", " ", 
       "numerical", " ", "stability", " ", "and", " ", 
       RowBox[{"efficiency", "."}]}], "Text"]}]},
   {"NumericalLieAlgebra", "False", 
    StyleBox[
     RowBox[{"Whether", " ", "or", " ", "not", " ", "to", " ", "call", " ", 
      TagBox[
       StyleBox["N", "Input",
        FontFamily->"Courier"],
       DisplayForm], " ", "on", " ", "inputs", " ", "in", " ", "order", " ", 
      "to", " ", "perform", " ", "all", " ", "calculations", " ", 
      RowBox[{"numerically", "."}]}], "Text"]},
   {"VSField", "Reals", 
    StyleBox[
     RowBox[{
      RowBox[{
      "The", " ", "field", " ", "of", " ", "the", " ", "Lie", " ", 
       "algebra"}], ",", " ", 
      RowBox[{"either", " ", 
       TagBox[
        StyleBox["Reals", "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", "or", " ", 
       RowBox[{
        TagBox[
         StyleBox["Complex", "Input",
          FontFamily->"Courier"],
         DisplayForm], ".", " ", "This"}], " ", "is", " ", "necessary", " ", 
       "as", " ", "two", " ", "linearly", " ", "independent", " ", "vectors", 
       " ", "in", " ", "one", " ", "vector", " ", "space", " ", "may", " ", 
       "not", " ", "be", " ", "linearly", " ", "independent", " ", "in", " ", 
       "the", " ", "same", " ", "space", " ", "but", " ", "with", " ", "a", 
       " ", "different", " ", "field"}], ",", " ", 
      RowBox[{
       RowBox[{
        RowBox[{"e", ".", "g", ".", " ", 
         TagBox[Cell[BoxData[
           FormBox[
            RowBox[{"(", 
             RowBox[{"1", ",", " ", "0"}], ")"}], TraditionalForm]]],
          DisplayForm]}], " ", "and", " ", 
        RowBox[{
         TagBox[Cell[BoxData[
           FormBox[
            RowBox[{"(", 
             RowBox[{"i", ",", " ", "0"}], ")"}], TraditionalForm]]],
          DisplayForm], ".", " ", "In"}], " ", "the", " ", "case", " ", "of", 
        " ", "anti"}], "-", 
       RowBox[{"hermitian", " ", "generators", " ", 
        RowBox[{"(", 
         RowBox[{"quantum", " ", "theory"}], ")"}]}]}], ",", " ", 
      RowBox[{"the", " ", "real", " ", "vector", " ", "space", " ", "is", " ", 
       RowBox[{"used", "."}]}]}], "Text"]},
   {"LinearIndependenceFun", "Automatic", 
    StyleBox[
     RowBox[{
      RowBox[{"A", " ", "boolean", " ", "function", " ", 
       TagBox[
        StyleBox[
         RowBox[{"Fun", "[", 
          RowBox[{"M", ",", "V"}], "]"}], "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", "which", " ", "decides", " ", "whether", " ", 
       "the", " ", "vector", " ", 
       TagBox[
        StyleBox["V", "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", "is", " ", "linearly", " ", "independent", " ", 
       "from", " ", "the", " ", "space", " ", "spanned", " ", "by", " ", 
       "the", " ", "list", " ", "of", " ", "vectors", " ", 
       RowBox[{
        TagBox[
         StyleBox["M", "Input",
          FontFamily->"Courier"],
         DisplayForm], ".", " ", "If"}], " ", 
       TagBox[
        StyleBox["Orthogonalize", "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", "is", " ", 
       TagBox[
        StyleBox["True", "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", "in", " ", 
       TagBox[
        StyleBox["GenerateLieAlgebra", "Input",
         FontFamily->"Courier"],
        DisplayForm]}], ",", " ", 
      RowBox[{
       TagBox[
        StyleBox["M", "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", "is", " ", "guaranteed", " ", "to", " ", "be", " ",
        "orthonormal"}], ",", " ", 
      RowBox[{"and", " ", 
       TagBox[
        StyleBox["v", "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", 
       RowBox[{"normalized", ".", " ", "The"}], " ", "default", " ", "value", 
       " ", "of", " ", 
       TagBox[
        StyleBox["Automatic", "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", "uses", " ", 
       TagBox[
        StyleBox[
         RowBox[{
          RowBox[{
           RowBox[{"Norm", "[", 
            RowBox[{"#1", ".", "#2"}], "]"}], ">", 
           RowBox[{"OptionValue", "[", "Threshold", "]"}]}], "&"}], "Input",
         FontFamily->"Courier"],
        DisplayForm], " ", "in", " ", "the", " ", "case", " ", 
       RowBox[{
        TagBox[
         StyleBox[
          RowBox[{"VSField", "\[Rule]", "Reals"}], "Input",
          FontFamily->"Courier"],
         DisplayForm], "."}]}]}], "Text"]},
   {"StoppingDepth", "\[Infinity]", 
    StyleBox[
     RowBox[{
     "The", " ", "commutation", " ", "depth", " ", "at", " ", "which", " ", 
      "to", " ", "stop", " ", "searching", " ", "for", " ", "new", " ", 
      "linearly", " ", "independent", " ", "operators", " ", "in", " ", 
      RowBox[{
       TagBox[
        StyleBox["GenerateLieAlgebra", "Input",
         FontFamily->"Courier"],
        DisplayForm], "."}]}], "Text"]},
   {"StoppingDimension", "Automatic", 
    StyleBox[
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{
        "The", " ", "dimension", " ", "at", " ", "which", " ", "to", " ", 
         "stop", " ", "searching", " ", "for", " ", "new", " ", "linearly", 
         " ", "independent", " ", "operators", " ", "in", " ", 
         RowBox[{
          TagBox[
           StyleBox["GenerateLieAlgebra", "Input",
            FontFamily->"Courier"],
           DisplayForm], ".", " ", "The"}], " ", "default", " ", "of", " ", 
         TagBox[
          StyleBox["Automatic", "Input",
           FontFamily->"Courier"],
          DisplayForm], " ", "results", " ", "in", " ", "a", " ", "value", 
         " ", 
         TagBox[Cell[BoxData[
           FormBox[
            RowBox[{
             SuperscriptBox["d", "2"], " ", "-", " ", "1"}], 
            TraditionalForm]]],
          DisplayForm]}], ",", " ", 
        RowBox[{"where", " ", 
         TagBox[Cell[BoxData[
           FormBox["d", TraditionalForm]]],
          DisplayForm], " ", "is", " ", "the", " ", "dimension", " ", "the", 
         " ", "Lie", " ", "algebra"}]}], "\[CloseCurlyQuote]"}], "s", " ", 
      "elements", " ", "act", " ", 
      RowBox[{"on", "."}]}], "Text"]},
   {"AssumeHamiltonians", "True", 
    StyleBox[
     RowBox[{
      RowBox[{"When", " ", 
       TagBox[
        StyleBox["True", "Input",
         FontFamily->"Courier"],
        DisplayForm]}], ",", " ", 
      RowBox[{
      "makes", " ", "the", " ", "assumption", " ", "that", " ", "inputs", " ",
        "are", " ", "hermitian"}], ",", " ", 
      RowBox[{
      "and", " ", "that", " ", "you", " ", "are", " ", "looking", " ", "for", 
       " ", "the", " ", "hermitian", " ", "algebra", " ", "\[ScriptCapitalA]",
        " ", "in", " ", 
       TagBox[Cell[BoxData[
         FormBox[
          RowBox[{
           RowBox[{"i", " ", "*", " ", "\[ScriptCapitalA]"}], " ", "=", " ", 
           RowBox[{"\[ScriptCapitalB]", " ", "=", " ", 
            RowBox[{"\[LeftAngleBracket]", " ", 
             RowBox[{"\[LeftAngleBracket]", " ", 
              RowBox[{
               RowBox[{"i", " ", "*", " ", 
                SubscriptBox["H", "1"]}], ",", " ", 
               RowBox[{"i", " ", "*", " ", 
                SubscriptBox["H", "2"]}], ",", " ", "..."}], " ", 
              "\[RightAngleBracket]"}], " ", "\[RightAngleBracket]"}]}]}], 
          TraditionalForm]]],
        DisplayForm], " ", "where", " ", "\[ScriptCapitalB]", " ", "is", " ", 
       "actually", " ", "a", " ", "skew", " ", "hermitian", " ", 
       RowBox[{"algebra", ".", " ", "This"}], " ", "is", " ", "a", " ", 
       "simple", " ", "though", " ", "subtle", " ", 
       RowBox[{"point", ".", " ", "Default"}], " ", "is", " ", 
       RowBox[{
        TagBox[
         StyleBox["True", "Input",
          FontFamily->"Courier"],
         DisplayForm], "."}]}]}], "Text"]},
   {"Threshold", 
    FractionBox["1", "10000000000"], 
    StyleBox[
     RowBox[{
      RowBox[{"Threshold", "[", 
       StyleBox["data", "TI"], "]"}], " ", "thresholds", " ", 
      StyleBox["data", "TI"], " ", "by", " ", "replacing", " ", "values", " ",
       "close", " ", "to", " ", "zero", " ", "by", " ", 
      RowBox[{"zero", ".", "\n", 
       RowBox[{"Threshold", "[", 
        RowBox[{
         StyleBox["data", "TI"], ",", 
         StyleBox["tspec", "TI"]}], "]"}]}], " ", "thresholds", " ", 
      StyleBox["data", "TI"], " ", "using", " ", "threshold", " ", 
      "specification", " ", 
      RowBox[{
       StyleBox["tspec", "TI"], ".", "\n", 
       RowBox[{"Threshold", "[", 
        RowBox[{
         StyleBox["image", "TI"], ",", 
         StyleBox["\[Ellipsis]", "TR"]}], "]"}]}], " ", "thresholds", " ", 
      "an", " ", 
      RowBox[{"image", ".", "\n", 
       RowBox[{"Threshold", "[", 
        RowBox[{
         StyleBox["sound", "TI"], ",", 
         StyleBox["\[Ellipsis]", "TR"]}], "]"}]}], " ", "thresholds", " ", 
      "a", " ", "sound", " ", 
      RowBox[{"object", "."}]}], "Text"]},
   {"GeneratorSymbols", "Automatic", 
    StyleBox[
     RowBox[{
     "A", " ", "list", " ", "of", " ", "symbols", " ", "or", " ", "strings", 
      " ", "to", " ", "use", " ", "as", " ", "human", " ", "readable", " ", 
      "labels", " ", "for", " ", "the", " ", "input", " ", 
      RowBox[{"generators", ".", " ", "The"}], " ", "default", " ", "value", 
      " ", "of", " ", 
      TagBox[
       StyleBox["Automatic", "Input",
        FontFamily->"Courier"],
       DisplayForm], " ", "results", " ", "in", " ", "generator", " ", 
      "labels", " ", 
      TagBox[
       StyleBox[
        RowBox[{"{", 
         RowBox[{
          RowBox[{"G", "[", "1", "]"}], ",", "...", ",", 
          RowBox[{"G", "[", "n", "]"}]}], "}"}], "Input",
        FontFamily->"Courier"],
       DisplayForm], "  ", "where", " ", 
      TagBox[
       StyleBox["G", "Input",
        FontFamily->"Courier"],
       DisplayForm], " ", "is", " ", "a", " ", 
      TagBox[
       StyleBox["Head", "Input",
        FontFamily->"Courier"],
       DisplayForm], " ", "reserved", " ", "for", " ", "this", " ", 
      RowBox[{"purpose", "."}]}], "Text"]}
  },
  GridBoxAlignment->{"Columns" -> {{Left}}},
  GridBoxDividers->{"Columns" -> {{False}}, "Rows" -> {{True}}},
  GridBoxItemSize->{"Columns" -> {Automatic, Automatic, {
       Scaled[0.5]}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}}, 
    "RowsIndexed" -> {}}]], "Text",
 CellChangeTimes->{
  3.634395120071556*^9, {3.634396430967825*^9, 3.6343964416653347`*^9}, {
   3.63439648529278*^9, 3.634396568950224*^9}}],

Cell[CellGroupData[{

Cell["LieOperation", "Subsubsection",
 CellChangeTimes->{{3.634330953849202*^9, 3.634330956983987*^9}}],

Cell["\<\
The Lie operation to be used between two Lie algebra elements in their matrix \
representation.\
\>", "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{{3.634330974025662*^9, 3.6343310173063097`*^9}},
 CellTags->"LieOperation::usage"]
}, Open  ]],

Cell[CellGroupData[{

Cell["VSField", "Subsubsection",
 CellChangeTimes->{{3.634331025257421*^9, 3.634331026456582*^9}}],

Cell[TextData[{
 "The field of the Lie algebra, either ",
 StyleBox["Reals", "Input"],
 " or ",
 StyleBox["Complex", "Input"],
 ". This is necessary as two linearly independent vectors in one vector space \
may not be linearly independent in the same space but with a different field, \
e.g. ",
 Cell[BoxData[
  FormBox[
   RowBox[{"(", 
    RowBox[{"1", ",", "0"}], ")"}], TraditionalForm]]],
 " and ",
 Cell[BoxData[
  FormBox[
   RowBox[{"(", 
    RowBox[{"i", ",", "0"}], ")"}], TraditionalForm]]],
 ". In the case of anti-hermitian generators (quantum theory), the real \
vector space is used."
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{{3.63433103372034*^9, 3.634331240041471*^9}, {
  3.634331621558449*^9, 3.63433162530208*^9}, {3.6343950363694477`*^9, 
  3.634395044143812*^9}},
 CellTags->"VSField::usage"]
}, Open  ]],

Cell[CellGroupData[{

Cell["LinearIndependenceFun", "Subsubsection",
 CellChangeTimes->{{3.634331245791628*^9, 3.634331251361903*^9}}],

Cell[TextData[{
 "A boolean function ",
 StyleBox["Fun[M,V]", "Input"],
 " which decides whether the vector ",
 StyleBox["V", "Input"],
 " is linearly independent from the space spanned by the list of vectors ",
 StyleBox["M", "Input"],
 ". If ",
 StyleBox["Orthogonalize", "Input"],
 " is ",
 StyleBox["True", "Input"],
 " in ",
 StyleBox["GenerateLieAlgebra", "Input"],
 ", ",
 StyleBox["M", "Input"],
 " is guaranteed to be orthonormal, and ",
 StyleBox["v", "Input"],
 " normalized. The default value of ",
 StyleBox["Automatic", "Input"],
 " uses ",
 StyleBox["Norm[#1.#2]>OptionValue[Threshold]&", "Input"],
 " in the case ",
 StyleBox["VSField\[Rule]Reals", "Input"],
 "."
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{
  3.634331260066362*^9, {3.634331314328472*^9, 3.6343313454959917`*^9}},
 CellTags->"LinearIndependenceFun::usage"]
}, Open  ]],

Cell[CellGroupData[{

Cell["StoppingDepth", "Subsubsection",
 CellChangeTimes->{{3.6343313560770483`*^9, 3.634331365477507*^9}}],

Cell[TextData[{
 "The commutation depth at which to stop searching for new linearly \
independent operators in ",
 StyleBox["GenerateLieAlgebra", "Input"],
 "."
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{{3.634331376680532*^9, 3.6343313843108463`*^9}},
 CellTags->"StoppingDepth::usage"]
}, Open  ]],

Cell[CellGroupData[{

Cell["StoppingDimension", "Subsubsection",
 CellChangeTimes->{{3.634331400438117*^9, 3.634331405599959*^9}}],

Cell[TextData[{
 "The dimension at which to stop searching for new linearly independent \
operators in ",
 StyleBox["GenerateLieAlgebra", "Input"],
 ". The default of ",
 StyleBox["Automatic", "Input"],
 " results in a value ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    SuperscriptBox["d", "2"], "-", "1"}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 ", where ",
 Cell[BoxData[
  FormBox["d", TraditionalForm]],
  FormatType->"TraditionalForm"],
 " is the dimension the Lie algebra\[CloseCurlyQuote]s elements act on."
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{{3.6343314163837757`*^9, 3.6343314894780903`*^9}},
 CellTags->"StoppingDimension::usage"]
}, Open  ]],

Cell[CellGroupData[{

Cell["AssumeHamiltonians", "Subsubsection",
 CellChangeTimes->{{3.634331500565523*^9, 3.63433150511078*^9}}],

Cell[TextData[{
 "When ",
 StyleBox["True", "Input"],
 ", makes the assumption that inputs are hermitian, and that you are looking \
for the hermitian algebra \[ScriptCapitalA] in ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    RowBox[{"i", "*", "\[ScriptCapitalA]"}], "=", 
    RowBox[{"\[ScriptCapitalB]", "=", 
     RowBox[{"\[LeftAngleBracket]", 
      RowBox[{"\[LeftAngleBracket]", 
       RowBox[{
        RowBox[{"i", "*", 
         SubscriptBox["H", "1"]}], ",", 
        RowBox[{"i", "*", 
         SubscriptBox["H", "2"]}], ",", "..."}], "\[RightAngleBracket]"}], 
      "\[RightAngleBracket]"}]}]}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 " where \[ScriptCapitalB] is actually a skew hermitian algebra. This is a \
simple though subtle point. Default is ",
 StyleBox["True", "Input"],
 "."
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{{3.63433156127283*^9, 3.634331606686656*^9}, {
   3.634331636764819*^9, 3.6343317601997833`*^9}, 3.63433179937297*^9},
 CellTags->"AssumeHamiltonians::usage"]
}, Open  ]],

Cell[CellGroupData[{

Cell["GeneratorSymbols", "Subsubsection",
 CellChangeTimes->{{3.634331507844613*^9, 3.634331510205758*^9}}],

Cell[TextData[{
 "A list of symbols or strings to use as human readable labels for the input \
generators. The default value of ",
 StyleBox["Automatic", "Input"],
 " results in generator labels ",
 StyleBox["{G[1],...,G[n]}", "Input"],
 "  where ",
 StyleBox["G", "Input"],
 " is a ",
 StyleBox["Head", "Input"],
 " reserved for this purpose."
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{
  3.634331570183009*^9, 3.6343318614205112`*^9, {3.634332353552425*^9, 
   3.63433243007194*^9}},
 CellTags->"GeneratorSymbols::usage"],

Cell[TextData[{
 StyleBox["G", "Input"],
 " is the default ",
 StyleBox["Head", "Input"],
 " used to label Lie algebra elements."
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{{3.634332029241769*^9, 3.6343320424330387`*^9}, {
  3.634332132656271*^9, 3.634332134759694*^9}},
 CellTags->"G::usage"]
}, Open  ]],

Cell[CellGroupData[{

Cell["NumericalLieAlgebra", "Subsubsection",
 CellChangeTimes->{{3.634331523652102*^9, 3.634331551259983*^9}}],

Cell[TextData[{
 "Whether or not to call ",
 StyleBox["N", "Input"],
 " on inputs in order to perform all calculations numerically."
}], "Text",
 CellDingbat->"\[FilledSquare]",
 CellChangeTimes->{3.6343315778878403`*^9},
 CellTags->"NumericalLieAlgebra::usage"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Qubit Example", "Subsection",
 CellChangeTimes->{{3.634332640572611*^9, 3.634332642366716*^9}, {
  3.63433353551357*^9, 3.6343335361766577`*^9}}],

Cell[TextData[{
 "We see that rotations of a qubit about the x and y axes indeed generate \
full control; we require a Lie algebra dimension of at least ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    RowBox[{
     SuperscriptBox["2", "2"], "-", "1"}], "=", "3"}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "."
}], "Text",
 CellChangeTimes->{{3.634333197907205*^9, 3.6343332289873123`*^9}, {
  3.6343335090187883`*^9, 3.6343335301133003`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"dim", ",", "labels", ",", "basis"}], "}"}], "=", 
   RowBox[{"GenerateLieAlgebra", "[", 
    RowBox[{
     RowBox[{"TP", "[", "X", "]"}], ",", 
     RowBox[{"TP", "[", "Y", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Dimension: \>\"", "<>", 
   RowBox[{"ToString", "[", "dim", "]"}]}], 
  "]"}], "\[IndentingNewLine]", "labels", "\[IndentingNewLine]", 
 RowBox[{"basis", "//", "MatrixListForm"}]}], "Input",
 CellChangeTimes->{{3.634332645544158*^9, 3.634332661426464*^9}, {
   3.634332887874597*^9, 3.634332889787505*^9}, {3.634333058956403*^9, 
   3.634333074456367*^9}, {3.6343331506282597`*^9, 3.634333187506667*^9}, 
   3.634395527956313*^9, {3.634395666867566*^9, 3.6343956800832663`*^9}, {
   3.634396608679371*^9, 3.634396612288747*^9}}],

Cell[BoxData["\<\"Dimension: 3\"\>"], "Print",
 CellChangeTimes->{{3.634333171294609*^9, 3.634333187815991*^9}, 
   3.6343337271314898`*^9, {3.634395524111643*^9, 3.634395550808226*^9}, 
   3.634395599861759*^9, {3.634395647230783*^9, 3.634395680544506*^9}, {
   3.6343965921593103`*^9, 3.634396612619556*^9}, 3.6343966491850853`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"G", "[", "1", "]"}], ",", 
   RowBox[{"G", "[", "2", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "1", "]"}], ",", 
     RowBox[{"G", "[", "2", "]"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.634333187850479*^9, 3.634333727192985*^9, {3.634395524142859*^9, 
   3.634395550837936*^9}, 3.634395599868321*^9, {3.634395647252955*^9, 
   3.634395680569879*^9}, {3.634396592202182*^9, 3.634396612660242*^9}, 
   3.634396649240135*^9}],

Cell[BoxData[
 TemplateBox[{TagBox[
    RowBox[{"(", "\[NoBreak]", 
      GridBox[{{"0", 
         FractionBox["1", 
          SqrtBox["2"]]}, {
         FractionBox["1", 
          SqrtBox["2"]], "0"}}, RowSpacings -> 1, ColumnSpacings -> 1, 
       RowAlignments -> Baseline, ColumnAlignments -> Center], "\[NoBreak]", 
      ")"}], 
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]],"\",\"",TagBox[
    RowBox[{"(", "\[NoBreak]", 
      GridBox[{{"0", 
         RowBox[{"-", 
           FractionBox["\[ImaginaryI]", 
            SqrtBox["2"]]}]}, {
         FractionBox["\[ImaginaryI]", 
          SqrtBox["2"]], "0"}}, RowSpacings -> 1, ColumnSpacings -> 1, 
       RowAlignments -> Baseline, ColumnAlignments -> Center], "\[NoBreak]", 
      ")"}], 
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]],"\",\"",TagBox[
    RowBox[{"(", "\[NoBreak]", 
      GridBox[{{
         RowBox[{"-", 
           FractionBox["1", 
            SqrtBox["2"]]}], "0"}, {"0", 
         FractionBox["1", 
          SqrtBox["2"]]}}, RowSpacings -> 1, ColumnSpacings -> 1, 
       RowAlignments -> Baseline, ColumnAlignments -> Center], "\[NoBreak]", 
      ")"}], 
    Function[BoxForm`e$, 
     MatrixForm[BoxForm`e$]]]},
  "RowDefault"]], "Output",
 CellChangeTimes->{
  3.634333187850479*^9, 3.634333727192985*^9, {3.634395524142859*^9, 
   3.634395550837936*^9}, 3.634395599868321*^9, {3.634395647252955*^9, 
   3.634395680569879*^9}, {3.634396592202182*^9, 3.634396612660242*^9}, 
   3.634396649244452*^9}]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Anisotropic hyperfine control", "Subsection",
 CellChangeTimes->{{3.6343335623771963`*^9, 3.634333571011238*^9}, 
   3.634334233277602*^9}],

Cell[TextData[{
 "We rederive a standard result (see references below): if two qubits are \
coupled with an anisotropic interaction, but we can only control one of them, \
we still have full control (we need dimension ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    RowBox[{
     SuperscriptBox["4", "2"], "-", "1"}], "=", "15"}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 ")."
}], "Text",
 CellChangeTimes->{{3.634333197907205*^9, 3.6343332289873123`*^9}, {
  3.6343335090187883`*^9, 3.6343335301133003`*^9}, {3.634333824364408*^9, 
  3.634333853894891*^9}, {3.634333892389818*^9, 3.634333945406625*^9}, {
  3.6343340929574347`*^9, 3.634334097462557*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"dim", ",", "labels", ",", "basis"}], "}"}], "=", 
   RowBox[{"GenerateLieAlgebra", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"TP", "[", 
       RowBox[{"ZX", "+", "ZZ"}], "]"}], "+", 
      RowBox[{"TP", "[", 
       RowBox[{"IZ", "+", "ZI"}], "]"}]}], ",", 
     RowBox[{"TP", "[", "XI", "]"}], ",", 
     RowBox[{"TP", "[", "YI", "]"}], ",", 
     RowBox[{"NumericalLieAlgebra", "\[Rule]", "True"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Dimension: \>\"", "<>", 
   RowBox[{"ToString", "[", "dim", "]"}]}], 
  "]"}], "\[IndentingNewLine]", "labels"}], "Input",
 CellChangeTimes->{{3.634332645544158*^9, 3.634332661426464*^9}, {
   3.634332887874597*^9, 3.634332889787505*^9}, {3.634333058956403*^9, 
   3.634333074456367*^9}, {3.6343331506282597`*^9, 3.634333187506667*^9}, {
   3.6343335888668013`*^9, 3.634333613339383*^9}, {3.634333659907202*^9, 
   3.634333660346468*^9}, {3.634333752263914*^9, 3.6343338111543093`*^9}, 
   3.6343338850349216`*^9, {3.634333966653171*^9, 3.6343340019694557`*^9}}],

Cell[BoxData["\<\"Dimension: 15\"\>"], "Print",
 CellChangeTimes->{{3.634333171294609*^9, 3.634333187815991*^9}, {
   3.634333599241817*^9, 3.634333613708527*^9}, 3.634333660817954*^9, {
   3.634333728781124*^9, 3.634333757404483*^9}, {3.6343337890023727`*^9, 
   3.6343338116722517`*^9}, 3.6343338854297657`*^9, {3.6343339708555517`*^9, 
   3.634334002263979*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"G", "[", "1", "]"}], ",", 
   RowBox[{"G", "[", "2", "]"}], ",", 
   RowBox[{"G", "[", "3", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "1", "]"}], ",", 
     RowBox[{"G", "[", "2", "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "1", "]"}], ",", 
     RowBox[{"G", "[", "3", "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "2", "]"}], ",", 
     RowBox[{"G", "[", "3", "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "1", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "1", "]"}], ",", 
       RowBox[{"G", "[", "2", "]"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "1", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "1", "]"}], ",", 
       RowBox[{"G", "[", "3", "]"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "2", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "1", "]"}], ",", 
       RowBox[{"G", "[", "2", "]"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "2", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "1", "]"}], ",", 
       RowBox[{"G", "[", "3", "]"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "2", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "2", "]"}], ",", 
       RowBox[{"G", "[", "3", "]"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "3", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "1", "]"}], ",", 
       RowBox[{"G", "[", "2", "]"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "3", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "1", "]"}], ",", 
       RowBox[{"G", "[", "3", "]"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "3", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "2", "]"}], ",", 
       RowBox[{"G", "[", "3", "]"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"G", "[", "1", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"G", "[", "1", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"G", "[", "1", "]"}], ",", 
         RowBox[{"G", "[", "2", "]"}]}], "}"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.634333187850479*^9, {3.634333599258515*^9, 3.634333613725987*^9}, 
   3.634333660837123*^9, {3.634333728794608*^9, 3.634333757416964*^9}, {
   3.6343337890141783`*^9, 3.6343338116863403`*^9}, 3.634333885454494*^9, {
   3.634333970852229*^9, 3.6343340022873287`*^9}}]
}, Open  ]],

Cell["\<\
The results of course scale to any number of qubits, where only the first \
qubit is controllable.\
\>", "Text",
 CellChangeTimes->{{3.6343340688302298`*^9, 3.634334087182949*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"{", 
    RowBox[{"dim", ",", "labels", ",", "basis"}], "}"}], "=", 
   RowBox[{"GenerateLieAlgebra", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"TP", "[", 
       RowBox[{"ZXI", "+", "ZZI"}], "]"}], "+", 
      RowBox[{"TP", "[", 
       RowBox[{"ZIX", "+", "ZIZ"}], "]"}], "+", 
      RowBox[{"TP", "[", 
       RowBox[{"ZII", "+", "IZI", "+", "IIZ"}], "]"}]}], ",", 
     RowBox[{"TP", "[", "XII", "]"}], ",", 
     RowBox[{"TP", "[", "YII", "]"}], ",", 
     RowBox[{"NumericalLieAlgebra", "\[Rule]", "True"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Dimension: \>\"", "<>", 
   RowBox[{"ToString", "[", "dim", "]"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.6343339568356123`*^9, 3.634333962216814*^9}, {
  3.6343340084699306`*^9, 3.63433405682554*^9}, {3.634334130248294*^9, 
  3.634334130407775*^9}}],

Cell[BoxData["\<\"Dimension: 63\"\>"], "Print",
 CellChangeTimes->{{3.634334038746757*^9, 3.634334057216703*^9}, 
   3.6343341308677273`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["References", "Subsubsection",
 CellChangeTimes->{{3.634334134618195*^9, 3.634334135731655*^9}}],

Cell["\<\
Hodges, J.S., Yang, J.C., Ramanathan, C., Cory, D.G., 2007. Universal Control \
of Nuclear Spins Via Anisotropic Hyperfine Interactions. arXiv:0707.2956. \
doi:10.1103/PhysRevA.78.010303\
\>", "Text",
 CellChangeTimes->{{3.6343341408819933`*^9, 3.634334142883691*^9}}]
}, Open  ]]
}, Closed]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1280, 972},
WindowMargins->{{87, Automatic}, {-112, Automatic}},
FrontEndVersion->"10.0 for Linux x86 (64-bit) (June 27, 2014)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{
 "GenerateLieAlgebra::usage"->{
  Cell[3846, 101, 2719, 74, 236, "Text",
   CellTags->"GenerateLieAlgebra::usage"]},
 "LieOperation::usage"->{
  Cell[19187, 514, 254, 6, 31, "Text",
   CellTags->"LieOperation::usage"]},
 "VSField::usage"->{
  Cell[19579, 528, 837, 24, 53, "Text",
   CellTags->"VSField::usage"]},
 "LinearIndependenceFun::usage"->{
  Cell[20568, 560, 860, 28, 74, "Text",
   CellTags->"LinearIndependenceFun::usage"]},
 "StoppingDepth::usage"->{
  Cell[21574, 596, 308, 8, 32, "Text",
   CellTags->"StoppingDepth::usage"]},
 "StoppingDimension::usage"->{
  Cell[22030, 612, 683, 20, 54, "Text",
   CellTags->"StoppingDimension::usage"]},
 "AssumeHamiltonians::usage"->{
  Cell[22861, 640, 1036, 27, 54, "Text",
   CellTags->"AssumeHamiltonians::usage"]},
 "GeneratorSymbols::usage"->{
  Cell[24044, 675, 545, 16, 53, "Text",
   CellTags->"GeneratorSymbols::usage"]},
 "G::usage"->{
  Cell[24592, 693, 314, 9, 32, "Text",
   CellTags->"G::usage"]},
 "NumericalLieAlgebra::usage"->{
  Cell[25056, 710, 262, 7, 32, "Text",
   CellTags->"NumericalLieAlgebra::usage"]}
 }
*)
(*CellTagsIndex
CellTagsIndex->{
 {"GenerateLieAlgebra::usage", 36448, 1042},
 {"LieOperation::usage", 36559, 1045},
 {"VSField::usage", 36657, 1048},
 {"LinearIndependenceFun::usage", 36765, 1051},
 {"StoppingDepth::usage", 36879, 1054},
 {"StoppingDimension::usage", 36988, 1057},
 {"AssumeHamiltonians::usage", 37103, 1060},
 {"GeneratorSymbols::usage", 37218, 1063},
 {"G::usage", 37315, 1066},
 {"NumericalLieAlgebra::usage", 37414, 1069}
 }
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 406, 5, 66, "Chapter"],
Cell[CellGroupData[{
Cell[1011, 31, 242, 3, 44, "Subsection"],
Cell[1256, 36, 488, 7, 32, "Input"],
Cell[1747, 45, 182, 4, 31, "Text"],
Cell[1932, 51, 512, 7, 55, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[2481, 63, 101, 1, 36, "Subsection"],
Cell[2585, 66, 531, 12, 39, "Output"]
}, Closed]],
Cell[CellGroupData[{
Cell[3153, 83, 110, 1, 51, "Section"],
Cell[3266, 86, 383, 6, 31, "Text"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3686, 97, 157, 2, 65, "Section"],
Cell[3846, 101, 2719, 74, 236, "Text",
 CellTags->"GenerateLieAlgebra::usage"],
Cell[CellGroupData[{
Cell[6590, 179, 100, 1, 35, "Subsubsection"],
Cell[6693, 182, 210, 4, 31, "Text"],
Cell[6906, 188, 154, 3, 31, "Text"]
}, Closed]],
Cell[CellGroupData[{
Cell[7097, 196, 97, 1, 37, "Subsection"],
Cell[7197, 199, 11859, 308, 748, "Text"],
Cell[CellGroupData[{
Cell[19081, 511, 103, 1, 35, "Subsubsection"],
Cell[19187, 514, 254, 6, 31, "Text",
 CellTags->"LieOperation::usage"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19478, 525, 98, 1, 35, "Subsubsection"],
Cell[19579, 528, 837, 24, 53, "Text",
 CellTags->"VSField::usage"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20453, 557, 112, 1, 35, "Subsubsection"],
Cell[20568, 560, 860, 28, 74, "Text",
 CellTags->"LinearIndependenceFun::usage"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21465, 593, 106, 1, 35, "Subsubsection"],
Cell[21574, 596, 308, 8, 32, "Text",
 CellTags->"StoppingDepth::usage"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21919, 609, 108, 1, 35, "Subsubsection"],
Cell[22030, 612, 683, 20, 54, "Text",
 CellTags->"StoppingDimension::usage"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22750, 637, 108, 1, 35, "Subsubsection"],
Cell[22861, 640, 1036, 27, 54, "Text",
 CellTags->"AssumeHamiltonians::usage"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23934, 672, 107, 1, 35, "Subsubsection"],
Cell[24044, 675, 545, 16, 53, "Text",
 CellTags->"GeneratorSymbols::usage"],
Cell[24592, 693, 314, 9, 32, "Text",
 CellTags->"G::usage"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24943, 707, 110, 1, 35, "Subsubsection"],
Cell[25056, 710, 262, 7, 32, "Text",
 CellTags->"NumericalLieAlgebra::usage"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[25367, 723, 151, 2, 37, "Subsection"],
Cell[25521, 727, 452, 12, 33, "Text"],
Cell[CellGroupData[{
Cell[25998, 743, 855, 18, 99, "Input"],
Cell[26856, 763, 334, 4, 23, "Print"],
Cell[27193, 769, 513, 13, 32, "Output"],
Cell[27709, 784, 1508, 39, 74, "Output"]
}, Open  ]]
}, Closed]],
Cell[CellGroupData[{
Cell[29266, 829, 145, 2, 37, "Subsection"],
Cell[29414, 833, 665, 15, 53, "Text"],
Cell[CellGroupData[{
Cell[30104, 852, 1119, 25, 77, "Input"],
Cell[31226, 879, 365, 5, 23, "Print"],
Cell[31594, 886, 2802, 89, 77, "Output"]
}, Open  ]],
Cell[34411, 978, 190, 4, 31, "Text"],
Cell[CellGroupData[{
Cell[34626, 986, 903, 23, 77, "Input"],
Cell[35532, 1011, 141, 2, 23, "Print"]
}, Open  ]],
Cell[CellGroupData[{
Cell[35710, 1018, 101, 1, 35, "Subsubsection"],
Cell[35814, 1021, 278, 5, 51, "Text"]
}, Open  ]]
}, Closed]]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)