Difference between revisions of "Q Evaluation routines"

From magneticfields
Jump to: navigation, search
Line 1: Line 1:
 
We have evaluated Q using two independent routines, one by Tashiro, second by Ferrer, and performed a triple check using Vachaspati's Mathematica notebook.
 
We have evaluated Q using two independent routines, one by Tashiro, second by Ferrer, and performed a triple check using Vachaspati's Mathematica notebook.
 +
 +
 +
== Tashiro's C routines ==
 +
  
 
Here are Tashiro's C routines.
 
Here are Tashiro's C routines.
  
Here are Ferrer's Python routines.
+
 
 +
== Ferrer's Python routines ==
 +
 
 +
The following Python script takes the files in nosourceoutput.tar.gz as input and calculates Q for each E1 and E2 combination [[Media:qcalc.py]]. A small modification of the script that calculates Q sepparately for the northern and southern hemispheres is given here [[Media:qcalcnorthsouth.py]]. Both scripts generate a figure of the output and write the value of Q and the standard error to file. For example, the first script produces the figure [[Media:jan14z100bgt80.pdf]] and the file arrays [[Media:qarrays.tar.gz]].
 +
 
 +
 
 +
 
 +
== Vachaspati's Mathematica routines ==
 +
 
  
 
Here is Vachaspati's Mathematica notebook: [[File:Cpgamma-northsouth.nb]]
 
Here is Vachaspati's Mathematica notebook: [[File:Cpgamma-northsouth.nb]]

Revision as of 20:09, 19 May 2014

We have evaluated Q using two independent routines, one by Tashiro, second by Ferrer, and performed a triple check using Vachaspati's Mathematica notebook.


Tashiro's C routines

Here are Tashiro's C routines.


Ferrer's Python routines

The following Python script takes the files in nosourceoutput.tar.gz as input and calculates Q for each E1 and E2 combination Media:qcalc.py. A small modification of the script that calculates Q sepparately for the northern and southern hemispheres is given here Media:qcalcnorthsouth.py. Both scripts generate a figure of the output and write the value of Q and the standard error to file. For example, the first script produces the figure Media:jan14z100bgt80.pdf and the file arrays Media:qarrays.tar.gz.


Vachaspati's Mathematica routines

Here is Vachaspati's Mathematica notebook: File:Cpgamma-northsouth.nb