Changeset 830 for trunk/yat/utility/vector.cc
- Timestamp:
- Mar 24, 2007, 3:54:42 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/yat/utility/vector.cc
r825 r830 55 55 throw utility::GSL_error("vector::vector failed to allocate memory"); 56 56 57 set_all(init_value);57 all(init_value); 58 58 } 59 59 … … 300 300 if (!v_) 301 301 throw utility::GSL_error("vector::vector failed to allocate memory"); 302 set_all(init_value);302 all(init_value); 303 303 } 304 304 … … 311 311 312 312 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) 322 314 { 323 315 assert(v_); … … 397 389 const vector& vector::operator=( const vector& other ) 398 390 { 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"); 400 394 return *this; 401 395 } … … 478 472 bool nan(const vector& templat, vector& flag) 479 473 { 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)); 483 477 else 484 flag. set_all(1.0);478 flag.all(1.0); 485 479 bool nan=false; 486 for (size_t i=0; i< rows; i++)480 for (size_t i=0; i<vsize; i++) 487 481 if (std::isnan(templat(i))) { 488 482 flag(i)=0;
Note: See TracChangeset
for help on using the changeset viewer.