Fix one index in distribute_local_to_global_function.

git-svn-id: https://svn.dealii.org/trunk@18426 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/include/dofs/dof_constraints.templates.h b/deal.II/deal.II/include/dofs/dof_constraints.templates.h
index a1ffb952c56801fad70c647b08e6b5da5b8aa9f2..d237fe415cb273801ed0a2707f323b8f26645770 100644 (file)
--- a/deal.II/deal.II/include/dofs/dof_constraints.templates.h
+++ b/deal.II/deal.II/include/dofs/dof_constraints.templates.h
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -986,7 +986,7 @@ distribute_local_to_global (const FullMatrix<double>        &local_matrix,
                                   // row of the inhomogeneous constraint in
                                   // the matrix with Gauss elimination
                  if (use_vectors == true)
-                   global_vector(local_dof_indices[i]) -= local_matrix(j,i) * 
+                   global_vector(local_dof_indices[i]) -= local_matrix(i,j) * 
                                                           position_j->inhomogeneity;
                 }
               else if ((is_constrained_i == true) &&
@@ -1011,7 +1011,7 @@ distribute_local_to_global (const FullMatrix<double>        &local_matrix,
 
                      if (use_vectors == true)
                        global_vector (position_i->entries[p].first) -=
-                         local_matrix(j,i) * position_i->entries[p].second *
+                         local_matrix(i,j) * position_i->entries[p].second *
                          position_j->inhomogeneity;
                    }