]> https://gitweb.dealii.org/ - dealii.git/commitdiff
find missing brackets back 7212/head
authorJiaxin Wang <jiaxin@Jiaxins-MacBook-Pro.local>
Tue, 30 Oct 2018 15:10:07 +0000 (16:10 +0100)
committerJiaxin Wang <jiaxin@Jiaxins-MacBook-Pro.local>
Tue, 30 Oct 2018 15:10:07 +0000 (16:10 +0100)
examples/step-30/step-30.cc

index 77c24d9e1b5227ebc0e7555aff9b69b882fd003c..7771341ce969a135e84498c5824a62adb909e57d 100644 (file)
@@ -571,47 +571,48 @@ namespace Step30
                     // cell with the higher level does the work. In case that
                     // both cells have the same level, the cell with lower index
                     // is selected.
-                    if (dim>1 && cell->neighbor_is_coarser(face_no) continue;
-                        if (neighbor->level() < cell->level() ||
-                              (neighbor->index() > cell->index() &&
-                               neighbor->level() == cell->level()))
+                    if (dim > 1 && cell->neighbor_is_coarser(face_no))
+                      continue;
+                    if (neighbor->level() < cell->level() ||
+                        (neighbor->index() > cell->index() &&
+                         neighbor->level() == cell->level()))
                       {
-                      // Here we know, that the neighbor is not coarser so we
-                      // can use the usual @p neighbor_of_neighbor
-                      // function. However, we could also use the more
-                      // general @p neighbor_face_no function.
-                      const unsigned int neighbor2 =
-                        cell->neighbor_of_neighbor(face_no);
-
-                      ue_vi_matrix = 0;
-                      ui_ve_matrix = 0;
-                      ue_ve_matrix = 0;
-
-                      fe_v_face.reinit(cell, face_no);
-                      fe_v_face_neighbor.reinit(neighbor, neighbor2);
-
-                      dg.assemble_face_term2(fe_v_face,
-                                             fe_v_face_neighbor,
-                                             ui_vi_matrix,
-                                             ue_vi_matrix,
-                                             ui_ve_matrix,
-                                             ue_ve_matrix);
-
-                      neighbor->get_dof_indices(dofs_neighbor);
-
-                      for (unsigned int i = 0; i < dofs_per_cell; ++i)
-                        for (unsigned int j = 0; j < dofs_per_cell; ++j)
-                          {
-                            system_matrix.add(dofs[i],
-                                              dofs_neighbor[j],
-                                              ue_vi_matrix(i, j));
-                            system_matrix.add(dofs_neighbor[i],
-                                              dofs[j],
-                                              ui_ve_matrix(i, j));
-                            system_matrix.add(dofs_neighbor[i],
-                                              dofs_neighbor[j],
-                                              ue_ve_matrix(i, j));
-                          }
+                        // Here we know, that the neighbor is not coarser so we
+                        // can use the usual @p neighbor_of_neighbor
+                        // function. However, we could also use the more
+                        // general @p neighbor_face_no function.
+                        const unsigned int neighbor2 =
+                          cell->neighbor_of_neighbor(face_no);
+
+                        ue_vi_matrix = 0;
+                        ui_ve_matrix = 0;
+                        ue_ve_matrix = 0;
+
+                        fe_v_face.reinit(cell, face_no);
+                        fe_v_face_neighbor.reinit(neighbor, neighbor2);
+
+                        dg.assemble_face_term2(fe_v_face,
+                                               fe_v_face_neighbor,
+                                               ui_vi_matrix,
+                                               ue_vi_matrix,
+                                               ui_ve_matrix,
+                                               ue_ve_matrix);
+
+                        neighbor->get_dof_indices(dofs_neighbor);
+
+                        for (unsigned int i = 0; i < dofs_per_cell; ++i)
+                          for (unsigned int j = 0; j < dofs_per_cell; ++j)
+                            {
+                              system_matrix.add(dofs[i],
+                                                dofs_neighbor[j],
+                                                ue_vi_matrix(i, j));
+                              system_matrix.add(dofs_neighbor[i],
+                                                dofs[j],
+                                                ui_ve_matrix(i, j));
+                              system_matrix.add(dofs_neighbor[i],
+                                                dofs_neighbor[j],
+                                                ue_ve_matrix(i, j));
+                            }
                       }
 
                     // We do not need to consider case d), as those faces are

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.