(* ::Package:: *)

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



(* ::Input::Initialization:: *)
SetDirectory["C:\\drorbn\\AcademicPensieve\\Projects\\APAI"];


(* ::Input::Initialization:: *)
Once[<<KnotTheory`; <<Rot.m];
Once[Get@"../Profile/Profile.m"];
BeginProfile[];


(* ::Input::Initialization:: *)
Subscript[R, 1][s_,i_,j_]:=s(Subscript[g, j\[InvisibleComma]i](Subscript[g, SuperPlus[j],j]+Subscript[g, j,SuperPlus[j]]-Subscript[g, i\[InvisibleComma]j])-Subscript[g, i\[InvisibleComma]i](Subscript[g, j,SuperPlus[j]]-1)-1/2);
\[Rho][K_]:=Subscript[PP, \[Rho]]@Module[{Cs,\[CurlyPhi],n,A,s,i,j,k,\[CapitalDelta],G,\[Rho]1,res},
Subscript[PP, Init]@({Cs,\[CurlyPhi]}=Rot[K]; n=Length[Cs];
A=IdentityMatrix[2n+1];
Cases[Cs,{s_,i_,j_}:>(A[[{i,j},{i+1,j+1}]]+=({
 {-T^s, T^s-1},
 {0, -1}
}))]);
\!\(\*SubscriptBox[\(PP\), \("\<\[CapitalDelta]\>"\)]\)@(\[CapitalDelta]=T^((-Total[\[CurlyPhi]]-Total[Cs[[All,1]]])/2) Det[A]);
\!\(\*SubscriptBox[\(PP\), \("\<G\>"\)]\)@(G=Inverse[A]);
\!\(\*SubscriptBox[\(PP\), \("\<\[Rho]1\>"\)]\)@(\[Rho]1=\!\(
\*SubsuperscriptBox[\(\[Sum]\), \(k = 1\), \(n\)]\(
\*SubscriptBox[\(R\), \(1\)] @@ Cs[[k]]\)\) - \!\(
\*SubsuperscriptBox[\(\[Sum]\), \(k = 1\), \(2  n\)]\(\[CurlyPhi][[k]] \((
\*SubscriptBox[\(g\), \(k\[InvisibleComma]k\)] - 1/2)\)\)\);
res=Factor@{\[CapitalDelta],\[CapitalDelta]^2 \[Rho]1/.SuperPlus[\[Alpha]_]:>\[Alpha]+1/.Subscript[g, \[Alpha]_,\[Beta]_]:>G[[\[Alpha],\[Beta]]]})
];


(* ::Input::Initialization:: *)
GST48=EPD[Subscript[X, 14,1],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 2,29],Subscript[X, 3,40],Subscript[X, 43,4],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 26,5],Subscript[X, 6,95],Subscript[X, 96,7],Subscript[X, 13,8],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 9,28],Subscript[X, 10,41],Subscript[X, 42,11],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 27,12],Subscript[X, 30,15],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 16,61],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 17,72],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 18,83],Subscript[X, 19,34],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 89,20],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 21,92],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 79,22],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 68,23],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 57,24],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 25,56],Subscript[X, 62,31],Subscript[X, 73,32],Subscript[X, 84,33],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 50,35],Subscript[X, 36,81],Subscript[X, 37,70],Subscript[X, 38,59],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 39,54],Subscript[X, 44,55],Subscript[X, 58,45],Subscript[X, 69,46],Subscript[X, 80,47],Subscript[X, 48,91],Subscript[X, 90,49],Subscript[X, 51,82],Subscript[X, 52,71],Subscript[X, 53,60],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 63,74],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 64,85],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 76,65],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 87,66],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 67,94],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 75,86],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 88,77],Subscript[\!\(\*OverscriptBox[\(X\), \(_\)]\), 78,93]];



