Index3ArrayPtr Index4Array::triangulate( ) const {
Index3ArrayPtr _index3List(new Index3Array(__A.size() * 2));
Index3Array::iterator _i3 = _index3List->begin();
for (const_iterator _i4 = __A.begin(); _i4 != __A.end(); _i4++) {
*(_i3++) = Index3(_i4->getAt(0),_i4->getAt(1),_i4->getAt(2));
*(_i3++) = Index3(_i4->getAt(0),_i4->getAt(2),_i4->getAt(3));
};
return _index3List;
};
Index3ArrayPtr IndexArray::triangulate( ) const {
uint_t _size = 0;
for (const_iterator _i4f = __A.begin(); _i4f != __A.end(); _i4f++)
if(_i4f->size() >= 3)_size += _i4f->size() - 2;
Index3ArrayPtr _index3List(new Index3Array(_size));
Index3Array::iterator _i3 = _index3List->begin();
for (const_iterator _i4 = __A.begin();
_i4 != __A.end();
_i4++) {
_size = _i4->size();
if(_size >= 3)
for (Index::const_iterator _i = _i4->begin() + 1;
_i != _i4->end() - 1;
_i++){
(*_i3) = Index3(_i4->getAt(0), (*_i), (*(_i + 1)));_i3++;
}
};
return _index3List;
};
double Diff (int period) { Get_Field (Index3 (period), &dvalue); return (dvalue); }
void Diff (int period, double value) { Put_Field (Index3 (period), value); }
inline double Rectangle::ErrorEstimate(int i,int j) {
return errorWeights[0]*std::fabs(f[Index3(i+1,j,1)]-f[Index3(i-1,j,1)])+errorWeights[1]*std::fabs(f[Index3(i,j+1,1)]-f[Index3(i,j-1,1)])+errorWeights[2]*std::fabs(f[Index3(i+1,j,1)]-2*f[Index3(i,j,1)]+f[Index3(i-1,j,1)])+errorWeights[3]*std::fabs(f[Index3(i,j+1,1)]-2*f[Index3(i,j,1)]+f[Index3(i,j-1,1)])+errorWeights[4]*std::fabs(f[Index3(i,j,1)]);
}
