
(* %%%%%%%%%%%%%%%%%%% NOTATION %%%%%%%%%%%%%%%%%% *)
(*                                                 *)
(*        savingP1[ a, b, a, spin][ block ]        *)
(*        savingP2[ a, b, a, spin][ block ]        *)
(*                                                 *)
(* %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  *)



(* %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% *)


savingP1[0, 2, 0, 2][1] = 1980*(-91728 - 45864*I - 5733*I^2 + 20896*T + 
      8450*I*T + 793*I^2*T - 2373*T^2 - 546*I*T^2 + 105*T^3)
 
savingP1[0, 2, 0, 2][2] = 3960*(-6534 - 5148*I - 1014*I^2 + 1728*T + 
      1560*I*T + 299*I^2*T - 504*T^2 - 182*I*T^2 + 35*T^3)
 
savingP1[0, 3, 0, 2][1] = (495*(-9871875 - 3159000*I - 252720*I^2 + 
       1194660*T + 312000*I*T + 19136*I^2*T - 76944*T^2 - 
       11648*I*T^2 + 2240*T^3))/8
 
savingP1[0, 3, 0, 2][2] = (495*(609375 + 195000*I + 15600*I^2 + 
       361132*T + 293696*I*T + 50752*I^2*T - 120624*T^2 - 
       34944*I*T^2 + 6720*T^3))/16
 
savingP1[0, 4, 0, 2][1] = 6600*(-227448 - 50544*I - 2808*I^2 + 
      18071*T + 3302*I*T + 143*I^2*T - 738*T^2 - 78*I*T^2 + 
      15*T^3)
 
savingP1[0, 4, 0, 2][2] = 2640*(-227448 - 50544*I - 2808*I^2 - 8785*T - 
      6058*I*T - 793*I^2*T + 14924*T^2 + 4784*I*T^2 + 
      416*I^2*T^2 - 1776*T^3 - 312*I*T^3 + 60*T^4)
 
savingP1[0, 4, 0, 2][3] = 660*(-44944 - 22048*I - 2704*I^2 + 17303*T + 
      9230*I*T + 1079*I^2*T - 3294*T^2 - 702*I*T^2 + 135*T^3)
 
savingP1[0, 5, 0, 2][1] = (165*(-596948625 - 97461000*I - 3978000*I^2 + 
       34506508*T + 4666688*I*T + 150592*I^2*T - 966960*T^2 - 
       74880*I*T^2 + 14400*T^3))/32
 
savingP1[0, 5, 0, 2][2] = (2475*(-128940903 - 21051576*I - 859248*I^2 - 
       19551888*T - 3740256*I*T - 194688*I^2*T + 4069728*T^2 + 
       783744*I*T^2 + 43264*I^2*T^2 - 272640*T^3 - 
       33280*I*T^3 + 6400*T^4))/128
 
savingP1[0, 5, 0, 2][3] = (165*(1591863 + 259896*I + 10608*I^2 + 
       2486476*T + 1082432*I*T + 110656*I^2*T - 431280*T^2 - 
       74880*I*T^2 + 14400*T^3))/16
 
savingP1[0, 6, 0, 2][1] = 210*(-26687232 - 3335904*I - 104247*I^2 + 
      1191672*T + 124254*I*T + 3107*I^2*T - 24255*T^2 - 
      1430*I*T^2 + 275*T^3)
 
savingP1[0, 6, 0, 2][2] = 120*(-47443968 - 5930496*I - 185328*I^2 - 
      12344778*T - 1627236*I*T - 55458*I^2*T + 1177856*T^2 + 
      153712*I*T^2 + 5681*I^2*T^2 - 48840*T^3 - 4290*I*T^3 + 
      825*T^4)
 
savingP1[0, 6, 0, 2][3] = 150*(-2965248 - 370656*I - 11583*I^2 - 
      709128*T - 184626*I*T - 13338*I^2*T + 426209*T^2 + 
      94900*I*T^2 + 5993*I^2*T^2 - 35970*T^3 - 4290*I*T^3 + 
      825*T^4)
 
