]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Minor adjustments.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 29 May 2000 11:17:09 +0000 (11:17 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 29 May 2000 11:17:09 +0000 (11:17 +0000)
git-svn-id: https://svn.dealii.org/trunk@2966 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_tools.cc

index f698b4f5e4e7bb40c5b15ca664e42aee79edc60b..63ad0ced20350f24e94585cbe798a398743818eb 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 1998, 1999, 2000 by the deal.II authors
+//    Copyright (C) 2000 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -67,30 +67,10 @@ void FETools::get_interpolation_matrix(const FiniteElement<dim> &fe1,
            interpolation_matrix(i,j)=0.;
        }  
     }
-  
-/*
-  vector<Point<dim> > unit_support_points (fe2.dofs_per_cell);
-  fe2.get_unit_support_points (unit_support_points);
-  
-  for (unsigned int i=0; i<fe2.dofs_per_cell; ++i)     
-    {
-      const unsigned int i1
-       = fe2.system_to_component_index(i).first;
-      for (unsigned int j=0; j<fe1.dofs_per_cell; ++j)
-       {
-         const unsigned int j1
-           = fe1.system_to_component_index(j).first;
-         if (i1==j1)
-           interpolation_matrix(i,j) =
-             fe1.shape_value (j, unit_support_points[i]);
-         else
-           interpolation_matrix(i,j)=0.;
-       }
-    }
-*/
 }
 
 
+
 template <int dim, typename number>
 void FETools::get_back_interpolation_matrix(const FiniteElement<dim> &fe1,
                                            const FiniteElement<dim> &fe2,
@@ -173,7 +153,7 @@ void FETools::interpolate(const DoFHandler<dim> &dof1,
 
   vector<unsigned int> index_multiplicity(dof2.n_dofs(),0);
   vector<unsigned int> dofs (dofs_per_cell2);
-  u2=0;
+  u2.clear ();
   
   for (; cell1!=endc1, cell2!=endc2; ++cell1, ++cell2) 
     {
@@ -189,7 +169,7 @@ void FETools::interpolate(const DoFHandler<dim> &dof1,
   for (unsigned int i=0; i<dof2.n_dofs(); ++i)
     {
       Assert(index_multiplicity[i]!=0, ExcInternalError());
-      u2(i)/=static_cast<double>(index_multiplicity[i]);
+      u2(i) /= index_multiplicity[i];
     }
 }
 
@@ -262,6 +242,7 @@ void FETools::interpolation_difference(const DoFHandler<dim> &dof1,
 }
 
 
+
 template <int dim, typename number>
 void FETools::extrapolate(const DoFHandler<dim> &dof1,
                          const Vector<number> &u1,

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.