datafarming 1.4.0 → 2.2.0

data/lib/datafarming/res_v_seqs.rb

@@ -0,0 +1,57 @@
+require 'fwt'
+
+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,
+           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]
+
+  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
+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.2.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.4.0
+  version: 2.2.0
 platform: ruby
 authors:
 - Paul J Sanchez
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2021-01-19 00:00:00.000000000 Z
+date: 2024-12-11 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,11 +52,38 @@ 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
@@ -64,26 +91,21 @@ executables:
 - 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: []
 extra_rdoc_files: []
 files:
-- DESCRIPTIONS.html
-- DESCRIPTIONS.md
 - LICENSE.md
 - README.md
 - datafarming.gemspec
-- exe/augment_design.rb
 - exe/batchrunner.rb
 - exe/blank2csv.rb
 - exe/cat.rb
@@ -91,23 +113,25 @@ files:
 - 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
+- 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
@@ -120,14 +144,14 @@ required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
-      version: '2.5'
+      version: '2.6'
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
   - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.2.3
+rubygems_version: 3.5.23
 signing_key:
 specification_version: 4
 summary: Useful scripts for data farming.