magick 0.1.1 → 0.1.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 4826f852b57d6f69bbf1d1d1c4671174e3b3d68dae520b5fd09f43a7c2ba2bac
-  data.tar.gz: 0d56f6b3d4c9a6ec5fe45759695375700aacdedc28dcd9b1dccc23c06e9bff37
+  metadata.gz: 5ae8b9cd9f9526ea3adac12fd63ab30cd38895c035493a3f4866c83d25df015b
+  data.tar.gz: d60e88dbedd1c8d9c7bedf778beeefed2c3012678ce5022a741c02471ba76abd
 SHA512:
-  metadata.gz: ebe917fea532e9cd78d94e247120278f9f155cb857ff0daaea8f9538342a33fefff6e8a3c813180763845fd822efababcfaa0d70a0101d5b665becded5815a46
-  data.tar.gz: 88573678ce6d983164280c784f14dd285f1136c0354cc1c40a3c92df1538637ea909b5f9ad4d1098a32eb067ec34a2b6a08a5058839d4d3a17aa0080115a09fb
+  metadata.gz: 9f033baf4e66b21b0c542d7f99ab95a717fb864acb48e5f51ba06e3c7eec9cd8e68f3d659272942fe4ffb487c618bb563e2b6c6765f4e016e196826cfd6cc45b
+  data.tar.gz: c967a3d6c062ee01b0f804cae609d8f281d4a33f66ba7449edaf2a08f8c505f52607090f769811c1db205745bf458bebdcd506cacfadc2f29b6295b3acf74bf9

data/lib/magick/loops.rb

@@ -31,6 +31,7 @@ module Magick
         }
       end
 
+      # sum and length both assume they will be given a list
       def sum = fold.(Magick::Maths::add).(0)
 
       def length = fold.(->(acc) { ->(elem) { acc + 1 }}).(0)

data/lib/magick/maths.rb

@@ -11,7 +11,7 @@ module Magick
       def subtract = ->(x){ ->(y) { add.(x).(negate.(y)) }}
 
       def multiply
-        power.(add)
+        power_2.(add)
       end
       alias :mult :multiply
 
@@ -48,7 +48,7 @@ module Magick
       end
 
       def exp
-        power.(multiply)
+        power_2.(multiply)
       end
 
       def power
@@ -73,6 +73,60 @@ module Magick
         }
       end
 
+      # power, as written, will have a stack trace error if we try to 
+      # do something like power.(mult).(2).(32)
+      # If we eliminate the recursion, we can avoid this, but it means we need
+      # to use some regular Ruby constructs (eg, `while`, instead of the Y combinator)
+      # The following is essentially how Stepanov & McJones define `power` in 
+      # Elements of Programming
+      def power_accumulate_positive
+        ->(op){
+          ->(r){
+            ->(a){
+              ->(n){
+                while (true)
+                  if n % 2 != 0
+                    r = op.(r).(a)
+                    return r if n == 1
+                  end
+                  a = op.(a).(a)
+                  n = n/2
+                end
+              }
+            }
+          }
+        }
+      end
+
+      def power_accumulate
+        ->(op){
+          ->(r){
+            ->(n){
+              ->(a){
+                return r if n == 0
+                power_accumulate_positive.(op).(r).(n).(a)
+              }
+            }
+          }
+        }
+      end
+
+      def power_2
+        ->(op){
+          ->(b){
+            ->(e){
+              while (e % 2 == 0)
+                b = op.(b).(b)
+                e = e/2
+              end
+              e = e / 2
+              return b if e == 0
+              power_accumulate.(op).(b).(op.(b).(b)).(e)
+            }
+          }
+        }
+      end
+
       def average
         Smullyan::Birds::Phi.
           (Magick::Maths::div_to_f).
@@ -80,19 +134,24 @@ module Magick
           (Magick::Loops::length)
       end
 
-      def powers_of_two = power.(multiply).(2)
+      def powers_of_two = exp.(2)
 
       # x and y are pairs of ints
       def fibonacci_matrix_multiply = ->(x){
         ->(y){
-          [x[0] * (y[1] + y[0]) + x[1] * y[0], x[0] * y[0] + x[1] * y[1]]
+          [
+            x[0] * (y[1] + y[0]) +
+              x[1] * y[0], 
+            x[0] * y[0] +
+              x[1] * y[1]
+          ]
         }
       }
 
       def fibonacci
         ->(n){
           n == 0 ? 0 :
-          power.(fibonacci_matrix_multiply).([1,0]).(n)[0]
+          power_2.(fibonacci_matrix_multiply).([1,0]).(n)[0]
         }
       end
     end

data/lib/magick/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module Magick
-  VERSION = "0.1.1"
+  VERSION = "0.1.2"
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: magick
 version: !ruby/object:Gem::Version
-  version: 0.1.1
+  version: 0.1.2
 platform: ruby
 authors:
 - Phil Crissman