Commits April 2012

commits@lists.einsteintoolkit.org

10 participants
167 discussions

[svn:einsteintoolkit] Workshop_Spring_2012/ (Rev. 49)
by bcmsma＠astro.rit.edu 04 Apr '12

04 Apr '12

User: bmundim Date: 2012/04/04 11:54 AM Added: /numerical_relativity/movies/ DiagShock2D_rho.mpg, DiagShock2D_rho_570.jpg Modified: /numerical_relativity/ numerical_relativity.tex Log: Shocktube test movie Directory Changes: Directory: /svn:mime-type/ ========================== + application/octet-stream File Changes: Directory: /numerical_relativity/movies/ ======================================== File [added]: DiagShock2D_rho.mpg Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: numerical_relativity/movies/DiagShock2D_rho.mpg ___________________________________________________________________ File [added]: DiagShock2D_rho_570.jpg Delta lines: +0 -0 =================================================================== (Binary files differ) Property changes on: numerical_relativity/movies/DiagShock2D_rho_570.jpg ___________________________________________________________________ Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +11 -0 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 16:45:17 UTC (rev 48) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 16:54:47 UTC (rev 49) @@ -685,7 +685,18 @@ \begin{figure} \includegraphics[width=10cm]{all_waves.pdf} \end{figure} +} +\frame{\frametitle{Shocktube Test} +\begin{figure} +% \hspace{1.0cm} + \href{run:movies/DiagShock2D_rho.mpg}{ + \includegraphics[width=6cm]{movies/DiagShock2D_rho_570.jpg} + } + \caption{Diagonal shock along xy-plane: final $\rho$ profile.} +\end{figure} + + }

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 48)
by bcmsma＠astro.rit.edu 04 Apr '12

04 Apr '12

User: bmundim Date: 2012/04/04 11:45 AM Modified: /numerical_relativity/ numerical_relativity.tex Log: algorithm overview. File Changes: Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +32 -0 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 16:29:55 UTC (rev 47) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 16:45:17 UTC (rev 48) @@ -658,4 +658,36 @@ \end{figure} } + +\frame{\frametitle{Algorithm Overview} +\begin{itemize} +\item Set boundary conditions or ghost zone values \pause +\item Calculate numerical flux at every cell interface by +solving an approximate Riemann problem \pause + \begin{itemize} + \item Reconstruct the primitive variables $p_l$ and $p_r$ at the +cell interfaces (by using slope limiters for example). \pause + \item Find the conservative variables $q_l$ and $q_r$ by using their +definition \pause + \item Calculate the fluxes $f(q_l)$ and $f(q_r)$ \pause + \item Calculate $A$, $r_{\alpha}$, $\lambda_{\alpha}$ and use them to +find $\omega_{\alpha}$ \pause + \item Calculate the Roe formula for the numerical flux: $F^{\rm Roe}_{i+1/2} +$ \pause + \end{itemize} + \item Evaluate the sources $\Psi(q)$. \pause + \item Update the average quantities $Q^n_i$ using method of lines. +\end{itemize} + +} + +\frame{\frametitle{Algorithm Overview} +\begin{figure} + \includegraphics[width=10cm]{all_waves.pdf} +\end{figure} + +} + + + \end{document}

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 47)
by bcmsma＠astro.rit.edu 04 Apr '12

04 Apr '12

User: bmundim Date: 2012/04/04 11:29 AM Added: /numerical_relativity/figs/ minmod.eps, minmod.pdf Modified: /numerical_relativity/ numerical_relativity.tex Log: slope limiters File Changes: Directory: /numerical_relativity/figs/ ====================================== File [added]: minmod.eps Delta lines: <Unavailable> File [added]: minmod.pdf Delta lines: <Unavailable> Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 46)
by bcmsma＠astro.rit.edu 04 Apr '12

04 Apr '12