savingP1[0, 6, 0, 2][4] = 240*(-196830 - 63180*I - 5070*I^2 + 89424*T + 
      28704*I*T + 2171*I^2*T - 10560*T^2 - 1430*I*T^2 + 
      275*T^3)
 
savingP1[0, 7, 0, 2][1] = (105*(-714670047 - 70584696*I - 1742832*I^2 + 
       25639124*T + 2125760*I*T + 42432*I^2*T - 395472*T^2 - 
       18304*I*T^2 + 3520*T^3))/8
 
savingP1[0, 7, 0, 2][2] = (105*(-9115689375 - 900315000*I - 
       22230000*I^2 - 4741623184*T - 477843808*I*T - 
       12215424*I^2*T + 280597088*T^2 + 26915200*I*T^2 + 
       708864*I^2*T^2 - 7825664*T^3 - 512512*I*T^3 + 98560*T^4))/
     128
 
savingP1[0, 7, 0, 2][3] = (15*(-10632540087 - 1050127416*I - 
       25929072*I^2 - 5412591760*T - 643388512*I*T - 
       22851712*I^2*T + 973786976*T^2 + 139308416*I*T^2 + 
       6166784*I^2*T^2 - 51980544*T^3 - 4612608*I*T^3 + 
       887040*T^4))/128
 
savingP1[0, 7, 0, 2][4] = (75*(1620567 + 160056*I + 3952*I^2 + 
       10918476*T + 2879552*I*T + 190528*I^2*T - 1132208*T^2 - 
       128128*I*T^2 + 24640*T^3))/16
 
savingP1[0, 8, 0, 2][1] = 2520*(-5850000 - 468000*I - 9360*I^2 + 
      173375*T + 11710*I*T + 191*I^2*T - 2094*T^2 - 
      78*I*T^2 + 15*T^3)
 
savingP1[0, 8, 0, 2][2] = 3360*T*(-2473497 - 198666*I - 4017*I^2 + 
      103356*T + 7680*I*T + 152*I^2*T - 2064*T^2 - 
      104*I*T^2 + 20*T^3)
 
savingP1[0, 8, 0, 2][3] = 840*T*(-2644668 - 223704*I - 5148*I^2 + 
      228329*T + 22690*I*T + 713*I^2*T - 8106*T^2 - 
      546*I*T^2 + 105*T^3)
 
savingP1[0, 8, 0, 2][4] = 240*T*(-1172925 - 169650*I - 7605*I^2 + 
      423020*T + 67840*I*T + 3128*I^2*T - 25872*T^2 - 
      2184*I*T^2 + 420*T^3)
 
savingP1[0, 8, 0, 2][5] = 1500*(-50544 - 11232*I - 624*I^2 + 25299*T + 
      5286*I*T + 275*I^2*T - 1974*T^2 - 182*I*T^2 + 35*T^3)



(* %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% *)


savingP2[0, 2, 0, 2][1] = -7920*(1742832 + 871416*I + 108927*I^2 - 
      339712*T - 139022*I*T - 13780*I^2*T + 90463*T^2 + 
      34736*I*T^2 + 3133*I^2*T^2 - 11508*T^3 - 2730*I*T^3 + 
      525*T^4)
 
savingP2[0, 2, 0, 2][2] = -15840*(-124146 - 97812*I - 19266*I^2 + 
      113886*T + 68484*I*T + 10205*I^2*T + 936*T^2 + 
      2938*I*T^2 + 871*I^2*T^2 - 2401*T^3 - 910*I*T^3 + 
      175*T^4)
 
savingP2[0, 3, 0, 2][1] = (-495*(1549884375 + 495963000*I + 
       39677040*I^2 - 130268880*T - 31493280*I*T - 
       1735552*I^2*T + 17148192*T^2 + 4021888*I*T^2 + 
       222976*I^2*T^2 - 1431808*T^3 - 232960*I*T^3 + 44800*T^4))/8
 
