forked from jedbrown/talks
-
Notifications
You must be signed in to change notification settings - Fork 0
/
20111205-AGU.tex
223 lines (193 loc) · 7.59 KB
/
20111205-AGU.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
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
% \documentclass[handout]{beamer}
\documentclass{beamer}
\mode<presentation>
{
\usetheme{default}
\usefonttheme[onlymath]{serif}
% \usetheme{Singapore}
% \usetheme{Warsaw}
% \usetheme{Malmoe}
% \useinnertheme{circles}
% \useoutertheme{infolines}
% \useinnertheme{rounded}
\setbeamercovered{transparent=5}
}
\usepackage[english]{babel}
\usepackage[latin1]{inputenc}
\usepackage{textpos,alltt,listings,multirow,ulem,siunitx}
\newcommand\hmmax{0}
\newcommand\bmmax{0}
\usepackage{bm}
% font definitions, try \usepackage{ae} instead of the following
% three lines if you don't like this look
\usepackage{mathptmx}
\usepackage[scaled=.90]{helvet}
% \usepackage{courier}
\usepackage[T1]{fontenc}
\usepackage{tikz}
\usetikzlibrary[shapes,shapes.arrows,arrows,shapes.misc,fit,positioning]
% \usepackage{pgfpages}
% \pgfpagesuselayout{4 on 1}[a4paper,landscape,border shrink=5mm]
\usepackage{JedMacros}
\title{Utilizing Emerging Hardware for Multiphysics Simulation Through Implicit High-Order Finite Element Methods With Tensor Product Structure}
\author{Jed Brown\inst{1}, Aron Ahmadia\inst{2}, Matt Knepley\inst{3}, Barry Smith\inst{1}}
% - Use the \inst command only if there are several affiliations.
% - Keep it simple, no one is interested in your street address.
\institute
{
\inst{1}{Mathematics and Computer Science Division, Argonne National Laboratory} \\
\inst{2}{King Abdullah University of Science and Technology} \\
\inst{3}{Computation Institute, University of Chicago}
}
\date{2011-12-05}
% This is only inserted into the PDF information catalog. Can be left
% out.
\subject{Talks}
% If you have a file called "university-logo-filename.xxx", where xxx
% is a graphic format that can be processed by latex or pdflatex,
% resp., then you can add a logo as follows:
% \pgfdeclareimage[height=0.5cm]{university-logo}{university-logo-filename}
% \logo{\pgfuseimage{university-logo}}
% Delete this, if you do not want the table of contents to pop up at
% the beginning of each subsection:
% \AtBeginSubsection[]
% {
% \begin{frame}<beamer>
% \frametitle{Outline}
% \tableofcontents[currentsection,currentsubsection]
% \end{frame}
% }
\AtBeginSection[]
{
\begin{frame}<beamer>
\frametitle{Outline}
\tableofcontents[currentsection]
\end{frame}
}
% If you wish to uncover everything in a step-wise fashion, uncomment
% the following command:
% \beamerdefaultoverlayspecification{<+->}
\begin{document}
\lstset{language=C}
\normalem
\begin{frame}
\titlepage
\end{frame}
\begin{frame}{The Roadmap}
\begin{block}{Hardware trends}
\begin{itemize}
\item More cores (keep hearing $\bigO(1000)$ per node)
\item Long vector registers (already 32 bytes for AVX and BG/Q)
\item Must use SMT to hide memory latency
\item Must use SMT for floating point performance (GPU, BG/Q)
\item Large penalty for non-contiguous memory access
\end{itemize}
\end{block}
\begin{block}{``Free flops'', but how can we use them?}
\begin{itemize}
\item High order methods good: better accuracy per storage
\item High order methods bad: work unit gets larger
\item GPU threads have very little memory, must keep work unit small
\item Need library composability, keep user contribution embarrassingly parallel
\end{itemize}
\end{block}
\end{frame}
\begin{frame}{How to program this beast?}
\begin{itemize}
\item Decouple physics from discretization
\begin{itemize*}
\item Expose small, embarrassingly parallel operations to user
\item Library schedules user threads for reuse between kernels
\item User provides physics in kernels run at each quadrature point
\item Continuous weak form: find $u \in \VV_D$
\[ v^T F(u) \sim \int_\Omega v \cdot {\color{green!70!black} f_0(u,\nabla u)}
+ \nabla v \tcolon {\color{green!70!black} f_1(u,\nabla u)} = 0, \qquad \forall v \in \VV_0 \]
\item Similar form at faces, but may involve Riemann solve
\end{itemize*}
\item Library manages reductions
\begin{itemize*}
\item Interpolation and differentiation on elements
\item Exploit tensor product structure to keep working set small
\item Assembly into solution/residual vector (sum over elements)
\end{itemize*}
\end{itemize}
\end{frame}
\input{slides/Dohp/TensorFEM.tex}
\input{slides/Dohp/RepresentationOfJacobians.tex}
\newcommand\smallterm[1]{{\color{gray} #1}}
\begin{frame}{Conservative (non-Boussinesq) two-phase ice flow}
Find momentum density $\rho\uu$, pressure $p$, and total energy density $E$:
\begin{gather*}
(\rho\uu)_t + \div (\smallterm{\rho\uu\otimes\uu} - \eta D\uu_i + p\bm 1) - \rho \bm g = 0 \\
\rho_t + \div \rho\uu = 0 \\
E_t + \div \big((E+p)\uu - k_T\nabla T - k_\omega\nabla\omega \big) - \eta D\uu_i\tcolon D\uu_i - \smallterm{\rho\uu\cdot\bm g} = 0
\end{gather*}
\begin{itemize}
\item Solve for density $\rho$, ice velocity $\uu_i$, temperature $T$, and melt fraction $\omega$ using constitutive relations.
\begin{itemize}
\item Simplified constitutive relations can be solved explicitly.
\item Temperature, moisture, and strain-rate dependent rheology $\eta$.
\item High order FEM, typically $Q_3$ momentum \& energy, SUPG (yuck).
\end{itemize}
\item DAEs solved implicitly after semidiscretizing in space.
\item Preconditioning using nested fieldsplit
\end{itemize}
\end{frame}
\begin{frame}[fragile,shrink=5]{Traversal code}
\begin{itemize*}
\item CPU traversal computes coefficients of test functions, \url{https://github.com/jedbrown/dohp/}
\begin{ccode}
while (IteratorHasPatch(iter)) {
IteratorGetPatchApplied(iter,&Q,&jw,
&x,&dx,NULL,NULL,
&u,&du,&u_,&du_, &p,&dp,&p_,NULL, &e,&de,&e_,&de_);
IteratorGetStash(iter,NULL,&stash);
for (dInt i=0; i<Q; i++) {
PointwiseFunction(context,x[i],dx[i],jw[i],
u[i],du[i],p[i],dp[i],e[i],de[i],
&stash[i], u_[i],du_[i],p_[i],e_[i],de_[i]);
}
IteratorCommitPatchApplied(iter,INSERT_VALUES, NULL,NULL,
u_,du_, p_,NULL, e_,de_);
IteratorNextPatch(iter);
}
\end{ccode}
\item GPU version calls \cfunc|PointwiseFunction| directly.
\item Unassembled Jacobian application reuses \cverb|stash|
\begin{ccode}
PointwiseJacobian(context,&stash[i],dx[i],jw[i],
u[i],du[i],p[i],dp[i],e[i],de[i],
u_[i],du_[i],p_[i],e_[i],de_[i]);
\end{ccode}
\end{itemize*}
\end{frame}
\begin{frame}
\includegraphics[width=\textwidth]{figures/VHT/TopViewStreamline}
\end{frame}
\begin{frame}
\includegraphics[width=\textwidth]{figures/VHT/VHTJakoContourStream}
\end{frame}
\input{slides/Dohp/TensorVsAssembly.tex}
\input{slides/HardwareArithmeticIntensity.tex}
\begin{frame}{Outlook}
\begin{itemize}
\item Sparse matrix assembly (for preconditioning) not shown
\begin{itemize}
\item $> 100$ GF/s for lowest order Stokes (Matt Knepley)
\item common physics code with CPU implementation
\item Dohp CPU version faster than libMesh and Deal.II for $Q_1$
\item $Q_1$ assembly embedded in higher order is 8\% slower than hand-rolled
\end{itemize}
\item Can't wait for OpenCL to implement indirect function calls
\item Symbolic differentiation too slow, tired of hand-differentiation
\item I want source-transformation AD with indirect function calls
\item Find correct amount of reuse between face and cell integration
\item Riemann solves harder to vectorize
\item Finer grained parallelism in GPU tensor product kernels
\item Hide dispatch to pointwise kernels inside library
\begin{itemize}
\item Easy, but scary. Library/framework becomes \alert{\bf F}ramework.
\end{itemize}
\end{itemize}
\end{frame}
\end{document}