User: bmundim Date: 2012/04/04 11:13 AM Modified: /numerical_relativity/ numerical_relativity.tex Log: approximate riemann solvers, an example. File Changes: Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +24 -0 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 15:55:15 UTC (rev 45) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 16:13:46 UTC (rev 46) @@ -609,7 +609,31 @@ \includegraphics[width=6cm]{rarefaction.pdf} \end{figure} +} +\frame{\frametitle{Approximate Riemann Solvers} +Usually the exact solution of the Riemann problem is computationally +very expensive. \pause + +Fortunately we can obtain very good approximation for the solutions +by approximating the conservation law as a quasi-linear system: +\begin{equation} +q_t + A q_x = 0 +\end{equation} \pause +where $A$ is a diagonalizable matrix given by: +\begin{equation} +A(q_l,q_r) = \left. \frac{\partial f}{\partial q} \right|_{q=1/2(q_l+q_R)} +\end{equation} \pause +Roe solver: +\begin{equation} +F^{\rm Roe}_{i+1/2}=\frac{1}{2} \left[f(q^r_{i+1/2})+f(q^l_{i-1/2}) +-\sum_{\alpha} |\lambda_{\alpha}| \omega_{\alpha} r_{\alpha} \right] +\end{equation} +where $\lambda_{\alpha}$ are the characteristics speeds, $\omega_{\alpha}$ +the jumps in the characteristics and $r_{\alpha}$ the right eigenvector +of $A$. } + + \end{document}

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 45)
by bcmsma＠astro.rit.edu 04 Apr '12

04 Apr '12

User: bmundim Date: 2012/04/04 10:55 AM Modified: /numerical_relativity/ numerical_relativity.tex Log: burgers eqn File Changes: Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +21 -0 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 15:40:22 UTC (rev 44) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 15:55:15 UTC (rev 45) @@ -591,4 +591,25 @@ } +\frame{\frametitle{Example: Burgers Equation} +\begin{equation} +u_t + f(u)_x = 0, \qquad f(u)=\frac{1}{2}u^2 +\end{equation} +with +\[ +u(x,0) = \left\{ \begin{array}{ll} + u_l & \mbox{if $x<0$}\\ + u_r & \mbox{if $x>0$} + \end{array} + \right. +\]\pause + +\begin{figure} + \includegraphics[width=6cm]{shock.pdf} + \includegraphics[width=6cm]{rarefaction.pdf} +\end{figure} + + +} + \end{document}

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 44)
by bcmsma＠astro.rit.edu 04 Apr '12

04 Apr '12

User: bmundim Date: 2012/04/04 10:40 AM Added: /numerical_relativity/figs/ riemann.eps, riemann.pdf Modified: /numerical_relativity/ Makefile, numerical_relativity.tex Log: update... not ready yet... File Changes: Directory: /numerical_relativity/figs/ ====================================== File [added]: riemann.eps Delta lines: <Unavailable> File [added]: riemann.pdf Delta lines: <Unavailable> Directory: /numerical_relativity/ ================================= File [modified]: Makefile Delta lines: <Unavailable> File [modified]: numerical_relativity.tex Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/example_wavetoy/Cactus/misc/ (Rev. 43)
by roland.haas＠physics.gatech.edu 04 Apr '12

04 Apr '12

User: rhaas Date: 2012/04/04 10:32 AM Modified: /example_wavetoy/Cactus/misc/ WaveToy.par Log: keep jpeg files around in wavetoy example par file File Changes: Directory: /example_wavetoy/Cactus/misc/ ======================================== File [modified]: WaveToy.par Delta lines: +1 -1 =================================================================== --- example_wavetoy/Cactus/misc/WaveToy.par 2012-04-04 15:24:28 UTC (rev 42) +++ example_wavetoy/Cactus/misc/WaveToy.par 2012-04-04 15:32:35 UTC (rev 43) @@ -26,7 +26,7 @@ IOJpeg::gridpoints = "interpolate" IOJpeg::out_every = 5 IOJpeg::refinement_factor = 10 -IOJpeg::mode = "remove" +IOJpeg::mode = "standard" IOJpeg::colormap = "custom" IOJpeg::colormap_min = -0.25 IOJpeg::colormap_max = +0.25

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/example_wavetoy/Cactus/misc/ (Rev. 42)
by roland.haas＠physics.gatech.edu 04 Apr '12

