matlabInterpolation/filterStructs.m at master · weedle/matlabInterpolation · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
function [ cellStructs, structs ] = filterStructs( structs, threshold )
    fields = fieldnames(structs);

    % Sort by slopes
    [ cellStructs, ~ ] = sortStructs( structs, 2 );


    sz = size(cellStructs);
    display( sprintf( '\tCurrently have %d mutants', sz(3) ) );

    % Remove mutants with higher slopes
    % Anything that's higher by an order or more can be rejected
    % Filtering by slopes first prevents us from locking on to a level
    % plane an order or more above the actual best slope
    cellStructs = filterByIndex( cellStructs, 2, mode( [ cellStructs{2,:} ] )+(threshold-1) );
    sz = size(cellStructs);
    display( sprintf( '\t\tFiltered by slope: now have %d mutants', sz(3) ) );


    % Remove all mutants with higher mean errors (by one order)
    cellStructs = filterByIndex( cellStructs, 3, 100*threshold*mode( [ cellStructs{3,:} ] ) );

    % Remove all mutants with higher maximum errors (by one order)
    cellStructs = filterByIndex( cellStructs, 4, 100*threshold*mode( [ cellStructs{4,:} ] ) );

    sz = size(cellStructs);
    display( sprintf( '\t\tFiltered by max: now have %d mutants', sz(3) ) );

    % Remove mutants with higher residuals
    cellStructs = filterByIndex( cellStructs, 5, threshold*mode( [ cellStructs{5,:} ] ) );
    sz = size(cellStructs);
    display( sprintf( '\t\tFiltered by residual: now have %d mutants', sz(3) ) );

    structs = cell2struct(cellStructs, fields, 1);
end

function cellStructs = filterByIndex( cellStructs, index, threshold )
    cellStructs = removeInvalid( cellStructs, index, threshold );
end

function cellStructs = removeInvalid( cellStructs, index, threshold )
    slope = 0;
    if( index == 2 )
        slope = 1;
    end
    values = [ cellStructs{ index,: } ];
    cellStructs( :,:,isnan( values ) ) = [];
    if( slope )
        cellStructs( :,:,( [ cellStructs{ index,1:length( cellStructs ) } ] > threshold ) ) = [];
    else
        cellStructs( :,:,( [ cellStructs{ index,1:length( cellStructs ) } ] > threshold ) ) = [];
    end
end

function cellStructs = toCells( structs ) %#ok<DEFNU>
    cellStructs = struct2cell(structs);
    s = size(cellStructs);
    cellStructs = reshape(cellStructs, s(1), []);
    cellStructs = reshape(cellStructs, s);
end