add Pulseq MATLAB calcAdcSeg as calc_adc_segments

[felixlandmeyer]

I have created a Python implementation of calcAdcSeg on https://iffgit.fz-juelich.de/f.landmeyer/pypulseq/-/tree/add_calc_adc_segments.
The code is a one-to-one translation of the Matlab code.
As far as I know, such an implementation is not yet published in pypulseq, so I can create a pull-request if desired.
I have tested the output of calc_adc_segments against the output of calcAdcSeg with the output of the attached Matlab code and it is identical for all parameters tested. The output of the Matlab code is saved in test/expected_output and read and tested by test_calc_adc_segments.

%% System Setup
sys = mr.opts('MaxGrad',30,'GradUnit','mT/m',...
    'MaxSlew',120,'SlewUnit','T/m/s',...
    'rfRingdownTime', 30e-6, 'rfDeadtime', 100e-6,...
    'adcDeadTime', 10e-6, 'adcSamplesLimit', 8192, ...
    'adcSamplesDivisor', 4, 'adcRasterTime', 1e-7,...
    'gradRasterTime', 1e-5);

%% Create test parameters
rng(42);
dwell_list = randi([10, 100], 1, 40) * 1e-7;  
num_samples_list = randi([200 12000], 1, 100);  
adc_limit_list = [2048, 4096, 8192];  
adc_divisor_list = [1, 2, 4];

%% Preallocate the output array
total_combinations = (...
    length(dwell_list) ...
    * length(num_samples_list) ...
    * length(adc_limit_list) ...
    * length(adc_divisor_list) ...
    * 2 ... % for modes
    );
output = zeros(total_combinations, 7);

%% Create segmented ADC for all combinations
counter = 1;
for dwell_idx = 1:length(dwell_list)
    dwell = dwell_list(dwell_idx);
    for num_samples_idx = 1:length(num_samples_list)
        num_samples = num_samples_list(num_samples_idx);
        for adc_limit_idx = 1:length(adc_limit_list)
            adc_limit = adc_limit_list(adc_limit_idx);
            sys.adcSamplesLimit = adc_limit;
            for adc_divisor_idx = 1:length(adc_divisor_list)
                adc_divisor = adc_divisor_list(adc_divisor_idx);
                for n_mode = 1:2
                    if n_mode == 1
                        mode = 'shorten';
                    else
                        mode = 'lenghen';
                    end

                   
                    sys.adcSamplesDivisor = adc_divisor;
                
                    [num_seg, num_samples_seg] = mr.calcAdcSeg(num_samples, dwell, sys, mode);
                    
                    output(counter, :) = [dwell, num_samples, adc_limit, adc_divisor, n_mode, num_seg, num_samples_seg];
                
                    if mod(counter, 1000) == 0 || counter == total_combinations
                        fprintf('Progress: %.2f%%\n', (counter / total_combinations) * 100);
                    end

                    counter = counter + 1;
                end
            end
        end
    end
end

%% Write to txt file
writematrix(output, 'pulseq_calcAdcSeg.txt', 'Delimiter', 'tab');

disp('Data written to pulseq_calcAdcSeg.txt');

[btasdelen]

Hi @felixlandmeyer, the GitLab link you provide requires a login (maybe repo is private?).

In any case, if you create a pull request, we can review the implementation.

[felixlandmeyer]

I'm sorry. It should be public now.