linmeric 0.1.0

data/lib/linmeric/CnGal_Matrix_class.rb

@@ -0,0 +1,358 @@
+#! /usr/bin/env ruby
+
+require_relative 'CnGal_new_classes.rb'
+require_relative 'CnGal_tools.rb'
+require_relative 'Listener.rb'
+# Definition of Matrix class
+#
+
+class InputError < ArgumentError; end
+class MyArgError < RuntimeError; end
+
+class Matrix
+  
+  def initialize(rws = 0,cls = 0)
+    listener = Listener.new
+#    @exception = "Matrix "
+    if rws == 0 or cls == 0 or rws.is_a? Float or cls.is_a? Float then 
+      raise MyArgError,"   Argument Error: invalid dimension #{rws}x#{cls} for Matrix object" 
+    elsif !(rws.is_a? Fixnum) or !(cls.is_a? Fixnum) then
+          e = rws unless e.is_a? Fixnum
+          e = cls unless cls.is_a? Fixnum
+          raise MyArgError,"  Argument Error: colums and rows in Integer format expected but #{e.class} found"
+    else
+      @MyRws = rws; @MyCls = cls
+      @mx=[]
+      if block_given? then
+        for i in 0...@MyRws
+          for j in 0...@MyCls
+            value = yield(i,j) if block_given? 
+            @mx << value
+          end
+        end       
+      else
+#        if @exception == "Matrix " then
+          puts "Insert the line values (separated by space) and press return to go on"
+          for i in 0...rws
+            listener.reset
+            cl = listener.gets.split
+            raise ArgumentError, "  #{cls} element(s) expected but #{cl.size} found. Retry" unless cl.size==cls
+            cl.map! { |e| 
+              cv = Tool.convert(e)
+              err = e if cv == nil
+              raise MyArgError, "  Argument Error: Invalid operation '#{err}' found in matrix rows" unless cv != nil
+              e = cv
+            }
+            @mx +=cl 
+#           end  
+#           puts 
+        end
+      end
+    end
+ 
+      
+    rescue ArgumentError => message
+      puts message
+      @mx=[]
+      retry
+      
+  end
+   
+  def Matrix.from_file(filename = "")
+#    @exception = "Matrix "
+    @mx = []
+    @MyCls = 0
+    @MyRws = 0
+    if !(File.exists? filename) or filename == "" then
+      raise MyArgError, "  Argument Error: Invalid directory or filename"
+    else
+      File.open(filename, "r").each_line do |ln|
+        @MyRws += 1
+        ln = ln.chomp.split(',')
+        ln.map! do |el| 
+          cv = Tool.convert(el)
+          err = el if cv == nil
+          raise InputError, "  Argument Error: Invalid operation '#{err}' found in matrix rows" unless cv != nil
+          el = cv
+        end
+        @mx += ln
+        if @mx.size % ln.size != 0 then
+          raise MyArgError, "  Argument Error: Irregular matrix rows"
+#          break
+        end
+        @MyCls = ln.length
+      end
+    end
+#    mat = Matrix.new(1,1) {1}
+#    mat.update(@mx,@MyRws,@MyCls)
+#    mat.inherit_exception(@exception)
+    return mat = Matrix.new(@MyRws,@MyCls) {|i,j| @mx[i*@MyCls+j]}
+    
+#    rescue InputError => error
+#      @exception += "#{error}"
+#      mat = Matrix.new(1,1) {1}
+#      mat.inherit_exception(@exception)
+#      return mat
+  end 
+  
+  def to_file(filename = "")
+    if !(Dir.exist? File.dirname(filename)) or filename == "" then
+      raise MyArgError, "  Argument Error: Invalid directory; Directory not found"
+    else
+      File.open(filename,"w") do |rw|
+        for i in 0...@MyRws do
+          for j in 0...(@MyCls - 1) do
+            rw.print "#{self[i,j]},"
+          end
+          rw.puts "#{self[i,j+1]}"
+        end
+      end
+    end
+  end
+  
+#  def raised_exception?
+#    return nil if @exception == "Matrix "
+#    return @exception
+#  end 
+  
+#  def got_exception()
+#    @exception = "Matrix "
+#  end
+  
+  def show()
+    max = 0
+    @mx.each{ |n| max = n.to_s.size if n.to_s.size > max}
+    for i in 0...@MyRws do
+      print '|'
+      for j in (i*(@MyCls))...(i*(@MyCls)+@MyCls) do
+        print "  #{" "*(max-@mx[j].to_s.size)}#{@mx[j] }"
+      end
+      puts '  |'
+    end
+    return nil
+  end
+  
+  def getCls()
+    return @MyCls
+  end
+  
+  def getRws()
+    return @MyRws
+  end
+  
+  def export()
+    return @mx
+  end
+  
+  def update(mat,rws,cls)
+    if !(mat.is_a? Array) then
+      raise MyArgError, "  Argument Error: invalid matrix array found"
+    else
+      [cls,rws].each do |e|
+        raise MyArgError, "  Argument Error: colums and rows in Integer format expected but #{e.class} found" unless e.is_a? Numeric
+      end
+    end
+#    if @exception == "Matrix " then
+     @mx = mat
+     @MyRws = rws
+     @MyCls = cls
+#    end
+  end
+  
+#  def inherit_exception(exception)
+#     @exception = exception
+#   end
+  
+  def [](x,y)
+    if x.is_a? Numeric and y.is_a? Numeric then
+      return @mx[x*@MyCls + y]
+    elsif x.is_a? Range and y.is_a? Numeric
+      nm = []
+      x.each do |i|
+        nm << @mx[i*@MyCls + y]
+      end
+      return Matrix.new(x.count,1){|i,j| nm[i + j]}
+    elsif x.is_a? Numeric and y.is_a? Range
+      nm = []
+      y.each do |j|
+        nm << @mx[x*@MyCls + j]
+      end
+      return Matrix.new(1,y.count){|i,j| nm[i + j]}
+    elsif x.is_a? Range and y.is_a? Range
+      nm = []
+      x.each do |i|
+        y.each do |j|
+          nm << @mx[i*@MyCls + j]
+        end
+      end
+      return Matrix.new(x.count,y.count){ |i,j| nm[i*y.count + j]}
+    end
+    return nil
+  end
+  
+  def []=(i,j,val)
+    @mx[i*@MyCls + j] = val
+  end
+  
+  def ==(obj)
+   if obj.is_a? Matrix then
+     if (self.getCls == obj.getCls && self.getRws == obj.getRws ) then
+       (@mx == obj.export) ? (return true) : (return false)
+     end
+   end
+   return false
+  end
+  
+  def !=(obj)
+    if obj.is_a? Matrix then
+      self == obj ? (return false) : (return true)
+    else
+      return true
+    end
+  end
+  
+  def similar_to?(obj)
+    if obj.is_a? Matrix then
+      (self.getCls == obj.getCls && 
+        self.getRws == obj.getRws ) ? (return true) : (return false)
+    end
+    return false
+  end
+  
+  def +(m)
+    if m.is_a? Matrix then
+      if m.similar_to? self then
+        return Matrix.new(@MyRws, @MyCls){ |i,j| self[i,j] + m[i,j]}
+      else
+        raise MyArgError, "  Argument Error: invalid matrix for + operation"
+      end
+    else
+      raise MyArgError, "  Argument Error: can't sum #{m.class} to Matrix"
+    end
+  end
+ 
+  def -(m)
+    if m.is_a? Matrix then
+      if m.similar_to? self then
+        return Matrix.new(@MyRws, @MyCls){ |i,j| self[i,j] - m[i,j]}
+      else
+        raise MyArgError, "  Argument Error: invalid matrix for - operation"
+      end
+    else
+      raise MyArgError, "  Argument Error: can't subtact #{m.class} to Matrix"
+    end
+  end
+  
+  def can_multiply?(obj)
+    if obj.is_a? Matrix then
+      (self.getCls == obj.getRws) ? (return true) : (return false)
+    elsif obj.is_a? Numeric
+      return true
+    end
+    return false
+  end
+  
+  def *(ob)
+    case 
+      when (ob.is_a? Numeric)
+        return Matrix.new(@MyRws,@MyCls){ |i,j| self[i,j] * ob}
+      when (ob.is_a? Matrix)
+        if self.can_multiply? ob then
+          temp = []
+          for i in 0...@MyRws
+            for j in 0...ob.getCls()
+              temp[i*@MyCls + j] = 0
+              for k in 0...@MyCls
+                temp[i*@MyCls + j] += @mx[i*@MyCls + k] * ob[k,j]
+              end
+            end
+          end
+          return Matrix.new(@MyRws,ob.getCls()){ |i,j| temp[i * @MyCls + j]}
+        else
+          raise MyArgError, "  Argument Error: Cannot multiply #{@MyRws}x#{@MyCls} Matrix with #{ob.getRws}x#{ob.getCls}"
+        end
+      else
+        raise MyArgError, "  Argument Error: Cannot multiply Matrix with #{ob.class}"
+    end
+  end
+  
+  
+  def can_divide?(obj)
+    return (obj.is_a? Numeric)
+  end
+  
+  
+  def /(obj)
+    return Matrix.new(@MyRws,@MyCls) { |i,j| self[i,j] / obj.to_f} if obj.is_a? Numeric
+    raise MyArgError, "  Argument Error: Invalid division between Matrix and #{obj.class} "
+  end
+  
+  def **(obj)
+    return Matrix.new(@MyRws,@MyCls) { |i,j| self[i,j] ** obj} if obj.is_a? Numeric
+    raise MyArgError, "  Argument Error: Invalid power Matrix-#{obj.class} "
+  end
+  
+  def tr()
+    return Matrix.new(@MyCls,@MyRws) { |i,j| @mx[j * @MyCls + i]}
+  end
+  
+  def norma()
+   sum=0
+   @mx.map do |e|
+     sum +=e**2
+   end
+   return Math.sqrt(sum)
+  end
+  
+  def is_squared?()
+    self.getCls == self.getRws ? (return true) : (return false)
+  end
+  
+  def laplace
+    if self.is_squared? then
+      return (self[0,0] * self[1,1]) - (self[0,1] * self[1,0]) if self.getRws == 2 
+      res = 0
+      for i in 0...@MyRws
+        res += (((-1) ** i) * self[i,0] * del_rwcl(self,i,0).laplace) if self[i,0] != 0
+      end
+      return res
+    else
+      raise MyArgError, '  Argument Error: Cannot calculate determinat on a non-squared matrix'
+      return nil
+    end
+  end
+  
+  def Matrix.identity(n)
+    return Matrix.new(n,n) { |i,j|  (i == j) ? 1 : 0 }
+  end
+  
+  def coerce(val)
+    return [self, val]
+  end
+  
+  private
+  
+   def del_rwcl(mx,rw,cl)
+     rd = mx.export.clone
+     cols = mx.getCls
+     for i in 0...cols
+       rd[i*cols+rw] = nil
+       rd[cl*(cols)+i] = nil
+     end
+     rd = rd.compact
+     rw = mx.getRws - 1
+     cl = mx.getCls - 1
+     return Matrix.new(rw,cl) { |i,j| rd[i*cl+j]}
+   end
+   
+  protected
+  
+   
+end
+
+
+
+
+
+
+

