[Commits] [svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 150)

[knarf at cct.lsu.edu]

User: knarf
Date: 2011/10/17 12:21 AM

Added:
 /examples/collapse/
  rho_max_conv.pdf, rho_max_conv.py

Modified:
 /examples/collapse/
  make_plots.sh, radii.pdf, radii.py, rho_max.pdf

Log:
 one more matplotlib plot

Directory Changes:

Directory: /svn:executable/
===========================

   + *
Directory: /svn:mime-type/
==========================

   + application/pdf

File Changes:

Directory: /examples/collapse/
==============================

File [modified]: make_plots.sh
Delta lines: +3 -3
===================================================================
--- examples/collapse/make_plots.sh	2011-10-17 04:58:40 UTC (rev 149)
+++ examples/collapse/make_plots.sh	2011-10-17 05:21:57 UTC (rev 150)
@@ -1,5 +1,5 @@
 #!/bin/bash
+#gnuplot *.plt
 
-# use find_edge.awk to find the outer edge of the star if needed
-
-gnuplot *.plt
+echo ./rho_max.py ./radii.py ./rho_max_conv.py \
+| xargs -n 1 -P 4 bash -c 

File [modified]: radii.pdf
Delta lines: +0 -0
===================================================================
(Binary files differ)

File [modified]: radii.py
Delta lines: +1 -2
===================================================================
--- examples/collapse/radii.py	2011-10-17 04:58:40 UTC (rev 149)
+++ examples/collapse/radii.py	2011-10-17 05:21:57 UTC (rev 150)
@@ -36,8 +36,7 @@
 ax2.xaxis.set_major_locator(mticker.MaxNLocator(7))
 ax2.xaxis.set_minor_locator(mticker.MaxNLocator(14))
 ax.set_ylabel(r'$R/R_G$')
-ax.yaxis.set_major_locator(mticker.MaxNLocator(5))
-ax.yaxis.set_minor_locator(mticker.MaxNLocator(10))
+ax.yaxis.set_major_locator(mticker.MaxNLocator(8))
 ax.yaxis.grid(False)
 set_tick_sizes(ax, 8, 4)
 

File [modified]: rho_max.pdf
Delta lines: +0 -0
===================================================================
(Binary files differ)

File [added]: rho_max_conv.pdf
Delta lines: +0 -0
===================================================================
(Binary files differ)



Property changes on: examples/collapse/rho_max_conv.pdf
___________________________________________________________________

File [added]: rho_max_conv.py
Delta lines: +46 -0
===================================================================
--- examples/collapse/rho_max_conv.py	                        (rev 0)
+++ examples/collapse/rho_max_conv.py	2011-10-17 05:21:57 UTC (rev 150)
@@ -0,0 +1,46 @@
+#!/usr/bin/python
+
+# Plot evolution of convergence factor of hamiltonian constraint violation at x=0
+from plot_defaults import *
+import math
+
+# stuff
+xlim = (0,1400)
+ylim = (1,5)
+
+# load data
+Fx, Fy3, Fy5, Fy6 = np.loadtxt("rho_maximum_convergence.asc", comments="#", usecols=(1,2,4,5), unpack=True)
+
+# Plot basics
+fig = plt.figure(figsize=(6, 3))
+fig.subplots_adjust(top=0.85, bottom=0.16, left=0.11,right=0.97)
+ax = fig.add_subplot(1,1,1)
+
+for i in range(len(Fy3)):
+  Fy3[i] = math.log(max(1.e-10, (Fy3[i]-Fy5[i])/(Fy5[i]-Fy6[i])),2)
+# plot
+ax.plot(Fx, Fy3, linestyle='-' , color='black')
+
+# plot properties
+#ax.set_xlim(xlim)
+ax.set_ylim(ylim) 
+
+#plt.title('PSD', fontsize=fontsize)
+
+ax.set_xlabel(r't [M]')
+ax.xaxis.set_major_locator(mticker.MaxNLocator(7))
+ax.xaxis.set_minor_locator(mticker.MaxNLocator(14))
+ax.xaxis.grid(False)
+ax2 = ax.twiny()
+ax2.set_xlabel(r't [ms]')
+ax2.set_xlim((ax.get_xlim()[0]/M_to_ms, ax.get_xlim()[1]/M_to_ms))
+ax2.xaxis.set_major_locator(mticker.MaxNLocator(7))
+ax2.xaxis.set_minor_locator(mticker.MaxNLocator(14))
+ax.set_ylabel(r'$C(\varrho_C)$')
+ax.yaxis.set_major_locator(mticker.MaxNLocator(8))
+ax.yaxis.grid(False)
+set_tick_sizes(ax, 8, 4)
+
+#plt.show()
+plt.savefig('rho_max_conv.pdf')
+



Property changes on: examples/collapse/rho_max_conv.py
___________________________________________________________________