Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Cheta with base-change: preliminary results
#9
I'm attaching SAGE variables, for those who use SAGE. The zip file contains two sobj files, which are SAGE objects. I'm also attaching a text file with the coefficients, one per line, for those who don't have SAGE. (Note: I'm using SAGE 4.0.1, a slightly outdated version.)

First, the file jay_sexp_vec_128terms_1024bits.sobj contains the coefficients in a vector(RealField(1024), 129). These are the coefficients of the Taylor series, not the derivatives directly (i.e., the factorial is already factored in).

The file jay_sexp_pol_128terms_1024bits.sobj is a polynomial over the ring PolynomialRing(ComplexField(1024), 'z'), with 129 terms (128 not including the constant). To use it, simply load it into a variable, e.g., sexp, then use it like a function:

Code:
sexp = load('jay_sexp_pol_128terms_1024bits.sobj');
APC = ComplexField(1024);
print sexp(APC(0.2, 0.1));

To convert the vector into a polynomial, you could do the following (there's probably an easier way, but the documentation is a bit sparse at times):

Code:
svec = load('jay_sexp_vec_128terms_1024bits.sobj');
APC = ComplexField(1024);
APCPol = PolynomialRing(APC, 'z');
z = APCPol('z');
sexp = sum(s[kk] * z**kk for kk in xrange(len(svec)));

The text file jay_sexp_coeffs_128terms_256bits.txt contains the coefficients in a text format, one per line. I truncated it to 256 bits, so that the decimal expansions would be a reasonable length.


Attached Files
.zip   jay_sexp_128terms_1024bits.zip (Size: 37.71 KB / Downloads: 282)
.txt   jay_sexp_coeffs_128terms_256bits.txt (Size: 10.31 KB / Downloads: 316)
~ Jay Daniel Fox
Reply


Messages In This Thread
RE: Cheta with base-change: preliminary results - by jaydfox - 08/12/2009, 12:24 AM



Users browsing this thread: 1 Guest(s)