pulo 0.1.1 → 0.1.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: c66281d85496053c2775de970706f5cdaae460d8
-  data.tar.gz: 0d6b04db8547de1a3c72d6cac992a59b6d09a2f5
+  metadata.gz: c3cd821eca19f9a4d1a19906109e39c39077b3d3
+  data.tar.gz: 7287a0e64e3f40f0f90b074a20bb30a4730ec0a3
 SHA512:
-  metadata.gz: c3ac909efb7efe59174d02a21e6a20944c6ce241f8f71627da31aca8eb3f68159c9eb5ad5fa4ae51bbdaafe1f1d49bd378b62a30a3d476b7f4587f9202dd2e28
-  data.tar.gz: 13b06f8726a6e088abb87fc7ae37ca05d114f92becb880adf2bed40da9219214b056ee9da3cf562b5b8c9c7dc44ae1c90933c1d010702c2d00279f4060f2c34a
+  metadata.gz: 6e9bfeff6d877d31741b238b137a86830e2b29568d86db2fcf07532b0c73007edbbd0d4063fd8535de318f0cb34a7d93093102e8a3caf94fb5a170cc6df242f3
+  data.tar.gz: a01f39e0ebd30eb5e60f887338e4601818415cdb528ef5f3f4f7e90513de0aec06f4d59d76254c694b00267f0df31799adc555f3258e71c844f8846decc9876b

data/lib/pulo.rb

@@ -0,0 +1,28 @@
+# encoding: utf-8
+
+require_relative 'pulo/version'
+
+require 'bigdecimal'
+require 'bigdecimal/util'
+require 'yaml'
+
+require_relative 'pulo/helpers'
+require_relative 'pulo/formatting'
+require_relative 'pulo/quantity_checker'
+require_relative 'pulo/exceptions'
+
+require_relative 'pulo/quantity/quantity'
+require_relative 'pulo/quantity/quantity_builder'
+require_relative 'pulo/quantity/unit'
+require_relative 'pulo/quantity/dimension'
+
+require_relative 'pulo/quantity/quantity_definitions'
+
+require_relative 'pulo/quantity/quantity_groups/quantity_groups'
+
+require_relative 'pulo/figure/figure2d'
+require_relative 'pulo/figure/figure3d'
+require_relative 'pulo/machine/machines'
+require_relative 'pulo/quantity/numeric_overloads'
+require_relative 'pulo/frames/frame'
+require_relative 'pulo/tables/tables'

data/lib/pulo/exceptions.rb

@@ -0,0 +1,6 @@
+# encoding: utf-8
+module Pulo
+  class QuantitiesException < StandardError
+
+  end
+end

data/lib/pulo/figure/figure2d.rb

