]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix accidental checkin of merge conflict marker :-(
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 30 Apr 2013 16:21:59 +0000 (16:21 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 30 Apr 2013 16:21:59 +0000 (16:21 +0000)
git-svn-id: https://svn.dealii.org/branches/branch_bigger_global_dof_indices_4@29411 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/include/deal.II/lac/full_matrix.templates.h

index a0bbe2ef1fbbc5ddca2cb9639bd89de87924df28..39a2d461fca3a6a4e7365f507715bb2188d9dd78 100644 (file)
@@ -744,27 +744,13 @@ void FullMatrix<number>::mTmult (FullMatrix<number2>       &dst,
 
   const size_type m = this->m(), n = src.m(), l = this->n();
 
-<<<<<<< .working
-  // arrange the loops in a way that
-  // we keep write operations low,
-  // (writing is usually more costly
-  // than reading).
-  for (size_type i=0; i<m; i++)
-    for (size_type j=0; j<n; j++)
-      {
-        number2 add_value = adding ? dst(i,j) : 0.;
-        for (size_type k=0; k<l; k++)
-          add_value += (number2)(*this)(i,k) * (number2)(src(j,k));
-        dst(i,j) = add_value;
-      }
-=======
   // symmetric matrix if the two matrices are the same
   if (PointerComparison::equal(this, &src))
-    for (unsigned int i=0; i<m; ++i)
-      for (unsigned int j=i; j<m; ++j)
+    for (size_type i=0; i<m; ++i)
+      for (size_type j=i; j<m; ++j)
         {
           number2 add_value = 0.;
-          for (unsigned int k=0; k<l; ++k)
+          for (size_type k=0; k<l; ++k)
             add_value += (number2)(*this)(i,k) * (number2)(*this)(j,k);
           if (adding)
             {
@@ -778,15 +764,14 @@ void FullMatrix<number>::mTmult (FullMatrix<number2>       &dst,
   else
     // arrange the loops in a way that we keep write operations low, (writing is
     // usually more costly than reading).
-    for (unsigned int i=0; i<m; i++)
-      for (unsigned int j=0; j<n; j++)
+    for (size_type i=0; i<m; i++)
+      for (size_type j=0; j<n; j++)
         {
           number2 add_value = adding ? dst(i,j) : 0.;
-          for (unsigned int k=0; k<l; k++)
+          for (size_type k=0; k<l; k++)
             add_value += (number2)(*this)(i,k) * (number2)(src(j,k));
           dst(i,j) = add_value;
         }
->>>>>>> .merge-right.r29408
 }
 
 

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.