Difference between revisions of "Q Evaluation routines"
(→Ferrer's Python routines) |
(→Vachaspati's Mathematica routines) |
||
Line 12: | Line 12: | ||
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 separately 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]]. | 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 separately 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 | + | == Vachaspati's Mathematica routine == |
Intended as a check of some of the numbers generated by other routines: [[File:Cpgamma-northsouth.nb]] | Intended as a check of some of the numbers generated by other routines: [[File:Cpgamma-northsouth.nb]] |
Revision as of 21:35, 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 separately 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 routine
Intended as a check of some of the numbers generated by other routines: File:Cpgamma-northsouth.nb