jfjoch_process: Fix

This commit is contained in:
2026-02-08 19:29:00 +01:00
parent 69ebbed8fb
commit 58f29409e6

View File

@@ -428,16 +428,11 @@ int main(int argc, char **argv) {
// --- French-Wilson: convert I → F ---
{
// Build d-spacings vector parallel to merged
std::vector<double> d_spacings(scale_result->merged.size());
for (size_t i = 0; i < scale_result->merged.size(); ++i)
d_spacings[i] = scale_result->merged[i].d;
FrenchWilsonOptions fw_opts;
fw_opts.acentric = true; // typical for MX
fw_opts.num_bins = 20;
fw_opts.num_shells = 20;
auto fw = FrenchWilson(scale_result->merged, d_spacings, fw_opts);
auto fw = FrenchWilson(scale_result->merged, fw_opts);
const std::string fw_path = output_prefix + "_amplitudes.hkl";
std::ofstream fw_file(fw_path);