datafarming 1.0.0 → 1.1.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 5aa6e5a57e29e45261e3322740ae3694c0c5eb7e32f841b71a15fb4da631d2d1
-  data.tar.gz: c7d299e10278175d1a87020f9e84bc91d697754ab981c1af07174d25c8467e9a
+  metadata.gz: 190ba17087fafc8f7308e5e4decbe18b412fb99342c4d55011ac28b01735d329
+  data.tar.gz: 2c9da6d6824f2899904f076a4d6184279a272b04fd3fc979538f8809dde83328
 SHA512:
-  metadata.gz: fffc9ed4734f4d1a8ee9ebd18bacbe7f9f679ade49b8910e61ac751de2adb8f7ec5bd2475b7572597d31307f8ab08d425fb2db176438527a4c68264d0e8fdbd6
-  data.tar.gz: 8bed8adb0b6d0e1f1ad19ad5b3e63517d590ce07e8128e18c4eb24d9545b3c03ac68fcd5bd26ec497f643995c4b31563ed310162c5df38da22669ea1405376d8
+  metadata.gz: 77347491591b5d18612832d102618a84897e4be99e96116b946e75b09a1b71e1bad19d2bf44fafa4afff122cce4a3146ad9113731147e27c6e3fb3e47173c8bb
+  data.tar.gz: c1a8f5d16b081d94731abb8e7c0cf6e6550a4e9929ddd6218e2b5e81cca25f9cf66789400173b510cd3c50143e8379d95966888e4d844c1651880c8aac0449d0

data/LICENSE.md

@@ -0,0 +1,22 @@
+
+The MIT License (MIT)
+
+Copyright (c) 2017-2018 Paul J. Sanchez
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

data/datafarming.gemspec

@@ -1,10 +1,10 @@
 # -*- ruby -*-
-_VERSION = "1.0.0"
+_VERSION = "1.1.0"
 
 Gem::Specification.new do |s|
   s.name = "datafarming"
   s.version = _VERSION
-  s.date = "2018-06-25"
+  s.date = "2018-07-15"
   s.summary = "Useful scripts for data farming."
   s.homepage = "https://gitlab.nps.edu/pjsanche/datafarmingrubyscripts.git"
   s.email = "pjs@alum.mit.edu"

data/exe/scaled_fde.rb

