16x32.nb 51 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179118011811182118311841185118611871188118911901191119211931194119511961197119811991200120112021203120412051206120712081209121012111212121312141215121612171218121912201221122212231224122512261227122812291230123112321233123412351236123712381239124012411242124312441245124612471248124912501251125212531254125512561257125812591260126112621263126412651266126712681269127012711272127312741275127612771278127912801281128212831284128512861287128812891290129112921293129412951296129712981299130013011302130313041305130613071308130913101311131213131314131513161317131813191320132113221323132413251326132713281329133013311332133313341335133613371338133913401341134213431344134513461347134813491350135113521353135413551356135713581359136013611362136313641365136613671368136913701371137213731374137513761377137813791380138113821383138413851386138713881389139013911392139313941395139613971398139914001401140214031404140514061407140814091410141114121413141414151416141714181419142014211422142314241425142614271428142914301431143214331434143514361437143814391440144114421443144414451446144714481449145014511452145314541455145614571458145914601461146214631464146514661467146814691470147114721473147414751476147714781479148014811482148314841485148614871488148914901491149214931494149514961497149814991500150115021503150415051506150715081509151015111512151315141515151615171518151915201521152215231524152515261527
  1. (* Content-type: application/mathematica *)
  2. (*** Wolfram Notebook File ***)
  3. (* http://www.wolfram.com/nb *)
  4. (* CreatedBy='Mathematica 7.0' *)
  5. (*CacheID: 234*)
  6. (* Internal cache information:
  7. NotebookFileLineBreakTest
  8. NotebookFileLineBreakTest
  9. NotebookDataPosition[ 145, 7]
  10. NotebookDataLength[ 51606, 1517]
  11. NotebookOptionsPosition[ 50273, 1475]
  12. NotebookOutlinePosition[ 50610, 1490]
  13. CellTagsIndexPosition[ 50567, 1487]
  14. WindowFrame->Normal*)
  15. (* Beginning of Notebook Content *)
  16. Notebook[{
  17. Cell[BoxData[
  18. RowBox[{"a", ":=",
  19. RowBox[{"ArrayFlatten", "[",
  20. RowBox[{"{",
  21. RowBox[{
  22. RowBox[{"{",
  23. RowBox[{
  24. RowBox[{"IdentityMatrix", "[", "8", "]"}], ",",
  25. RowBox[{"Reverse", "[",
  26. RowBox[{"IdentityMatrix", "[", "8", "]"}], "]"}]}], "}"}], ",",
  27. RowBox[{"{",
  28. RowBox[{
  29. RowBox[{"Reverse", "[",
  30. RowBox[{"IdentityMatrix", "[", "8", "]"}], "]"}], ",",
  31. RowBox[{
  32. RowBox[{"-", "1"}], "*",
  33. RowBox[{"IdentityMatrix", "[", "8", "]"}]}]}], "}"}]}], "}"}],
  34. "]"}]}]], "Input",
  35. CellChangeTimes->{{3.545829278882098*^9, 3.545829319074397*^9}}],
  36. Cell[BoxData[
  37. RowBox[{"v", ":=",
  38. RowBox[{
  39. RowBox[{"{",
  40. RowBox[{
  41. RowBox[{"{",
  42. RowBox[{
  43. "1", ",", "q0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  44. "0"}], "}"}], ",",
  45. RowBox[{"{",
  46. RowBox[{
  47. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  48. "0"}], "}"}], ",",
  49. RowBox[{"{",
  50. RowBox[{
  51. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  52. "0"}], "}"}], ",",
  53. RowBox[{"{",
  54. RowBox[{
  55. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  56. "0"}], "}"}], ",",
  57. RowBox[{"{",
  58. RowBox[{
  59. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  60. "0"}], "}"}], ",",
  61. RowBox[{"{",
  62. RowBox[{
  63. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  64. "0"}], "}"}], ",",
  65. RowBox[{"{",
  66. RowBox[{
  67. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  68. "0"}], "}"}], ",",
  69. RowBox[{"{",
  70. RowBox[{
  71. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  72. "1"}], "}"}]}], "}"}], ".",
  73. RowBox[{"{",
  74. RowBox[{
  75. RowBox[{"{",
  76. RowBox[{
  77. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  78. "0"}], "}"}], ",",
  79. RowBox[{"{",
  80. RowBox[{
  81. "0", ",", "1", ",", "q1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  82. "0"}], "}"}], ",",
  83. RowBox[{"{",
  84. RowBox[{
  85. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  86. "0"}], "}"}], ",",
  87. RowBox[{"{",
  88. RowBox[{
  89. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  90. "0"}], "}"}], ",",
  91. RowBox[{"{",
  92. RowBox[{
  93. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  94. "0"}], "}"}], ",",
  95. RowBox[{"{",
  96. RowBox[{
  97. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  98. "0"}], "}"}], ",",
  99. RowBox[{"{",
  100. RowBox[{
  101. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  102. "0"}], "}"}], ",",
  103. RowBox[{"{",
  104. RowBox[{
  105. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  106. "1"}], "}"}]}], "}"}], ".",
  107. RowBox[{"{",
  108. RowBox[{
  109. RowBox[{"{",
  110. RowBox[{
  111. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  112. "0"}], "}"}], ",",
  113. RowBox[{"{",
  114. RowBox[{
  115. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  116. "0"}], "}"}], ",",
  117. RowBox[{"{",
  118. RowBox[{
  119. "0", ",", "0", ",", "1", ",", "q2", ",", "0", ",", "0", ",", "0", ",",
  120. "0"}], "}"}], ",",
  121. RowBox[{"{",
  122. RowBox[{
  123. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  124. "0"}], "}"}], ",",
  125. RowBox[{"{",
  126. RowBox[{
  127. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  128. "0"}], "}"}], ",",
  129. RowBox[{"{",
  130. RowBox[{
  131. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  132. "0"}], "}"}], ",",
  133. RowBox[{"{",
  134. RowBox[{
  135. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  136. "0"}], "}"}], ",",
  137. RowBox[{"{",
  138. RowBox[{
  139. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  140. "1"}], "}"}]}], "}"}], ".",
  141. RowBox[{"{",
  142. RowBox[{
  143. RowBox[{"{",
  144. RowBox[{
  145. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  146. "0"}], "}"}], ",",
  147. RowBox[{"{",
  148. RowBox[{
  149. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  150. "0"}], "}"}], ",",
  151. RowBox[{"{",
  152. RowBox[{
  153. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  154. "0"}], "}"}], ",",
  155. RowBox[{"{",
  156. RowBox[{
  157. "0", ",", "0", ",", "0", ",", "1", ",", "q3", ",", "0", ",", "0", ",",
  158. "0"}], "}"}], ",",
  159. RowBox[{"{",
  160. RowBox[{
  161. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  162. "0"}], "}"}], ",",
  163. RowBox[{"{",
  164. RowBox[{
  165. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  166. "0"}], "}"}], ",",
  167. RowBox[{"{",
  168. RowBox[{
  169. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  170. "0"}], "}"}], ",",
  171. RowBox[{"{",
  172. RowBox[{
  173. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  174. "1"}], "}"}]}], "}"}], ".",
  175. RowBox[{"{",
  176. RowBox[{
  177. RowBox[{"{",
  178. RowBox[{
  179. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  180. "0"}], "}"}], ",",
  181. RowBox[{"{",
  182. RowBox[{
  183. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  184. "0"}], "}"}], ",",
  185. RowBox[{"{",
  186. RowBox[{
  187. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  188. "0"}], "}"}], ",",
  189. RowBox[{"{",
  190. RowBox[{
  191. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  192. "0"}], "}"}], ",",
  193. RowBox[{"{",
  194. RowBox[{
  195. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "q4", ",", "0", ",",
  196. "0"}], "}"}], ",",
  197. RowBox[{"{",
  198. RowBox[{
  199. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  200. "0"}], "}"}], ",",
  201. RowBox[{"{",
  202. RowBox[{
  203. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  204. "0"}], "}"}], ",",
  205. RowBox[{"{",
  206. RowBox[{
  207. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  208. "1"}], "}"}]}], "}"}], ".",
  209. RowBox[{"{",
  210. RowBox[{
  211. RowBox[{"{",
  212. RowBox[{
  213. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  214. "0"}], "}"}], ",",
  215. RowBox[{"{",
  216. RowBox[{
  217. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  218. "0"}], "}"}], ",",
  219. RowBox[{"{",
  220. RowBox[{
  221. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  222. "0"}], "}"}], ",",
  223. RowBox[{"{",
  224. RowBox[{
  225. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  226. "0"}], "}"}], ",",
  227. RowBox[{"{",
  228. RowBox[{
  229. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  230. "0"}], "}"}], ",",
  231. RowBox[{"{",
  232. RowBox[{
  233. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "q5", ",",
  234. "0"}], "}"}], ",",
  235. RowBox[{"{",
  236. RowBox[{
  237. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  238. "0"}], "}"}], ",",
  239. RowBox[{"{",
  240. RowBox[{
  241. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  242. "1"}], "}"}]}], "}"}], ".",
  243. RowBox[{"{",
  244. RowBox[{
  245. RowBox[{"{",
  246. RowBox[{
  247. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  248. "0"}], "}"}], ",",
  249. RowBox[{"{",
  250. RowBox[{
  251. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  252. "0"}], "}"}], ",",
  253. RowBox[{"{",
  254. RowBox[{
  255. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  256. "0"}], "}"}], ",",
  257. RowBox[{"{",
  258. RowBox[{
  259. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  260. "0"}], "}"}], ",",
  261. RowBox[{"{",
  262. RowBox[{
  263. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  264. "0"}], "}"}], ",",
  265. RowBox[{"{",
  266. RowBox[{
  267. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  268. "0"}], "}"}], ",",
  269. RowBox[{"{",
  270. RowBox[{
  271. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  272. "q6"}], "}"}], ",",
  273. RowBox[{"{",
  274. RowBox[{
  275. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  276. "1"}], "}"}]}], "}"}], ".",
  277. RowBox[{"{",
  278. RowBox[{
  279. RowBox[{"{",
  280. RowBox[{
  281. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  282. "0"}], "}"}], ",",
  283. RowBox[{"{",
  284. RowBox[{
  285. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  286. "0"}], "}"}], ",",
  287. RowBox[{"{",
  288. RowBox[{
  289. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  290. "0"}], "}"}], ",",
  291. RowBox[{"{",
  292. RowBox[{
  293. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  294. "0"}], "}"}], ",",
  295. RowBox[{"{",
  296. RowBox[{
  297. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  298. "0"}], "}"}], ",",
  299. RowBox[{"{",
  300. RowBox[{
  301. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  302. "0"}], "}"}], ",",
  303. RowBox[{"{",
  304. RowBox[{
  305. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  306. "0"}], "}"}], ",",
  307. RowBox[{"{",
  308. RowBox[{
  309. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "p6", ",",
  310. "1"}], "}"}]}], "}"}], ".",
  311. RowBox[{"{",
  312. RowBox[{
  313. RowBox[{"{",
  314. RowBox[{
  315. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  316. "0"}], "}"}], ",",
  317. RowBox[{"{",
  318. RowBox[{
  319. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  320. "0"}], "}"}], ",",
  321. RowBox[{"{",
  322. RowBox[{
  323. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  324. "0"}], "}"}], ",",
  325. RowBox[{"{",
  326. RowBox[{
  327. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  328. "0"}], "}"}], ",",
  329. RowBox[{"{",
  330. RowBox[{
  331. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  332. "0"}], "}"}], ",",
  333. RowBox[{"{",
  334. RowBox[{
  335. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  336. "0"}], "}"}], ",",
  337. RowBox[{"{",
  338. RowBox[{
  339. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "p5", ",", "1", ",",
  340. "0"}], "}"}], ",",
  341. RowBox[{"{",
  342. RowBox[{
  343. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  344. "1"}], "}"}]}], "}"}], ".",
  345. RowBox[{"{",
  346. RowBox[{
  347. RowBox[{"{",
  348. RowBox[{
  349. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  350. "0"}], "}"}], ",",
  351. RowBox[{"{",
  352. RowBox[{
  353. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  354. "0"}], "}"}], ",",
  355. RowBox[{"{",
  356. RowBox[{
  357. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  358. "0"}], "}"}], ",",
  359. RowBox[{"{",
  360. RowBox[{
  361. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  362. "0"}], "}"}], ",",
  363. RowBox[{"{",
  364. RowBox[{
  365. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  366. "0"}], "}"}], ",",
  367. RowBox[{"{",
  368. RowBox[{
  369. "0", ",", "0", ",", "0", ",", "0", ",", "p4", ",", "1", ",", "0", ",",
  370. "0"}], "}"}], ",",
  371. RowBox[{"{",
  372. RowBox[{
  373. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  374. "0"}], "}"}], ",",
  375. RowBox[{"{",
  376. RowBox[{
  377. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  378. "1"}], "}"}]}], "}"}], ".",
  379. RowBox[{"{",
  380. RowBox[{
  381. RowBox[{"{",
  382. RowBox[{
  383. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  384. "0"}], "}"}], ",",
  385. RowBox[{"{",
  386. RowBox[{
  387. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  388. "0"}], "}"}], ",",
  389. RowBox[{"{",
  390. RowBox[{
  391. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  392. "0"}], "}"}], ",",
  393. RowBox[{"{",
  394. RowBox[{
  395. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  396. "0"}], "}"}], ",",
  397. RowBox[{"{",
  398. RowBox[{
  399. "0", ",", "0", ",", "0", ",", "p3", ",", "1", ",", "0", ",", "0", ",",
  400. "0"}], "}"}], ",",
  401. RowBox[{"{",
  402. RowBox[{
  403. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  404. "0"}], "}"}], ",",
  405. RowBox[{"{",
  406. RowBox[{
  407. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  408. "0"}], "}"}], ",",
  409. RowBox[{"{",
  410. RowBox[{
  411. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  412. "1"}], "}"}]}], "}"}], ".",
  413. RowBox[{"{",
  414. RowBox[{
  415. RowBox[{"{",
  416. RowBox[{
  417. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  418. "0"}], "}"}], ",",
  419. RowBox[{"{",
  420. RowBox[{
  421. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  422. "0"}], "}"}], ",",
  423. RowBox[{"{",
  424. RowBox[{
  425. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  426. "0"}], "}"}], ",",
  427. RowBox[{"{",
  428. RowBox[{
  429. "0", ",", "0", ",", "p2", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  430. "0"}], "}"}], ",",
  431. RowBox[{"{",
  432. RowBox[{
  433. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  434. "0"}], "}"}], ",",
  435. RowBox[{"{",
  436. RowBox[{
  437. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  438. "0"}], "}"}], ",",
  439. RowBox[{"{",
  440. RowBox[{
  441. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  442. "0"}], "}"}], ",",
  443. RowBox[{"{",
  444. RowBox[{
  445. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  446. "1"}], "}"}]}], "}"}], ".",
  447. RowBox[{"{",
  448. RowBox[{
  449. RowBox[{"{",
  450. RowBox[{
  451. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  452. "0"}], "}"}], ",",
  453. RowBox[{"{",
  454. RowBox[{
  455. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  456. "0"}], "}"}], ",",
  457. RowBox[{"{",
  458. RowBox[{
  459. "0", ",", "p1", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  460. "0"}], "}"}], ",",
  461. RowBox[{"{",
  462. RowBox[{
  463. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  464. "0"}], "}"}], ",",
  465. RowBox[{"{",
  466. RowBox[{
  467. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  468. "0"}], "}"}], ",",
  469. RowBox[{"{",
  470. RowBox[{
  471. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  472. "0"}], "}"}], ",",
  473. RowBox[{"{",
  474. RowBox[{
  475. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  476. "0"}], "}"}], ",",
  477. RowBox[{"{",
  478. RowBox[{
  479. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  480. "1"}], "}"}]}], "}"}], ".",
  481. RowBox[{"{",
  482. RowBox[{
  483. RowBox[{"{",
  484. RowBox[{
  485. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  486. "0"}], "}"}], ",",
  487. RowBox[{"{",
  488. RowBox[{
  489. "p0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  490. "0"}], "}"}], ",",
  491. RowBox[{"{",
  492. RowBox[{
  493. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  494. "0"}], "}"}], ",",
  495. RowBox[{"{",
  496. RowBox[{
  497. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  498. "0"}], "}"}], ",",
  499. RowBox[{"{",
  500. RowBox[{
  501. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  502. "0"}], "}"}], ",",
  503. RowBox[{"{",
  504. RowBox[{
  505. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  506. "0"}], "}"}], ",",
  507. RowBox[{"{",
  508. RowBox[{
  509. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  510. "0"}], "}"}], ",",
  511. RowBox[{"{",
  512. RowBox[{
  513. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  514. "1"}], "}"}]}], "}"}], ".",
  515. RowBox[{"{",
  516. RowBox[{
  517. RowBox[{"{",
  518. RowBox[{
  519. "s0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  520. "0"}], "}"}], ",",
  521. RowBox[{"{",
  522. RowBox[{
  523. "0", ",", "s1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  524. "0"}], "}"}], ",",
  525. RowBox[{"{",
  526. RowBox[{
  527. "0", ",", "0", ",", "s2", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  528. "0"}], "}"}], ",",
  529. RowBox[{"{",
  530. RowBox[{
  531. "0", ",", "0", ",", "0", ",", "s3", ",", "0", ",", "0", ",", "0", ",",
  532. "0"}], "}"}], ",",
  533. RowBox[{"{",
  534. RowBox[{
  535. "0", ",", "0", ",", "0", ",", "0", ",", "s4", ",", "0", ",", "0", ",",
  536. "0"}], "}"}], ",",
  537. RowBox[{"{",
  538. RowBox[{
  539. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "s5", ",", "0", ",",
  540. "0"}], "}"}], ",",
  541. RowBox[{"{",
  542. RowBox[{
  543. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "s6", ",",
  544. "0"}], "}"}], ",",
  545. RowBox[{"{",
  546. RowBox[{
  547. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  548. "s7"}], "}"}]}], "}"}]}]}]], "Input",
  549. CellChangeTimes->{{3.5458293327981825`*^9, 3.5458293359303617`*^9}, {
  550. 3.5458294146008615`*^9, 3.5458296302641964`*^9}, {3.545829660343917*^9,
  551. 3.545829713685968*^9}, {3.5458297688811245`*^9, 3.545829986622579*^9}, {
  552. 3.545830471926337*^9, 3.545830493855591*^9}}],
  553. Cell[BoxData[
  554. RowBox[{"b", ":=",
  555. RowBox[{"ArrayFlatten", "[",
  556. RowBox[{"{",
  557. RowBox[{
  558. RowBox[{"{",
  559. RowBox[{
  560. RowBox[{"IdentityMatrix", "[", "8", "]"}], ",", "0"}], "}"}], ",",
  561. RowBox[{"{",
  562. RowBox[{"0", ",", "v"}], "}"}]}], "}"}], "]"}]}]], "Input",
  563. CellChangeTimes->{{3.545830273509988*^9, 3.545830273802005*^9}, {
  564. 3.5458310447871027`*^9, 3.5458310471452374`*^9}, {3.5458310910837507`*^9,
  565. 3.545831121041464*^9}}],
  566. Cell[BoxData[
  567. RowBox[{"p", ":=",
  568. RowBox[{
  569. RowBox[{"1", "/", "2"}], "*",
  570. RowBox[{"a", ".", "b", ".", "a"}]}]}]], "Input",
  571. CellChangeTimes->{{3.545831029669238*^9, 3.5458310431910114`*^9}, {
  572. 3.5713513889224052`*^9, 3.5713513894743233`*^9}}],
  573. Cell[BoxData[
  574. RowBox[{"vi", ":=",
  575. RowBox[{
  576. RowBox[{"{",
  577. RowBox[{
  578. RowBox[{"{",
  579. RowBox[{
  580. RowBox[{"1", "/", "s0"}], ",", "0", ",", "0", ",", "0", ",", "0", ",",
  581. "0", ",", "0", ",", "0"}], "}"}], ",",
  582. RowBox[{"{",
  583. RowBox[{"0", ",",
  584. RowBox[{"1", "/", "s1"}], ",", "0", ",", "0", ",", "0", ",", "0", ",",
  585. "0", ",", "0"}], "}"}], ",",
  586. RowBox[{"{",
  587. RowBox[{"0", ",", "0", ",",
  588. RowBox[{"1", "/", "s2"}], ",", "0", ",", "0", ",", "0", ",", "0", ",",
  589. "0"}], "}"}], ",",
  590. RowBox[{"{",
  591. RowBox[{"0", ",", "0", ",", "0", ",",
  592. RowBox[{"1", "/", "s3"}], ",", "0", ",", "0", ",", "0", ",", "0"}],
  593. "}"}], ",",
  594. RowBox[{"{",
  595. RowBox[{"0", ",", "0", ",", "0", ",", "0", ",",
  596. RowBox[{"1", "/", "s4"}], ",", "0", ",", "0", ",", "0"}], "}"}], ",",
  597. RowBox[{"{",
  598. RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  599. RowBox[{"1", "/", "s5"}], ",", "0", ",", "0"}], "}"}], ",",
  600. RowBox[{"{",
  601. RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  602. RowBox[{"1", "/", "s6"}], ",", "0"}], "}"}], ",",
  603. RowBox[{"{",
  604. RowBox[{
  605. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  606. RowBox[{"1", "/", "s7"}]}], "}"}]}], "}"}], ".", "\[IndentingNewLine]",
  607. RowBox[{"{",
  608. RowBox[{
  609. RowBox[{"{",
  610. RowBox[{
  611. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  612. "0"}], "}"}], ",",
  613. RowBox[{"{",
  614. RowBox[{
  615. RowBox[{"-", "p0"}], ",", "1", ",", "0", ",", "0", ",", "0", ",", "0",
  616. ",", "0", ",", "0"}], "}"}], ",",
  617. RowBox[{"{",
  618. RowBox[{
  619. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  620. "0"}], "}"}], ",",
  621. RowBox[{"{",
  622. RowBox[{
  623. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  624. "0"}], "}"}], ",",
  625. RowBox[{"{",
  626. RowBox[{
  627. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  628. "0"}], "}"}], ",",
  629. RowBox[{"{",
  630. RowBox[{
  631. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  632. "0"}], "}"}], ",",
  633. RowBox[{"{",
  634. RowBox[{
  635. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  636. "0"}], "}"}], ",",
  637. RowBox[{"{",
  638. RowBox[{
  639. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  640. "1"}], "}"}]}], "}"}], ".", "\[IndentingNewLine]",
  641. RowBox[{"{",
  642. RowBox[{
  643. RowBox[{"{",
  644. RowBox[{
  645. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  646. "0"}], "}"}], ",",
  647. RowBox[{"{",
  648. RowBox[{
  649. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  650. "0"}], "}"}], ",",
  651. RowBox[{"{",
  652. RowBox[{"0", ",",
  653. RowBox[{"-", "p1"}], ",", "1", ",", "0", ",", "0", ",", "0", ",", "0",
  654. ",", "0"}], "}"}], ",",
  655. RowBox[{"{",
  656. RowBox[{
  657. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  658. "0"}], "}"}], ",",
  659. RowBox[{"{",
  660. RowBox[{
  661. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  662. "0"}], "}"}], ",",
  663. RowBox[{"{",
  664. RowBox[{
  665. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  666. "0"}], "}"}], ",",
  667. RowBox[{"{",
  668. RowBox[{
  669. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  670. "0"}], "}"}], ",",
  671. RowBox[{"{",
  672. RowBox[{
  673. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  674. "1"}], "}"}]}], "}"}], ".", "\[IndentingNewLine]",
  675. RowBox[{"{",
  676. RowBox[{
  677. RowBox[{"{",
  678. RowBox[{
  679. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  680. "0"}], "}"}], ",",
  681. RowBox[{"{",
  682. RowBox[{
  683. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  684. "0"}], "}"}], ",",
  685. RowBox[{"{",
  686. RowBox[{
  687. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  688. "0"}], "}"}], ",",
  689. RowBox[{"{",
  690. RowBox[{"0", ",", "0", ",",
  691. RowBox[{"-", "p2"}], ",", "1", ",", "0", ",", "0", ",", "0", ",",
  692. "0"}], "}"}], ",",
  693. RowBox[{"{",
  694. RowBox[{
  695. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  696. "0"}], "}"}], ",",
  697. RowBox[{"{",
  698. RowBox[{
  699. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  700. "0"}], "}"}], ",",
  701. RowBox[{"{",
  702. RowBox[{
  703. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  704. "0"}], "}"}], ",",
  705. RowBox[{"{",
  706. RowBox[{
  707. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  708. "1"}], "}"}]}], "}"}], ".", "\[IndentingNewLine]",
  709. RowBox[{"{",
  710. RowBox[{
  711. RowBox[{"{",
  712. RowBox[{
  713. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  714. "0"}], "}"}], ",",
  715. RowBox[{"{",
  716. RowBox[{
  717. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  718. "0"}], "}"}], ",",
  719. RowBox[{"{",
  720. RowBox[{
  721. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  722. "0"}], "}"}], ",",
  723. RowBox[{"{",
  724. RowBox[{
  725. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  726. "0"}], "}"}], ",",
  727. RowBox[{"{",
  728. RowBox[{"0", ",", "0", ",", "0", ",",
  729. RowBox[{"-", "p3"}], ",", "1", ",", "0", ",", "0", ",", "0"}], "}"}],
  730. ",",
  731. RowBox[{"{",
  732. RowBox[{
  733. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  734. "0"}], "}"}], ",",
  735. RowBox[{"{",
  736. RowBox[{
  737. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  738. "0"}], "}"}], ",",
  739. RowBox[{"{",
  740. RowBox[{
  741. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  742. "1"}], "}"}]}], "}"}], ".", "\[IndentingNewLine]",
  743. RowBox[{"{",
  744. RowBox[{
  745. RowBox[{"{",
  746. RowBox[{
  747. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  748. "0"}], "}"}], ",",
  749. RowBox[{"{",
  750. RowBox[{
  751. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  752. "0"}], "}"}], ",",
  753. RowBox[{"{",
  754. RowBox[{
  755. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  756. "0"}], "}"}], ",",
  757. RowBox[{"{",
  758. RowBox[{
  759. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  760. "0"}], "}"}], ",",
  761. RowBox[{"{",
  762. RowBox[{
  763. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  764. "0"}], "}"}], ",",
  765. RowBox[{"{",
  766. RowBox[{"0", ",", "0", ",", "0", ",", "0", ",",
  767. RowBox[{"-", "p4"}], ",", "1", ",", "0", ",", "0"}], "}"}], ",",
  768. RowBox[{"{",
  769. RowBox[{
  770. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  771. "0"}], "}"}], ",",
  772. RowBox[{"{",
  773. RowBox[{
  774. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  775. "1"}], "}"}]}], "}"}], ".", "\[IndentingNewLine]",
  776. RowBox[{"{",
  777. RowBox[{
  778. RowBox[{"{",
  779. RowBox[{
  780. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  781. "0"}], "}"}], ",",
  782. RowBox[{"{",
  783. RowBox[{
  784. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  785. "0"}], "}"}], ",",
  786. RowBox[{"{",
  787. RowBox[{
  788. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  789. "0"}], "}"}], ",",
  790. RowBox[{"{",
  791. RowBox[{
  792. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  793. "0"}], "}"}], ",",
  794. RowBox[{"{",
  795. RowBox[{
  796. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  797. "0"}], "}"}], ",",
  798. RowBox[{"{",
  799. RowBox[{
  800. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  801. "0"}], "}"}], ",",
  802. RowBox[{"{",
  803. RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  804. RowBox[{"-", "p5"}], ",", "1", ",", "0"}], "}"}], ",",
  805. RowBox[{"{",
  806. RowBox[{
  807. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  808. "1"}], "}"}]}], "}"}], ".", "\[IndentingNewLine]",
  809. RowBox[{"{",
  810. RowBox[{
  811. RowBox[{"{",
  812. RowBox[{
  813. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  814. "0"}], "}"}], ",",
  815. RowBox[{"{",
  816. RowBox[{
  817. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  818. "0"}], "}"}], ",",
  819. RowBox[{"{",
  820. RowBox[{
  821. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  822. "0"}], "}"}], ",",
  823. RowBox[{"{",
  824. RowBox[{
  825. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  826. "0"}], "}"}], ",",
  827. RowBox[{"{",
  828. RowBox[{
  829. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  830. "0"}], "}"}], ",",
  831. RowBox[{"{",
  832. RowBox[{
  833. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  834. "0"}], "}"}], ",",
  835. RowBox[{"{",
  836. RowBox[{
  837. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  838. "0"}], "}"}], ",",
  839. RowBox[{"{",
  840. RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  841. RowBox[{"-", "p6"}], ",", "1"}], "}"}]}], "}"}], ".",
  842. "\[IndentingNewLine]",
  843. RowBox[{"{",
  844. RowBox[{
  845. RowBox[{"{",
  846. RowBox[{
  847. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  848. "0"}], "}"}], ",",
  849. RowBox[{"{",
  850. RowBox[{
  851. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  852. "0"}], "}"}], ",",
  853. RowBox[{"{",
  854. RowBox[{
  855. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  856. "0"}], "}"}], ",",
  857. RowBox[{"{",
  858. RowBox[{
  859. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  860. "0"}], "}"}], ",",
  861. RowBox[{"{",
  862. RowBox[{
  863. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  864. "0"}], "}"}], ",",
  865. RowBox[{"{",
  866. RowBox[{
  867. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  868. "0"}], "}"}], ",",
  869. RowBox[{"{",
  870. RowBox[{
  871. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  872. RowBox[{"-", "q6"}]}], "}"}], ",",
  873. RowBox[{"{",
  874. RowBox[{
  875. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  876. "1"}], "}"}]}], "}"}], ".",
  877. RowBox[{"{",
  878. RowBox[{
  879. RowBox[{"{",
  880. RowBox[{
  881. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  882. "0"}], "}"}], ",",
  883. RowBox[{"{",
  884. RowBox[{
  885. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  886. "0"}], "}"}], ",",
  887. RowBox[{"{",
  888. RowBox[{
  889. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  890. "0"}], "}"}], ",",
  891. RowBox[{"{",
  892. RowBox[{
  893. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  894. "0"}], "}"}], ",",
  895. RowBox[{"{",
  896. RowBox[{
  897. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  898. "0"}], "}"}], ",",
  899. RowBox[{"{",
  900. RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  901. RowBox[{"-", "q5"}], ",", "0"}], "}"}], ",",
  902. RowBox[{"{",
  903. RowBox[{
  904. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  905. "0"}], "}"}], ",",
  906. RowBox[{"{",
  907. RowBox[{
  908. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  909. "1"}], "}"}]}], "}"}], ".",
  910. RowBox[{"{",
  911. RowBox[{
  912. RowBox[{"{",
  913. RowBox[{
  914. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  915. "0"}], "}"}], ",",
  916. RowBox[{"{",
  917. RowBox[{
  918. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  919. "0"}], "}"}], ",",
  920. RowBox[{"{",
  921. RowBox[{
  922. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  923. "0"}], "}"}], ",",
  924. RowBox[{"{",
  925. RowBox[{
  926. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  927. "0"}], "}"}], ",",
  928. RowBox[{"{",
  929. RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  930. RowBox[{"-", "q4"}], ",", "0", ",", "0"}], "}"}], ",",
  931. RowBox[{"{",
  932. RowBox[{
  933. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  934. "0"}], "}"}], ",",
  935. RowBox[{"{",
  936. RowBox[{
  937. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  938. "0"}], "}"}], ",",
  939. RowBox[{"{",
  940. RowBox[{
  941. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  942. "1"}], "}"}]}], "}"}], ".",
  943. RowBox[{"{",
  944. RowBox[{
  945. RowBox[{"{",
  946. RowBox[{
  947. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  948. "0"}], "}"}], ",",
  949. RowBox[{"{",
  950. RowBox[{
  951. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  952. "0"}], "}"}], ",",
  953. RowBox[{"{",
  954. RowBox[{
  955. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  956. "0"}], "}"}], ",",
  957. RowBox[{"{",
  958. RowBox[{"0", ",", "0", ",", "0", ",", "1", ",",
  959. RowBox[{"-", "q3"}], ",", "0", ",", "0", ",", "0"}], "}"}], ",",
  960. RowBox[{"{",
  961. RowBox[{
  962. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  963. "0"}], "}"}], ",",
  964. RowBox[{"{",
  965. RowBox[{
  966. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  967. "0"}], "}"}], ",",
  968. RowBox[{"{",
  969. RowBox[{
  970. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  971. "0"}], "}"}], ",",
  972. RowBox[{"{",
  973. RowBox[{
  974. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  975. "1"}], "}"}]}], "}"}], ".",
  976. RowBox[{"{",
  977. RowBox[{
  978. RowBox[{"{",
  979. RowBox[{
  980. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  981. "0"}], "}"}], ",",
  982. RowBox[{"{",
  983. RowBox[{
  984. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  985. "0"}], "}"}], ",",
  986. RowBox[{"{",
  987. RowBox[{"0", ",", "0", ",", "1", ",",
  988. RowBox[{"-", "q2"}], ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}],
  989. ",",
  990. RowBox[{"{",
  991. RowBox[{
  992. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  993. "0"}], "}"}], ",",
  994. RowBox[{"{",
  995. RowBox[{
  996. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  997. "0"}], "}"}], ",",
  998. RowBox[{"{",
  999. RowBox[{
  1000. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  1001. "0"}], "}"}], ",",
  1002. RowBox[{"{",
  1003. RowBox[{
  1004. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  1005. "0"}], "}"}], ",",
  1006. RowBox[{"{",
  1007. RowBox[{
  1008. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1009. "1"}], "}"}]}], "}"}], ".",
  1010. RowBox[{"{",
  1011. RowBox[{
  1012. RowBox[{"{",
  1013. RowBox[{
  1014. "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1015. "0"}], "}"}], ",",
  1016. RowBox[{"{",
  1017. RowBox[{"0", ",", "1", ",",
  1018. RowBox[{"-", "q1"}], ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1019. "0"}], "}"}], ",",
  1020. RowBox[{"{",
  1021. RowBox[{
  1022. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1023. "0"}], "}"}], ",",
  1024. RowBox[{"{",
  1025. RowBox[{
  1026. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  1027. "0"}], "}"}], ",",
  1028. RowBox[{"{",
  1029. RowBox[{
  1030. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  1031. "0"}], "}"}], ",",
  1032. RowBox[{"{",
  1033. RowBox[{
  1034. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  1035. "0"}], "}"}], ",",
  1036. RowBox[{"{",
  1037. RowBox[{
  1038. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  1039. "0"}], "}"}], ",",
  1040. RowBox[{"{",
  1041. RowBox[{
  1042. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1043. "1"}], "}"}]}], "}"}], ".",
  1044. RowBox[{"{",
  1045. RowBox[{
  1046. RowBox[{"{",
  1047. RowBox[{"1", ",",
  1048. RowBox[{"-", "q0"}], ",", "0", ",", "0", ",", "0", ",", "0", ",", "0",
  1049. ",", "0"}], "}"}], ",",
  1050. RowBox[{"{",
  1051. RowBox[{
  1052. "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1053. "0"}], "}"}], ",",
  1054. RowBox[{"{",
  1055. RowBox[{
  1056. "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1057. "0"}], "}"}], ",",
  1058. RowBox[{"{",
  1059. RowBox[{
  1060. "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",", "0", ",",
  1061. "0"}], "}"}], ",",
  1062. RowBox[{"{",
  1063. RowBox[{
  1064. "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",", "0", ",",
  1065. "0"}], "}"}], ",",
  1066. RowBox[{"{",
  1067. RowBox[{
  1068. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0", ",",
  1069. "0"}], "}"}], ",",
  1070. RowBox[{"{",
  1071. RowBox[{
  1072. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",",
  1073. "0"}], "}"}], ",",
  1074. RowBox[{"{",
  1075. RowBox[{
  1076. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1077. "1"}], "}"}]}], "}"}]}]}]], "Input",
  1078. CellChangeTimes->{{3.5458303069829025`*^9, 3.545830307743946*^9},
  1079. 3.5458303484782763`*^9, {3.54583039506394*^9, 3.5458304538903055`*^9}, {
  1080. 3.545830517654952*^9, 3.545830670162675*^9}}],
  1081. Cell[BoxData[
  1082. RowBox[{"c", ":=",
  1083. RowBox[{"ArrayFlatten", "[",
  1084. RowBox[{"{",
  1085. RowBox[{
  1086. RowBox[{"{",
  1087. RowBox[{
  1088. RowBox[{"IdentityMatrix", "[", "8", "]"}], ",", "0"}], "}"}], ",",
  1089. RowBox[{"{",
  1090. RowBox[{"0", ",", "vi"}], "}"}]}], "}"}], "]"}]}]], "Input",
  1091. CellChangeTimes->{{3.5457683897333775`*^9, 3.5457683916024837`*^9}, {
  1092. 3.545830784950241*^9, 3.5458307852102556`*^9}}],
  1093. Cell[BoxData[
  1094. RowBox[{"pi", ":=",
  1095. RowBox[{
  1096. RowBox[{"1", "/", "2"}], "*",
  1097. RowBox[{"a", ".", "c", ".", "a"}]}]}]], "Input",
  1098. CellChangeTimes->{{3.5713513967514887`*^9, 3.57135139718403*^9}}],
  1099. Cell[BoxData[
  1100. RowBox[{
  1101. RowBox[{"dct", "[", "k_", "]"}], ":=",
  1102. RowBox[{"Table", "[",
  1103. RowBox[{
  1104. RowBox[{
  1105. RowBox[{"Sqrt", "[",
  1106. RowBox[{"2", "/", "k"}], "]"}], "*",
  1107. RowBox[{"Cos", "[",
  1108. RowBox[{
  1109. RowBox[{"(",
  1110. RowBox[{"j", "+",
  1111. RowBox[{"1", "/", "2"}]}], ")"}], "*", "i", "*",
  1112. RowBox[{"\[Pi]", "/", "k"}]}], "]"}], "*",
  1113. RowBox[{"If", "[",
  1114. RowBox[{
  1115. RowBox[{"i", "\[Equal]", "0"}], ",",
  1116. RowBox[{"Sqrt", "[",
  1117. RowBox[{"1", "/", "2"}], "]"}], ",", "1"}], "]"}]}], ",",
  1118. RowBox[{"{",
  1119. RowBox[{"i", ",", "0", ",",
  1120. RowBox[{"k", "-", "1"}]}], "}"}], ",",
  1121. RowBox[{"{",
  1122. RowBox[{"j", ",", "0", ",",
  1123. RowBox[{"k", "-", "1"}]}], "}"}]}], "]"}]}]], "Input",
  1124. CellChangeTimes->{{3.545830925593285*^9, 3.5458309314066176`*^9}}],
  1125. Cell[BoxData[
  1126. RowBox[{"d", ":=",
  1127. RowBox[{"ArrayFlatten", "[",
  1128. RowBox[{"{",
  1129. RowBox[{"{",
  1130. RowBox[{
  1131. "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",",
  1132. "0", ",",
  1133. RowBox[{"dct", "[", "16", "]"}], ",", "0", ",", "0", ",", "0", ",", "0",
  1134. ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], "}"}],
  1135. "]"}]}]], "Input",
  1136. CellChangeTimes->{{3.5713514036194563`*^9, 3.571351403734996*^9}}],
  1137. Cell[BoxData[
  1138. RowBox[{"e", ":=",
  1139. RowBox[{"ArrayFlatten", "[",
  1140. RowBox[{"{",
  1141. RowBox[{
  1142. RowBox[{"{",
  1143. RowBox[{"p", ",", "0"}], "}"}], ",",
  1144. RowBox[{"{",
  1145. RowBox[{"0", ",", "p"}], "}"}]}], "}"}], "]"}]}]], "Input",
  1146. CellChangeTimes->{{3.5457274768888016`*^9, 3.545727477021809*^9}, {
  1147. 3.545764597801342*^9, 3.5457645985633855`*^9}, {3.545830810389696*^9,
  1148. 3.545830853981189*^9}, {3.545830919428932*^9, 3.5458309561620336`*^9}}],
  1149. Cell[BoxData[
  1150. RowBox[{"g", ":=",
  1151. RowBox[{"d", ".", "e"}]}]], "Input"],
  1152. Cell[BoxData[
  1153. RowBox[{"g", ":=", "d"}]], "Input",
  1154. CellChangeTimes->{{3.571351408452491*^9, 3.571351410637926*^9}}],
  1155. Cell[BoxData[
  1156. RowBox[{"f", ":=",
  1157. RowBox[{"ArrayFlatten", "[",
  1158. RowBox[{"{",
  1159. RowBox[{
  1160. RowBox[{"{",
  1161. RowBox[{"pi", ",", "0"}], "}"}], ",",
  1162. RowBox[{"{",
  1163. RowBox[{"0", ",", "pi"}], "}"}]}], "}"}], "]"}]}]], "Input"],
  1164. Cell[BoxData[
  1165. RowBox[{"h", ":=",
  1166. RowBox[{"f", ".",
  1167. RowBox[{"Transpose", "[", "d", "]"}]}]}]], "Input"],
  1168. Cell[BoxData[
  1169. RowBox[{"h", ":=",
  1170. RowBox[{"Transpose", "[", "d", "]"}]}]], "Input",
  1171. CellChangeTimes->{{3.57135141550872*^9, 3.571351418102202*^9}}],
  1172. Cell[BoxData[
  1173. RowBox[{
  1174. RowBox[{"ar", "[",
  1175. RowBox[{"k_", ",", "l_", ",", "r_"}], "]"}], ":=",
  1176. RowBox[{"Table", "[",
  1177. RowBox[{
  1178. RowBox[{"r", "^",
  1179. RowBox[{"(",
  1180. RowBox[{
  1181. RowBox[{"Abs", "[",
  1182. RowBox[{"i", "-", "j"}], "]"}], "+", "l"}], ")"}]}], ",",
  1183. RowBox[{"{",
  1184. RowBox[{"i", ",", "k"}], "}"}], ",",
  1185. RowBox[{"{",
  1186. RowBox[{"j", ",", "k"}], "}"}]}], "]"}]}]], "Input",
  1187. CellChangeTimes->{{3.54583097742725*^9, 3.5458309780892878`*^9}, {
  1188. 3.571351431671892*^9, 3.5713514383357077`*^9}}],
  1189. Cell[BoxData[
  1190. RowBox[{"obj", ":=",
  1191. RowBox[{"GeometricMean", "[",
  1192. RowBox[{
  1193. RowBox[{"Diagonal", "[",
  1194. RowBox[{"g", ".",
  1195. RowBox[{"ar", "[",
  1196. RowBox[{"32", ",", "0", ",", "0.95"}], "]"}], ".",
  1197. RowBox[{"Transpose", "[", "g", "]"}]}], "]"}], "*",
  1198. RowBox[{"Diagonal", "[",
  1199. RowBox[{
  1200. RowBox[{"Transpose", "[", "h", "]"}], ".", "h"}], "]"}]}],
  1201. "]"}]}]], "Input",
  1202. CellChangeTimes->{{3.5458309882078667`*^9, 3.5458309885658865`*^9},
  1203. 3.5466157425283413`*^9, 3.5466160069624662`*^9, {3.546616296586032*^9,
  1204. 3.546616299230183*^9}, {3.5466163309889994`*^9, 3.54661634025953*^9}, {
  1205. 3.546616817198809*^9, 3.546616821172036*^9}, {3.571351466766985*^9,
  1206. 3.571351467025209*^9}}],
  1207. Cell[CellGroupData[{
  1208. Cell[BoxData[
  1209. RowBox[{"cg", "=",
  1210. RowBox[{
  1211. RowBox[{"10", "*",
  1212. RowBox[{"Log", "[",
  1213. RowBox[{"10", ",",
  1214. RowBox[{"1", "/", "obj"}]}], "]"}]}], "//.",
  1215. RowBox[{"{",
  1216. RowBox[{
  1217. RowBox[{"s0", "\[Rule]",
  1218. RowBox[{"90.0", "/", "64"}]}], ",",
  1219. RowBox[{"s1", "\[Rule]",
  1220. RowBox[{"74.0", "/", "64"}]}], ",",
  1221. RowBox[{"s2", "\[Rule]",
  1222. RowBox[{"73.0", "/", "64"}]}], ",",
  1223. RowBox[{"s3", "\[Rule]",
  1224. RowBox[{"71.0", "/", "64"}]}], ",",
  1225. RowBox[{"s4", "\[Rule]",
  1226. RowBox[{"67.0", "/", "64"}]}], ",",
  1227. RowBox[{"s5", "\[Rule]",
  1228. RowBox[{"67.0", "/", "64"}]}], ",",
  1229. RowBox[{"s6", "\[Rule]",
  1230. RowBox[{"67.0", "/", "64"}]}], ",",
  1231. RowBox[{"s7", "\[Rule]",
  1232. RowBox[{"72.0", "/", "64"}]}], ",",
  1233. RowBox[{"p0", "\[Rule]",
  1234. RowBox[{
  1235. RowBox[{"-", "24.0"}], "/", "64"}]}], ",",
  1236. RowBox[{"p1", "\[Rule]",
  1237. RowBox[{
  1238. RowBox[{"-", "23.0"}], "/", "64"}]}], ",",
  1239. RowBox[{"p2", "\[Rule]",
  1240. RowBox[{
  1241. RowBox[{"-", "17.0"}], "/", "64"}]}], ",",
  1242. RowBox[{"p3", "\[Rule]",
  1243. RowBox[{
  1244. RowBox[{"-", "12.0"}], "/", "64"}]}], ",",
  1245. RowBox[{"p4", "\[Rule]",
  1246. RowBox[{
  1247. RowBox[{"-", "14.0"}], "/", "64"}]}], ",",
  1248. RowBox[{"p5", "\[Rule]",
  1249. RowBox[{
  1250. RowBox[{"-", "13.0"}], "/", "64"}]}], ",",
  1251. RowBox[{"p6", "\[Rule]",
  1252. RowBox[{
  1253. RowBox[{"-", "7.0"}], "/", "64"}]}], ",",
  1254. RowBox[{"q0", "\[Rule]",
  1255. RowBox[{"50.0", "/", "64"}]}], ",",
  1256. RowBox[{"q1", "\[Rule]",
  1257. RowBox[{"40.0", "/", "64"}]}], ",",
  1258. RowBox[{"q2", "\[Rule]",
  1259. RowBox[{"31.0", "/", "64"}]}], ",",
  1260. RowBox[{"q3", "\[Rule]",
  1261. RowBox[{"22.0", "/", "64"}]}], ",",
  1262. RowBox[{"q4", "\[Rule]",
  1263. RowBox[{"18.0", "/", "64"}]}], ",",
  1264. RowBox[{"q5", "\[Rule]",
  1265. RowBox[{"16.0", "/", "64"}]}], ",",
  1266. RowBox[{"q6", "\[Rule]",
  1267. RowBox[{"11.0", "/", "64"}]}]}], "}"}]}]}]], "Input",
  1268. CellChangeTimes->{{3.546615331788849*^9, 3.546615552381466*^9}, {
  1269. 3.546615615810094*^9, 3.5466157264434214`*^9}, 3.5466160329119506`*^9,
  1270. 3.546616834253785*^9, {3.570877445159778*^9, 3.570877606901031*^9},
  1271. 3.571351456738694*^9}],
  1272. Cell[BoxData["9.89338176232053`"], "Output",
  1273. CellChangeTimes->{
  1274. 3.5466157337758408`*^9, 3.546616038879292*^9, 3.5466168404231377`*^9,
  1275. 3.570876126258851*^9, 3.5708774092622013`*^9, 3.570877609912533*^9, {
  1276. 3.5713514531814833`*^9, 3.5713515034340487`*^9}, 3.579069516179351*^9}]
  1277. }, Open ]],
  1278. Cell[BoxData[
  1279. RowBox[{
  1280. RowBox[{"ar2d", "[",
  1281. RowBox[{"k_", ",", "r_"}], "]"}], ":=",
  1282. RowBox[{"ArrayFlatten", "[",
  1283. RowBox[{"Table", "[",
  1284. RowBox[{
  1285. RowBox[{"ar", "[",
  1286. RowBox[{"k", ",",
  1287. RowBox[{"Abs", "[",
  1288. RowBox[{"i", "-", "j"}], "]"}], ",", "r"}], "]"}], ",",
  1289. RowBox[{"{",
  1290. RowBox[{"i", ",", "k"}], "}"}], ",",
  1291. RowBox[{"{",
  1292. RowBox[{"j", ",", "k"}], "}"}]}], "]"}], "]"}]}]], "Input",
  1293. CellChangeTimes->{{3.571215344262658*^9, 3.571215350241576*^9},
  1294. 3.571215796960052*^9, 3.5712164912715063`*^9, {3.571217517951624*^9,
  1295. 3.57121752072649*^9}, 3.571217741460986*^9}],
  1296. Cell[BoxData[
  1297. RowBox[{
  1298. RowBox[{"g2d", "[", "k_", "]"}], ":=",
  1299. RowBox[{"ArrayFlatten", "[",
  1300. RowBox[{"Table", "[",
  1301. RowBox[{
  1302. RowBox[{"If", "[",
  1303. RowBox[{
  1304. RowBox[{"i", "\[Equal]", "j"}], ",", "g", ",", "0"}], "]"}], ",",
  1305. RowBox[{"{",
  1306. RowBox[{"i", ",", "k"}], "}"}], ",",
  1307. RowBox[{"{",
  1308. RowBox[{"j", ",", "k"}], "}"}]}], "]"}], "]"}]}]], "Input",
  1309. CellChangeTimes->{{3.571217377033023*^9, 3.5712174099145327`*^9}, {
  1310. 3.5712175508007708`*^9, 3.571217558924345*^9}}],
  1311. Cell[BoxData[
  1312. RowBox[{
  1313. RowBox[{"h2d", "[", "k_", "]"}], ":=",
  1314. RowBox[{"ArrayFlatten", "[",
  1315. RowBox[{"Table", "[",
  1316. RowBox[{
  1317. RowBox[{"If", "[",
  1318. RowBox[{
  1319. RowBox[{"i", "\[Equal]", "j"}], ",", "h", ",", "0"}], "]"}], ",",
  1320. RowBox[{"{",
  1321. RowBox[{"i", ",", "k"}], "}"}], ",",
  1322. RowBox[{"{",
  1323. RowBox[{"j", ",", "k"}], "}"}]}], "]"}], "]"}]}]], "Input",
  1324. CellChangeTimes->{{3.57121768285324*^9, 3.571217702211115*^9}}],
  1325. Cell[BoxData[
  1326. RowBox[{
  1327. RowBox[{"t2d", "[", "k_", "]"}], ":=",
  1328. RowBox[{
  1329. RowBox[{"IdentityMatrix", "[",
  1330. RowBox[{"k", "*",
  1331. RowBox[{"k", "/", "2"}]}], "]"}], "[",
  1332. RowBox[{"[",
  1333. RowBox[{"Flatten", "[",
  1334. RowBox[{"Transpose", "[",
  1335. RowBox[{"Partition", "[",
  1336. RowBox[{
  1337. RowBox[{"Range", "[",
  1338. RowBox[{"k", "*",
  1339. RowBox[{"k", "/", "2"}]}], "]"}], ",",
  1340. RowBox[{"k", "/", "2"}]}], "]"}], "]"}], "]"}], "]"}],
  1341. "]"}]}]], "Input",
  1342. CellChangeTimes->{{3.571216514685655*^9, 3.571216535449212*^9}, {
  1343. 3.571217837719702*^9, 3.571217852699953*^9}}],
  1344. Cell[BoxData[
  1345. RowBox[{"obj2d", ":=",
  1346. RowBox[{"GeometricMean", "[",
  1347. RowBox[{
  1348. RowBox[{"Diagonal", "[",
  1349. RowBox[{
  1350. RowBox[{"g2d", "[", "16", "]"}], ".",
  1351. RowBox[{"t2d", "[", "32", "]"}], ".",
  1352. RowBox[{"g2d", "[", "32", "]"}], ".",
  1353. RowBox[{"ar2d", "[",
  1354. RowBox[{"32", ",", "0.95"}], "]"}], ".",
  1355. RowBox[{"Transpose", "[",
  1356. RowBox[{
  1357. RowBox[{"g2d", "[", "16", "]"}], ".",
  1358. RowBox[{"t2d", "[", "32", "]"}], ".",
  1359. RowBox[{"g2d", "[", "32", "]"}]}], "]"}]}], "]"}], "*",
  1360. RowBox[{"Diagonal", "[",
  1361. RowBox[{
  1362. RowBox[{"Transpose", "[",
  1363. RowBox[{
  1364. RowBox[{"h2d", "[", "32", "]"}], ".",
  1365. RowBox[{"t2d", "[", "32", "]"}], ".",
  1366. RowBox[{"h2d", "[", "16", "]"}]}], "]"}], ".",
  1367. RowBox[{"h2d", "[", "32", "]"}], ".",
  1368. RowBox[{"t2d", "[", "32", "]"}], ".",
  1369. RowBox[{"h2d", "[", "16", "]"}]}], "]"}]}], "]"}]}]], "Input",
  1370. CellChangeTimes->{{3.571217495222146*^9, 3.57121749695989*^9}, {
  1371. 3.571217529685446*^9, 3.5712175379810658`*^9}, {3.571217568436331*^9,
  1372. 3.571217574381446*^9}, {3.571217660893468*^9, 3.5712176794285927`*^9}, {
  1373. 3.5712178183248453`*^9, 3.5712178224715137`*^9}, {3.5712178578153143`*^9,
  1374. 3.571217899033683*^9}, {3.5712180481933937`*^9, 3.571218065866927*^9}, {
  1375. 3.5712189138225317`*^9, 3.571218932628921*^9}, {3.571220575284101*^9,
  1376. 3.571220575519759*^9}, {3.571262758598935*^9, 3.571262768217821*^9}, {
  1377. 3.5713515506821203`*^9, 3.571351607969413*^9}}],
  1378. Cell[CellGroupData[{
  1379. Cell[BoxData[
  1380. RowBox[{"cg2d", "=",
  1381. RowBox[{
  1382. RowBox[{"10", "*",
  1383. RowBox[{"Log", "[",
  1384. RowBox[{"10", ",",
  1385. RowBox[{"1", "/", "obj2d"}]}], "]"}]}], "//.",
  1386. RowBox[{"{",
  1387. RowBox[{
  1388. RowBox[{"s0", "\[Rule]",
  1389. RowBox[{"90.0", "/", "64"}]}], ",",
  1390. RowBox[{"s1", "\[Rule]",
  1391. RowBox[{"74.0", "/", "64"}]}], ",",
  1392. RowBox[{"s2", "\[Rule]",
  1393. RowBox[{"73.0", "/", "64"}]}], ",",
  1394. RowBox[{"s3", "\[Rule]",
  1395. RowBox[{"71.0", "/", "64"}]}], ",",
  1396. RowBox[{"s4", "\[Rule]",
  1397. RowBox[{"67.0", "/", "64"}]}], ",",
  1398. RowBox[{"s5", "\[Rule]",
  1399. RowBox[{"67.0", "/", "64"}]}], ",",
  1400. RowBox[{"s6", "\[Rule]",
  1401. RowBox[{"67.0", "/", "64"}]}], ",",
  1402. RowBox[{"s7", "\[Rule]",
  1403. RowBox[{"72.0", "/", "64"}]}], ",",
  1404. RowBox[{"p0", "\[Rule]",
  1405. RowBox[{
  1406. RowBox[{"-", "24.0"}], "/", "64"}]}], ",",
  1407. RowBox[{"p1", "\[Rule]",
  1408. RowBox[{
  1409. RowBox[{"-", "23.0"}], "/", "64"}]}], ",",
  1410. RowBox[{"p2", "\[Rule]",
  1411. RowBox[{
  1412. RowBox[{"-", "17.0"}], "/", "64"}]}], ",",
  1413. RowBox[{"p3", "\[Rule]",
  1414. RowBox[{
  1415. RowBox[{"-", "12.0"}], "/", "64"}]}], ",",
  1416. RowBox[{"p4", "\[Rule]",
  1417. RowBox[{
  1418. RowBox[{"-", "14.0"}], "/", "64"}]}], ",",
  1419. RowBox[{"p5", "\[Rule]",
  1420. RowBox[{
  1421. RowBox[{"-", "13.0"}], "/", "64"}]}], ",",
  1422. RowBox[{"p6", "\[Rule]",
  1423. RowBox[{
  1424. RowBox[{"-", "7.0"}], "/", "64"}]}], ",",
  1425. RowBox[{"q0", "\[Rule]",
  1426. RowBox[{"50.0", "/", "64"}]}], ",",
  1427. RowBox[{"q1", "\[Rule]",
  1428. RowBox[{"40.0", "/", "64"}]}], ",",
  1429. RowBox[{"q2", "\[Rule]",
  1430. RowBox[{"31.0", "/", "64"}]}], ",",
  1431. RowBox[{"q3", "\[Rule]",
  1432. RowBox[{"22.0", "/", "64"}]}], ",",
  1433. RowBox[{"q4", "\[Rule]",
  1434. RowBox[{"18.0", "/", "64"}]}], ",",
  1435. RowBox[{"q5", "\[Rule]",
  1436. RowBox[{"16.0", "/", "64"}]}], ",",
  1437. RowBox[{"q6", "\[Rule]",
  1438. RowBox[{"11.0", "/", "64"}]}]}], "}"}]}]}]], "Input",
  1439. CellChangeTimes->{{3.571218095000925*^9, 3.5712180969589567`*^9},
  1440. 3.571358649990808*^9, {3.5713587744026117`*^9, 3.5713587751039753`*^9}}],
  1441. Cell[BoxData["$Aborted"], "Output",
  1442. CellChangeTimes->{3.579069718407742*^9}]
  1443. }, Open ]]
  1444. },
  1445. WindowSize->{1424, 775},
  1446. WindowMargins->{{0, Automatic}, {Automatic, 33}},
  1447. FrontEndVersion->"7.0 for Linux x86 (64-bit) (February 25, 2009)",
  1448. StyleDefinitions->"Default.nb"
  1449. ]
  1450. (* End of Notebook Content *)
  1451. (* Internal cache information *)
  1452. (*CellTagsOutline
  1453. CellTagsIndex->{}
  1454. *)
  1455. (*CellTagsIndex
  1456. CellTagsIndex->{}
  1457. *)
  1458. (*NotebookFileOutline
  1459. Notebook[{
  1460. Cell[545, 20, 631, 18, 32, "Input"],
  1461. Cell[1179, 40, 17489, 516, 517, "Input"],
  1462. Cell[18671, 558, 460, 12, 32, "Input"],
  1463. Cell[19134, 572, 252, 6, 32, "Input"],
  1464. Cell[19389, 580, 17915, 507, 605, "Input"],
  1465. Cell[37307, 1089, 412, 11, 32, "Input"],
  1466. Cell[37722, 1102, 199, 5, 32, "Input"],
  1467. Cell[37924, 1109, 845, 25, 32, "Input"],
  1468. Cell[38772, 1136, 441, 11, 32, "Input"],
  1469. Cell[39216, 1149, 460, 11, 32, "Input"],
  1470. Cell[39679, 1162, 73, 2, 32, "Input"],
  1471. Cell[39755, 1166, 116, 2, 32, "Input"],
  1472. Cell[39874, 1170, 246, 8, 32, "Input"],
  1473. Cell[40123, 1180, 110, 3, 32, "Input"],
  1474. Cell[40236, 1185, 151, 3, 32, "Input"],
  1475. Cell[40390, 1190, 544, 16, 32, "Input"],
  1476. Cell[40937, 1208, 730, 17, 32, "Input"],
  1477. Cell[CellGroupData[{
  1478. Cell[41692, 1229, 2245, 63, 77, "Input"],
  1479. Cell[43940, 1294, 285, 4, 31, "Output"]
  1480. }, Open ]],
  1481. Cell[44240, 1301, 640, 17, 32, "Input"],
  1482. Cell[44883, 1320, 517, 14, 32, "Input"],
  1483. Cell[45403, 1336, 463, 13, 32, "Input"],
  1484. Cell[45869, 1351, 616, 18, 32, "Input"],
  1485. Cell[46488, 1371, 1511, 33, 55, "Input"],
  1486. Cell[CellGroupData[{
  1487. Cell[48024, 1408, 2153, 61, 77, "Input"],
  1488. Cell[50180, 1471, 77, 1, 31, "Output"]
  1489. }, Open ]]
  1490. }
  1491. ]
  1492. *)
  1493. (* End of internal cache information *)