]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fixed MatrixCreator::create_mass_matrix and create_laplace_matrix
authorallmaras <allmaras@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 8 Feb 2007 23:50:27 +0000 (23:50 +0000)
committerallmaras <allmaras@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 8 Feb 2007 23:50:27 +0000 (23:50 +0000)
git-svn-id: https://svn.dealii.org/trunk@14428 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/numerics/matrices.cc
deal.II/doc/news/changes.html

index 6959df271a464e85751528056321a559c69f0629..0cfa1ccbd43ac892423e7f50b3af501e9373604b 100644 (file)
@@ -821,9 +821,9 @@ MatrixCreator::create_mass_matrix_2 (const hp::MappingCollection<dim>       &map
                                                 fe_values.shape_value(j,point) *
                                                 weight *
                                                 coefficient_values[point]);
-                         local_rhs(i) += fe_values.shape_value(i,point) *
-                                         rhs_values[point] * weight;
                        }
+                     local_rhs(i) += fe_values.shape_value(i,point) *
+                                     rhs_values[point] * weight;
                    }
                }
            }
@@ -846,9 +846,9 @@ MatrixCreator::create_mass_matrix_2 (const hp::MappingCollection<dim>       &map
                                                 fe_values.shape_value(j,point) *
                                                 weight *
                                                 coefficient_vector_values[point](component_i));
-                         local_rhs(i) += fe_values.shape_value(i,point) *
-                                         rhs_values[point] * weight;
-                       }
+                       }                         
+                     local_rhs(i) += fe_values.shape_value(i,point) *
+                                     rhs_values[point] * weight;
                    }
                }
            }
@@ -868,8 +868,8 @@ MatrixCreator::create_mass_matrix_2 (const hp::MappingCollection<dim>       &map
                      cell_matrix(i,j) += (fe_values.shape_value(i,point) *
                                           fe_values.shape_value(j,point) *
                                           weight);
-                   local_rhs(i) += v * rhs_values[point] * weight;
                  }
+               local_rhs(i) += v * rhs_values[point] * weight;
              }
          }
 
@@ -2115,8 +2115,8 @@ MatrixCreator::create_laplace_matrix_2 (const Mapping<dim>       &mapping,
                                 fe.system_to_component_index(j).first))
                              cell_matrix(i,j) += (Du * Dv * weight *
                                                   coefficient_values[point]);
-                           local_rhs(i) += v * rhs_values[point] * weight;
                          }
+                       local_rhs(i) += v * rhs_values[point] * weight;  
                      }
                  }
            }
@@ -2140,8 +2140,8 @@ MatrixCreator::create_laplace_matrix_2 (const Mapping<dim>       &mapping,
                                (fe.system_to_component_index(j).first == component_i))
                              cell_matrix(i,j) += (Du * Dv * weight *
                                                   coefficient_vector_values[point](component_i));
-                           local_rhs(i) += v * rhs_values[point] * weight;
                          }
+                       local_rhs(i) += v * rhs_values[point] * weight;
                      }
                  }
            }
@@ -2161,8 +2161,8 @@ MatrixCreator::create_laplace_matrix_2 (const Mapping<dim>       &mapping,
                        (fe.system_to_component_index(i).first ==
                         fe.system_to_component_index(j).first))
                      cell_matrix(i,j) += (Du * Dv * weight);
-                   local_rhs(i) += v * rhs_values[point] * weight;
                  }
+               local_rhs(i) += v * rhs_values[point] * weight;
              }
          }
 
@@ -2514,8 +2514,8 @@ MatrixCreator::create_laplace_matrix_2 (const hp::MappingCollection<dim>       &
                                 fe.system_to_component_index(j).first))
                              cell_matrix(i,j) += (Du * Dv * weight *
                                                   coefficient_values[point]);
-                           local_rhs(i) += v * rhs_values[point] * weight;
                          }
+                       local_rhs(i) += v * rhs_values[point] * weight;
                      }
                  }
            }
@@ -2539,8 +2539,8 @@ MatrixCreator::create_laplace_matrix_2 (const hp::MappingCollection<dim>       &
                                (fe.system_to_component_index(j).first == component_i))
                              cell_matrix(i,j) += (Du * Dv * weight *
                                                   coefficient_vector_values[point](component_i));
-                           local_rhs(i) += v * rhs_values[point] * weight;
                          }
+                       local_rhs(i) += v * rhs_values[point] * weight;
                      }
                  }
            }
@@ -2560,8 +2560,8 @@ MatrixCreator::create_laplace_matrix_2 (const hp::MappingCollection<dim>       &
                        (fe.system_to_component_index(i).first ==
                         fe.system_to_component_index(j).first))
                      cell_matrix(i,j) += (Du * Dv * weight);
-                   local_rhs(i) += v * rhs_values[point] * weight;
                  }
+               local_rhs(i) += v * rhs_values[point] * weight;
              }
          }
 
index b9434e91c6206d8773472648cf7fa98be0e93f62..4dd40005e69e22475248d5f0363a001c90b2dd5b 100644 (file)
@@ -975,6 +975,14 @@ inconvenience this causes.
 <h3>deal.II</h3>
 
 <ol>
+   <li> <p>Fixed: <code class="class">MatrixCreator</code>::<code 
+       class="member">create_mass_matrix</code> and <code class="member">
+       create_laplace_matrix</code> computed wrong values for the right 
+       hand sides. This has been fixed. 
+       <br>
+       (Moritz Allmaras 2007/02/08)
+       </p>
+
    <li> <p>Extended: <code>DataOutBase::Patch</code> has been extended by
        a new boolean flag <tt>points_are_available</tt>, which
        defaults to <tt>false</tt>. It is set to <tt>true</tt> if the

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.