data/lib/linmeric/CnGal_new_classes.rb

@@ -0,0 +1,231 @@
+#! /usr/bin/env ruby
+
+require_relative 'Function_class.rb'
+
+# Adding some method to Numeric class
+# to work properly with CnGal project
+
+class Numeric 
+
+  def is_similar?(obj)
+    (obj.is_a? Numeric) ? (return true) : (return false)
+  end
+  
+  def can_multiply?(obj)
+    (obj.is_a? Numeric) ? (return true) : ((obj.is_a? Matrix) ? (return true) : (return false))
+  end
+  
+  def can_divide?(obj)
+    (obj.is_a? Numeric) ? (return true) : ((obj.is_a? Matrix) ? (return true) : (return false))
+  end
+
+  def show()
+    puts self
+  end
+
+end
+
+# Adding some method to Fixnum class
+# to work properly with CnGal project
+
+class Fixnum 
+  def similar_to?(obj)
+    (obj.is_a? Numeric) ? (return true) : (return false)
+  end
+  
+  def can_multiply?(obj)
+    (obj.is_a? Numeric) ? (return true) : ((obj.is_a? Matrix) ? (return true) : (return false))
+  end
+  
+  def can_divide?(obj)
+    (obj.is_a? Numeric) ? (return true) : ((obj.is_a? Matrix) ? (return true) : (return false))
+  end
+  
+  def show()
+    puts self
+  end
+end
+
+
+# Adding some method to Float class
+# to work properly with CnGal project
+
+class Float 
+  def similar_to?(obj)
+    (obj.is_a? Numeric) ? (return true) : (return false)
+  end
+  
+  def can_multiply?(obj)
+    (obj.is_a? Numeric) ? (return true) : ((obj.is_a? Matrix) ? (return true) : (return false))
+  end
+  
+  def can_divide?(obj)
+    (obj.is_a? Numeric) ? (return true) : ((obj.is_a? Matrix) ? (return true) : (return false))
+  end
+
+  def show()
+    puts self
+  end
+  
+end
+
+class NilClass
+  def show
+    puts "nil"
+  end
+  
+  def size
+    return 0
+  end
+  
+end
+
+# Adding some method to String class:
+# contain? returns true whether the current string contains one of the
+# chars in str
+# contain_all? returns true whether the current string contains all the 
+# chars in str
+# compact returns a string without spaces
+
+class String
+
+  def contain?(str)
+    raise ArgumentError, '' unless str.is_a? String
+    str.each_char do |ch|
+      return true if self.include? ch
+    end
+    return false
+    
+    rescue ArgumentError
+      str = str.to_s
+      retry
+  end
+  
+  def contain_all?(str)
+    raise ArgumentError, '' unless str.is_a? String
+    str.each_char do |ch|
+      return false if not self.include? ch
+    end
+    return true
+    rescue ArgumentError
+      str = str.to_s
+      retry
+  end
+  
+  def compact()
+    nstr = ''
+    self.each_char do |ch|
+      nstr += ch if ch != ' '
+    end
+    return nstr
+  end
+  
+  def -(str)
+    myStr = str.to_s
+    temp = ''
+    i = 0
+    if self.include? myStr then
+      while i < self.length do
+        if self[i...(i + myStr.size)] == myStr then
+          i += myStr.size
+        else
+          temp += self[i]
+          i += 1
+        end
+      end
+    else
+      return self
+    end
+    return temp
+  end
+  
+  def integer?
+    [                          # In descending order of likeliness:
+      /^[-+]?[1-9]([0-9]*)?$/, # decimal
+      /^0[0-7]+$/,             # octal
+      /^0x[0-9A-Fa-f]+$/,      # hexadecimal
+      /^0b[01]+$/,             # binary
+      /[0]/                     # zero
+    ].each do |match_pattern|
+      return true if self =~ match_pattern
+    end
+    return false
+  end
+
+  def float?
+    pat = /^[-+]?[1-9]([0-9]*)?\.[0-9]+$/
+    return true if self=~pat 
+    return false 
+  end
+  
+  def number?
+    (self.integer? or self.float?) ? (return true) : (return false)
+  end
+
+  def to_n
+    return self.to_f if self.float?
+    return self.to_i if self.integer? 
+  end
+  
+  def remove(index)
+    return self if index < 0
+    n_str = ''
+    for i in 0...self.size
+      n_str += self[i] unless i == index
+    end
+    return n_str
+  end
+end
+
+class Filename < String
+  def show
+    puts self
+  end
+  
+  def name
+    i = self.size - 1
+    nm = ''
+    while self[i] != '/' do
+      nm += self[i]
+      i -= 1
+    end
+    return nm.reverse
+  end
+  
+  def is_a?(obj)
+    self.class == obj.class ? (return true) : (return false)
+  end
+  
+end
+
+class Dim
+  def initialize(str)
+    str = str - '"'
+    str = str.split(',')
+    if str.size == 2 and str[0] != "" and str[1] != "" then
+      @sx = str[0].to_n if str[0].integer? or str[0].float?
+      @dx = str[1].to_n if str[1].integer? or str[1].float? 
+    else
+      puts '  Argument Error: wrong range format'#; puts
+    end
+  end
+  
+  def sx()
+    return @sx
+  end
+  
+  def dx()
+    return @dx
+  end
+end
+
+
+
+
+
+
+
+
+
+
+

