parameter_vector.C

Go to the documentation of this file.

// The libMesh Finite Element Library.
// Copyright (C) 2002-2014 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner

// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.

// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.

// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA



// C++ Includes   -----------------------------------

// Local Includes -----------------------------------
#include "libmesh/parameter_vector.h"

namespace libMesh
{

// ------------------------------------------------------------
// ParameterVector implementation



void ParameterVector::deep_copy(ParameterVector &target) const
{
  const unsigned int Np = cast_int<unsigned int>
    (this->_params.size());
  target.clear();
  target._params.resize(Np);
  target._my_data.resize(Np);
  for (unsigned int i=0; i != Np; ++i)
    {
      target._params[i] =
        new ParameterPointer<Number>(&target._my_data[i]);
      target._my_data[i] = *(*this)[i];
    }
}



void ParameterVector::shallow_copy(ParameterVector &target) const
{
  target._my_data.clear();
  target._params = this->_params;
  target._is_shallow_copy = true;
}



void ParameterVector::value_copy(ParameterVector &target) const
{
  const unsigned int Np = cast_int<unsigned int>
    (this->_params.size());
  libmesh_assert_equal_to (target._params.size(), Np);

  for (unsigned int i=0; i != Np; ++i)
    *target[i] = *(*this)[i];
}



void ParameterVector::resize(unsigned int s)
{
  libmesh_assert(!_is_shallow_copy);

  const std::size_t old_size = this->_params.size();

  // If we're shrinking the vector, we don't want to leak memory.
  // Note that we're using < in these for loops, not !=
  // We don't know a priori if we're shrinking or growing
  for (unsigned int i=s; i < old_size; ++i)
    delete _params[i];

  this->_params.resize(s);

  for (unsigned int i=old_size; i < s; ++i)
    this->_params[i] =
      new ParameterPointer<Number>(NULL);
}



void ParameterVector::deep_resize(unsigned int s)
{
  libmesh_assert(!_is_shallow_copy);

  this->_params.resize(s);
  this->_my_data.resize(s);
  for (unsigned int i=0; i != s; ++i)
    this->_params[i] =
      new ParameterPointer<Number>(&this->_my_data[i]);
}



ParameterVector& ParameterVector::operator *= (const Number a)
{
  const unsigned int Np = cast_int<unsigned int>
    (this->_params.size());
  for (unsigned int i=0; i != Np; ++i)
    *(*this)[i] *= a;
  return *this;
}



ParameterVector& ParameterVector::operator += (const ParameterVector& a)
{
  const unsigned int Np = cast_int<unsigned int>
    (this->_params.size());
  libmesh_assert_equal_to (a._params.size(), Np);
  for (unsigned int i=0; i != Np; ++i)
    *(*this)[i] += *a[i];
  return *this;
}


} // namespace libMesh