savingP2[0, 3, 0, 2][2] = (-495*(-95671875 - 30615000*I - 2449200*I^2 + 
       122250224*T + 65776672*I*T + 9436544*I^2*T - 545312*T^2 + 
       1988480*I*T^2 + 802048*I^2*T^2 - 2091264*T^3 - 
       698880*I*T^3 + 134400*T^4))/16
 
savingP2[0, 4, 0, 2][1] = -26400*(14556672 + 3234816*I + 179712*I^2 - 
      764128*T - 125632*I*T - 4576*I^2*T + 56435*T^2 + 
      8606*I*T^2 + 299*I^2*T^2 - 3354*T^3 - 390*I*T^3 + 
      75*T^4)
 
savingP2[0, 4, 0, 2][2] = -10560*(14556672 + 3234816*I + 179712*I^2 - 
      663328*T - 320320*I*T - 42016*I^2*T - 74891*T^2 + 
      67522*I*T^2 + 12493*I^2*T^2 + 50708*T^3 + 15392*I*T^3 + 
      1976*I^2*T^3 - 7536*T^4 - 1560*I*T^4 + 300*T^5)
 
savingP2[0, 4, 0, 2][3] = -2640*(-2876416 - 1411072*I - 173056*I^2 + 
      2052256*T + 779584*I*T + 76960*I^2*T - 29869*T^2 + 
      5798*I*T^2 + 3731*I^2*T^2 - 13446*T^3 - 3510*I*T^3 + 
      675*T^4)
 
savingP2[0, 5, 0, 2][1] = (-165*(222661837125 + 36352953000*I + 
       1483794000*I^2 - 8445067216*T - 1023447776*I*T - 
       27424384*I^2*T + 385861600*T^2 + 40327040*I*T^2 + 
       881920*I^2*T^2 - 17337600*T^3 - 1497600*I*T^3 + 
       288000*T^4))/32
 
savingP2[0, 5, 0, 2][2] = (-2475*(48094956819 + 7852237848*I + 
       320499504*I^2 + 7151760324*T + 1117853568*I*T + 
       36264384*I^2*T - 775183392*T^2 - 65544960*I*T^2 + 
       2119936*I^2*T^2 + 62666880*T^3 + 11048960*I*T^3 + 
       865280*I^2*T^3 - 4563200*T^4 - 665600*I*T^4 + 128000*T^5))/
     128
 
savingP2[0, 5, 0, 2][3] = (-165*(-593764899 - 96941208*I - 3956784*I^2 + 
       1240128272*T + 452886304*I*T + 43202432*I^2*T - 
       20335520*T^2 + 590720*I*T^2 + 1680640*I^2*T^2 - 
       6624000*T^3 - 1497600*I*T^3 + 288000*T^4))/16
 
savingP2[0, 6, 0, 2][1] = -840*(3389278464 + 423659808*I + 13239369*I^2 - 
      99962424*T - 9378018*I*T - 196144*I^2*T + 3038265*T^2 + 
      226720*I*T^2 + 3055*I^2*T^2 - 107800*T^3 - 7150*I*T^3 + 
      1375*T^4)
 
savingP2[0, 6, 0, 2][2] = -480*(6025383936 + 753172992*I + 23536656*I^2 + 
      1607005386*T + 201309732*I*T + 6053346*I^2*T - 
      97143182*T^2 - 8689564*I*T^2 - 79157*I^2*T^2 + 
      4452920*T^3 + 531310*I*T^3 + 27365*I^2*T^3 - 203775*T^4 - 
      21450*I*T^4 + 4125*T^5)
 
savingP2[0, 6, 0, 2][3] = -600*(376586496 + 47073312*I + 1471041*I^2 - 
      25070904*T - 14628978*I*T - 1396629*I^2*T + 8271049*T^2 + 
      6186362*I*T^2 + 564967*I^2*T^2 + 798035*T^3 + 
      156130*I*T^3 + 25805*I^2*T^3 - 139425*T^4 - 21450*I*T^4 + 
      4125*T^5)
 
