Advertisement
MatsGranvik

Sum of eigenvalues of matrices A and B equal C

Jan 2nd, 2020
234
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.43 KB | None | 0 0
  1. (*start*)
  2. (*Mathematica 8.0.1*)
  3. Clear[C1, C2, A, B, s, nn, i];
  4. mm = 32;
  5.  
  6. "Counts to 32"
  7. t = 1;
  8. Monitor[Table[
  9. A = Table[
  10. Table[If[Mod[n, t*k] == 0, Exp[Exp[Exp[1]]], 0], {k, 1, nn}], {n,
  11. 1, nn}].Table[
  12. Table[If[Mod[k, t*n] == 0, MoebiusMu[n], 0], {k, 1, nn}], {n, 1,
  13. nn}];
  14. a = Eigenvalues[A];
  15. N[Sum[Sign[a[[i]]] If[a[[i]] == 0, 0,
  16. Log[Log[Log[Abs[a[[i]]]]]]], {i, 1, nn}], 6], {nn, 1, mm}], nn]
  17. sumA = Round[%]
  18.  
  19. "Counts to 32"
  20. t = 2;
  21. Monitor[Table[
  22. B = Table[
  23. Table[If[Mod[n, t*k] == 0, Exp[Exp[Exp[1]]], 0], {k, 1, nn}], {n,
  24. 1, nn}].Table[
  25. Table[If[Mod[k, t*n] == 0, MoebiusMu[n], 0], {k, 1, nn}], {n, 1,
  26. nn}];
  27. b = Eigenvalues[B];
  28. N[Sum[Sign[b[[i]]] If[b[[i]] == 0, 0,
  29. Log[Log[Log[Abs[b[[i]]]]]]], {i, 1, nn}], 6], {nn, 1, mm}], nn]
  30. sumB = Round[%]
  31.  
  32. "Counts to 32"
  33. C1 = A + B;
  34. Monitor[Table[C2 = Table[Table[C1[[n, k]], {k, 1, nn}], {n, 1, nn}];
  35. c = Eigenvalues[C2];
  36. N[Sum[Sign[c[[i]]] If[c[[i]] == 0, 0,
  37. Log[Log[Log[Abs[c[[i]]]]]]], {i, 1, nn}], 6], {nn, 1, mm}], nn]
  38. sumC1 = Round[%]
  39.  
  40. sumA + sumB
  41. Sum[Table[
  42. Sum[If[Mod[n, k] == 0, MoebiusMu[n/k], 0], {n, 1, nn}], {nn, 1,
  43. mm}], {k, 1, t}]
  44. %%% - %%
  45. %%% - %%
  46.  
  47. t = 1;
  48. muA = Table[
  49. Sum[If[Mod[n, t] == 0, MoebiusMu[n/t], 0], {n, 1, nn}], {nn, 1, mm}]
  50. t = 2;
  51. muB = Table[
  52. Sum[If[Mod[n, t] == 0, MoebiusMu[n/t], 0], {n, 1, nn}], {nn, 1, mm}]
  53. muC = muA + muB
  54. (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement