• 0 Vote(s) - 0 Average
• 1
• 2
• 3
• 4
• 5
 Matrix Operator Method Gottfried Ultimate Fellow Posts: 766 Threads: 119 Joined: Aug 2007 09/24/2008, 08:22 PM (This post was last modified: 09/25/2008, 07:38 AM by Gottfried.) Just derived a method to compute exact entries for powers of the (square) matrix-operator for T-tetration. It is applicable to positive integer powers only, but for any base. The restriction to positive integer powers lets look such solutions useless, since integer iteration-height can easily computed just using the scalar values. But I'll use this for further analysis of powerseries, series of powertowers and hopefully one time for the fractional iteration... Let's use the following notational conventions: Code:```´ b^^h - the powertower of height h using base b V(x) - the vandermonde-column-vector containing consecutive powers of its parameter x:         V(x) = column(1,x,x^2,x^3,...) dV(x) - used as diagonal-matrix T - the matrix which performs T-tetration to base b (in our forum:"exp_b°t(x)" ):       V(x)~ * T = V(b^x)~ U - the matrix which performs U-tetration to base b (in our forum:"dxp_b°t(x)" ):       V(x)~ * U = V(b^x - 1) ~       Note, that U is lower triangular.       The triangularity allows to compute exact entries for the integer matrix-powers.``` Then the entries for positive integer powers of T can be finitely computed and are "exact", as far as we assume scalar logarithms and exponentials as exact: Code:```´ T^2 = U*dV(b^^0) * T*dV(b^^1) T^3 = U*dV(b^^0) * U*dV(b^^1) * T*dV(b^^2) ... T^h = prod_{k=0}^{h-2} (U * dV(b^^k))            * (T * dV(b^^(h-1)))```This finding is interesting, because in my matrix-method I had to use fixpoint-shift to get exact entries even for integer powers, and since the fixpoints for T-tetration are real only for a small range of bases we had to deal with complex-valued U-matrices when considering the general case. Here we do not need a fixpoint-shift. I did not check how this computation is related to Ioannis Galidakis' method for exact entries yet, but I think, this is interesting too. Here is the top left of the symbolic T^2, where lambda=log(b). Each row has to be multiplied by the entry in the most left column and each column must also be multiplied by the entry in the first row.     Here is the top left of the symbolic T^3. (b^^2 means b^b). Legend as before     The difference between the symbolic computation and the simple matrixpower is interesting. I used dim=64x64, base b=sqrt(2), which provides a good approximation when the simple matrix-power is computed. Here are two (zoomed) images: very good aproximation in the leading 12 columns (abs differences to the exact values <1e-20 ), but in the columns 52 to 63 the differences grow up to absolute values greater than 1e10. Surely I "knew" that differences should occur, but I hadn't guessed, that they are so large - I just didn't investigate this in detail. The leading first twelve columns of the matrix of differences:     The twelve rightmost columns:     The large errors are actually still relatively small for that base b=sqrt(2). A measure for the quality of approximation is, whether the resulting vector of V(x)~*T^3 = Y~ is actually vandermonde and thus Y = V(y) . This means, that the ratios of logarithms of its entries : log(Y[k])/log(Y[1]), k=0..63, should give the exact sequence [0,1,2,3,...], because this means, that Y contains indeed the consecutive powers of Y[1]. Here is a table of that ratios. ( Remember: we check the col-sums of the third power of T, using x=1) Code:```´ column    symbolic              "naive"             difference ------------------------------------------------------------   0  -1.13831366798E-19        0.E-201       -1.13831366798E-19   1       1.00000000000  1.00000000000                  0.E-202   2       2.00000000000  2.00000000000        1.11189479752E-19   3       3.00000000000  3.00000000000        2.20208219209E-19   4       4.00000000000  4.00000000000        3.27448139227E-19 .... 42      42.0000000000  42.0000000000        1.19448154596E-11 43      43.0000000000  43.0000000000        3.09972987348E-11 44      44.0000000000  43.9999999999        7.77154497620E-11 45      45.0000000000  44.9999999998        1.88541576938E-10 46      46.0000000000  45.9999999996  0.000000000443257385765 47      47.0000000000  46.9999999990   0.00000000101122118773 48      48.0000000000  47.9999999978   0.00000000224146913405 49      49.0000000000  48.9999999952   0.00000000483322217091 50      50.0000000000  49.9999999899    0.0000000101495726916 51      51.0000000000  50.9999999792    0.0000000207790885254 52      52.0000000000  51.9999999585    0.0000000415151935602 53      53.0000000000  52.9999999190    0.0000000810212593226 54      54.0000000000  53.9999998454     0.000000154592714129 55      55.0000000000  54.9999997114     0.000000288630924056 56      56.0000000000  55.9999994723     0.000000527723771045 57      57.0000000000  56.9999990544     0.000000945602091679 58      58.0000000000  57.9999983383      0.00000166172556358 59      59.0000000000  58.9999971341      0.00000286585838319 60      60.0000000000  59.9999951463      0.00000485372880576 61      61.0000000000  60.9999919223      0.00000807772027087 62      62.0000000000  61.9999867825       0.0000132174922318 63      63.0000000000  62.9999787236       0.0000212764324046``` For base b=2 this looks already catastrophic for the "naive"-computation: Code:```´ column    symbolic              "naive"             difference ------------------------------------------------------------   0  -4.36636233681E-20        0.E-202      -4.36636233681E-20   1       1.00000000000  1.00000000000                 0.E-202   2       2.00000000000  2.00000000000       2.46662212171E-14   3       3.00000000000  2.99999999998       2.31650095614E-11   4       4.00000000000  3.99999999669  0.00000000331064284896   5       5.00000000000  4.99999985867    0.000000141325875791   6       6.00000000000  5.99999733964     0.00000266036059669 .... 50      50.0000000000  37.9413247398           12.0586752602 51      51.0000000000  38.3796009369           12.6203990631 52      52.0000000000  38.8098393554           13.1901606446 53      53.0000000000  39.2323121795           13.7676878205 54      54.0000000000  39.6472796473           14.3527203527 55      55.0000000000  40.0549905348           14.9450094652 56      56.0000000000  40.4556826502           15.5443173498 57      57.0000000000  40.8495833279           16.1504166721 58      58.0000000000  41.2369099170           16.7630900830 59      59.0000000000  41.6178702587           17.3821297413 60      60.0000000000  41.9926631496           18.0073368504 61      61.0000000000  42.3614787893           18.6385212107 62      62.0000000000  42.7244992089           19.2755007911 63      63.0000000000  43.0818986820           19.9181013180``` It is obvious, that we should use the "exact" (symbolic) description, if we ever explicitely consider powers of the tetration-matrix T. Gottfried Helms, Kassel « Next Oldest | Next Newest »

 Messages In This Thread Matrix Operator Method - by Gottfried - 08/12/2007, 08:08 PM RE: Matrix Operator Method - by bo198214 - 08/13/2007, 04:15 AM RE: Matrix Operator Method - by jaydfox - 08/13/2007, 05:40 AM RE: Matrix Operator Method - by Gottfried - 08/13/2007, 09:22 AM RE: Matrix Operator Method - by bo198214 - 08/14/2007, 03:43 PM RE: Matrix Operator Method - by Gottfried - 08/14/2007, 04:15 PM RE: Matrix Operator Method - by bo198214 - 08/26/2007, 12:18 AM RE: Matrix Operator Method - by Gottfried - 08/26/2007, 11:24 AM RE: Matrix Operator Method - by bo198214 - 08/26/2007, 11:39 AM RE: Matrix Operator Method - by Gottfried - 08/26/2007, 04:22 PM RE: Matrix Operator Method - by Gottfried - 08/26/2007, 10:54 PM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 08:29 AM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 11:04 AM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 11:35 AM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 11:58 AM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 12:13 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 01:19 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 02:29 PM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 02:36 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 03:09 PM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 07:15 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 08:15 PM RE: Matrix Operator Method - by bo198214 - 08/29/2007, 05:28 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 12:43 PM RE: Matrix Operator Method - by Gottfried - 10/08/2007, 12:11 PM RE: Matrix Operator Method - by Gottfried - 10/14/2007, 09:32 PM RE: Matrix Operator Method - by Gottfried - 04/04/2008, 09:41 AM RE: Matrix Operator Method - by Gottfried - 04/17/2008, 09:21 PM RE: Matrix Operator Method - by bo198214 - 04/25/2008, 03:39 PM RE: Matrix Operator Method - by Gottfried - 04/26/2008, 06:09 PM RE: Matrix Operator Method - by bo198214 - 04/26/2008, 06:47 PM RE: Matrix Operator Method - by Gottfried - 04/18/2008, 01:55 PM RE: Matrix Operator Method - by Gottfried - 07/08/2008, 06:46 AM Diagonalization for dxp/basic facts/Pari-routine - by Gottfried - 08/08/2008, 01:12 PM Exact entries for T-tetration Bell-matrix - by Gottfried - 09/24/2008, 08:22 PM RE: Exact entries for T-tetration Bell-matrix - by bo198214 - 09/26/2008, 07:30 AM RE: Exact entries for T-tetration Bell-matrix - by Gottfried - 09/26/2008, 09:56 AM

 Possibly Related Threads... Thread Author Replies Views Last Post My interpolation method [2020] tommy1729 1 315 02/20/2020, 08:40 PM Last Post: tommy1729 Kneser method question tommy1729 9 1,419 02/11/2020, 01:26 AM Last Post: sheldonison Half-iterates and periodic stuff , my mod method [2019] tommy1729 0 613 09/09/2019, 10:55 PM Last Post: tommy1729 A fundamental flaw of an operator who's super operator is addition JmsNxn 4 7,526 06/23/2019, 08:19 PM Last Post: Chenjesu Tetration and Sign operator tetration101 0 603 05/15/2019, 07:55 PM Last Post: tetration101 2 fixpoints , 1 period --> method of iteration series tommy1729 0 1,888 12/21/2016, 01:27 PM Last Post: tommy1729 Tommy's matrix method for superlogarithm. tommy1729 0 1,879 05/07/2016, 12:28 PM Last Post: tommy1729 [split] Understanding Kneser Riemann method andydude 7 8,764 01/13/2016, 10:58 PM Last Post: sheldonison [2015] New zeration and matrix log ? tommy1729 1 3,531 03/24/2015, 07:07 AM Last Post: marraco Kouznetsov-Tommy-Cauchy method tommy1729 0 2,215 02/18/2015, 07:05 PM Last Post: tommy1729

Users browsing this thread: 1 Guest(s)