| 146 |
|
//error |
| 147 |
|
cerr << "size does not match"<< endl; |
| 148 |
|
} |
| 149 |
+ |
|
| 150 |
+ |
if (storageLayout_ & dslSitePotential && sitePotential.size() != size_) { |
| 151 |
+ |
//error |
| 152 |
+ |
cerr << "size does not match"<< endl; |
| 153 |
+ |
} |
| 154 |
|
|
| 155 |
|
return size_; |
| 156 |
|
|
| 226 |
|
internalResize(flucQFrc, newSize); |
| 227 |
|
} |
| 228 |
|
|
| 229 |
+ |
if (storageLayout_ & dslSitePotential) { |
| 230 |
+ |
internalResize(sitePotential, newSize); |
| 231 |
+ |
} |
| 232 |
+ |
|
| 233 |
|
size_ = newSize; |
| 234 |
|
} |
| 235 |
|
|
| 301 |
|
if (storageLayout_ & dslFlucQForce) { |
| 302 |
|
flucQFrc.reserve(size); |
| 303 |
|
} |
| 304 |
+ |
|
| 305 |
+ |
if (storageLayout_ & dslSitePotential) { |
| 306 |
+ |
sitePotential.reserve(size); |
| 307 |
+ |
} |
| 308 |
|
} |
| 309 |
|
|
| 310 |
|
void DataStorage::copy(int source, int num, int target) { |
| 378 |
|
if (storageLayout_ & dslFlucQForce) { |
| 379 |
|
internalCopy(flucQFrc, source, num, target); |
| 380 |
|
} |
| 381 |
+ |
|
| 382 |
+ |
if (storageLayout_ & dslSitePotential) { |
| 383 |
+ |
internalCopy(sitePotential, source, num, target); |
| 384 |
+ |
} |
| 385 |
|
} |
| 386 |
|
|
| 387 |
|
int DataStorage::getStorageLayout() { |
| 447 |
|
case dslFlucQForce: |
| 448 |
|
return internalGetArrayPointer(flucQFrc); |
| 449 |
|
|
| 450 |
+ |
case dslSitePotential: |
| 451 |
+ |
return internalGetArrayPointer(sitePotential); |
| 452 |
+ |
|
| 453 |
|
default: |
| 454 |
|
//error message |
| 455 |
|
return NULL; |
| 569 |
|
if (layout & dslFlucQForce) { |
| 570 |
|
bytes += sizeof(RealType); |
| 571 |
|
} |
| 572 |
+ |
if (layout & dslSitePotential) { |
| 573 |
+ |
bytes += sizeof(RealType); |
| 574 |
+ |
} |
| 575 |
|
|
| 576 |
|
return bytes; |
| 577 |
|
} |