Ignore:
Timestamp:
Mar 24, 2007, 3:54:42 PM (16 years ago)
Author:
Jari Häkkinen
Message:

Fixes #206.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/yat/utility/vector.cc

    r825 r830  
    5555      throw utility::GSL_error("vector::vector failed to allocate memory");
    5656
    57     set_all(init_value);
     57    all(init_value);
    5858  }
    5959
     
    300300    if (!v_)
    301301      throw utility::GSL_error("vector::vector failed to allocate memory");
    302     set_all(init_value);
     302    all(init_value);
    303303  }
    304304
     
    311311
    312312
    313   void vector::set(const vector& vec)
    314   {
    315     assert(v_);
    316     if (gsl_vector_memcpy(v_,vec.gsl_vector_p()))
    317       throw utility::GSL_error("vector::set dimension mis-match");
    318   }
    319 
    320 
    321   void vector::set_all(const double& value)
     313  void vector::all(const double& value)
    322314  {
    323315    assert(v_);
     
    397389  const vector& vector::operator=( const vector& other )
    398390  {
    399     set(other);
     391    assert(v_);
     392    if (gsl_vector_memcpy(v_,other.gsl_vector_p()))
     393      throw utility::GSL_error("vector::set dimension mis-match");
    400394    return *this;
    401395  }
     
    478472  bool nan(const vector& templat, vector& flag)
    479473  {
    480     size_t rows=templat.size();
    481     if (rows!=flag.size())
    482       flag.clone(vector(rows,1.0));
     474    size_t vsize=templat.size();
     475    if (vsize!=flag.size())
     476      flag.clone(vector(vsize,1.0));
    483477    else
    484       flag.set_all(1.0);
     478      flag.all(1.0);
    485479    bool nan=false;
    486     for (size_t i=0; i<rows; i++)
     480    for (size_t i=0; i<vsize; i++)
    487481      if (std::isnan(templat(i))) {
    488482        flag(i)=0;
Note: See TracChangeset for help on using the changeset viewer.