@@ -0,0 +1,169 @@
+#!/usr/bin/env ruby -w
+
+require 'colorize'
+String.disable_colorization false
+
+require 'datafarming/error_handling'
+require 'datafarming/freq_sets'
+
+help_msg = [
+  'Generate scaled FDE designs with shifting and stacking. The resulting',
+  'design is written to ' +
+    'stdout'.light_blue + '.', '',
+  'Syntax:',
+  "\n\t#{$PROGRAM_NAME.split(%r{/|\\})[-1]} [--help]".yellow +
+    " [-c] [-d #] [-e] [-k #] [-s #] [file_name]\n".yellow,
+  "Arguments in square brackets are optional.  A vertical bar '|'",
+  'indicates valid alternatives for invoking the option.  Prefix',
+  'the command with "' + 'ruby'.yellow +
+    '" if it is not on your PATH.', '',
+  '  --help | -h | -? | ?'.green,
+  "\tProduce this help message.  Supersedes any other choices.",
+  '  --csv | -c'.green,
+  "\tGenerate output as comma-separated-values.",
+  '  --num_factors # | -k #'.green,
+  "\tSpecify number of factors using standard (+/-1) scaling.",
+  '  --design # | -d #'.green,
+  "\tSpecify which design number you want.  Default is highest",
+  "\tfrequency set available for the target number of factors.",
+  "\tAn input larger than the available number will produce an",
+  "\terror message that identifies the allowable range.",
+  '  --stack # | -s #'.green,
+  "\t# specifies the number of stackings. The specified value cannot",
+  "\texceed the number of columns in the design being used. If this",
+  "\toption is omitted then only the base design is generated.",
+  '  --excel-style-input | -e'.green,
+  "\tSpecify factor ranges similarly to the NOLH spreadsheet, i.e.,",
+  "\tthe first line is the set of minimum range values for all factors,",
+  "\tand the second line is maximum range values.  Without this option,",
+  "\tthe default input format is one line per factor, comprised of the",
+  "\tmin and max separated by commas or whitespace.",
+  '  file_name'.green,
+  "\tThe name of a file containing the factor specifications.  If no",
+  "\tfilename is given, the user can enter the values interactively in",
+  "\tthe desired form or use file redirection with '<'.", '',
+  'Options may be given in any order, but must come before the file name',
+  'if one is provided.'
+]
+
+# Scaler objects will rescale a FDE design from standard units
+# to a range specified by min and max
+class Scaler
+  def initialize(min = -1, max = 1)
+    @mid = 0.5 * (max + min)
+    @amplitude = 0.5 * (max - min)
+  end
+
+  def scale(value)
+    @mid + @amplitude * value
+  end
+end
+
+def make_design(number_of_factors, design_num = -1)
+  design = Array.new(number_of_factors)
+  ds = FDE::DESIGN_SETS[number_of_factors]
+  freqs = ds.freqs[design_num].map { |f| Rational(f, ds.nyq) }
+  design.each_index do |column|
+    omega = freqs[column]
+    design[column] = Array.new(ds.nyq) do |i|
+      Math.sin(((i * omega) % 1) * FDE::TWO_PI)
+    end
+  end
+  design.transpose
+end
+
+excel_style_inputs = false
+csv_out = false
+design_num = -1
+num_stackings = 1
+n = nil
+while ARGV[0] && (ARGV[0][0] == '-' || ARGV[0][0] == 45 || ARGV[0][0] == '?')
+  current_value = ARGV.shift
+  case current_value
+  when '--csv', '-c'
+    csv_out = true
+  when '--design', '-d'
+    design_num = ARGV.shift.to_i - 1
+  when '--num_factors', '-k'
+    n = ARGV.shift.to_i
+  when '--stack', '-s'
+    num_stackings = ARGV.shift.to_i
+  when '--excel-style-input', '-e'
+    excel_style_inputs = true
+  when '--help', '-h', '-help', '-?', '?'
+    ErrorHandling.clean_abort help_msg
+  else
+    ErrorHandling.message ['Unknown argument: '.red + current_value.yellow]
+    ErrorHandling.clean_abort help_msg
+  end
+end
+
+if n
+  factor = Array.new(n) { Scaler.new }
+else
+  begin
+    if excel_style_inputs
+      if ARGF.filename == '-'
+        STDERR.puts  'Enter one line of min values, and one of max values.'.green
+      end
+      min_values = ARGF.gets.strip.split(/\s*[,;:]\s*|\s+/).map(&:to_f)
+      max_values = ARGF.gets.strip.split(/\s*[,;:]\s*|\s+/).map(&:to_f)
+    else
+      if ARGF.filename == '-'
+        STDERR.puts  'To terminate input enter '.green + 'ctrl-d'.cyan +
+                    ' (Mac/Unix/Linux)'.green + ' or '.green + 'ctrl-z'.cyan +
+                    ' (Windows).'.green
+        STDERR.puts  'Enter ranges for each factor on a separate line.'.green
+        STDERR.puts  "\nMIN\tMAX".cyan
+      end
+      min_values = []
+      max_values = []
+      while line = ARGF.gets
+        values = line.strip.split(/\s*[,;:]\s*|\s+/).map(&:to_f)
+        min_values << values.shift
+        max_values << values.shift
+      end
+    end
+  rescue StandardError => e
+    ErrorHandling.message [e.message.red]
+    ErrorHandling.clean_abort help_msg
+  end
+
+  n = min_values.size
+  if max_values.size != n
+    ErrorHandling.message ['Unequal counts for min, max'.red]
+    ErrorHandling.clean_abort help_msg
+  end
+
+  factor = Array.new(n) do |i|
+    Scaler.new(min_values[i], max_values[i])
+  end
+end
+
+if design_num >= FDE::DESIGN_SETS[n].freqs.length
+  ErrorHandling.clean_abort [
+    "Design number cannot exceed #{FDE::DESIGN_SETS[n].freqs.length} for #{n} factors".red
+  ]
+end
+
+design = make_design(n, design_num)
+
+num_columns = design[0].length
+if num_stackings > num_columns
+  ErrorHandling.clean_abort [
+    'Requested stacking exceeds number of columns in FDE '.red +
+    "(#{num_columns})".red
+  ]
+end
+
+num_stackings.times do |_stack_num|
+  design.each_with_index do |dp, i|
+    scaled_dp = dp.slice(0, n).map.with_index { |x, k| factor[k].scale(x) }
+    if csv_out
+      puts scaled_dp.join ','
+    else
+      puts scaled_dp.join "\t"
+    end
+    design[i] = dp.rotate
+  end
+end

