]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Move the dVector class closer to the STL's vector<>.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 6 May 1998 16:03:56 +0000 (16:03 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 6 May 1998 16:03:56 +0000 (16:03 +0000)
git-svn-id: https://svn.dealii.org/trunk@263 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/dofs/dof_accessor.cc
deal.II/deal.II/source/dofs/dof_constraints.cc
deal.II/deal.II/source/dofs/dof_handler.cc
deal.II/deal.II/source/grid/tria.cc
deal.II/deal.II/source/numerics/assembler.cc
deal.II/deal.II/source/numerics/base.cc
deal.II/lac/include/lac/dvector.h
deal.II/lac/source/dfmatrix.cc
deal.II/lac/source/dsmatrix.cc
deal.II/lac/source/dvector.cc

index e2e620dc0a60d15966bda44842b4c76a8f1ecf20..18b5c5416f31161e3251db014b5e2b21d7ae71fc 100644 (file)
@@ -316,7 +316,7 @@ DoFCellAccessor<1>::get_dof_values (const dVector  &values,
   Assert (&dof_handler->get_selected_fe() != 0, ExcInvalidObject());
   Assert (dof_values.size() == dof_handler->get_selected_fe().total_dofs,
          ExcVectorDoesNotMatch());
-  Assert (values.n() == dof_handler->n_dofs(),
+  Assert (values.size() == dof_handler->n_dofs(),
          ExcVectorDoesNotMatch());
 
   vector<double>::iterator next_dof_value=dof_values.begin();
@@ -337,7 +337,7 @@ DoFCellAccessor<2>::get_dof_values (const dVector  &values,
   Assert (&dof_handler->get_selected_fe() != 0, ExcInvalidObject());
   Assert (dof_values.size() == dof_handler->get_selected_fe().total_dofs,
          ExcVectorDoesNotMatch());
-  Assert (values.n() == dof_handler->n_dofs(),
+  Assert (values.size() == dof_handler->n_dofs(),
          ExcVectorDoesNotMatch());
 
   vector<double>::iterator next_dof_value=dof_values.begin();
index c418c3cfcec1f0a9646f0bf5445812b4b4c31422..eeab1152af900e83b8e56e5faf484ed2bdbcbab2 100644 (file)
@@ -505,7 +505,7 @@ void ConstraintMatrix::condense (dSMatrix &uncondensed) const {
 void ConstraintMatrix::condense (const dVector &uncondensed,
                                 dVector       &condensed) const {
   Assert (sorted == true, ExcMatrixNotClosed());
-  Assert (condensed.n()+n_constraints() == uncondensed.n(),
+  Assert (condensed.size()+n_constraints() == uncondensed.size(),
          ExcWrongDimension());
   
                                   // store for each line of the vector
@@ -514,11 +514,11 @@ void ConstraintMatrix::condense (const dVector &uncondensed,
                                   // -1, this line will be condensed away
   vector<int> new_line;
 
-  new_line.reserve (uncondensed.n());
+  new_line.reserve (uncondensed.size());
 
   vector<ConstraintLine>::const_iterator next_constraint = lines.begin();
   unsigned int                           shift           = 0;
-  unsigned int n_rows = uncondensed.n();
+  unsigned int n_rows = uncondensed.size();
 
   if (next_constraint == lines.end()) 
                                     // if no constraint is to be handled
@@ -553,7 +553,7 @@ void ConstraintMatrix::condense (const dVector &uncondensed,
                                   // only evaluated so often as there are
                                   // entries in new_line[*] which tells us
                                   // which constraints exist
-  for (unsigned int row=0; row<uncondensed.n(); ++row)
+  for (unsigned int row=0; row<uncondensed.size(); ++row)
     if (new_line[row] != -1)
                                       // line not constrained
                                       // copy entry
@@ -581,7 +581,7 @@ void ConstraintMatrix::condense (dVector &vec) const {
     return;
   
   vector<ConstraintLine>::const_iterator next_constraint = lines.begin();
-  for (unsigned int row=0; row<vec.n(); ++row)
+  for (unsigned int row=0; row<vec.size(); ++row)
     if (row == next_constraint->line)
                                       // line must be distributed
       {
@@ -605,7 +605,7 @@ void ConstraintMatrix::condense (dVector &vec) const {
 void ConstraintMatrix::distribute (const dVector &condensed,
                                   dVector       &uncondensed) const {
   Assert (sorted == true, ExcMatrixNotClosed());
-  Assert (condensed.n()+n_constraints() == uncondensed.n(),
+  Assert (condensed.size()+n_constraints() == uncondensed.size(),
          ExcWrongDimension());
 
                                   // store for each line of the new vector
@@ -614,11 +614,11 @@ void ConstraintMatrix::distribute (const dVector &condensed,
                                   // -1, this line was condensed away
   vector<int> old_line;
 
-  old_line.reserve (uncondensed.n());
+  old_line.reserve (uncondensed.size());
 
   vector<ConstraintLine>::const_iterator next_constraint = lines.begin();
   unsigned int                           shift           = 0;
-  unsigned int n_rows = uncondensed.n();
+  unsigned int n_rows = uncondensed.size();
 
   if (next_constraint == lines.end()) 
                                     // if no constraint is to be handled
@@ -653,7 +653,7 @@ void ConstraintMatrix::distribute (const dVector &condensed,
                                   // only evaluated so often as there are
                                   // entries in new_line[*] which tells us
                                   // which constraints exist
-  for (unsigned int line=0; line<uncondensed.n(); ++line) 
+  for (unsigned int line=0; line<uncondensed.size(); ++line) 
     if (old_line[line] != -1)
                                       // line was not condensed away
       uncondensed(line) = condensed(old_line[line]);
index 1f325e92709d2dacef0e0c71a79010cce72448f6..c5cdcbc89cd8dd312fb989ac730435da504e9972 100644 (file)
@@ -1470,8 +1470,8 @@ unsigned int DoFHandler<2>::max_transfer_entries (const unsigned int max_level_d
 template <int dim>
 void DoFHandler<dim>::distribute_cell_to_dof_vector (const dVector &cell_data,
                                                     dVector       &dof_data) const {
-  Assert (cell_data.n()==tria->n_active_cells(),
-         ExcWrongSize (cell_data.n(), tria->n_active_cells()));
+  Assert (cell_data.size()==tria->n_active_cells(),
+         ExcWrongSize (cell_data.size(), tria->n_active_cells()));
 
                                   // assign the right size to the output vector
   dof_data.reinit (n_dofs());
index 0a2daf6ca7154bef72eeb3078567d0ee57f7a749..70019f75d240bd8098a3200f72f8ddf1b608ea68 100644 (file)
@@ -1731,11 +1731,11 @@ void Triangulation<2>::print_gnuplot (ostream &out,
 template <int dim>
 void Triangulation<dim>::refine (const dVector &criteria,
                                 const double   threshold) {
-  Assert (criteria.n() == n_active_cells(),
-         ExcInvalidVectorSize(criteria.n(), n_active_cells()));
+  Assert (criteria.size() == n_active_cells(),
+         ExcInvalidVectorSize(criteria.size(), n_active_cells()));
 
   active_cell_iterator cell = begin_active();
-  const unsigned int n_cells = criteria.n();
+  const unsigned int n_cells = criteria.size();
   
   for (unsigned int index=0; index<n_cells; ++cell, ++index)
     if (criteria(index) > threshold)
@@ -1749,12 +1749,12 @@ void Triangulation<dim>::refine_fixed_number (const dVector &criteria,
                                              const double   fraction) {
   dVector tmp(criteria);
   nth_element (tmp.begin(),
-              tmp.begin()+static_cast<int>(fraction*tmp.n()),
+              tmp.begin()+static_cast<int>(fraction*tmp.size()),
               tmp.end(),
               greater<double>());
 
   refine (criteria, *(tmp.begin() +
-                     static_cast<int>(fraction*tmp.n())));
+                     static_cast<int>(fraction*tmp.size())));
 };
 
 
index 475545b464b0bfaca294f1fe6dad3c36396867be..d016ae11cbcede3602da4c06b161b0562212bdb0 100644 (file)
@@ -90,7 +90,7 @@ Assembler<dim>::Assembler (Triangulation<dim> *tria,
   Assert (matrix.n() == dof_handler->n_dofs(), ExcInvalidData());
   Assert (((AssemblerData<dim>*)local_data)->fe == dof_handler->get_selected_fe(),
          ExcInvalidData());
-  Assert (rhs_vector.n() == dof_handler->n_dofs(),
+  Assert (rhs_vector.size() == dof_handler->n_dofs(),
          ExcInvalidData());
 };
 
index b411e92bcd9e01686b1578843c81f0b16920d977..13b5f823f31b11da4a9b44a1a27169525f348b13 100644 (file)
@@ -145,7 +145,7 @@ void ProblemBase<dim>::solve () {
   double tolerance = 1.e-16;
   
   Control                          control1(max_iter,tolerance);
-  PrimitiveVectorMemory<dVector>   memory(right_hand_side.n());
+  PrimitiveVectorMemory<dVector>   memory(right_hand_side.size());
   CG<dSMatrix,dVector>             cg(control1,memory);
 
                                   // solve
index 6eb44dc1d677303f38ad76ed7068b5535f440cc3..692986fc0672d702c80cc5b80fa37bdfaa4598a4 100644 (file)
@@ -151,7 +151,7 @@ class dVector : public VectorBase
                                      *  Inquire Dimension. returns Dimension , 
                                      *             INLINE
                                      */
-    unsigned int n () const;
+    unsigned int size () const;
 
                                     /**
                                      * Make the #dVector# class a bit like the
@@ -388,7 +388,7 @@ class dVector : public VectorBase
 
 
 
-inline unsigned int dVector::n () const
+inline unsigned int dVector::size () const
 {
   return dim;
 }
index dd42cbabeda527293874e8d1ef76c256fe2289d2..b2b3340f76ce88d47c1092f579343cd7dad55acf 100644 (file)
@@ -49,8 +49,8 @@ void dFMatrix::reinit (const unsigned int mm, const unsigned int nn)
 void dFMatrix::vmult (dVector& dst, const dVector& src,
                      const bool adding) const
 {
-  Assert(dst.n() == m(), ExcDimensionMismatch(dst.n(), m()));
-  Assert(src.n() == n(), ExcDimensionMismatch(src.n(), n()));
+  Assert(dst.size() == m(), ExcDimensionMismatch(dst.size(), m()));
+  Assert(src.size() == n(), ExcDimensionMismatch(src.size(), n()));
 
   double s;
   if ((n()==3) && (m()==3))
@@ -173,8 +173,8 @@ void dFMatrix::vmult (dVector& dst, const dVector& src,
 void dFMatrix::gsmult (dVector& dst, const dVector& src, const iVector& gl) const
 {
   Assert(n() == m(), ExcNotQuadratic());
-  Assert(dst.n() == n(), ExcDimensionMismatch(dst.n(), n()));
-  Assert(src.n() == n(), ExcDimensionMismatch(src.n(), n()));
+  Assert(dst.size() == n(), ExcDimensionMismatch(dst.size(), n()));
+  Assert(src.size() == n(), ExcDimensionMismatch(src.size(), n()));
   Assert(gl.n() == n(), ExcDimensionMismatch(gl.n(), n()));
 
   double s;
@@ -303,8 +303,8 @@ void dFMatrix::gsmult (dVector& dst, const dVector& src, const iVector& gl) cons
 
 void dFMatrix::Tvmult (dVector& dst, const dVector& src, const bool adding) const
 {
-  Assert(dst.n() == n(), ExcDimensionMismatch(dst.n(), n()));
-  Assert(src.n() == m(), ExcDimensionMismatch(src.n(), m()));
+  Assert(dst.size() == n(), ExcDimensionMismatch(dst.size(), n()));
+  Assert(src.size() == m(), ExcDimensionMismatch(src.size(), m()));
 
   unsigned int i,j;
   double s;
@@ -323,9 +323,9 @@ void dFMatrix::Tvmult (dVector& dst, const dVector& src, const bool adding) cons
 double dFMatrix::residual (dVector& dst, const dVector& src,
                          const dVector& right) const
 {
-  Assert(dst.n() == m(), ExcDimensionMismatch(dst.n(), m()));
-  Assert(src.n() == n(), ExcDimensionMismatch(src.n(), n()));
-  Assert(right.n() == m(), ExcDimensionMismatch(right.n(), m()));
+  Assert(dst.size() == m(), ExcDimensionMismatch(dst.size(), m()));
+  Assert(src.size() == n(), ExcDimensionMismatch(src.size(), n()));
+  Assert(right.size() == m(), ExcDimensionMismatch(right.size(), m()));
 
   unsigned int i,j;
   double s, res = 0.;
@@ -345,8 +345,8 @@ double dFMatrix::residual (dVector& dst, const dVector& src,
 void dFMatrix::forward (dVector& dst, const dVector& src) const
 {
   Assert(n() == m(), ExcNotQuadratic());
-  Assert(dst.n() == n(), ExcDimensionMismatch(dst.n(), n()));
-  Assert(src.n() == n(), ExcDimensionMismatch(src.n(), n()));
+  Assert(dst.size() == n(), ExcDimensionMismatch(dst.size(), n()));
+  Assert(src.size() == n(), ExcDimensionMismatch(src.size(), n()));
 
   unsigned int i,j;
   unsigned int nu = MIN(m(),n());
@@ -362,8 +362,8 @@ void dFMatrix::forward (dVector& dst, const dVector& src) const
 void dFMatrix::backward (dVector& dst, const dVector& src) const
 {
   Assert(n() == m(), ExcNotQuadratic());
-  Assert(dst.n() == n(), ExcDimensionMismatch(dst.n(), n()));
-  Assert(src.n() == n(), ExcDimensionMismatch(src.n(), n()));
+  Assert(dst.size() == n(), ExcDimensionMismatch(dst.size(), n()));
+  Assert(src.size() == n(), ExcDimensionMismatch(src.size(), n()));
 
   unsigned int j;
   unsigned int nu = MIN(m(),n());
index 6ca5d49372e4cef8e8ad5f0a1e15048aad8a4755..a0ac457eb9cad70d5e60e54c8890cc1a05da64d9 100644 (file)
@@ -386,8 +386,8 @@ void
 dSMatrix::vmult (dVector& dst, const dVector& src) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
-  Assert(m() == dst.n(), ExcDimensionsDontMatch(m(),dst.n()));
-  Assert(n() == src.n(), ExcDimensionsDontMatch(n(),src.n()));
+  Assert(m() == dst.size(), ExcDimensionsDontMatch(m(),dst.size()));
+  Assert(n() == src.size(), ExcDimensionsDontMatch(n(),src.size()));
 
   for (unsigned int i=0; i<m(); ++i)
     {
@@ -403,8 +403,8 @@ void
 dSMatrix::Tvmult (dVector& dst, const dVector& src) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
-  Assert(n() == dst.n(), ExcDimensionsDontMatch(n(),dst.n()));
-  Assert(m() == src.n(), ExcDimensionsDontMatch(m(),src.n()));
+  Assert(n() == dst.size(), ExcDimensionsDontMatch(n(),dst.size()));
+  Assert(m() == src.size(), ExcDimensionsDontMatch(m(),src.size()));
 
   unsigned int i;
   
@@ -426,8 +426,8 @@ double
 dSMatrix::matrix_norm (const dVector& v) const
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
-  Assert(m() == v.n(), ExcDimensionsDontMatch(m(),v.n()));
-  Assert(n() == v.n(), ExcDimensionsDontMatch(n(),v.n()));
+  Assert(m() == v.size(), ExcDimensionsDontMatch(m(),v.size()));
+  Assert(n() == v.size(), ExcDimensionsDontMatch(n(),v.size()));
 
   double sum = 0.;
   for (unsigned int i=0;i<m();i++)
@@ -447,9 +447,9 @@ double
 dSMatrix::residual (dVector& dst, const dVector& u, const dVector& b)
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
-  Assert(m() == dst.n(), ExcDimensionsDontMatch(m(),dst.n()));
-  Assert(m() == b.n(), ExcDimensionsDontMatch(m(),b.n()));
-  Assert(n() == u.n(), ExcDimensionsDontMatch(n(),u.n()));
+  Assert(m() == dst.size(), ExcDimensionsDontMatch(m(),dst.size()));
+  Assert(m() == b.size(), ExcDimensionsDontMatch(m(),b.size()));
+  Assert(n() == u.size(), ExcDimensionsDontMatch(n(),u.size()));
 
   double s,norm=0.;   
   
@@ -472,7 +472,7 @@ dSMatrix::Jacobi_precond (dVector& dst, const dVector& src, const double om)
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
 
-  const unsigned int n = src.n();
+  const unsigned int n = src.size();
 
   for (unsigned int i=0;i<n;++i)
     {
@@ -486,7 +486,7 @@ dSMatrix::SSOR_precond (dVector& dst, const dVector& src, const double om)
   Assert (cols != 0, ExcMatrixNotInitialized());
 
   int p;
-  unsigned int  n = src.n();
+  const unsigned int  n = src.size();
   unsigned int  j;
   
   for (unsigned i=0; i<n; i++)
@@ -527,7 +527,7 @@ dSMatrix::SOR (dVector& dst, const double om)
 {
   Assert (cols != 0, ExcMatrixNotInitialized());
   Assert(n() == m(), ExcDimensionsDontMatch(n(),m()));
-  Assert(m() == dst.n(), ExcDimensionsDontMatch(m(),dst.n()));
+  Assert(m() == dst.size(), ExcDimensionsDontMatch(m(),dst.size()));
 
   for (unsigned int i=0;i<m();i++)
     {
@@ -548,7 +548,7 @@ dSMatrix::SSOR (dVector& dst, const double om)
   Assert (cols != 0, ExcMatrixNotInitialized());
 
   int p;
-  unsigned int  n = dst.n();
+  const unsigned int  n = dst.size();
   unsigned int  j;
   double s;
   
index c65fb20e2b67a43768c1f4b68917dd0b987cee15..b04e9b74bc20eaa8e9c5acd71f6a17ae6bbd5bda 100644 (file)
@@ -33,8 +33,8 @@ dVector::dVector (const unsigned int n) :
 
 dVector::dVector (const dVector& v) :
                VectorBase(v),
-               dim(v.n()),
-               maxdim(v.n()),
+               dim(v.size()),
+               maxdim(v.size()),
                val(0)
 {
   if (dim)
@@ -67,7 +67,7 @@ void dVector::reinit (const unsigned int n, const bool fast)
 
 void dVector::reinit (const dVector& v, const bool fast)
 {
-  const unsigned int n = v.n();
+  const unsigned int n = v.size();
   if (n>maxdim)
   {
     if (val) delete[] val;
@@ -170,7 +170,7 @@ double dVector::mean_value () const
   for (unsigned int i=(dim/4)*4; i<dim; ++i)
     sum0 += val[i];
   
-  return (sum0+sum1+sum2+sum3)/n();
+  return (sum0+sum1+sum2+sum3)/size();
 };
 
 
@@ -510,7 +510,7 @@ const char* dVector::name () const
 void dVector::print (FILE* f, const char* format) const
 {
   if (!format) format = " %5.2f";
-  for (unsigned int j=0;j<n();j++)
+  for (unsigned int j=0;j<size();j++)
     fprintf(f, format, val[j]);
   fputc('\n',f);
 }
@@ -520,7 +520,7 @@ void dVector::print (FILE* f, const char* format) const
 void dVector::print (const char* format) const
 {
   if (!format) format = " %5.2f";
-  for (unsigned int j=0;j<n();j++)
+  for (unsigned int j=0;j<size();j++)
     printf (format, val[j]);
   printf ("\n");
 }
@@ -528,6 +528,6 @@ void dVector::print (const char* format) const
 
 
 void dVector::print (ostream &out) const {
-  for (unsigned int i=0; i<n(); ++i)
+  for (unsigned int i=0; i<size(); ++i)
     out << val[i] << endl;
 };

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.