]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Responses to review comments
authorLei Qiao <qiaol618@gmail.com>
Thu, 26 Feb 2015 18:35:10 +0000 (12:35 -0600)
committerLei Qiao <qiaol618@gmail.com>
Thu, 26 Feb 2015 18:35:10 +0000 (12:35 -0600)
examples/step-33/doc/intro.dox
examples/step-33/doc/results.dox
examples/step-33/step-33.cc

index b534e51c8ad442e7c353f421d3c61a67aca3fc27..c0bf755195fd9bda27bcc3232f4f8256b081f0ef 100644 (file)
@@ -128,7 +128,7 @@ Hartmann's PhD thesis ("Adaptive Finite Element Methods for the
 Compressible Euler Equations", PhD thesis, University of Heidelberg, 2002).
 
 We use a time stepping scheme to substitute the time derivative in the
-above equations. For simplicity, we define $ \mathbf{B}({\mathbf{w}_{n})(\mathbf z) $ as the spacial residual at time step $n$ :
+above equations. For simplicity, we define $ \mathbf{B}({\mathbf{w}_{n})(\mathbf z) $ as the spatial residual at time step $n$ :
 
 @f{eqnarray*}
  \mathbf{B}({\mathbf{w}_{n})(\mathbf z)  &=&
index 8348e431a04d0fe4d9d89de19e7c6b41d9bfe83f..03b5889565bf2640cf38d15a139fcd54a7c0f630 100644 (file)
@@ -261,15 +261,15 @@ faster.
 <h4>Cache the explicit part of residual</h4>
 
 The residual calulated in ConservationLaw::assemble_cell_term function
-read
+read
    $R_i = \left(\frac{\mathbf{w}^{k}_{n+1} - \mathbf{w}_n}{\delta t}
     , \mathbf{z}_i \right)_K  +
       \theta \mathbf{B}({\mathbf{w}^{k}_{n+1})(\mathbf{z}_i)_K +
       (1-\theta) \mathbf{B}({\mathbf{w}_{n}) (\mathbf{z}_i)_K $
-means that we calculate the spacial residual twice at one Newton
+This means that we calculate the spacial residual twice at one Newton
 iteration step: once respect to the current solution ${\mathbf{w}^{k}_{n+1}$
-and another respect to the last time step solution $\mathbf{w}_{n}$ which
-remains the same during all Newton interations through one timestep.
+and once more respect to the last time step solution $\mathbf{w}_{n}$ which
+remains the same during all Newton iterations through one timestep.
 Cache up the explicit part of residual
  $ \mathbf{B}({\mathbf{w}_{n}) (\mathbf{z}_i)_K}$
 during Newton iteration will save lots of labor.
index afd3e12820287e3c12a36f646b663624caa8d616..8f9a6d29fd21e5a8e000908687cb25956ba9acc3 100644 (file)
@@ -185,12 +185,12 @@ namespace Step33
     // to the <code>i</code>th element, and then dereference it. This works
     // for both kinds of vectors -- not the prettiest solution, but one that
     // works.
-    template <typename number, typename InputVector>
+    template <typename Number, typename InputVector>
     static
-    number
+    Number
     compute_kinetic_energy (const InputVector &W)
     {
-      number kinetic_energy = 0;
+      Number kinetic_energy = 0;
       for (unsigned int d=0; d<dim; ++d)
         kinetic_energy += *(W.begin()+first_momentum_component+d) *
                           *(W.begin()+first_momentum_component+d);
@@ -200,14 +200,14 @@ namespace Step33
     }
 
 
-    template <typename number, typename InputVector>
+    template <typename Number, typename InputVector>
     static
-    number
+    Number
     compute_pressure (const InputVector &W)
     {
       return ((gas_gamma-1.0) *
               (*(W.begin() + energy_component) -
-               compute_kinetic_energy<number>(W)));
+               compute_kinetic_energy<Number>(W)));
     }
 
 
@@ -227,15 +227,15 @@ namespace Step33
     // use the automatic differentiation type here.  Similarly, we will call
     // the function with different input vector data types, so we templatize
     // on it as well:
-    template <typename InputVector, typename number>
+    template <typename InputVector, typename Number>
     static
     void compute_flux_matrix (const InputVector &W,
-                              number (&flux)[n_components][dim])
+                              Number (&flux)[n_components][dim])
     {
       // First compute the pressure that appears in the flux matrix, and then
       // compute the first <code>dim</code> columns of the matrix that
       // correspond to the momentum terms:
-      const number pressure = compute_pressure<number> (W);
+      const Number pressure = compute_pressure<number> (W);
 
       for (unsigned int d=0; d<dim; ++d)
         {
@@ -266,16 +266,16 @@ namespace Step33
     // numerical flux function to enforce boundary conditions.  This routine
     // is the basic Lax-Friedrich's flux with a stabilization parameter
     // $\alpha$. It's form has also been given already in the introduction:
-    template <typename InputVector, typename number>
+    template <typename InputVector, typename Number>
     static
     void numerical_normal_flux (const Point<dim>          &normal,
                                 const InputVector         &Wplus,
                                 const InputVector         &Wminus,
                                 const double               alpha,
-                                number (&normal_flux)[n_components])
+                                Number (&normal_flux)[n_components])
     {
-      number iflux[n_components][dim];
-      number oflux[n_components][dim];
+      Number iflux[n_components][dim];
+      Number oflux[n_components][dim];
 
       compute_flux_matrix (Wplus, iflux);
       compute_flux_matrix (Wminus, oflux);
@@ -300,10 +300,10 @@ namespace Step33
     // \right)^T$, shown here for the 3d case. More specifically, we will
     // consider only $\mathbf g=(0,0,-1)^T$ in 3d, or $\mathbf g=(0,-1)^T$ in
     // 2d. This naturally leads to the following function:
-    template <typename InputVector, typename number>
+    template <typename InputVector, typename Number>
     static
     void compute_forcing_vector (const InputVector &W,
-                                 number (&forcing)[n_components])
+                                 Number (&forcing)[n_components])
     {
       const double gravity = -1.0;
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.