CoordUtils.cpp

Go to the documentation of this file.

#include "CoordUtils.hpp"
#include "MayaFlux/Kakshya/Region/Region.hpp"
 
namespace MayaFlux::Kakshya {
 
uint64_t coordinates_to_linear(const std::vector<uint64_t>& coords, const std::vector<DataDimension>& dimensions)
{
    uint64_t index = 0;
    uint64_t multiplier = 1;
    for (int i = dimensions.size() - 1; i >= 0; --i) {
        if (i < coords.size()) {
            index += coords[i] * multiplier;
        }
        multiplier *= dimensions[i].size;
    }
    return index;
}
 
std::vector<uint64_t> linear_to_coordinates(uint64_t index, const std::vector<DataDimension>& dimensions)
{
    std::vector<uint64_t> coords(dimensions.size());
    for (int i = dimensions.size() - 1; i >= 0; --i) {
        coords[i] = index % dimensions[i].size;
        index /= dimensions[i].size;
    }
    return coords;
}
 
std::vector<uint64_t> calculate_strides(const std::vector<DataDimension>& dimensions)
{
    std::vector<uint64_t> strides(dimensions.size());
    uint64_t stride = 1;
    for (int i = dimensions.size() - 1; i >= 0; --i) {
        strides[i] = stride;
        stride *= dimensions[i].size;
    }
    return strides;
}
 
bool validate_region_bounds(const Region& region, const std::vector<DataDimension>& dimensions)
{
    if (region.start_coordinates.size() != dimensions.size() || region.end_coordinates.size() != dimensions.size()) {
        return false;
    }
 
    for (size_t i = 0; i < dimensions.size(); ++i) {
        if (region.start_coordinates[i] > region.end_coordinates[i] || region.end_coordinates[i] >= dimensions[i].size) {
            return false;
        }
    }
 
    return true;
}
 
bool validate_slice_bounds(const std::vector<uint64_t>& slice_start,
    const std::vector<uint64_t>& slice_end,
    const std::vector<DataDimension>& dimensions)
{
    if (slice_start.size() != dimensions.size() || slice_end.size() != dimensions.size()) {
        return false;
    }
 
    for (size_t i = 0; i < dimensions.size(); ++i) {
        if (slice_start[i] > slice_end[i] || slice_end[i] >= dimensions[i].size) {
            return false;
        }
    }
 
    return true;
}
 
void clamp_coordinates_to_bounds(std::vector<uint64_t>& coords,
    const std::vector<DataDimension>& dimensions)
{
    for (size_t i = 0; i < coords.size() && i < dimensions.size(); ++i) {
        coords[i] = std::min(coords[i], dimensions[i].size - 1);
    }
}
 
std::vector<uint64_t> transform_coordinates(const std::vector<uint64_t>& coords,
    const std::vector<double>& scale_factors,
    const std::vector<int64_t>& offset_values,
    const std::unordered_map<std::string, std::any>& rotation_params)
{
    std::vector<uint64_t> transformed = coords;
 
    if (!scale_factors.empty()) {
        for (size_t i = 0; i < coords.size() && i < scale_factors.size(); ++i) {
            transformed[i] = static_cast<uint64_t>(scale_factors[i] * static_cast<double>(coords[i]));
        }
    }
 
    if (!offset_values.empty()) {
        for (size_t i = 0; i < transformed.size() && i < offset_values.size(); ++i) {
            int64_t new_coord = static_cast<int64_t>(transformed[i]) + offset_values[i];
            transformed[i] = static_cast<uint64_t>(std::max(int64_t(0), new_coord));
        }
    }
 
    if (!rotation_params.empty() && coords.size() >= 2) {
        auto angle_it = rotation_params.find("angle_radians");
        if (angle_it != rotation_params.end()) {
            try {
                auto angle = safe_any_cast_or_throw<double>(angle_it->second);
                double cos_a = std::cos(angle);
                double sin_a = std::sin(angle);
 
                auto x = static_cast<double>(transformed[0]);
                auto y = static_cast<double>(transformed[1]);
 
                transformed[0] = static_cast<uint64_t>(x * cos_a - y * sin_a);
                transformed[1] = static_cast<uint64_t>(x * sin_a + y * cos_a);
            } catch (const std::bad_any_cast&) {
                // TODO: DO better than ignore invalid
            }
        }
    }
 
    return transformed;
}
 
std::vector<uint64_t> wrap_position_with_loop(
    const std::vector<uint64_t>& positions,
    const Region& loop_region,
    bool looping_enabled)
{
    if (!looping_enabled || loop_region.start_coordinates.size() < 2 || loop_region.end_coordinates.size() < 2) {
        return positions;
    }
 
    uint64_t loop_start_frame = loop_region.start_coordinates[0];
    uint64_t loop_start_channel = loop_region.start_coordinates[1];
    uint64_t loop_end_frame = loop_region.end_coordinates[0];
    uint64_t loop_end_channel = loop_region.end_coordinates[1];
 
    if (loop_end_frame <= loop_start_frame) {
        return positions;
    }
 
    std::vector<uint64_t> wrapped_positions = positions;
    uint64_t loop_length = loop_end_frame - loop_start_frame + 1;
 
    for (size_t ch = loop_start_channel; ch <= loop_end_channel && ch < positions.size(); ++ch) {
        if (positions[ch] > loop_end_frame) {
            uint64_t overflow = positions[ch] - loop_end_frame;
            wrapped_positions[ch] = loop_start_frame + (overflow % loop_length);
        }
    }
 
    return wrapped_positions;
}
 
std::vector<uint64_t> advance_position(
    const std::vector<uint64_t>& current_positions,
    uint64_t frames_to_advance,
    const ContainerDataStructure& structure,
    bool looping_enabled,
    const Region& loop_region)
{
    uint64_t num_channels = structure.get_channel_count();
    uint64_t total_frames = structure.get_samples_count_per_channel();
 
    if (current_positions.size() != num_channels) {
        throw std::invalid_argument(
            "Position vector size " + std::to_string(current_positions.size()) + " must match channel count " + std::to_string(num_channels));
    }
 
    std::vector<uint64_t> new_positions;
    new_positions.reserve(num_channels);
 
    for (size_t ch = 0; ch < num_channels; ++ch) {
        uint64_t current_frame = current_positions[ch];
        uint64_t new_frame = current_frame + frames_to_advance;
 
        if (new_frame >= total_frames) {
            if (looping_enabled && !loop_region.start_coordinates.empty()) {
                uint64_t loop_start = (ch < loop_region.start_coordinates.size())
                    ? loop_region.start_coordinates[ch]
                    : loop_region.start_coordinates[0];
                uint64_t loop_end = (ch < loop_region.end_coordinates.size())
                    ? loop_region.end_coordinates[ch]
                    : loop_region.end_coordinates[0];
 
                if (loop_end > loop_start) {
                    uint64_t loop_length = loop_end - loop_start + 1;
                    uint64_t overflow = new_frame - total_frames;
                    new_frame = loop_start + (overflow % loop_length);
                } else {
                    new_frame = loop_start;
                }
            } else {
                new_frame = total_frames - 1;
            }
        }
 
        new_positions.push_back(new_frame);
    }
 
    return new_positions;
}
 
std::vector<uint64_t> advance_position(
    const std::vector<uint64_t>& current_positions,
    const std::vector<uint64_t>& frames_per_channel,
    const ContainerDataStructure& structure,
    bool looping_enabled,
    const Region& loop_region)
{
    uint64_t num_channels = structure.get_channel_count();
    uint64_t total_frames = structure.get_samples_count_per_channel();
 
    if (current_positions.size() != num_channels || frames_per_channel.size() != num_channels) {
        throw std::invalid_argument("All vectors must match channel count");
    }
 
    std::vector<uint64_t> new_positions;
    new_positions.reserve(num_channels);
 
    for (size_t ch = 0; ch < num_channels; ++ch) {
        uint64_t current_frame = current_positions[ch];
        uint64_t frames_to_advance = frames_per_channel[ch];
        uint64_t new_frame = current_frame + frames_to_advance;
 
        if (new_frame >= total_frames) {
            if (looping_enabled && !loop_region.start_coordinates.empty()) {
                uint64_t loop_start = (ch < loop_region.start_coordinates.size())
                    ? loop_region.start_coordinates[ch]
                    : loop_region.start_coordinates[0];
                uint64_t loop_end = (ch < loop_region.end_coordinates.size())
                    ? loop_region.end_coordinates[ch]
                    : loop_region.end_coordinates[0];
 
                if (loop_end > loop_start) {
                    uint64_t loop_length = loop_end - loop_start + 1;
                    uint64_t overflow = new_frame - total_frames;
                    new_frame = loop_start + (overflow % loop_length);
                } else {
                    new_frame = loop_start;
                }
            } else {
                new_frame = total_frames - 1;
            }
        }
 
        new_positions.push_back(new_frame);
    }
 
    return new_positions;
}
 
uint64_t time_to_position(double time, double sample_rate)
{
    return static_cast<uint64_t>(time * sample_rate);
}
 
double position_to_time(uint64_t position, double sample_rate)
{
    return static_cast<double>(position) / sample_rate;
}
 
uint64_t calculate_frame_size_for_dimension(const std::vector<DataDimension>& dimensions, size_t primary_dim)
{
    if (dimensions.empty() || primary_dim >= dimensions.size()) {
        return 0;
    }
    uint64_t frame_size = 1;
    for (size_t i = 0; i < dimensions.size(); ++i) {
        if (i != primary_dim) {
            frame_size *= dimensions[i].size;
        }
    }
    return frame_size;
}
 
std::unordered_map<std::string, std::any> create_coordinate_mapping(const std::shared_ptr<SignalSourceContainer>& container)
{
    if (!container) {
        throw std::invalid_argument("Container is null");
    }
 
    std::unordered_map<std::string, std::any> mapping_info;
    const auto dimensions = container->get_dimensions();
 
    std::vector<std::unordered_map<std::string, std::any>> dim_mappings;
    dim_mappings.reserve(dimensions.size());
 
    for (size_t i = 0; i < dimensions.size(); ++i) {
        std::unordered_map<std::string, std::any> dim_map;
        dim_map["index"] = i;
        dim_map["name"] = dimensions[i].name;
        dim_map["size"] = dimensions[i].size;
        dim_map["stride"] = dimensions[i].stride;
        dim_map["role"] = static_cast<int>(dimensions[i].role);
 
        uint64_t offset = (i == 0) ? 0ULL : std::accumulate(dimensions.begin(), dimensions.begin() + i, 1ULL, [](uint64_t acc, const auto& d) { return acc * d.size; });
        dim_map["offset"] = offset;
 
        dim_mappings.push_back(std::move(dim_map));
    }
 
    mapping_info["dimensions"] = std::move(dim_mappings);
    mapping_info["total_elements"] = container->get_total_elements();
    mapping_info["memory_layout"] = static_cast<int>(container->get_memory_layout());
 
    auto strides = calculate_strides(dimensions);
    mapping_info["calculated_strides"] = strides;
 
    return mapping_info;
}
 
std::vector<int> extract_dimension_roles(const std::vector<DataDimension>& dimensions)
{
    std::vector<int> roles;
    roles.reserve(dimensions.size());
 
    for (const auto& dim : dimensions) {
        roles.push_back(static_cast<int>(dim.role));
    }
 
    return roles;
}
 
std::vector<uint64_t> extract_dimension_sizes(const std::vector<DataDimension>& dimensions)
{
    std::vector<uint64_t> sizes;
    sizes.reserve(dimensions.size());
 
    for (const auto& dim : dimensions) {
        sizes.push_back(dim.size);
    }
 
    return sizes;
}
 
std::vector<std::unordered_map<std::string, std::any>> create_dimension_info(const std::vector<DataDimension>& dimensions)
{
    std::vector<std::unordered_map<std::string, std::any>> dim_info;
    dim_info.reserve(dimensions.size());
 
    for (const auto& dim : dimensions) {
        std::unordered_map<std::string, std::any> info;
        info["name"] = dim.name;
        info["size"] = dim.size;
        info["stride"] = dim.stride;
        info["role"] = static_cast<int>(dim.role);
        dim_info.push_back(std::move(info));
    }
 
    return dim_info;
}
 
std::pair<size_t, uint64_t> coordinates_to_planar_indices(
    const std::vector<uint64_t>& coords,
    const std::vector<DataDimension>& dimensions)
{
    size_t channel_dim_idx = 0;
    size_t time_dim_idx = 0;
 
    for (size_t i = 0; i < dimensions.size(); ++i) {
        if (dimensions[i].role == DataDimension::Role::CHANNEL) {
            channel_dim_idx = i;
        } else if (dimensions[i].role == DataDimension::Role::TIME) {
            time_dim_idx = i;
        }
    }
 
    return { coords[channel_dim_idx], coords[time_dim_idx] };
}
 
}