data/exe/scaled_rf_cubed.rb

@@ -0,0 +1,227 @@
+#!/usr/bin/env ruby -w
+
+require 'colorize'
+String.disable_colorization false
+
+require 'datafarming/error_handling'
+require 'fwt'
+
+# Generate a Resolution V Fractional Factorial design for the specified
+# number of factors.  The design is guaranteed to yield unconfounded
+# interactions for all pairs of factors.  The design uses standardized
+# notation, i.e., -1 represents a low setting and 1 represents a high
+# setting for each factor.
+#
+# *Arguments*::
+#   - +number_of_factors+ -> the number of factors in your design. Limit is 120.
+# *Returns*::
+#   - a two-dimensional array specifying the design, where each column
+#     corresponds to a factor and each row is a design point.
+#
+# Author:: Paul J Sanchez (mailto:pjs@alum.mit.edu)
+# Copyright:: Copyright (c) Paul J Sanchez
+# License:: LGPL
+#
+def make_design(number_of_factors)
+  index = [1, 2, 4, 8, 15, 16, 32, 51, 64, 85, 106, 128, 150, 171,
+           219, 237, 247, 256, 279, 297, 455, 512, 537, 557, 594, 643, 803,
+           863, 998, 1024, 1051, 1070, 1112, 1169, 1333, 1345, 1620, 1866,
+           2048, 2076, 2085, 2185, 2372, 2456, 2618, 2800, 2873, 3127, 3284,
+           3483, 3557, 3763, 4096, 4125, 4135, 4174, 4435, 4459, 4469, 4497,
+           4752, 5255, 5732, 5804, 5915, 6100, 6369, 6907, 7069, 8192, 8263,
+           8351, 8422, 8458, 8571, 8750, 8858, 9124, 9314, 9500, 10_026,
+           10_455, 10_556, 11_778, 11_885, 11_984, 13_548, 14_007, 14_514,
+           14_965, 15_125, 15_554, 16_384, 16_457, 16_517, 16_609, 16_771,
+           16_853, 17_022, 17_453, 17_891, 18_073, 18_562, 18_980, 19_030,
+           19_932, 20_075, 20_745, 21_544, 22_633, 23_200, 24_167, 25_700,
+           26_360, 26_591, 26_776, 28_443, 28_905, 29_577, 32_705
+         ]
+  power = [1, 2, 4, 8, 16, 16, 32, 64, 64, 128, 128, 128, 256, 256,
+           256, 256, 256, 256, 512, 512, 512, 512, 1024, 1024, 1024, 1024,
+           1024, 1024, 1024, 1024, 2048, 2048, 2048, 2048, 2048, 2048, 2048,
+           2048, 2048, 4096, 4096, 4096, 4096, 4096, 4096, 4096, 4096, 4096,
+           4096, 4096, 4096, 4096, 4096, 8192, 8192, 8192, 8192, 8192, 8192,
+           8192, 8192, 8192, 8192, 8192, 8192, 8192, 8192, 8192, 8192, 8192,
+           16_384, 16_384, 16_384, 16_384, 16_384, 16_384, 16_384, 16_384,
+           16_384, 16_384, 16_384, 16_384, 16_384, 16_384, 16_384, 16_384,
+           16_384, 16_384, 16_384, 16_384, 16_384, 16_384, 16_384, 32_768,
+           32_768, 32_768, 32_768, 32_768, 32_768, 32_768, 32_768, 32_768,
+           32_768, 32_768, 32_768, 32_768, 32_768, 32_768, 32_768, 32_768,
+           32_768, 32_768, 32_768, 32_768, 32_768, 32_768, 32_768, 32_768,
+           32_768, 32_768, 32_768
+         ]
+
+  design = Array.new(number_of_factors)
+  design.each_index do |i|
+    design[i] = Array.new(power[number_of_factors], 0)
+    design[i][index[i]] = 1
+    design[i].hadamard
+  end
+  design.transpose
+end
+
+def print_descriptor
+  description = [
+    "\nCreate a Resolution V fractional factorial design, either in standard +/-1",
+    'notation or scaled as per user-provided input. The design can optionally ',
+    'be supplemented with CCD center and star points.  Output is written to',
+    'stdout'.blue + ' with formatting controlled by command-line options.'
+  ]
+  STDERR.puts description
+end
+
+def print_directions # :nodoc:
+  prog = $PROGRAM_NAME.split(/\//)[-1]
+  help_msg = [
+    '', 'Syntax:', '',
+    "  ruby #{prog} [Options] #factors | low/high-settings".yellow, '',
+    "The vertical bar '|' indicates valid alternatives.", '',
+    "Command-line args must specify either:",
+    '  - the number of factors for the design; or',
+    '  - low and hi values for all factors.', '',
+    "\nOptions:",
+    '  --help | -h | -? | ?'.green,
+    "\tProduce this help message.",
+    "  --design-only | -d".green,
+    "\tprint the design without labeling",
+    "  --csv | -c".green,
+    "\tprint the design as comma separated",
+    "  --center | -a".green,
+    "\tadd global center pt to check linearity",
+    "  --star | -s".green,
+    "\taugment with Central Composite Design pts (includes --center)",
+    "  --rotatable | -r".green,
+    "\tscale a CCD design to be rotatable (includes --star)",
+    ''
+  ]
+  STDERR.puts help_msg
+end
+
+print_descriptor && print_directions && return if ARGV.length == 0
+begin
+  labels = true
+  csv = false
+  scaled = false
+  center = false
+  star = false
+  rotatable = false
+  fmt_size = '5'
+  fmt_type = 'd'
+
+  while ARGV[0] =~ /^(-\D|\?)/
+    arg = ARGV.shift
+    case arg
+    when '-c', '--csv'
+      csv = true
+    when '-d', '--design-only'
+      labels = false
+    when '-a', '--center'
+      center = true
+    when '-s', '--star'
+      star = true
+      center = true
+    when '-r', '--rotatable'
+      rotatable = true
+      star = true
+      center = true
+      fmt_size = '11.8'
+      fmt_type = 'f'
+    when '-h', '--help', '-?', '?'
+      print_descriptor
+      print_directions
+      return
+    else
+      raise "Unknown option #{arg}"
+    end
+  end
+
+  if ARGV.length == 0
+    raise 'Must specify a design size'
+  elsif ARGV.length == 1
+    n = ARGV.shift.to_i
+  elsif (ARGV.length & 1).nonzero? # Odd number of args remaining!
+    raise 'Number of levels is odd: low/hi values must occur in pairs'
+  else
+    n = ARGV.length / 2
+    scaled = true
+    unless rotatable
+      fmt_type = 's'
+      fmt_size = [ARGV.map(&:length).max + 1, fmt_size.to_i].max
+    end
+  end
+
+  my_design = make_design n
+  if rotatable && star
+    inverse_len = 1.0 / Math.sqrt(n)
+    my_design.each_with_index do |row, i|
+      my_design[i] = row.map { |elt| elt * inverse_len }
+    end
+  end
+  my_design += [Array.new(n, 0)] if center
+
+  if star
+    star_design = Array.new(2 * n)
+    star_design.each_index { |row| star_design[row] = Array.new(n, 0) }
+    n.times do |col|
+      row = 2 * col
+      star_design[row][col] = -1
+      star_design[row + 1][col] = 1
+    end
+    my_design += star_design
+  end
+
+  if scaled
+    temp = my_design.transpose
+    temp.each_with_index do |factor, idx|
+      factor.map! do |value|
+        case value
+        when -1
+          ARGV[2 * idx + 1]
+        when 1
+          ARGV[2 * idx]
+        when 0
+          begin
+            0.5 * (Float(ARGV[2 * idx + 1]) + Float(ARGV[2 * idx]))
+          rescue
+            raise 'All factors must be numeric for augmented designs.'
+          end
+        else
+          begin
+            lo = Float(ARGV[2 * idx])
+            hi = Float(ARGV[2 * idx + 1])
+            range = hi - lo
+            0.5 * (hi + lo + value * range)
+          rescue
+            raise 'All factors must be numeric for augmented designs.'
+          end
+        end
+      end
+    end
+    my_design = temp.transpose
+  end
+
+  if labels
+    if csv
+      print ' DP#'
+      n.times { |i| print ",X#{i + 1}" }
+    else
+      print '  DP#'
+      n.times { |i| printf " %#{fmt_size}s", "X#{i + 1}" }
+    end
+    puts
+  end
+
+  my_design.each_with_index do |design_point, row_number|
+    if csv
+      printf('%d,', row_number + 1) if labels
+      puts design_point.join(',')
+    else
+      printf('%5d', row_number + 1) if labels
+      design_point.each { |value| printf " %#{fmt_size}#{fmt_type}", value }
+      puts
+    end
+  end
+rescue Exception => e
+  STDERR.print "\n", e.message, "\n"
+  print_directions
+end

data/lib/datafarming/freq_sets.rb

@@ -0,0 +1,318 @@
+#!/usr/bin/env ruby -w
+
+module FDE
+  TWO_PI = 2.0 * Math::PI
+
+  DesignSet = Struct.new(:nyq, :freqs)
+
+  DESIGN_SETS = {
+    2 => DesignSet.new(
+      nyq = 13, freqs = [
+        [1, 5],
+        [2, 3],
+        [4, 6]
+      ]
+    ),
+    3 => DesignSet.new(
+      nyq = 27, freqs = [
+        [1, 5, 8],
+        [2, 5, 13],
+        [2, 10, 11],
+        [4, 11, 13],
+        [8, 10, 13]
+      ]
+    ),
+    4 => DesignSet.new(
+      nyq = 46, freqs = [
+        [1, 4, 10, 17],
+        [1, 6, 16, 19],
+        [2, 3, 11, 18],
+        [2, 9, 10, 15],
+        [3, 5, 12, 16],
+        [4, 5, 7, 20],
+        [6, 8, 9, 13],
+        [7, 18, 19, 22],
+        [8, 11, 20, 21],
+        [12, 14, 15, 21],
+        [13, 14, 17, 22]
+      ]
+    ),
+    5 => DesignSet.new(
+      nyq = 69, freqs = [
+        [1, 4, 13, 19, 29],
+        [1, 10, 14, 17, 22],
+        [2, 5, 11, 25, 26],
+        [2, 8, 11, 26, 31],
+        [2, 10, 13, 31, 32],
+        [2, 11, 14, 31, 32],
+        [2, 20, 25, 28, 34],
+        [4, 5, 7, 20, 26],
+        [4, 7, 16, 17, 22],
+        [4, 10, 17, 19, 22],
+        [5, 7, 11, 26, 34],
+        [8, 10, 14, 17, 29],
+        [8, 14, 25, 32, 34],
+        [8, 20, 25, 31, 34],
+        [11, 16, 20, 28, 34],
+        [13, 17, 29, 31, 32]
+      ]
+    ),
+    6 => DesignSet.new(
+      nyq = 103, freqs = [
+        [1, 5, 22, 35, 42, 50],
+        [1, 7, 10, 25, 29, 41],
+        [1, 11, 18, 24, 27, 32],
+        [1, 11, 28, 31, 35, 49],
+        [1, 14, 21, 25, 30, 33],
+        [1, 28, 33, 40, 44, 50],
+        [2, 22, 36, 39, 48, 49],
+        [3, 4, 13, 28, 40, 42],
+        [3, 4, 17, 19, 29, 47],
+        [3, 16, 20, 21, 28, 30],
+        [4, 5, 7, 25, 31, 44],
+        [4, 9, 10, 21, 37, 44],
+        [5, 6, 14, 37, 41, 44],
+        [5, 13, 17, 32, 46, 48],
+        [5, 31, 37, 39, 48, 51],
+        [6, 8, 9, 34, 38, 45],
+        [6, 32, 40, 42, 43, 47],
+        [7, 15, 35, 39, 41, 51],
+        [8, 10, 14, 15, 41, 50],
+        [8, 11, 12, 18, 43, 45],
+        [8, 15, 18, 20, 29, 42],
+        [9, 12, 16, 38, 46, 51],
+        [9, 12, 17, 19, 23, 39],
+        [10, 12, 15, 21, 28, 29],
+        [11, 21, 25, 34, 49, 51],
+        [11, 27, 35, 36, 48, 50],
+        [12, 13, 16, 18, 27, 35],
+        [13, 23, 31, 32, 38, 43],
+        [14, 20, 22, 25, 35, 51],
+        [14, 27, 34, 36, 46, 51],
+        [16, 23, 24, 43, 45, 49],
+        [20, 24, 30, 42, 45, 47]
+      ]
+    ),
+    7 => DesignSet.new(
+      nyq = 130, freqs = [
+        [1, 4, 19, 31, 44, 53, 60],
+        [1, 7, 18, 22, 27, 57, 60],
+        [1, 9, 14, 21, 40, 46, 57],
+        [1, 10, 16, 29, 34, 37, 41],
+        [2, 3, 12, 29, 37, 50, 57],
+        [2, 17, 22, 23, 30, 33, 59],
+        [3, 7, 19, 28, 30, 43, 48],
+        [3, 8, 10, 27, 41, 42, 63],
+        [3, 21, 41, 49, 50, 54, 64],
+        [4, 7, 9, 24, 30, 49, 59],
+        [6, 9, 19, 20, 36, 41, 43],
+        [6, 31, 40, 47, 51, 61, 64],
+        [7, 9, 17, 20, 32, 62, 63],
+        [8, 33, 47, 48, 51, 53, 60],
+        [10, 11, 18, 23, 37, 53, 62],
+        [10, 11, 36, 44, 49, 51, 63],
+        [11, 14, 23, 24, 29, 31, 50],
+        [11, 16, 17, 20, 46, 59, 61],
+        [12, 17, 21, 27, 40, 47, 58],
+        [19, 29, 30, 33, 54, 56, 61],
+        [20, 33, 37, 42, 43, 58, 61],
+        [21, 27, 34, 51, 56, 59, 60],
+        [23, 28, 38, 47, 49, 50, 63],
+        [31, 32, 38, 40, 43, 53, 57]
+      ]
+    ),
+    8 => DesignSet.new(
+      nyq = 168, freqs = [
+        [1, 9, 13, 16, 40, 46, 51, 74],
+        [2, 8, 11, 25, 26, 57, 64, 69],
+        [2, 8, 15, 26, 27, 59, 64, 73],
+        [3, 19, 32, 34, 39, 62, 79, 80],
+        [3, 22, 23, 32, 37, 39, 50, 80],
+        [5, 32, 34, 45, 62, 65, 80, 81],
+        [8, 15, 17, 27, 53, 58, 64, 82],
+        [8, 31, 57, 58, 64, 67, 69, 82],
+        [9, 10, 16, 38, 40, 43, 51, 55],
+        [10, 16, 29, 33, 38, 40, 41, 75],
+        [16, 33, 40, 46, 71, 74, 75, 83],
+        [22, 32, 45, 47, 50, 61, 80, 81]
+      ]
+    ),
+    9 => DesignSet.new(
+      nyq = 209, freqs = [
+        [1, 5, 19, 22, 31, 56, 64, 71, 99],
+        [1, 22, 29, 38, 42, 48, 53, 56, 88],
+        [2, 19, 54, 61, 66, 72, 88, 97, 98],
+        [3, 13, 41, 62, 76, 77, 81, 99, 101],
+        [5, 25, 54, 63, 71, 77, 95, 98, 99],
+        [6, 7, 11, 26, 47, 55, 57, 82, 85],
+        [7, 20, 22, 26, 38, 43, 75, 99, 100],
+        [8, 30, 33, 39, 40, 44, 57, 59, 94],
+        [9, 11, 12, 38, 45, 51, 55, 70, 86],
+        [9, 11, 14, 40, 44, 52, 59, 76, 86],
+        [11, 17, 43, 44, 47, 85, 93, 95, 100],
+        [11, 19, 21, 24, 28, 44, 78, 90, 104],
+        [11, 28, 32, 55, 69, 89, 95, 102, 104],
+        [12, 14, 22, 39, 45, 52, 94, 95, 99],
+        [14, 16, 51, 52, 57, 60, 70, 77, 99],
+        [16, 21, 60, 66, 78, 80, 88, 91, 95],
+        [17, 19, 20, 33, 44, 65, 75, 87, 93],
+        [18, 22, 24, 37, 69, 76, 90, 99, 102],
+        [18, 22, 28, 37, 57, 80, 88, 91, 104],
+        [19, 27, 32, 33, 42, 49, 53, 77, 89],
+        [23, 34, 35, 38, 40, 59, 66, 79, 88],
+        [27, 49, 55, 57, 66, 69, 73, 92, 102],
+        [31, 38, 54, 55, 64, 66, 84, 98, 103],
+        [33, 46, 51, 68, 70, 76, 77, 80, 91],
+        [36, 42, 53, 66, 73, 74, 76, 88, 92]
+      ]
+    ),
+    10 => DesignSet.new(
+      nyq = 268, freqs = [
+        [1, 5, 13, 24, 51, 54, 71, 86, 93, 114],
+        [1, 5, 18, 21, 32, 62, 70, 77, 105, 117],
+        [1, 5, 23, 34, 42, 49, 74, 87, 101, 104],
+        [1, 14, 21, 50, 60, 66, 69, 77, 103, 107],
+        [1, 35, 51, 90, 93, 97, 107, 112, 118, 130],
+        [1, 49, 62, 71, 87, 107, 114, 117, 122, 128],
+        [2, 3, 11, 23, 53, 68, 86, 105, 115, 122],
+        [2, 5, 13, 25, 34, 65, 71, 87, 106, 120],
+        [2, 13, 22, 51, 56, 59, 87, 101, 119, 126],
+        [2, 15, 22, 33, 47, 56, 75, 83, 121, 126],
+        [2, 19, 29, 68, 81, 86, 111, 122, 123, 131],
+        [2, 34, 47, 63, 69, 106, 109, 120, 121, 129],
+        [3, 7, 15, 35, 44, 46, 69, 102, 121, 126],
+        [3, 7, 53, 55, 74, 82, 83, 98, 116, 121],
+        [3, 10, 11, 15, 39, 55, 72, 74, 106, 115],
+        [3, 15, 37, 47, 54, 58, 63, 82, 83, 96],
+        [3, 37, 41, 42, 53, 61, 63, 70, 88, 118],
+        [4, 9, 58, 79, 85, 86, 109, 111, 123, 126],
+        [4, 11, 23, 25, 49, 55, 58, 86, 125, 126],
+        [5, 11, 25, 42, 49, 82, 90, 105, 108, 117],
+        [5, 16, 23, 25, 31, 58, 98, 101, 102, 115],
+        [6, 7, 15, 39, 50, 55, 73, 75, 92, 102],
+        [6, 7, 35, 39, 66, 89, 91, 100, 110, 115],
+        [6, 9, 10, 33, 47, 64, 69, 77, 98, 109],
+        [6, 10, 25, 57, 64, 65, 87, 98, 101, 125],
+        [6, 19, 43, 45, 66, 95, 99, 100, 110, 127],
+        [6, 50, 59, 61, 79, 92, 95, 102, 119, 127],
+        [7, 18, 26, 30, 35, 73, 75, 76, 97, 107],
+        [7, 20, 22, 23, 55, 89, 94, 106, 115, 125],
+        [8, 17, 18, 31, 38, 85, 89, 91, 113, 118],
+        [8, 18, 21, 38, 43, 45, 49, 103, 117, 118]
+      ]
+    ),
+    11 => DesignSet.new(
+      nyq = 335, freqs = [
+        [1, 4, 14, 26, 33, 50, 71, 94, 110, 152, 163]
+      ]
+    ),
+    12 => DesignSet.new(
+      nyq = 440, freqs = [
+        [1, 4, 14, 21, 72, 84, 113, 137, 146, 173, 192, 203]
+      ]
+    ),
+    13 => DesignSet.new(
+      nyq = 579, freqs = [
+        [1, 4, 10, 17, 29, 52, 95, 117, 132, 172, 196, 243, 269]
+      ]
+    ),
+    14 => DesignSet.new(
+      nyq = 684, freqs = [
+        [1, 4, 10, 17, 29, 52, 83, 107, 150, 172, 221, 236, 268, 306]
+      ]
+    ),
+    15 => DesignSet.new(
+      nyq = 811, freqs = [
+        [1, 4, 10, 17, 29, 52, 78, 93, 143, 165, 202, 245, 277, 343, 383]
+      ]
+    ),
+    16 => DesignSet.new(
+      nyq = 811, freqs = [
+        [1, 4, 10, 17, 29, 52, 88, 131, 155, 196, 246, 312, 359, 408, 423, 445]
+      ]
+    ),
+    17 => DesignSet.new(
+      nyq = 1089, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 116, 181, 224, 264, 318, 359, 403, 450, 481, 526]
+      ]
+    ),
+    18 => DesignSet.new(
+      nyq = 1283, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 165, 210, 265, 335, 400, 436, 480, 523, 563, 587]
+      ]
+    ),
+    19 => DesignSet.new(
+      nyq = 1458, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 176, 202, 299, 339, 413, 437, 468, 513, 595, 659]
+      ]
+    ),
+    20 => DesignSet.new(
+      nyq = 1673, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 291, 335, 411, 437, 511, 535, 594, 648, 688, 749]
+      ]
+    ),
+    21 => DesignSet.new(
+      nyq = 1871, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 218, 242, 312, 343, 425, 498, 598, 700, 755, 795, 882]
+      ]
+    ),
+    22 => DesignSet.new(
+      nyq = 2070, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 335, 399, 444, 613, 691, 715, 841, 902, 942, 989]
+      ]
+    ),
+    23 => DesignSet.new(
+      nyq = 2306, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 314, 438, 540, 564, 684, 758, 784, 924, 971, 1035, 1096]
+      ]
+    ),
+    24 => DesignSet.new(
+      nyq = 2625, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 303, 350, 470, 529, 659, 737, 903, 977, 1017, 1103, 1127, 1285]
+      ]
+    ),
+    25 => DesignSet.new(
+      nyq = 2886, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 303, 350, 405, 505, 618, 642, 728, 885, 1005, 1036, 1194, 1230, 1402]
+      ]
+    ),
+    26 => DesignSet.new(
+      nyq = 3227, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 303, 350, 405, 505, 529, 634, 748, 824, 1071, 1130, 1238, 1321, 1465, 1529]
+      ]
+    ),
+    27 => DesignSet.new(
+      nyq = 3549, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 303, 350, 405, 505, 529, 702, 941, 1092, 1168, 1251, 1277, 1359, 1542, 1582, 1722]
+      ]
+    ),
+    28 => DesignSet.new(
+      nyq = 3931, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 303, 350, 405, 505, 529, 588, 757, 1046, 1072, 1270, 1377, 1490, 1535, 1702, 1738, 1915]
+      ]
+    ),
+    29 => DesignSet.new(
+      nyq = 4298, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 303, 350, 405, 505, 529, 588, 680, 903, 1016, 1283, 1328, 1406, 1645, 1793, 1869, 1951, 2058]
+      ]
+    ),
+    30 => DesignSet.new(
+      nyq = 4692, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 303, 350, 405, 505, 529, 588, 680, 903, 1061, 1131, 1171, 1358, 1502, 1721, 1898, 1934, 2012, 2106]
+      ]
+    ),
+    50 => DesignSet.new(
+      nyq = 17_853, freqs = [
+        [1, 4, 10, 17, 29, 52, 67, 89, 132, 164, 205, 259, 303, 350, 405, 505, 529, 588, 680, 903, 1016, 1061, 1248, 1358, 1445, 1838, 1878, 2086, 2117, 2195, 2613, 2840, 3060, 3314, 3422, 3816, 4097, 4247, 4565, 4705, 5326, 5511, 6620, 6857, 6983, 7235, 7964, 8058, 8287, 8655]
+      ]
+    )
+  }.freeze
+end
+
+
+if __FILE__ == $PROGRAM_NAME
+  require 'colorize'
+  String.disable_colorization false
+  require_relative 'error_handling'
+  ErrorHandling.message ['This module is not intended to be executable.'.yellow]
+end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: datafarming
 version: !ruby/object:Gem::Version