savingP2[0, 6, 0, 2][4] = -960*(-24997410 - 8023860*I - 643890*I^2 + 
      14479398*T + 3961308*I*T + 279617*I^2*T - 220080*T^2 - 
      7670*I*T^2 + 7735*I^2*T^2 - 39325*T^3 - 7150*I*T^3 + 
      1375*T^4)
 
savingP2[0, 7, 0, 2][1] = (-105*(472396901067 + 46656484056*I + 
       1152011952*I^2 - 11315871952*T - 849520672*I*T - 
       14368640*I^2*T + 241693088*T^2 + 13293696*I*T^2 + 
       103168*I^2*T^2 - 6992128*T^3 - 366080*I*T^3 + 70400*T^4))/8
 
savingP2[0, 7, 0, 2][2] = (-105*(6025470676875 + 595108215000*I + 
       14694030000*I^2 + 3171510885476*T + 313779619712*I*T + 
       7641190336*I^2*T - 129909629216*T^2 - 9914854144*I*T^2 - 
       130973440*I^2*T^2 + 4023520384*T^3 + 348854272*I*T^3 + 
       12473344*I^2*T^3 - 130828544*T^4 - 10250240*I*T^4 + 
       1971200*T^5))/128
 
savingP2[0, 7, 0, 2][3] = (-15*(7028108997507 + 694134221976*I + 
       17139116592*I^2 + 2845768075652*T + 217125234560*I*T + 
       954060224*I^2*T - 198512377376*T^2 + 11620740352*I*T^2 + 
       2123999488*I^2*T^2 + 12279582848*T^3 + 1360566272*I*T^3 + 
       116519936*I^2*T^3 - 808448256*T^4 - 92252160*I*T^4 + 
       17740800*T^5))/128
 
savingP2[0, 7, 0, 2][4] = (-75*(-1071194787 - 105797016*I - 2612272*I^2 + 
       7060598928*T + 1886708512*I*T + 125497216*I^2*T - 
       107163040*T^2 - 7236736*I*T^2 + 2852096*I^2*T^2 - 
       16222976*T^3 - 2562560*I*T^3 + 492800*T^4))/16
 
savingP2[0, 8, 0, 2][1] = -10080*(1216800000 + 97344000*I + 1946880*I^2 - 
      24362000*T - 1499680*I*T - 21008*I^2*T + 381563*T^2 + 
      15862*I*T^2 + 59*I^2*T^2 - 9222*T^3 - 390*I*T^3 + 
      75*T^4)
 
savingP2[0, 8, 0, 2][2] = -13440*T*(512187624 + 40811472*I + 
      807144*I^2 - 15521571*T - 986622*I*T - 12571*I^2*T + 
      349188*T^2 + 22960*I*T^2 + 592*I^2*T^2 - 8656*T^3 - 
      520*I*T^3 + 100*T^4)
 
savingP2[0, 8, 0, 2][3] = -3360*T*(514609056 + 38645568*I + 
      632736*I^2 - 25182404*T - 725608*I*T + 45916*I^2*T + 
      835757*T^2 + 66730*I*T^2 + 3533*I^2*T^2 - 31794*T^3 - 
      2730*I*T^3 + 525*T^4)
 
savingP2[0, 8, 0, 2][4] = -960*T*(22206600 - 13993200*I - 1155960*I^2 + 
      22727185*T + 8836490*I*T + 538681*I^2*T + 310484*T^2 + 
      40816*I*T^2 + 12752*I^2*T^2 - 94416*T^3 - 10920*I*T^3 + 
      2100*T^4)
 
savingP2[0, 8, 0, 2][5] = -6000*(-10513152 - 2336256*I - 129792*I^2 + 
      5161104*T + 1077024*I*T + 55952*I^2*T - 67233*T^2 - 
      4226*I*T^2 + 991*I^2*T^2 - 6958*T^3 - 910*I*T^3 + 
      175*T^4)
