uint64_t IndexSetRegionSequence::computeNumberOfIndices(LocalArray &lower,
        LocalArray &upper) {
    unsigned numberOfDimensions = lower.getElementCount(1);
    _numberOfRegions = lower.getElementCount(0);
    if (numberOfDimensions == 0 || _numberOfRegions == 0)
        throw EInvalidArgument();
    if (upper.getElementCount(1) != numberOfDimensions
            || upper.getElementCount(0) != _numberOfRegions)
        throw EInvalidArgument();
    _lowerIndexList = new uint64_t[numberOfDimensions * _numberOfRegions];
    _upperIndexList = new uint64_t[numberOfDimensions * _numberOfRegions];
    if (_lowerIndexList == NULL || _upperIndexList == NULL)
        throw ENotEnoughMemory();
    if (lower.getElementType() == ArrayShape::UINT64)
        memcpy(_lowerIndexList, lower.getData(),
                sizeof (uint64_t) * numberOfDimensions * _numberOfRegions);
    else {
        LocalArray properTyped(ArrayShape::UINT64, lower);
        memcpy(_lowerIndexList, properTyped.getData(),
                sizeof (uint64_t) * numberOfDimensions * _numberOfRegions);
    }
    if (upper.getElementType() == ArrayShape::UINT64)
        memcpy(_upperIndexList, upper.getData(),
                sizeof (uint64_t) * numberOfDimensions * _numberOfRegions);
    else {
        LocalArray properTyped(ArrayShape::UINT64, upper);
        memcpy(_upperIndexList, properTyped.getData(),
                sizeof (uint64_t) * numberOfDimensions * _numberOfRegions);
    }

    uint64_t result = 0;
    uint64_t k = 0;
    for (uint64_t iRegion = 0; iRegion < _numberOfRegions; iRegion++) {
        uint64_t currentRegionSize = 1;
        for (unsigned iDim = 0; iDim < numberOfDimensions; iDim++) {
            uint64_t currentLower = _lowerIndexList[k];
            uint64_t currentUpper = _upperIndexList[k];
            k++;
            if (currentUpper < currentLower)
                throw EInvalidRegionDescriptor(iDim, iRegion, 
                        currentLower, currentUpper, 1);
            currentRegionSize *= currentUpper - currentLower + 1;
        }
        result += currentRegionSize;
    }
    return result;
}
Exemplo n.º 2
0
IndexSetPointSequence::IndexSetPointSequence(LocalArray &points) :
IndexSet(points.getElementCount(1), points.getElementCount(0)) {
    _numberOfRegions = 0;
    uint64_t indexLength = _numberOfIndices * _numberOfDimensions;
    _indexList = new uint64_t[indexLength];
    if (_indexList == NULL)
        throw ENotEnoughMemory();
    memcpy(_indexList, points.getData(), indexLength * sizeof (uint64_t));
    this->initConstantIterators();
    if (_begin == NULL || _end == NULL || _curr == NULL)
        throw ENotEnoughMemory();
}