@@ -0,0 +1,248 @@
+
+module Pulo
+
+  module Figure2D
+    include Quantity_Checker
+    attr_reader :area,:perimeter
+    def * (other)
+      if other.is_a?(Length)
+        raise "No extrusion figure known for #{self.class}" unless self.respond_to?(:extrusion_figure)
+        self.extrusion_figure.new(face: self, length: other)
+      else
+        raise "Dont know how to multiply a #{self.class} by a #{other.class}"
+      end
+    end
+  end
+
+  class Circle
+    include Figure2D
+    attr_reader :radius,:diameter
+    def extrusion_figure
+      Cylinder
+    end
+    def initialize(radius: nil, diameter: nil, area: nil)
+      raise 'Circle needs area or diameter or radius.' unless (area || diameter || radius)
+      quantity_check [radius,Length] ,[diameter,Length] , [area,Area]
+      if area
+        @area=area
+        @radius=(@area/Math::PI).rt(2)
+        @diameter=@radius*2
+      else
+        if diameter
+          @diameter=diameter
+          @radius=@diameter/2
+        else
+          @radius=radius
+          @diameter=@radius*2
+        end
+        @area=(@radius**2)*Math::PI
+      end
+      @perimeter=@radius*2*Math::PI
+    end
+  end
+
+  class Square
+    include Figure2D
+
+    attr_reader :width
+    def initialize(width: nil, area: nil)
+      raise "Square needs area or width." unless (area || width)
+
+      quantity_check [area,Area] ,[width,Length]
+      if area
+        @area=area; @width=@area.rt(2)
+      else
+        @width=width; @area=@width**2
+      end
+      @perimeter=@width*4
+    end
+  end
+
+  class Rectangle
+    include Figure2D
+
+    attr_reader :width, :height
+
+    def extrusion_figure
+      Cuboid
+    end
+    def initialize(area: nil, width: nil, height: nil)
+      quantity_check [area,Area] ,[width,Length] , [height,Length]
+      raise 'Rectangle needs width and height or area and width or height.' unless (width && height) || (area && (width || height))
+      if area
+        @area=area
+        if width
+          @width=width; @height=@area/@width
+        else
+          @height=height; @width=@area/@height
+        end
+      else
+        @width=width; @height=height
+        @area=@width*@height
+      end
+      @perimeter=@width*2+@height*2
+    end
+  end
+
+  class Triangle
+    include Figure2D
+    attr_reader :angles,:lengths
+
+    def extrusion_figure
+      Prism
+    end
+    def initialize(angles: [], lengths: [], area: nil)
+      quantity_check [area,Area] ,[angles,Array] , [lengths,Array]
+
+      raise 'Need more arguments for triangle' unless (
+            ((not area) && angles.count+lengths.count>=3 && lengths.count>0) ||
+            (area && (angles.count>=2 || lengths.count>=2 || (lengths.count==1 && angles.count>0)))
+      )
+      while lengths.count<3; lengths<<nil; end
+      while angles.count<3; angles<<nil; end
+
+      lengths_count=lengths.compact.count
+      angles_count=angles.compact.count
+      if area
+        #area=area.value if area.is_a?(Area)
+        if angles_count==2
+          while angles[0].nil?
+            angles=angles.unshift(angles.last).take(3)
+            lengths=lengths.unshift(lengths.last).take(3)
+          end
+          angles[2]=Angle.pi-angles[0]-angles[1]
+          k=((area*2)/(Math.sin(angles[0])*Math.sin(angles[1])*Math.sin(angles[2]))).rt(2)
+          lengths[0]=k*Math.sin(angles[0])
+          lengths[1]=k*Math.sin(angles[1])
+          lengths[2]=k*Math.sin(angles[2])
+        else
+          if lengths_count==2
+            while lengths[0].nil? or lengths[1].nil?
+              lengths=lengths.unshift(lengths.last).take(3)
+              angles=angles.unshift(angles.last).take(3)
+            end
+            angles[0]=angle_from_area(area,lengths[0],lengths[1])
+            lengths[2]=length_from_cosine(lengths[0],lengths[1],angles[0])
+            angles[1]=angle_from_sine(lengths[0],lengths[2],angles[0])
+            angles[2]=Angle.pi-angles[0]-angles[1]
+          else
+            while lengths[0].nil?
+              lengths=lengths.unshift(lengths.last).take(3)
+              angles=angles.unshift(angles.last).take(3)
+            end
+            if not angles[0].nil?
+              lengths[1]=length_from_area(area,lengths[0],angles[0])
+              lengths[2]=length_from_cosine(lengths[0],lengths[1],angles[0])
+              angles[1]=angle_from_sine(lengths[0],lengths[2],angles[0])
+              angles[2]=Angle.pi-angles[0]-angles[1]
+            elsif not angles[1].nil?
+              #assume isosceles
+              angles[0]=angles[2]=(Angle.pi-angles[1])/2
+              lengths[1]=length_from_area(area,lengths[0],angles[0])
+              lengths[2]=length_from_cosine(lengths[0],lengths[1],angles[0])
+            else
+              lengths[2]=length_from_area(area,lengths[0],angles[2])
+              lengths[1]=length_from_cosine(lengths[0],lengths[2],angles[2])
+              angles[0]=angle_from_sine(lengths[2],lengths[1],angles[2])
+              angles[1]=Angle.pi-angles[0]-angles[2]
+            end
+          end
+        end
+
+        @area=area
+      else
+        if lengths_count==3
+          angles[0]=angle_from_cosine(lengths[2],lengths[0],lengths[1])
+          angles[1]=angle_from_sine(lengths[0],lengths[2],angles[0])
+          angles[2]=Angle.pi-angles[0]-angles[1]
+        elsif lengths_count==2
+          while lengths[0].nil? or lengths[1].nil?
+            lengths=lengths.unshift(lengths.last).take(3)
+            angles=angles.unshift(angles.last).take(3)
+          end
+          if not angles[0].nil?
+            lengths[2]=length_from_cosine(lengths[0],lengths[1],angles[0])
+            angles[1]=angle_from_sine(lengths[0],lengths[2],angles[0])
+            angles[2]=Angle.pi-angles[0]-angles[1]
+          elsif not angles[1].nil?
+            angles[2]=angle_from_sine(lengths[1],lengths[0],angles[1])
+            angles[0]=Angle.pi-angles[1]-angles[2]
+            lengths[2]=length_from_sine(lengths[0],angles[0],angles[1])
+          else
+            angles[1]=angle_from_sine(lengths[0],lengths[1],angles[2])
+            angles[0]=Angle.pi-angles[1]-angles[2]
+            lengths[2]=length_from_sine(lengths[0],angles[0],angles[1])
+          end
+        else
+          while angles[0].nil? or angles[1].nil?
+            lengths=lengths.unshift(lengths.last).take(3)
+            angles=angles.unshift(angles.last).take(3)
+          end
+          angles[2]=Angle.pi-angles[0]-angles[1]
+          if not lengths[0].nil?
+            lengths[1]=length_from_sine(lengths[0],angles[2],angles[1])
+            lengths[2]=length_from_sine(lengths[0],angles[0],angles[1])
+          elsif not lengths[1].nil?
+            lengths[0]=length_from_sine(lengths[1],angles[1],angles[2])
+            lengths[2]=length_from_sine(lengths[0],angles[0],angles[1])
+          elsif not lengths[2].nil?
+            lengths[0]=length_from_sine(lengths[2],angles[1],angles[0])
+            lengths[1]=length_from_sine(lengths[0],angles[2],angles[1])
+          end
+        end
+        @area=0.5*lengths[2]*lengths[0]*Math.sin(angles[2])
+      end
+      @angles=angles
+      @lengths=lengths
+      @perimeter=@lengths[0]+@lengths[1]+@lengths[2]
+    end
+    def length_from_cosine(b,c,ang_a); Math.sqrt(b**2+c**2-(2*b*c*Math.cos(ang_a))); end
+    def angle_from_cosine(a,b,c); Math.acos((b**2+c**2-a**2)/(2*b*c)); end
+    def length_from_sine(b,ang_a,ang_b); b/Math.sin(ang_b)*Math.sin(ang_a); end
+    def angle_from_sine(a,b,ang_b); Math.asin(a*Math.sin(ang_b)/b); end
+    def length_from_area(area,b,ang_c); (2*area)/(b*Math.sin(ang_c)); end
+    def angle_from_area(area,a,b); Math.asin(2*area/(a*b)); end
+  end
+
+  class Trapezoid #isoscelese only
+    include Figure2D
+    attr_reader :base,:top,:height,:angle, :side_triangle, :side
+    def extrusion_figure
+      TrapezoidalPrism
+    end
+    def initialize(angle: nil,height: nil,base: nil, area: nil,top: nil)
+      quantity_check [angle,Angle] ,[height,Length] , [base,Length] , [top,Length], [area,Area]
+      raise 'Need more arguments for trapezoid' unless ((angle and height and (base or top) or (base and top and height) or (area and base and top)))
+      if angle and height and (base or top)
+        @angle=angle; @height=height
+        @side_triangle=Triangle.new(angles: [@angle,nil,Angle.pi/2], lengths: [@height,nil,nil])
+        @side=@side_triangle.lengths[2]
+        if base
+          @base=base; @top=@base-@side_triangle.lengths[0]*2
+        else #top
+          @top=top; @base=@top+@side_triangle.lengths[0]*2
+        end
+        @area=(@base+@top)/2*@height
+      elsif base and top and height
+        @height=height;@base=base;@top=top
+        @side_triangle=Triangle.new(lengths: [@height,(@base-@top)/2,nil],angles: [Angle.pi/2])
+        @side=@side_triangle.lengths[2]
+        @area=(@base+@top)/2*@height
+      elsif area and base and top
+        @area=area; @base=base; @top=top
+        @height=(@area*2)/(@base+@top)
+        @side_triangle=Triangle.new(angles: [Angle.pi/2], lengths: [@height,(@base-@top)/2])
+        @side=@side_triangle.lengths[2]
+      elsif area and angle and base
+        #TODO: area and angle and base
+
+        #@angle=angle; @area=area; @base=base
+        #@height=(Math.sqrt((@base**2)+(4*Math.tan(@angle)*@area))-@base)/(2*Math.tan(@angle))
+        #@side_triangle=Triangle.new(angles: [@angle,nil,Angle.pi/2], lengths: [@height,nil,nil])
+        #@side=@side_triangle.lengths[2]
+        #@top=@base+@side_triangle.lengths[0]*2
+      end
+      @perimeter=@base+@top+(@side*2)
+    end
+  end
+end

