-
Notifications
You must be signed in to change notification settings - Fork 0
/
paper.tex
202 lines (159 loc) · 6.77 KB
/
paper.tex
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
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
%\documentclass[12pt, preprint]{aastex}
\documentclass[preprint,times,tighten]{aastex63}
\usepackage{epsf,color,amsmath}
\usepackage{multirow}
\usepackage{cancel}
\newcommand{\sfrac}[2]{\mathchoice%
{\kern0em\raise.5ex\hbox{\the\scriptfont0 #1}\kern-.15em/
\kern-.15em\lower.25ex\hbox{\the\scriptfont0 #2}}
{\kern0em\raise.5ex\hbox{\the\scriptfont0 #1}\kern-.15em/
\kern-.15em\lower.25ex\hbox{\the\scriptfont0 #2}}
{\kern0em\raise.5ex\hbox{\the\scriptscriptfont0 #1}\kern-.2em/
\kern-.15em\lower.25ex\hbox{\the\scriptscriptfont0 #2}} {#1\!/#2}}
\newcommand{\myhalf}{\sfrac{1}{2}}
\newcommand{\nph}{{n+\myhalf}}
\newcommand{\nmh}{{n-\myhalf}}
\newcommand{\inp}{\mathrm{in}}
\newcommand{\outp}{\mathrm{out}}
% boldsymbol means bold italic
\newcommand{\eb}{{\bf{e}}}
\newcommand{\Ub}{{\bf{U}}}
\newcommand{\xb}{{\bf{x}}}
\newcommand{\kb}{{\bf{k}}}
\newcommand{\Vb}{{\bf{V}}_n}
\newcommand{\Vbhat}{{\bf{\widehat{V}}}_n}
\newcommand{\Omegab}{{\bf{\Omega}}}
\newcommand{\gb}{{\bf{g}}}
\newcommand{\rb}{{\bf{r}}}
\newcommand{\pb}{p_\mathrm{base}}
\newcommand{\epsdot}{\dot{\epsilon}}
\newcommand{\qburn}{q_\mathrm{burn}}
\newcommand{\rt}{\tilde{r}_0}
\newcommand{\nablab}{\mathbf{\nabla}}
\newcommand{\dt}{\Delta\ t}
\newcommand{\omegadot}{\dot{\omega}}
\newcommand{\Hext}{H_{\rm ext}}
\newcommand{\Hnuc}{H_{\rm nuc}}
\newcommand{\kth}{k_{\rm th}}
\newcommand{\Gammaonebar}{\overline{\Gamma}_1}
\newcommand{\Sbar}{\overline{S}}
\newcommand{\etarho}{\eta_\rho}
\newcommand{\etarhoh}{\eta_{\rho~h}}
\newcommand{\Ubt}{\widetilde{\Ub}}
\newcommand{\wt}{\widetilde{w}}
\newcommand{\He}{$^4$He}
\newcommand{\C}{$^{12}$C}
\newcommand{\Fe}{$^{56}$Fe}
\newcommand{\isot}[2]{$^{#2}\mathrm{#1}$}
\newcommand{\isotm}[2]{{}^{#2}\mathrm{#1}}
\newcommand{\maestro}{{\sf MAESTRO}}
\newcommand{\castro}{{\sf Castro}}
\newcommand{\amrex}{{\sf AMReX}}
\newcommand{\pynucastro}{{\sf pynucastro}}
\newcommand{\avg}[1]{\overline{#1}}
\newcommand{\avgtwod}[1]{\langle~#1 \rangle}
\newcommand{\rms}[2]{\left(\delta#1\right)_{r_{#2}}}
\newcommand{\mymax}[1]{\left(#1\right)_{\rm max}}
\newcommand{\Tb}{\ensuremath{T_\mathrm{base}}}
\newcommand{\gcc}{\mathrm{g~cm^{-3} }}
\newcommand{\cms}{\mathrm{cm~s^{-1} }}
\newcommand{\half}{\frac{1}{2}}
\setlength{\marginparwidth}{0.5in}
\newcommand{\MarginPar}[1]{
\marginpar{\vskip-\baselineskip%
\raggedright%
\tiny\sffamily%
{\color{red}\hrule%
\smallskip%
#1\par%
\smallskip%
\hrule}}%
}
\newcommand{\AssignTo}[1]{
\marginpar{\vskip-\baselineskip%
\raggedright%
\tiny\sffamily%
{\color{blue}\hrule%
\smallskip%
#1\par%
\smallskip%
\hrule}}%
}
\begin{document}
%======================================================================
% Title
%======================================================================
\title{Pynucastro Methods and Verification}
\shorttitle{Pynucastro}
\shortauthors{Pynucastro Development Team}
\author[0000-0000-0000-0000]{Pynucastro Development Team}
\affiliation{Planet Earth}
\correspondingauthor{Pynucastro Development Team}
\email{weshouldprobablyhaveanemail@somewhere}
%======================================================================
% Abstract and Keywords
%======================================================================
\begin{abstract}
We present pynucastro, a software tool for interacting with nuclear
reaction rate databases for pedagogical and research purposes. We
also present verification tests for pynucastro as well as code
comparisons with other community codes for integrating general
reaction networks.
\end{abstract}
\keywords{X-ray bursts (1814), Nucleosynthesis (1131), Open source software (1866), Computational methods (1965)}
%======================================================================
% Introduction
%======================================================================
\section{Introduction}\label{Sec:Introduction}
\begin{itemize}
\item Outline existing need for a code like this
\item Review previous community codes for reaction networks
\item Describe what we are doing differently (vis, code generation, ...)
\item Introduce the verification and code comparison problems we will talk about later
\item Discuss the big picture: what data we turn into what ODEs for what other codes (Castro, MAESTROeX)
\item Sources of data for pynucastro (Reaclib, tabulated weak rates, CODATA, Atomic Mass Evaluation, ...)
\end{itemize}
%======================================================================
% Software Design
%======================================================================
\section{Software Design}\label{Sec:Design}
\begin{itemize}
\item Overview of the structure of the python code comprising pynucastro
\item Present visualization capabilities with Jupyter notebooks
\item Discuss ODE construction with Sympy
\item Discuss code template processing system for code generation
\item Discuss different code outputs (python networks, StarKiller Microphysics networks)
\item Discuss CI approach and what parts of pynucastro we test
\end{itemize}
%======================================================================
% Verification
%======================================================================
\section{Simple Verification Tests}\label{Sec:Verification}
\begin{itemize}
\item Find a way to verify sub-components of pynucastro
\item Present simple tests where we could get the answer analytically?
\item At least one test for each Reaclib chapter, plus tabulated rates?
\end{itemize}
%======================================================================
% Code Comparisons
%======================================================================
\section{Code Comparisons}\label{Sec:Comparisons}
\begin{itemize}
\item Describe other codes for the comparison (XNet, Skynet, MESA?)
\item Outline similarities and differences in the microphysics and integration strategy for the codes
\item Outline real problems we will compare with other codes using a range of networks
\item XRBs? (one zone burns, thermodynamic trajectories from multi-D flames)
\item Discuss performance?
\item Discuss ease of use/python vs other APIs?
\end{itemize}
%======================================================================
% Conclusions/Outlook
%======================================================================
\section{Conclusions and Outlook}\label{Sec:Conclusions}
Conclude something, invite people to create issues/submit PRs.
%======================================================================
% References
%======================================================================
\bibliographystyle{aasjournal}
\bibliography{ws}
\end{document}