datafarming 1.3.0 → 2.1.0

data/lib/datafarming/{factorial_generator.rb → res_v_seqs.rb}

@@ -1,23 +1,23 @@
 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) 2018 Paul J Sanchez
-# License:: MIT
-#
-def make_design(number_of_factors)
-  index = [1, 2, 4, 8, 15, 16, 32, 51, 64, 85, 106, 128, 150, 171,
+module RESV
+  # 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) 2023 Paul J Sanchez
+  # License:: MIT
+  #
+  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,
@@ -28,9 +28,8 @@ def make_design(number_of_factors)
            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,
+           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,
@@ -42,14 +41,17 @@ def make_design(number_of_factors)
            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, 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
+  def RESV.make_design(number_of_factors)
+    Array.new(number_of_factors) do |i|
+      Array.new(POWER[number_of_factors], 0).tap { |a| a[INDEX[i]] = 1 }.hadamard
+    end.transpose
+  end
+
+  def RESV.star_pts(k)
+    Array.new(2*k+1) do |i|
+      Array.new(k) { |j| (j+1) == (i+1) / 2 ? (i.odd? ? -1 : 1) : 0 }
+    end
   end
-  design.transpose
 end

data/lib/datafarming/scale_design.rb

@@ -0,0 +1,29 @@
+#!/usr/bin/env ruby -w
+
+# A Scaler object will rescale a design from standard range [-1, 1] to a
+# range specified by min & max, with the desired number of decimals
+class Scaler
+  def initialize(min: -1, max: 1, decimals: nil)
+    unless decimals.nil? || ((decimals >= 0) && (decimals == decimals.to_i))
+      fail "Invalid decimals argument: #{decimals}. Must be a non-negative integer"
+    end
+    @range = (max - min).to_r / 2
+    @mid = min + @range
+    @scale_factor = (decimals.nil? || decimals.zero?) ? decimals : 10r**decimals
+  end
+
+  def scale(value)
+    new_value = @mid + @range * value
+    if @scale_factor.nil?
+      if new_value.to_r.denominator == 1
+        new_value.to_i
+      else
+        new_value.to_f
+      end
+    elsif @scale_factor.zero?
+      new_value.round
+    else
+      ((@scale_factor * new_value).round / @scale_factor).to_f
+    end
+  end
+end

data/lib/datafarming/screen_freq_sets.rb

@@ -0,0 +1,143 @@
+# frozen_string_literal: true
+
+module FBS
+  TWO_PI = 2.0 * Math::PI
+
+  DesignSet = Struct.new(:nyq, :freqs, keyword_init: true)
+
+  DESIGN_SETS = {
+    2 => DesignSet.new(
+      nyq: 11, freqs: [
+        [1, 3],
+        [1, 4],
+        [1, 5]
+      ]
+    ),
+    3 => DesignSet.new(
+      nyq: 19, freqs: [
+        [1, 3, 7],
+        [1, 7, 8],
+        [2, 3, 5],
+        [4, 6, 9]
+      ]
+    ),
+    4 => DesignSet.new(
+      nyq: 29, freqs: [
+        [1, 4, 5, 11]
+      ]
+    ),
+    5 => DesignSet.new(
+      nyq: 43, freqs: [
+        [1, 5, 6, 8, 20],
+        [1, 8, 9, 11, 14],
+        [3, 4, 7, 9, 19]
+      ]
+    ),
+    6 => DesignSet.new(
+      nyq: 59, freqs: [
+        [1, 6, 10, 16, 24, 28],
+        [1, 7, 11, 16, 19, 23],
+        [2, 8, 20, 23, 25, 29]
+      ]
+    ),
+    7 => DesignSet.new(
+      nyq: 79, freqs: [
+        [1, 3, 8, 11, 26, 29, 33],
+        [1, 4, 11, 12, 29, 30, 35],
+        [1, 5, 6, 14, 15, 17, 38],
+        [1, 10, 11, 16, 17, 24, 25],
+        [2, 3, 5, 11, 20, 21, 25],
+        [2, 3, 17, 22, 25, 33, 36],
+        [2, 5, 7, 13, 15, 16, 39],
+        [2, 7, 9, 10, 15, 19, 33],
+        [2, 20, 22, 29, 31, 32, 34],
+        [3, 8, 15, 17, 18, 24, 28]
+      ]
+    ),
+    8 => DesignSet.new(
+      nyq: 101, freqs: [
+        [1, 4, 5, 17, 21, 27, 41, 45],
+        [1, 5, 6, 14, 15, 17, 38, 41],
+        [2, 8, 10, 11, 19, 34, 42, 47],
+        [2, 14, 29, 31, 38, 40, 48, 50],
+        [3, 5, 8, 24, 28, 33, 42, 47],
+        [3, 13, 15, 18, 22, 42, 45, 50],
+        [4, 7, 16, 17, 20, 22, 26, 45],
+        [5, 7, 13, 25, 29, 31, 40, 42],
+        [5, 8, 12, 20, 26, 34, 43, 45]
+      ]
+    ),
+    9 => DesignSet.new(
+      nyq: 127, freqs: [
+        [1, 3, 10, 11, 15, 33, 38, 49, 57],
+        [1, 7, 8, 18, 19, 24, 39, 52, 53],
+        [1, 7, 8, 18, 19, 24, 52, 53, 62],
+        [1, 9, 13, 22, 28, 37, 42, 51, 58],
+        [1, 13, 14, 22, 23, 54, 55, 61, 62],
+        [1, 29, 38, 39, 41, 54, 55, 61, 62],
+        [2, 3, 5, 11, 15, 31, 38, 50, 52],
+        [2, 5, 20, 26, 28, 39, 41, 55, 63],
+        [3, 4, 7, 12, 27, 29, 38, 47, 52],
+        [4, 5, 9, 11, 15, 23, 36, 48, 62],
+        [5, 7, 16, 20, 22, 28, 45, 46, 48],
+        [5, 7, 16, 20, 22, 28, 46, 48, 62],
+        [7, 13, 16, 18, 24, 28, 34, 51, 54],
+        [7, 17, 20, 24, 26, 30, 38, 42, 53]
+      ]
+    ),
+    10 => DesignSet.new(
+      nyq: 149, freqs: [
+        [3, 15, 18, 25, 40, 41, 49, 56, 60, 61]
+      ]
+    ),
+    11 => DesignSet.new(
+      nyq: 191, freqs: [
+        [1, 8, 9, 12, 19, 34, 41, 46, 62, 85, 93],
+        [1, 9, 22, 29, 41, 46, 52, 56, 59, 61, 94],
+        [1, 15, 16, 19, 20, 26, 27, 62, 75, 79, 84],
+        [1, 19, 25, 26, 35, 36, 40, 41, 48, 64, 94],
+        [2, 3, 10, 13, 19, 33, 44, 54, 58, 72, 82],
+        [2, 3, 11, 12, 31, 34, 44, 47, 61, 82, 87],
+        [2, 5, 7, 13, 18, 35, 38, 48, 67, 82, 85],
+        [2, 7, 15, 20, 31, 36, 48, 54, 57, 69, 73],
+        [3, 16, 19, 36, 39, 44, 47, 62, 65, 92, 95],
+        [4, 7, 11, 16, 20, 37, 41, 49, 55, 83, 84],
+        [4, 18, 19, 21, 22, 28, 31, 51, 60, 65, 76],
+        [5, 7, 13, 15, 16, 34, 35, 59, 67, 71, 79],
+        [6, 10, 16, 19, 25, 33, 39, 40, 46, 51, 93],
+        [8, 9, 17, 21, 28, 29, 31, 32, 70, 76, 88],
+        [9, 35, 38, 39, 47, 59, 65, 75, 84, 90, 92],
+        [12, 15, 16, 22, 34, 51, 53, 54, 65, 74, 79],
+        [16, 25, 30, 31, 33, 39, 43, 49, 77, 84, 85]
+      ]
+    ),
+    12 => DesignSet.new(
+      nyq: 223, freqs: [
+        [1, 15, 16, 19, 20, 27, 28, 61, 62, 64, 86, 87],
+        [2, 8, 17, 19, 41, 43, 47, 48, 50, 61, 70, 73],
+        [2, 11, 13, 19, 21, 31, 33, 37, 82, 84, 109, 111],
+        [2, 19, 21, 39, 41, 47, 49, 50, 52, 64, 74, 108]
+      ]
+    ),
+    13 => DesignSet.new(
+      nyq: 263, freqs: [
+        [1, 4, 5, 22, 23, 35, 36, 55, 56, 84, 98, 99, 126]
+      ]
+    ),
+    14 => DesignSet.new(
+      nyq: 331, freqs: [
+        [1, 3, 7, 19, 26, 36, 47, 60, 70, 79, 128, 148, 153, 161]
+      ]
+    ),
+    15 => DesignSet.new(
+      nyq: 389, freqs: [
+        [1, 3, 7, 12, 29, 40, 49, 65, 91, 100, 125, 131, 148, 154, 172]
+      ]
+    ),
+    20 => DesignSet.new(
+      nyq: 853, freqs: [
+        [1, 3, 7, 12, 15, 25, 41, 58, 60, 80, 107, 154, 161, 193, 232, 249, 284, 291, 377, 412]
+      ]
+    )
+  }.freeze
+end

data/lib/datafarming/version.rb

@@ -0,0 +1,5 @@
+# frozen_string_literal: true
+
+module DataFarming
+  VERSION = "2.1.0"
+end

data/lib/datafarming.rb

@@ -0,0 +1,8 @@
+# frozen_string_literal: true
+
+require_relative "datafarming/version"
+
+module DataFarming
+  class Error < StandardError; end
+  # Your code goes here...
+end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: datafarming
 version: !ruby/object:Gem::Version
-  version: 1.3.0
+  version: 2.1.0
 platform: ruby
 authors:
 - Paul J Sanchez
-autorequire: 
+autorequire:
 bindir: exe
 cert_chain: []
-date: 2018-08-17 00:00:00.000000000 Z
+date: 2023-11-22 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: fwt
@@ -30,14 +30,14 @@ dependencies:
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '0.8'
+        version: '1'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - "~>"
       - !ruby/object:Gem::Version
-        version: '0.8'
+        version: '1'
 - !ruby/object:Gem::Dependency
   name: quickstats
   requirement: !ruby/object:Gem::Requirement
@@ -52,26 +52,52 @@ dependencies:
     - - "~>"
       - !ruby/object:Gem::Version
         version: '2'
-description: Ruby scripts for data farming, including pre- and post-processing, design
+- !ruby/object:Gem::Dependency
+  name: optparse
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.4'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.4'
+- !ruby/object:Gem::Dependency
+  name: yaml
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.3'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.3'
+description: Ruby scripts for data farming, including pre- andpost-processing, design
   generation and scaling, and run control.
 email: pjs@alum.mit.edu
 executables:
-- augment_design.rb
 - batchrunner.rb
 - blank2csv.rb
 - cat.rb
 - convert_line_endings.rb
 - cross.rb
 - csv2blank.rb
+- datafarming.rb
+- generate_design.rb
 - moving_average.rb
 - mser.rb
 - mser_nolbm.rb
 - mser_olbm.rb
 - pool_files.rb
 - rundesign_general.rb
-- scaled_fde.rb
-- scaled_rf_cubed.rb
-- stack_nolhs.rb
 - stripheaderdups.rb
 - stripheaders.rb
 extensions: []
@@ -80,35 +106,37 @@ files:
 - LICENSE.md
 - README.md
 - datafarming.gemspec
-- exe/augment_design.rb
 - exe/batchrunner.rb
 - exe/blank2csv.rb
 - exe/cat.rb
 - exe/convert_line_endings.rb
 - exe/cross.rb
 - exe/csv2blank.rb
+- exe/datafarming.rb
+- exe/generate_design.rb
 - exe/moving_average.rb
 - exe/mser.rb
 - exe/mser_nolbm.rb
 - exe/mser_olbm.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.rb
 - lib/datafarming/cross.rb
 - lib/datafarming/error_handling.rb
-- lib/datafarming/factorial_generator.rb
 - lib/datafarming/freq_sets.rb
 - lib/datafarming/moving_average.rb
 - lib/datafarming/nolh_designs.rb
-homepage: https://gitlab.nps.edu/pjsanche/datafarmingrubyscripts.git
+- lib/datafarming/res_v_seqs.rb
+- lib/datafarming/scale_design.rb
+- lib/datafarming/screen_freq_sets.rb
+- lib/datafarming/version.rb
+homepage: https://bitbucket.org/paul_j_sanchez/datafarmingrubyscripts
 licenses:
 - MIT
 metadata: {}
-post_install_message: 
+post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -116,16 +144,15 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
-      version: '2.3'
+      version: '2.6'
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.7.7
-signing_key: 
+rubygems_version: 3.4.20
+signing_key:
 specification_version: 4
 summary: Useful scripts for data farming.
 test_files: []

data/exe/augment_design.rb

@@ -1,66 +0,0 @@
-#!/usr/bin/env ruby -w
-
-# This script generates star points to augment a fractional factorial
-# so you can check for quadratic effects.
-
-require 'colorize'
-
-String.disable_colorization false
-
-require 'optparse'
-require 'datafarming/error_handling'
-
-help_msg = [
-  'Generate star points to augment a fractional factorial ' \
-  'with quadratic effects.',
-  'Results are white-space delimited data written to ' +
-    'stdout'.light_blue + ', and can be redirected', 'as desired.', '',
-  'Syntax:',
-  "\n\t#{ErrorHandling.prog_name} [--help] FACTOR_INFO".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.",
-  '  FACTOR_INFO'.green,
-  "\tEITHER the number of factors (produces standardized +/-1 design),",
-  "\tOR pairs of values " + 'low1 hi1 low2 hi2...lowN hiN'.green +
-    ' for each of the',
-  "\tN factors.", ''
-]
-
-OptionParser.new do |opts|
-  opts.banner = "Usage: #{$PROGRAM_NAME} [-h|--help] [number of factors]"
-  opts.on('-h', '-?', '--help') { ErrorHandling.clean_abort help_msg }
-end.parse!
-
-if ARGV.length == 0 || (ARGV[0] == '?') || (ARGV.length > 1 && ARGV.length.odd?)
-  ErrorHandling.clean_abort help_msg
-else
-  num_factors = 0
-  if ARGV.length == 1
-    num_factors = ARGV.shift.to_i
-    x = Array.new(num_factors, ' 0')
-    ranges = Array.new(2 * num_factors)
-    ranges.each_index { |i| ranges[i] = (i.even? ? '-1' : ' 1') }
-  else
-    num_factors = ARGV.length / 2
-    x = Array.new(num_factors)
-    ranges = ARGV
-    num_factors.times do |i|
-      x[i] = 0.5 * (ranges[2 * i].to_f + ranges[2 * i + 1].to_f)
-    end
-  end
-  # create and print an array at the global center point (0,0,...,0)
-  puts x.join(' ')
-
-  # now generate and print the star points at +/-1 along each factor axis
-  num_factors.times do |i|
-    y = x.clone
-    y[i] = ranges[2 * i]
-    puts y.join(' ')
-    y[i] = ranges[2 * i + 1]
-    puts y.join(' ')
-  end
-end

data/exe/scaled_fde.rb

@@ -1,169 +0,0 @@
-#!/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