data/lib/linmeric/CnGal_tools.rb

@@ -0,0 +1,80 @@
+#! /usr/bin/env ruby
+
+require_relative 'CnGal_new_classes.rb'
+
+module Tool
+  def self.keys 
+    return ["mx:","t:", "shw:","shwvar:","det:","solve:","from:","as:","f:","norm:","id_mx:","integ:"]
+  end
+  
+  def self.operators
+    return ["=","-","+","*","/","^",">"]
+  end
+  
+  def self.letters
+    return "abcdefghijklmnopqrstuvwxyz"
+  end
+  
+  def self.f
+    return ["log","sin","cos","PI","tan","exp"]
+  end
+  
+  def self.is_keyword?(ob)   
+    (self.keys.include? ob) ? (return true) : (return false)
+  end
+  
+  def self.is_exp?(str)
+    (str.contain? "=+*-/^") ? (return true) : (('0123456789.'.contain_all? str) ? 
+    (return true):((self.letters.contain? str) ? (return true):(return false)))
+  end
+  
+  def self.print_stack(stack)
+    stack.keys.each do |k|
+      puts "#{k} = " if stack[k].is_a? Matrix
+      print "#{k} = " unless stack[k].is_a? Matrix
+      stack[k].show
+    end
+  end
+  
+  def self.convert(string)
+     raise ArgumentError, '' unless string.is_a? String
+     if string.contain? '*/-+^' then
+       case
+         when (string.include? '^')
+           string = string.split('^')
+           string << '**'
+         when (string.include? '*')
+           string = string.split('*')
+           string << '*'
+         when (string.include? '/')
+           string = string.split('/')
+           string << '/'
+         when (string.include? '-')
+           string = string.split('-')
+           string << '-'
+         when (string.include? '+')
+           string = string.split('+')
+           string << '+'
+       end
+       if string.size == 3 then
+         return eval("convert(string[0]) #{string[2]} convert(string[1])")
+       elsif string.size == 2 then
+         return eval("#{string[2]} convert(string[0])")
+       else
+         return nil
+       end
+     elsif '0123456789.'.contain_all? string then
+         return string.to_f
+     else
+       return nil
+     end
+     
+     rescue ArgumentError
+       string = string.to_s
+       retry
+  end
+
+  
+end
+
+