]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Use [][] instead of (.,.) for shape_values
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 8 Mar 2002 10:46:37 +0000 (10:46 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 8 Mar 2002 10:46:37 +0000 (10:46 +0000)
git-svn-id: https://svn.dealii.org/trunk@5553 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-12/step-12.cc

index b757801ff4252fa1064863c7d9a5f6737ef24dc1..399e4fbc4e6af58acaadf1edce3f572da69716bc 100644 (file)
@@ -295,10 +295,10 @@ void DGTransportEquation<dim>::assemble_cell_term(
       {
        for (unsigned int j=0; j<fe_v.dofs_per_cell; ++j)
          u_v_matrix(i,j) -= beta[point]*grad_v[i][point]*
-                             v(j,point) *
+                             v[j][point] *
                              JxW[point];
        
-       cell_vector(i) += rhs[point] *v(i,point) * JxW[point];
+       cell_vector(i) += rhs[point] * v[i][point] * JxW[point];
       }
 }
 
@@ -356,8 +356,8 @@ void DGTransportEquation<dim>::assemble_boundary_term(
        for (unsigned int i=0; i<fe_v.dofs_per_cell; ++i)
          for (unsigned int j=0; j<fe_v.dofs_per_cell; ++j)
            u_v_matrix(i,j) += beta_n *
-                              v(j,point) *
-                              v(i,point) *
+                              v[j][point] *
+                              v[i][point] *
                               JxW[point];
       else
                                         // and the term $(\beta\cdot
@@ -366,7 +366,7 @@ void DGTransportEquation<dim>::assemble_boundary_term(
        for (unsigned int i=0; i<fe_v.dofs_per_cell; ++i)
          cell_vector(i) -= beta_n *
                            g[point] *
-                           v(i,point) *
+                           v[i][point] *
                            JxW[point];
     }
 }
@@ -444,8 +444,8 @@ void DGTransportEquation<dim>::assemble_face_term1(
        for (unsigned int i=0; i<fe_v.dofs_per_cell; ++i)
          for (unsigned int j=0; j<fe_v.dofs_per_cell; ++j)
            u_v_matrix(i,j) += beta_n *
-                              v(j,point) *
-                              v(i,point) *
+                              v[j][point] *
+                              v[i][point] *
                               JxW[point];
       else
                                         // and the
@@ -455,8 +455,8 @@ void DGTransportEquation<dim>::assemble_face_term1(
        for (unsigned int i=0; i<fe_v.dofs_per_cell; ++i)
          for (unsigned int k=0; k<fe_v_neighbor.dofs_per_cell; ++k)
            un_v_matrix(i,k) += beta_n *
-                               v_neighbor(k,point) *
-                               v(i,point) *
+                               v_neighbor[k][point] *
+                               v[i][point] *
                                JxW[point];
     }
 }
@@ -506,8 +506,8 @@ void DGTransportEquation<dim>::assemble_face_term2(
          for (unsigned int i=0; i<fe_v.dofs_per_cell; ++i)
            for (unsigned int j=0; j<fe_v.dofs_per_cell; ++j)
              u_v_matrix(i,j) += beta_n *
-                                v(j,point) *
-                                v(i,point) *
+                                v[j][point] *
+                                v[i][point] *
                                 JxW[point];
 
                                           // We additionally assemble
@@ -517,8 +517,8 @@ void DGTransportEquation<dim>::assemble_face_term2(
          for (unsigned int k=0; k<fe_v_neighbor.dofs_per_cell; ++k)
            for (unsigned int j=0; j<fe_v.dofs_per_cell; ++j)
              u_vn_matrix(k,j) -= beta_n *
-                                 v(j,point) *
-                                 v_neighbor(k,point) *
+                                 v[j][point] *
+                                 v_neighbor[k][point] *
                                  JxW[point];
        }
       else
@@ -528,8 +528,8 @@ void DGTransportEquation<dim>::assemble_face_term2(
          for (unsigned int i=0; i<fe_v.dofs_per_cell; ++i)
            for (unsigned int l=0; l<fe_v_neighbor.dofs_per_cell; ++l)
              un_v_matrix(i,l) += beta_n *
-                                 v_neighbor(l,point) *
-                                 v(i,point) *
+                                 v_neighbor[l][point] *
+                                 v[i][point] *
                                  JxW[point];
 
                                           // And this is another new
@@ -539,8 +539,8 @@ void DGTransportEquation<dim>::assemble_face_term2(
          for (unsigned int k=0; k<fe_v_neighbor.dofs_per_cell; ++k)
            for (unsigned int l=0; l<fe_v_neighbor.dofs_per_cell; ++l)
              un_vn_matrix(k,l) -= beta_n *
-                                  v_neighbor(l,point) *
-                                  v_neighbor(k,point) *
+                                  v_neighbor[l][point] *
+                                  v_neighbor[k][point] *
                                   JxW[point];
        }
     }

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.