remove unused code

Author: Cam-Sim

+nodes/QkdPassSimulation.m

@@ -122,200 +122,6 @@
         end
 end
 
-
-function needs_env = map_receivers_and_environments(receiver, options)
-    arguments
-        receiver { ...
-            nodes.mustBeReceiverOrTransmitter(receiver), ...
-            nodes.mustHaveDetector(receiver) }
-        options.Environment environment.Environment
-    end
-
-    n_rx = numel(receiver);
-    n_env = 0;
-    have_environment = ismember(fieldnames(options), "Environment");
-    if have_environment
-        n_env = numel(options.Environment);
-    end
-
-    needs_env = arrayfun(@(maybe_has_env) isa(maybe_has_env, 'nodes.Ground_Station'), receiver);
-
-    if n_rx >= n_env && n_env > 1 && have_environment
-        % convert logical values to indices
-        needs_env = needs_env .* cumsum(needs_env);
-
-        % so we can catch them later set "0" values to -1
-        needs_env(needs_env == 0) = -1;
-
-        assert(max(needs_env) == numel(options.Environment), ...
-            [ ...
-            'Number of receivers and environments do not match, either supply ', ...
-            'a single environment for all receivers or an environment per receiver' ...
-            ])
-    end
-end
-
-
-function [min_val, max_val] = extrema(array)
-    min_val = min(array);
-    max_val = max(array);
-end
-
-
-function loc_time = location_time(tx_idx, rx_idx, direction, heading, elevation, range, elev_limit, elev_mask, time)
-    arguments
-        tx_idx {mustBeNumeric}
-        rx_idx {mustBeNumeric}
-        direction nodes.LinkDirection
-        heading (1, :) {mustBeNumeric}
-        elevation (1, :) {mustBeNumeric}
-        range (1, :) {mustBeNumeric}
-        elev_limit {mustBeScalarOrEmpty, mustBeNumeric}
-        elev_mask (1, :) {mustBeNumericOrLogical}
-        time (1, :) {mustBeA(time, "datetime")}
-    end
-
-    loc_time = struct( ...
-        "tx_idx", tx_idx, ...
-        "rx_idx", rx_idx, ...
-        "direction", direction, ...
-        "heading", heading, ...
-        "elevation", elevation, ...
-        "range", range, ...
-        "elevation_limit", elev_limit, ...
-        "elevation_mask", elev_mask, ...
-        "time", time);
-end
-
-
-function tt = timetable_for_elevation_window(rel_loc, elev_limit)
-    elev_mask = rel_loc.elevation > elev_limit;
-    if sum(elev_mask) < 1
-        tt = [];
-        return
-    end
-    [t_min, t_max] = extrema(rel_loc.time(elev_mask));
-    tt = timetable([t_min; t_max]);
-end
-
-
-function [single_link_idx, double_link_idx] = find_links(receivers, transmitters)
-
-    arguments
-        receivers { ...
-            nodes.mustBeReceiverOrTransmitter(receivers), ...
-            nodes.mustHaveDetector(receivers) }
-        transmitters { ...
-            nodes.mustBeReceiverOrTransmitter(transmitters), ...
-            nodes.mustHaveSource(transmitters) }
-    end
-
-    single_link_idx = {};
-    double_link_idx = {};
-
-    n_receivers = numel(receivers);
-
-    loc_time = struct( ...
-        "tx_idx", 0, ...
-        "rx_idx", 0, ...
-        "direction", nodes.LinkDirection.Downlink, ...
-        "heading", [], ...
-        "elevation", [], ...
-        "range", [], ...
-        "elevation_limit", 0, ...
-        "elevation_mask", [], ...
-        "time", []);
-
-    % relative_locations = createArray(1, n_receivers, Like=loc_time(0, 0, [], [], [], [], []));
-    % relative_locations = {};
-
-    for T = 1:numel(transmitters)
-        tx = transmitters(T);
-        i = 1;
-        valid_links = zeros(1, n_receivers); % for current transmitter
-        elevation_limits = zeros(1, n_receivers);
-
-        % relative_locations = createArray(1, n_receivers, ...
-        %     Like=loc_time(0, 0, 0, [], [], [], 0, [], []));
-
-        relative_locations = createArray(1, n_receivers, Like=loc_time);
-        for R = 1:numel(receivers)
-            rx = receivers(R);
-            direction = nodes.LinkDirection.DetermineLinkDirection(rx, tx);
-
-            switch direction
-                case nodes.LinkDirection.Downlink
-                    [h, e, r] = tx.RelativeHeadingAndElevation(rx);
-                    limit = rx.Elevation_Limit;
-                    elevation_limits(i) = limit;
-                    elevation_limit_mask = e > elevation_limits(i);
-                    rel_loc = location_time(T, R, direction, h, e, r, limit, elevation_limit_mask, tx.Times);
-
-                case nodes.LinkDirection.Uplink
-                    [h, e, r] = rx.RelativeHeadingAndElevation(tx);
-                    limit = tx.Elevation_Limit;
-                    elevation_limits(i) = limit;
-                    elevation_limit_mask = e > elevation_limits(i);
-                    rel_loc = location_time(T, R, direction, h, e, r, limit, elevation_limit_mask, rx.Times);
-            end
-
-            % elevation_limit_mask = relative_locations(i).elevation > elevation_limits(i);
-            % if sum(elevation_limit_mask) > 0
-            if sum(rel_loc.elevation_mask) > 0
-                valid_links(i) = 1;
-                single_link_count = numel(single_link_idx) + 1;
-                single_link_idx{single_link_count} = rel_loc;
-                %relative_locations_count = numel(relative_locations) + 1;
-                % relative_locations{relative_locations_count} = rel_loc;
-                relative_locations(R) = rel_loc;
-            end
-
-            i = i + 1;
-        end
-
-        relative_locations = relative_locations(arrayfun(@(rl) ~isempty(rl.elevation), relative_locations));
-
-        %if sum(valid_links) < 2
-        if numel(relative_locations) < 2
-            % skip to next transmitter if there is one
-            continue
-        end
-
-        %n = sum(valid_links);
-        n = numel(relative_locations);
-        for rx_pairs = nchoosek(linspace(1, n, n), 2)'
-
-            rel_loc_i = relative_locations(rx_pairs(1));
-            e_lim = elevation_limits(rx_pairs(1));
-            tt_i = timetable_for_elevation_window(rel_loc_i, e_lim);
-
-            rel_loc_j = relative_locations(rx_pairs(2));
-            e_lim = elevation_limits(rx_pairs(2));
-            tt_j = timetable_for_elevation_window(rel_loc_j, e_lim);
-
-            if any([sum(size(tt_i)) < 1, sum(size(tt_j)) < 1])
-                % one of the timetables is empty so it mustn't have seen the satellite
-                continue
-            end
-
-            if ~overlapsrange(tt_i, tt_j)
-                %{
-                    skip to next pair
-                    eventhough both rx, tx pairs could see each other they
-                    couldn't see each other at the same time
-                %}
-                continue
-            end
-
-            dl_idx = numel(double_link_idx) + 1;
-            % double_link_idx{dl_idx} = {T, rx_pairs(1), rx_pairs(2)}; % {transmitter, receiver_i, receivers_j}
-            double_link_idx{dl_idx} = {[rel_loc_i, rel_loc_j]};
-        end
-
-    end
-end
-
-
 function [loss_results, noise] = loss_and_noise_for_channel(transmitter, receiver, qkd_protocol)
     arguments
         transmitter (1,1) { ...