04 Apr '12

User: rhaas Date: 2012/04/04 10:24 AM Modified: /example_wavetoy/Cactus/misc/ WaveToy.th Log: update WaveToy thorn list File Changes: Directory: /example_wavetoy/Cactus/misc/ ======================================== File [modified]: WaveToy.th Delta lines: +5 -0 =================================================================== --- example_wavetoy/Cactus/misc/WaveToy.th 2012-04-04 14:46:16 UTC (rev 41) +++ example_wavetoy/Cactus/misc/WaveToy.th 2012-04-04 15:24:28 UTC (rev 42) @@ -117,6 +117,11 @@ !CHECKOUT = Introduction/WaveToy +!TARGET = $ROOT +!TYPE = svn +!URL = https://svn.einsteintoolkit.org/documents/Workshop_Spring_2012/example_wave… +!CHECKOUT = wavetoy-misc + # Private thorns !TARGET = $ARR !TYPE = ignore

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/numerical_relativity/ (Rev. 41)
by diener＠cct.lsu.edu 04 Apr '12

04 Apr '12

User: diener Date: 2012/04/04 09:46 AM Modified: /numerical_relativity/ numerical_relativity.tex Log: Add command /mypause to easily turn on or off all pauses. Change the date. Remove the Dr's. File Changes: Directory: /numerical_relativity/ ================================= File [modified]: numerical_relativity.tex Delta lines: +51 -49 =================================================================== --- numerical_relativity/numerical_relativity.tex 2012-04-04 14:42:26 UTC (rev 40) +++ numerical_relativity/numerical_relativity.tex 2012-04-04 14:46:16 UTC (rev 41) @@ -3,12 +3,14 @@ \usecolortheme[RGB={200,200,200}]{structure} \subtitle[Introduction to numerical methods \& numerical Relativity]{{\large Einstein Toolkit Workshop}\\*[0.3em] Numerical methods for spacetime and matter evolution} -\author[\mbox{Peter Diener \& Bruno Mundim}]{Dr Peter Diener \& Dr Bruno Mundim} -\date{Apr 05 2012} +\author[\mbox{Peter Diener \& Bruno Mundim}]{Peter Diener \& Bruno Mundim} +\date{Apr 4 2012} \begin{document} \newcommand{\code}[1]{\texttt{#1}} +%\newcommand{\mypause}{\pause} +\newcommand{\mypause}{} \frame{\titlepage} @@ -34,11 +36,11 @@ \begin{align} {\cal H} & \equiv R + K^2 - K_{ij} K^{ij} = 0, \\ {\cal M}^i & \equiv D_j (K^{ij} - \gamma^{ij} K) = 0. -\end{align}\pause -These equations were used a lot in the early years of numerical relativity.\pause +\end{align}\mypause +These equations were used a lot in the early years of numerical relativity.\mypause However, it was later discovered that the ADM evolution equations are only -weakly hyperbolic.\pause +weakly hyperbolic.\mypause In the mid to late 90's a new formulation was introduced that proved to be much more robust and stable: The BSSN formulation. @@ -49,18 +51,18 @@ Introduce a conformal rescaling of the three metric\vspace{\vs} \begin{equation} \gamma_{ij} = \psi^4 \tg_{ij}.\vspace{\vs} -\end{equation}\pause +\end{equation}\mypause We choose $\psi = \gamma^{1/12}$ such that the determinant of $\tg_{ij}$ is 1%\vspace{\vs} \begin{equation} \mathrm{det}(\tg_{ij}) = \mathrm{det}(\psi^{-4}\gamma_{ij}) = \mathrm{det}(\gamma^{-1/3}\gamma_{ij}) = \gamma^{-1} \mathrm{det}(\gamma_{ij}) = 1.\vspace{\vs} -\end{equation}\pause +\end{equation}\mypause In addition we introduce a trace decomposition of the extrinsic curvature.\vspace{\vs} \begin{align} K & = \gamma^{ij} K_{ij}, \\ A_{ij} & = K_{ij} - \frac{1}{3} \gamma_{ij} K.\vspace{\vs} -\end{align}\pause +\end{align}\mypause We then promote the following variables to evolution variables\vspace{\vs} \begin{align} \phi &= \ln \psi = \frac{1}{12} \ln \gamma, @@ -78,17 +80,17 @@ \begin{equation} \tG^i = \tg^{jk} \tG^i{}_{jk} = - \partial_j \tg^{ij}, \end{equation} -to evolved variables as well.\pause +to evolved variables as well.\mypause The final set of evolution variables are $\phi$, $K$, $\tg_{ij}$, $\tA_{ij}$ -and $\tG^i$.\pause +and $\tG^i$.\mypause The BSSN evolution equations can be derived from the ADM equations. As an example take the equation for $\phi$. \begin{equation} \dt \phi = \partial_0\phi = \partial_0 \left (\frac{1}{12}\ln\gamma\right) = \frac{1}{12}\frac{1}{\gamma}\partial_0\gamma. -\end{equation}\pause +\end{equation}\mypause Using the expression for the derivative of the determinant of the metric in terms of the derivatives of the metric ($\partial_0\gamma = \gamma \gamma^{ij} \partial_0\gamma_{ij}$) we find @@ -144,18 +146,18 @@ \tilde{{\cal G}} & \equiv \tg-1 = 0, \\ \tilde{{\cal A}} & \equiv \tg^{ij}\tA_{ij} = 0, \\ \tilde{{\cal L}}^{i} & \equiv \tG^{i}+\partial_j\tg^{ij} = 0. -\end{align}\pause +\end{align}\mypause The constraints $\tilde{{\cal G}}$ and $\tilde{{\cal A}}$ are enforced actively -at each timestep.\pause +at each timestep.\mypause The other constraints ($\tilde{{\cal H}}$, $\tilde{{\cal M}}^i$ and -$\tilde{{\cal L}}^{i}$) are not enforced.\pause +$\tilde{{\cal L}}^{i}$) are not enforced.\mypause To improve stability and to help maintain $\tilde{{\cal L}}^{i}$ at a low level the following rule is employed in an implementation \begin{itemize} \item Where derivatives of $\tG^{i}$ are needed the evolved $\tG^{i}$ are used - directly.\pause + directly.\mypause \item Where $\tG^{i}$ are needed without taking derivatives $\tg^{jk}\tG^{i}{}_{jk}$ are used instead. \end{itemize} @@ -163,7 +165,7 @@ \frame{\frametitle{The BSSN formulation (continued)} In order to evolve a spacetime with the BSSN equations, you have to specify -the gauges $\alpha$ and $\beta^{i}$.\pause +the gauges $\alpha$ and $\beta^{i}$.\mypause Most codes use the ``moving puncture'' gauges \begin{align} @@ -181,12 +183,12 @@ \gamma^{\mathrm{ph}}_{ij} = \psi^4\gamma_{ab}, \mbox{\hspace{2em}} K^{\mathrm{ph}}_{ij} = \psi^{-2} K_{ij}, \end{equation} where $\gamma_{ij}$ is chosen to be the flat metric and $K_{ij}$ is assumed -to be tracefree.\pause \\ +to be tracefree.\mypause \\ The constraint equations become\vspace{\vs} \begin{align} 0 & = \Delta\psi +\frac{1}{8}K^{ij}K_{ij}\psi^{-7} \\ 0 & = D_j K^{ij}. -\end{align}\pause +\end{align}\mypause The momentum constraint has an analytic solution\vspace{\vs} \begin{align} K^{ij}_{\mathrm{BY}} = {} & \frac{3}{2r^2}\left ( P^i n^j + P^j n^i - @@ -203,29 +205,29 @@ black holes. \begin{equation} \psi = 1+\sum_{i=1}^N\frac{m_i}{2|\vec{r}-\vec{r}_i|}. -\end{equation}\pause +\end{equation}\mypause Inspired by this, for $K_{ij}\ne 0$ we make the ansatz \begin{equation} \psi=\frac{1}{\alpha}+u, \mbox{\hspace{2em}} \frac{1}{\alpha} = \sum_{i=1}^N\frac{m_i}{2|\vec{r}-\vec{r}_i|}, -\end{equation}\pause +\end{equation}\mypause In which case the Hamiltonian constraint becomes an equation for $u$ \begin{equation} \Delta u+\frac{1}{8}\alpha^7 K^{ij}K_{ij}(1+\alpha u)^{-7} = 0. -\end{equation}\pause +\end{equation}\mypause It can be shown that $u$ will be $C^2$ at the location of the `punctures' and that $\psi$ has a unique solution. } \frame{\frametitle{Finite differencing} With finite differencing we discretize a function by sampling it at a -collection of grid points.\pause +collection of grid points.\mypause -The grid points are usually (but not necessarily) equally spaced. \pause +The grid points are usually (but not necessarily) equally spaced. \mypause We can then approximate derivatives of a function at a grid point by a weigthed sum of function values at grid points in the neighbourhood (the -stencil) of the grid point. \pause +stencil) of the grid point. \mypause As an example consider a stencil containing the grid point ($f_i$) and it's two nearest neighbors ($f_{i-1}$ and $f_{i+1}$) with @@ -245,7 +247,7 @@ f(x_{i+1}) & = & f(x_i) + \left. \frac{df}{dx}\right |_{x_i} \Delta x + \frac{1}{2} \left. \frac{d^2 f}{dx^2}\right |_{x_i} \Delta x^2 + O((\Delta x)^3) \nonumber -\end{eqnarray} \pause +\end{eqnarray} \mypause Requiring that the weighted sum approximates the derivative yields the following equations for $a_{-1}$, $a_0$ and $a_1$ \begin{eqnarray} @@ -261,20 +263,20 @@ \begin{equation} \left. \frac{df}{dx} \right |_{x_i} = \frac{f_{i+1}-f_{i-1}}{2\Delta x} + O((\Delta x)^2). -\end{equation}\pause +\end{equation}\mypause Similarly we find for the second derivative that \begin{equation} \left. \frac{d^2f}{dx^2} \right |_{x_i} = \frac{f_{i+1}-2 f_i+f_{i-1}}{(\Delta x)^2} + O((\Delta x)^2). -\end{equation}\pause -These finite difference operators are second order accurate.\pause +\end{equation}\mypause +These finite difference operators are second order accurate.\mypause -Higher order accuracy or higher order derivatives require larger stencils.\pause +Higher order accuracy or higher order derivatives require larger stencils.\mypause Another way of looking at finite differencing operators is through -interpolating polynomials.\pause +interpolating polynomials.\mypause -Either approach gives the same coefficients for the same stencil.\pause +Either approach gives the same coefficients for the same stencil.\mypause It is also clear from either approach that the error estimates are only correct if $f$ is smooth enough. @@ -285,20 +287,20 @@ \begin{equation} \partial_t \mathbf{q}+\mathbf{A}^i(\mathbf{q})\partial_i \mathbf{B}(\mathbf{q}) = \mathbf{S}(\mathbf{q}). -\end{equation}\pause +\end{equation}\mypause The idea then is to discretize in space first, i.e.\ write the equations as \begin{equation} \partial_t \mathbf{q} = \mathbf{L}(\mathbf{q}), \end{equation} where $\mathbf{L}(\mathbf{q})$ is a discrete approximation to the equations -(e.g.\ using finite differencing).\pause +(e.g.\ using finite differencing).\mypause This then turns the equations into a set of coupled ODE's with respect to -time.\pause +time.\mypause If the spatial discretization (including boundary conditions) is stable we we can then evolve the system of equations using any stable ODE time -integrator.\pause +integrator.\mypause Often Runge-Kutta integration schemes are the scheme of choice. } @@ -307,37 +309,37 @@ Advantages of using the method of lines: \begin{itemize} \item It is easy to change the time integration scheme (e.g.\ going to higher - order).\pause + order).\mypause \item It is easy to couple different evolution equations maintaining - high order coupling.\pause + high order coupling.\mypause \end{itemize} The method of lines is implemented in Cactus in the thorn -\code{CactusNumerical/MoL}.\pause +\code{CactusNumerical/MoL}.\mypause To use: \begin{itemize} \item Schedule a routine to register the evolution variables and RHS - variables with \code{MoL}.\pause - \item Schedule routines to set the RHS variables.\pause + variables with \code{MoL}.\mypause + \item Schedule routines to set the RHS variables.\mypause \item Set parameters when launching job to choose the time integration scheme. \end{itemize} } \frame{\frametitle{McLachlan} -\code{McLachlan} is the Einstein Toolkit implementation of BSSN.\pause +\code{McLachlan} is the Einstein Toolkit implementation of BSSN.\mypause \code{McLachlan} is named in honor of the Canadian Singer/Song writer Sarah -McLachlan.\pause +McLachlan.\mypause -It is based on finite differencing and the method of lines.\pause +It is based on finite differencing and the method of lines.\mypause -It supports high order finite differencing (8th order).\pause +It supports high order finite differencing (8th order).\mypause Since it is generated from the tensor equations by Kranc it is easy to -maintain and modify.\pause +maintain and modify.\mypause It is highly optimized (supports vectorization and OpenMP) and an effort is -ongoing to make it be able to run on GPU's as well.\pause +ongoing to make it be able to run on GPU's as well.\mypause In the future it may be necessary to generate different versions optimized for specific computer architectures. @@ -348,7 +350,7 @@ Once we have established the spacetime evolution and initial data equations, we need to obtain the evolution equations for the matter fields and the magnetic field evolution equation in -ideal MHD case. \pause +ideal MHD case. \mypause These equations can be expressed as the local conservation laws of baryon number and energy momentum. For baryon number we have: @@ -357,7 +359,7 @@ \nabla_\nu J^\nu=0, \end{equation} where $J^\mu=\rho u^\mu$ is the rest-mass current, $\rho$ the rest-mass -density and $u^\mu$ is the four-velocity of a fluid comoving observer.\pause +density and $u^\mu$ is the four-velocity of a fluid comoving observer.\mypause The conservation of energy-momentum is given by: % @@ -385,7 +387,7 @@ \end{equation} where we define $p^*=p+b^2/2$, $h^*=h+b^2/\rho$, $\varepsilon^*=\varepsilon+b^2/(2\rho)$, that results into -$h^*=1+\varepsilon^*+p^*/\rho$. \pause +$h^*=1+\varepsilon^*+p^*/\rho$. \mypause Note that in the expression above we have used the magnetic field $b^{\mu}$ as measured by the comoving observer. It can expressed

1 0

[svn:einsteintoolkit] Workshop_Spring_2012/ (Rev. 40)
by knarf＠cct.lsu.edu 04 Apr '12

04 Apr '12

User: knarf Date: 2012/04/04 09:42 AM Modified: /cactuswriting/ cactuswriting.pdf /genrel/ genrel.pdf Log: change mime type of pdfs Directory Changes: Directory: /svn:mime-type/ ========================== - application/octet-stream + application/pdf File Changes: Directory: /cactuswriting/ ========================== File [modified]: cactuswriting.pdf Delta lines: +0 -0 =================================================================== Property changes on: cactuswriting/cactuswriting.pdf ___________________________________________________________________ Directory: /genrel/ =================== File [modified]: genrel.pdf Delta lines: +0 -0 =================================================================== Property changes on: genrel/genrel.pdf ___________________________________________________________________

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits April 2012