(* ::Package:: *)

(************************************************************************)\012(* This file was generated automatically by the Mathematica front end.  *)\012(* It contains Initialization cells from a Notebook file, which         *)\012(* typically will have the same name as this file except ending in      *)\012(* ".nb" instead of ".m".                                               *)\012(*                                                                      *)\012(* This file is intended to be loaded into the Mathematica kernel using *)\012(* the package loading commands Get or Needs.  Doing so is equivalent   *)\012(* to using the Evaluate Initialization Cells menu command in the front *)\012(* end.                                                                 *)\012(*                                                                      *)\012(* DO NOT EDIT THIS FILE.  This entire file is regenerated              *)\012(* automatically each time the parent Notebook file is saved in the     *)\012(* Mathematica front end.  Any changes you make to this file will be    *)\012(* overwritten.                                                         *)\012(************************************************************************)\012\012

(* ::Input::Initialization:: *)
SetDirectory["C:\\drorbn\\AcademicPensieve\\Talks\\Bonn-2505"];


(* ::Input::Initialization:: *)
Once[<<KnotTheory`; <<Rot.m];


(* ::Input::Initialization:: *)
CF[\[Omega]_. \[ScriptCapitalE]_\[DoubleStruckCapitalE]]:=CF[\[Omega]]CF/@\[ScriptCapitalE];
CF[\[ScriptCapitalE]_List]:=CF/@\[ScriptCapitalE];
CF[\[ScriptCapitalE]_] := Module[{vs,ps,c},
vs=Cases[\[ScriptCapitalE],Subscript[(x|p|\[Xi]|\[Pi]|g), __],\[Infinity]]\[Union]{\[Epsilon]};
Total[CoefficientRules[Expand[\[ScriptCapitalE]], vs] /. (ps_->c_) :> Factor[c](Times@@(vs^ps))] ];


(* ::Input::Initialization:: *)
\[DoubleStruckCapitalE]/:\[DoubleStruckCapitalE][A_]\[DoubleStruckCapitalE][B_]:=\[DoubleStruckCapitalE][A+B];


(* ::Input::Initialization:: *)
$\[Pi]=Identity; (* The Wisdom Projection *)


(* ::Input::Initialization:: *)
Unprotect[Integrate];
\[Integral]\[Omega]_. \[DoubleStruckCapitalE][L_]\[DifferentialD](vs_List):=Module[{n,L0,Q,\[CapitalDelta],G,Z0,Z,\[Lambda],DZ,DDZ,FZ,a,b},
n=Length@vs; L0=L/.\[Epsilon]->0;
Q=Table[-D[L0,vs[[a]],vs[[b]]]/.Thread[vs->0]/.Subscript[(p|x), __]->0,{a,n},{b,n}];
If[(\[CapitalDelta]=Det[Q])==0,Return@"Degenerate Q!"];
Z=Z0=CF@$\[Pi][L+vs . Q . vs/2]; G=Inverse[Q];
FixedPoint[(DZ=Table[\!\(
\*SubscriptBox[\(\[PartialD]\), \(v\)]Z\),{v,vs}];
DDZ=Table[\!\(
\*SubscriptBox[\(\[PartialD]\), \(u\)]DZ\),{u,vs}];
FZ=Sum[G[[a,b]](DDZ[[a,b]]+DZ[[a]]DZ[[b]]),{a,n},{b,n}]/2;
Z=CF[Z0+\!\(
\*SubsuperscriptBox[\(\[Integral]\), \(0\), \(\[Lambda]\)]\($\[Pi][FZ]\ \[DifferentialD]\[Lambda]\)\)])&,Z];
PowerExpand@Factor[\[Omega] \[CapitalDelta]^(-1/2)]\[DoubleStruckCapitalE][CF[Z/.\[Lambda]->1/.Thread[vs->0]]]];
Protect[Integrate];


(* ::Input::Initialization:: *)
\[ScriptCapitalL][Subscript[X, i_,j_][s_]]:=T^(s/2) \[DoubleStruckCapitalE][
Subscript[x, i](Subscript[p, i+1]-Subscript[p, i])+Subscript[x, j](Subscript[p, j+1]-Subscript[p, j])+(T^s-1)Subscript[x, i](Subscript[p, i+1]-Subscript[p, j+1])+(\[Epsilon] s/2)*(Subscript[x, i](Subscript[p, i]-Subscript[p, j])((T^s-1)Subscript[x, i] Subscript[p, j]+2(1-Subscript[x, j] Subscript[p, j]))-1)]
\[ScriptCapitalL][Subscript[C, i_][\[CurlyPhi]_]]:=T^(\[CurlyPhi]/2) \[DoubleStruckCapitalE][Subscript[x, i](Subscript[p, i+1]-Subscript[p, i])+\[Epsilon] \[CurlyPhi](1/2-Subscript[x, i] Subscript[p, i])]
\[ScriptCapitalL][K_]:=CF[\[ScriptCapitalL]/@Features[K][[2]]]
vs[K_]:=Join@@Table[{Subscript[p, i],Subscript[x, i]},{i,Features[K][[1]]}]



