(* ::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.                                                         *)
(************************************************************************)



ca[0,__] :=0;
ca /: ca[f_,is__] + ca[g_, is__] := ca[f+g, is];
ca /: -ca[f_, is__] := ca[Expand[-f], is];
B[a[f_,j_,k_], c[l_]] /; DistinctQ[j,k,l]:= 0;
B[a[f_,j_,k_], c[j_]] := \[Gamma][-f,j,k];
B[a[f_,j_,k_], c[k_]] := \[Gamma][f,j,k];
B[x_, ca[f_, i_,j_,k_]] := Plus[
B[x, c[i]] /. \[Gamma][g_,l_,m_] :> \[Gamma]a[f g, l,m,j,k],
B[x, a[f,j,k]] /. {a[g_,l_,m_] :> ca[g, i,l,m], _\[Gamma] -> 0, _\[Gamma]a -> 0}
];
\[Delta]aa[f_, i_,j_,k_,l_] /; OrderedQ[{l,j}]\[And]DistinctQ[i,j,k,l] := \[Delta]aa[f, k,l,i,j];
\[Delta]aa[f_, i_,j_,i_,l_] /; OrderedQ[{l,j}]\[And]DistinctQ[j,l]:= \[Delta]aa[f,i,l,i,j];
(* missing a 4T! *)
E\[Gamma] = {
\[Gamma][f_,j_,k_,l_] :> ca[f, l,j,k]+ca[Expand[-f], k,j,l],
\[Gamma]a[f_, i_,j_,k_,l_] :> \[Delta]aa[f, i,j,k,l] + ca[Expand[-Subscript[b, i]f], j,k,l]
};
Check3[y_\[Gamma]] := {t1=B[a[f,j,k], y] /. E\[Gamma], t2=B[a[f,j,k], y /. E\[Gamma]] /. E\[Gamma], t1==t2 // Simplify}



