aboutsummaryrefslogtreecommitdiffstats
path: root/other/ipython_log.py
blob: fa1fd2973b7147dddbf6cb8a16f5df044c9da68d (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
#log# Automatic Logger file. *** THIS MUST BE THE FIRST LINE ***
#log# DO NOT CHANGE THIS LINE OR THE TWO BELOW
#log# opts = Struct({'__allownew': True,
 'c': '\nimport sage.misc.misc; print sage.misc.misc.branch_current_hg_notice(sage.misc.misc.branch_current_hg()); from sage.misc.interpreter import preparser; preparser(True);import sage.all_cmdline; sage.all_cmdline._init_cmdline(globals());from sage.all import Integer, RealNumber;import os; os.chdir("/home/bnewbold/thesis/other");import sage.misc.interpreter;from sage.misc.interpreter import attached_files;from sage.all_cmdline import *;_=sage.misc.interpreter.load_startup_file("/home/bnewbold/.sage//init.sage");',
 'logfile': 'ipython_log.py'})
#log# args = []
#log# It is safe to make manual edits below here.
#log#-----------------------------------------------------------------------
import sage.misc.misc; print sage.misc.misc.branch_current_hg_notice(sage.misc.misc.branch_current_hg()); from sage.misc.interpreter import preparser; preparser(True);import sage.all_cmdline; sage.all_cmdline._init_cmdline(globals());from sage.all import Integer, RealNumber;import os; os.chdir("/home/bnewbold/thesis/other");import sage.misc.interpreter;from sage.misc.interpreter import attached_files;from sage.all_cmdline import *;_=sage.misc.interpreter.load_startup_file("/home/bnewbold/.sage//init.sage");


from pylab import *

d = loadtxt('./huge-output.data')
plot(d)
d
d = d[Integer(0)]
d
d = loadtxt('./huge-output.data')
transpose(d)
d = transpose(d)
d = transpose(d)
#?plot

d
d = transpose(d)
d
p = plot(d[Integer(0)],d[Integer(1)],'b.')
p
show(p)
p
p = p[Integer(0)]
p
show(p)
_ip.system("rm huge-output.png")
savefig('huge-output.png')
d
max(d[Integer(0)])
xrange([-Integer(1),Integer(1)])
#?xrange
p
p.xaxis()
f = figure()
gca()
a = gca()
a.xaxis
a.xaxis()
#?a.set_xlim
#?a.set_xlim
a.set_xlim((-Integer(1),Integer(1)))
a.set_ylim((-Integer(1),Integer(1)))
savefig('huge-output.png')
p = plot(d[Integer(0)],d[Integer(1)],'b.')
savefig('huge-output.png')
a = gca()
a.set_ylim((-Integer(1),Integer(1)))
a.set_xlim((-Integer(1),Integer(1)))
min(d[Integer(0)])
max(d[Integer(0)])
min(d[Integer(1)])
max(d[Integer(1)])
size(d)
d
mean(d)
mean(d[Integer(0)])
median(d[Integer(0)])
median(d[Integer(1)])
p = plot(d[Integer(0)][:Integer(10)],d[Integer(1)][:Integer(10)],'k.')
savefig('huge-output.png')
a = gca()
a.set_xlim((-Integer(1),Integer(1)))
a.set_ylim((-Integer(1),Integer(1)))
savefig('huge-output.png')
#?plot
figure()
plot(d[Integer(0)],d[Integer(1)],'k',markersize=RealNumber('.0001'))
def doaxes(): a = gca(); a.set_xlim((-Integer(1),Integer(1))); a.set_ylim((-Integer(1),Integer(1)));

doaxes()
savefig('huge-output.png')
#?savefig
savefig('huge-output.png', dpi=Integer(1000))
savefig('huge-output.png', dpi=Integer(200))
savefig('huge-output.png', dpi=Integer(200), format='png')
#?savefig
savefig('huge-output.png')
a = gca()
a.set_xlim[-RealNumber('.001'),RealNumber('.001')]
a.set_xlim((-RealNumber('.001'),RealNumber('.001')))
savefig('huge-output.png')
figure()
plot(d[Integer(0)][Integer(5000)],d[Integer(1)][Integer(5000)],'k',markersize=Integer(1))
plot(d[Integer(0)][Integer(5000)],d[Integer(1)][Integer(5000)],'k',markersize=Integer(1))
d
size(d)
size(d[Integer(0)])
size(d[Integer(1)])
plot(d[Integer(0)][Integer(5000)],d[Integer(1)][Integer(5000)],'k',markersize=RealNumber('.1'))
_ip.magic("logon ")
_ip.magic("logstart ")