data/lib/pulo/figure/figure3d.rb

@@ -0,0 +1,166 @@
+
+module Pulo
+
+  module Figure3D
+    include Quantity_Checker
+    attr_reader :volume, :surfacearea
+  end
+
+  class Sphere
+    include Figure3D
+    attr_reader :radius,:diameter
+    def initialize(radius: nil, diameter: nil, volume: nil)
+      quantity_check [radius,Length] ,[diameter,Length], [volume, Volume]
+      raise "Need more arguments for Sphere" unless volume || diameter || radius
+      if volume
+        @volume=volume
+        @radius=(@volume*3/(Angle.pi*4)).rt(3)
+        @diameter=@radius*2
+      else
+        if diameter
+          @diameter=diameter
+          @radius=@diameter/2
+        else
+          @radius=radius
+          @diameter=@radius*2
+        end
+        @volume=(@radius**3)*(4.0/3.0)*Angle.pi
+      end
+      @surfacearea=@radius**2*(4*Angle.pi)
+    end
+  end
+
+  class Cube
+    include Figure3D
+    attr_reader :length, :face
+
+    def initialize(length: nil, face: nil, volume: nil)
+      quantity_check [length,Length] ,[face,Square], [volume, Volume]
+      raise 'Cube needs volume or face or length.' unless (length || volume || face)
+
+      if volume
+        @volume=volume
+        @length=volume.rt(3)
+        @face=Square.new(width: @length)
+      elsif face
+        @face=face
+        @length=@face.width
+        @volume=@face.area*@length
+      elsif length
+        @length=length
+        @face=Square.new(width: @length)
+        @volume=@face.area*@length
+      end
+      @surfacearea=@face.area*6
+    end
+  end
+
+  class Cuboid
+    include Figure3D
+    attr_reader :width, :length, :height, :faces
+    def initialize(width: nil, height: nil, length: nil, volume: nil, face: nil)
+      raise 'Cuboid needs width, length and height or volume and face or length and face.' unless ((width && length && height) || (volume && face) || (length && face))
+
+      @faces=[]
+      if width and length and height
+        @width=width; @length=length; @height=height
+        @faces[0]=Rectangle.new(width: width, height: height)
+        @faces[1]=Rectangle.new(width: width, height: length)
+        @faces[2]=Rectangle.new(width: length, height: height)
+        @volume=@faces[0].area*@length
+      end
+      if volume and face
+        @volume=volume
+        @faces[0]=face
+        @width=@faces[0].width
+        @height=@faces[0].height
+        @length=@volume/@faces[0].area
+        @faces[1]=Rectangle.new(width: @width, height: @length)
+        @faces[2]=Rectangle.new(width: @length, height: @height)
+      end
+      if face and length
+        @faces[0]=face
+        @length=length
+        @volume=face.area * length
+        @width=@faces[0].width
+        @height=@faces[0].height
+        @faces[1]=Rectangle.new(width: @width, height: @length)
+        @faces[2]=Rectangle.new(width: @length, height: @height)
+      end
+      @surfacearea=(@faces[0].area+@faces[1].area+@faces[2].area)*2
+    end
+  end
+
+  class Cylinder
+    include Figure3D
+
+    attr_reader :face, :length
+    def initialize (face: nil, length: nil, volume: nil, radius: nil, diameter: nil)
+      quantity_check [face,Circle] ,[length,Length] , [volume,Volume], [radius, Length], [diameter, Length]
+      raise 'Cylinder needs length and face or volume and length or volume and radius or diameter.' unless
+          (face && length) || (volume && length) || (volume && (radius || diameter)) || (length && (radius || diameter))
+
+      if face and length
+        @face=face
+        @length=length
+        @volume=@face.area*@length
+      else
+        if volume and length
+          @length=length
+          @volume=volume
+          @face=Circle.new(area: @volume/@length)
+        else
+          if (radius or diameter) and volume
+            @face=Circle.new(radius: radius, diameter: diameter)
+            @volume=volume
+            @length=volume/@face.area
+          else
+            @face=Circle.new(radius: radius, diameter: diameter)
+            @length=length
+            @volume=@face.area*@length
+          end
+        end
+      end
+
+      @surfacearea=@face.area*2+@face.perimeter*@length
+      @area=@face.area
+      @radius=@face.radius
+    end
+  end
+
+  class Prism
+    include Figure3D
+    attr_reader :face, :length
+    attr_reader :side_areas
+    def initialize(face: nil, length: nil)
+      quantity_check [face,Triangle] ,[length,Length]
+      raise 'Prism needs a face and length.' unless face && length
+      @face=face
+      @length=length
+      @volume=@face.area*@length
+      @side_areas=[]
+      @side_areas[0]=@face.lengths[0]*@length
+      @side_areas[1]=@face.lengths[1]*@length
+      @side_areas[2]=@face.lengths[2]*@length
+      @surfacearea=@face.area*2+@side_areas[0]+@side_areas[1]+@side_areas[2]
+    end
+  end
+
+  class TrapezoidalPrism
+    include Figure3D
+    attr_reader :face, :length
+    attr_reader :base_area,:side_area,:top_area
+    def initialize(face: nil, length: nil)
+      quantity_check [face,Trapezoid] ,[length,Length]
+      raise 'TrapezoidalPrism needs a face and length.' unless face && length
+      @face=face
+      @length=length
+      @volume=@face.area*@length
+      @base_area=@face.base*@length
+      @side_area=@face.side*@length
+      @top_area=@face.top*@length
+      @surfacearea=@face.area*2+@base_area+@top_area+@side_area*2
+    end
+  end
+
+end