-  version: 1.0.0
+  version: 1.1.0
 platform: ruby
 authors:
 - Paul J Sanchez
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2018-06-25 00:00:00.000000000 Z
+date: 2018-07-15 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: fwt
@@ -66,12 +66,15 @@ executables:
 - mser.rb
 - pool_files.rb
 - rundesign_general.rb
+- scaled_fde.rb
+- scaled_rf_cubed.rb
 - stack_nolhs.rb
 - stripheaderdups.rb
 - stripheaders.rb
 extensions: []
 extra_rdoc_files: []
 files:
+- LICENSE.md
 - README.md
 - datafarming.gemspec
 - exe/augment_design.rb
@@ -84,11 +87,14 @@ files:
 - exe/mser.rb
 - exe/pool_files.rb
 - exe/rundesign_general.rb
+- exe/scaled_fde.rb
+- exe/scaled_rf_cubed.rb
 - exe/stack_nolhs.rb
 - exe/stripheaderdups.rb
 - exe/stripheaders.rb
 - lib/datafarming/cross.rb
 - lib/datafarming/error_handling.rb
+- lib/datafarming/freq_sets.rb
 - lib/datafarming/nolh_designs.rb
 homepage: https://gitlab.nps.edu/pjsanche/datafarmingrubyscripts.git
 licenses: