31 lines
1004 B
C++
31 lines
1004 B
C++
// SPDX-FileCopyrightText: 2024 Filip Leonarski, Paul Scherrer Institute <filip.leonarski@psi.ch>
|
|
// SPDX-License-Identifier: GPL-3.0-only
|
|
|
|
|
|
#pragma once
|
|
|
|
#include <algorithm>
|
|
#include <string>
|
|
#include <vector>
|
|
#include <unordered_map>
|
|
#include <cstdint>
|
|
|
|
inline uint64_t hkl_key_16(int64_t h, int64_t k, int64_t l) {
|
|
const uint16_t bias = 512; // maps -512..512 -> 0..1024
|
|
uint64_t uh = static_cast<uint16_t>(std::clamp(h + bias, int64_t(0), int64_t(1024)));
|
|
uint64_t uk = static_cast<uint16_t>(std::clamp(k + bias, int64_t(0), int64_t(1024)));
|
|
uint64_t ul = static_cast<uint16_t>(std::clamp(l + bias, int64_t(0), int64_t(1024)));
|
|
return uh | (uk << 16) | (ul << 32);
|
|
}
|
|
|
|
struct HKLData {
|
|
int64_t h ,k,l;
|
|
double I = 0.0;
|
|
double sigma = 0.0;
|
|
float last_image = 0;
|
|
std::vector<double> tail; // any remaining numeric columns
|
|
};
|
|
|
|
using IntegrateMap = std::unordered_map<uint64_t, std::vector<HKLData>>;
|
|
|
|
IntegrateMap ParseXdsIntegrateHkl(const std::string& filename); |