math_expansion 0.0.1

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 02ca15e77c69bf92a0c178965cda06d7c3882d07
+  data.tar.gz: 65b6b03f99356b64cea1e9f442b49d7a778efe8b
+SHA512:
+  metadata.gz: 53f7eed81c242e10dd77602a6128e3a42ae6c815325d568f6419e8583cb0376782ab63ef67ddfad47ecc696101871ce1a3589a183969af3b86aa1f2db1e10fa8
+  data.tar.gz: 5ab60e8a0fb97c1210fa33a23364d1845e057d4757072bd8daee0e903a4d25f7f4382429f57c7e1c88a03151ab31b90368be6c661b35e67fce24959ecde08ede

data/.gitignore

@@ -0,0 +1,18 @@
+*.gem
+*.rbc
+.bundle
+.config
+.yardoc
+Gemfile.lock
+InstalledFiles
+_yardoc
+coverage
+doc/
+lib/bundler/man
+pkg
+rdoc
+spec/reports
+test/tmp
+test/version_tmp
+tmp
+*~

data/.travis.yml

@@ -0,0 +1,8 @@
+language: ruby
+rvm:
+  - 1.9.3
+  - jruby-18mode # JRuby in 1.8 mode
+  - jruby-19mode # JRuby in 1.9 mode
+  - rbx-18mode
+  - rbx-19mode
+  - 1.8.7

data/Gemfile

@@ -0,0 +1,7 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in math_expansion.gemspec
+gemspec
+gem 'guard'
+gem 'guard-rspec'
+gem 'guard-bundler'

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2013 Eliasib Garcia, Daniel Herzog
+
+MIT License
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+"Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/README.md

@@ -0,0 +1,26 @@
+# MathExpansion
+
+Útil gema de Ruby para representar matrices densas y dispersas.
+
+## Instalación
+
+Añada esta línea al fichero Gemfile de su aplicación:
+
+    gem 'math_expansion'
+
+Y a continuación ejecute:
+
+    $ bundle
+
+O instale manualmente:
+
+    $ gem install math_expansion
+
+## Planteamiento
+
+La clase Matriz contendrá el número de filas, número de columnas y un vector de contenido.
+
+La matriz densa se representará como en la práctica anterior.
+La matriz dispersa se representará como un Array que contendrá un Hash por cada fila. Cada par clave-valor del hash representará la fila de un elemento no nulo y el valor de ese elemento, respectivamente.
+
+Las operaciones entre matrices se hará de manera natural, operando con las matrices dadas independientemente de que sean densas o dispersas. Dado que en cualquier operación puede cambiar el número de valores nulos de manera impredecible, será en la matriz resultado (en principio densa) donde se hará una comprobación de los valores nulos de la misma, y si éste es superior al 60%, la matriz resultado se convertirá a dispersa; en caso contrario, la matriz resultado seguirá siendo densa.

data/Rakefile

@@ -0,0 +1,24 @@
+$:.unshift File.dirname(__FILE__) + 'lib'
+$:.unshift './lib', './spec'
+
+require "bundler/gem_tasks"
+
+require 'rspec/core/rake_task'
+RSpec::Core::RakeTask.new
+
+task :default => :spec
+
+desc "Espectativas de la clase Matriz" 
+task :test do
+	sh "rspec -I. spec/matriz_spec.rb --format documentation"
+end
+
+desc "Espectativas de la clase Matriz, con documentacion HTML"
+task :thtml do
+	sh "rspec -I. spec/matriz_spec.rb --format html"
+end
+
+desc "Pruebas unitarias de las clases Matriz_Densa y Matriz_Dispersa" 
+task :tc do
+	sh "ruby -I. test/tc_mathexpansion.rb"
+end

data/lib/math_expansion.rb

@@ -0,0 +1,12 @@
+require "./lib/math_expansion/version"
+require "./lib/math_expansion/matriz.rb"
+require "./lib/math_expansion/matriz_densa.rb"
+require "./lib/math_expansion/matriz_dispersa.rb"
+
+require "./lib/math_expansion/racional.rb"
+require "./lib/math_expansion/gcd.rb"
+
+# Vacío
+module MathExpansion
+
+end

data/lib/math_expansion/gcd.rb

@@ -0,0 +1,9 @@
+module MathExpansion
+  def self.gcd(a, b)
+    a, b = a.abs, b.abs
+    while b != 0
+      a, b = b, a % b
+    end
+    a
+  end
+end

data/lib/math_expansion/matriz.rb

@@ -0,0 +1,15 @@
+require './lib/math_expansion/racional.rb'
+
+module MathExpansion
+  class Matriz
+    attr_reader :N, :M
+    attr_accessor :contenido
+    
+    def initialize(n, m)
+      raise ArgumentError, 'Indice no valido' unless n.is_a? Fixnum and n > 0 and m.is_a? Fixnum and m > 0
+    
+      @N, @M = n, m
+    end
+
+  end
+end

data/lib/math_expansion/matriz_densa.rb

@@ -0,0 +1,191 @@
+require "./lib/math_expansion/matriz.rb"
+
+module MathExpansion
+  class Matriz_Densa < Matriz
+    def initialize(n, m)
+      super
+      
+      @contenido = Array.new(@N,0)
+      i = 0
+      while i < @N
+        @contenido[i] = Array.new(@M,0)
+        i += 1  
+      end
+    end
+      
+    def get(i, j)
+      if( i < 0 or i >=@N or j < 0 or j >= @M)
+        return nil
+      end
+
+      @contenido[i][j]
+    end
+    
+   def null_percent
+      total = @N*@M
+      no_nulos = 0
+      
+      i = 0
+      while(i < @N)
+        j = 0
+        while(j < @M)
+          if(@contenido[i][j] != @contenido[i][j].class.null)
+            no_nulos += 1
+	  end
+          j += 1
+        end
+        i += 1
+      end
+      
+      nulos = total - no_nulos
+      nulos.to_f/total.to_f
+    end #endmethod null_percent      
+
+    def set(i, j, value)
+      if( i < 0 or i >=@N or j < 0 or j >= @M)
+       return nil
+      end
+        
+      if(!(value.class.respond_to? :null))
+        puts "Se debe definir el metodo \"null\" que devuelva un elemento nulo para la clase #{value.class}"
+        return nil
+      end
+
+      # Contar elementos nulos y comprobar si se hace una matriz dispersa
+      # De momento, no dejamos añadir elementos nulos
+      # ¿o si?
+      #if(value != nil and value != value.class.null) # Y se puede comprobar para todos los tipos si es necesario. (con un método zero, por ejemplo)
+        @contenido[i][j] = value
+      #end
+    end
+    
+    def to_s
+      s = ""
+      i = 0
+      while(i < @N)
+        j = 0
+        while(j < @M)
+          s += "#{@contenido[i][j].to_s}\t"
+          j += 1
+        end
+        s += "\n"
+        i += 1
+      end
+      s
+    end
+    
+    def +(other)
+      		raise ArgumentError , 'Tipo invalido' unless other.is_a? Matriz
+      		raise ArgumentError , 'Matriz no compatible' unless @N == other.N and @M == other.M
+  
+      		c = Matriz_Densa.new(@N, @M)
+      		i = 0
+		while(i < @N)
+        		j = 0
+        		while(j < @M)
+	      			c.set(i, j, get(i,j) + other.get(i,j))
+    	  			j += 1
+    			end 
+    			i += 1
+      		end
+		if (c.null_percent > 0.6)
+			c = MathExpansion::Matriz_Dispersa.copy(c)
+		end
+      		c
+    end
+	
+    def -(other)
+      		raise ArgumentError , 'Tipo invalido' unless other.is_a? Matriz
+      		raise ArgumentError , 'Matriz no compatible' unless @N == other.N and @M == other.M
+  
+      		c = Matriz_Densa.new(@N, @M)
+      		i = 0
+      		while(i < @N)
+        		j = 0
+        		while(j < @M)
+	      			c.set(i, j, get(i,j) - other.get(i,j))
+    	  			j += 1
+    			end
+    			i += 1
+      		end
+		if (c.null_percent > 0.6)
+			c = MathExpansion::Matriz_Dispersa.copy(c)
+		end
+      		c
+    end
+	
+    def *(other)
+    		raise ArgumentError , 'Parametro invalido' unless other.is_a? Numeric or other.is_a? Matriz
+
+    		if(other.is_a? Numeric) # Matriz * numero
+      			c = Matriz_Densa.new(@N, @M)
+	  		i = 0
+      			while(i < @N)
+        			j = 0
+        			while(j < @M)
+         				c.set(i, j, get(i,j)*other)
+         				j += 1
+	    			end # while j
+	    			i += 1
+      			end # while i
+    		else # Matriz * Matriz
+      			raise ArgumentError , 'Matriz no compatible (A.N == B.M)' unless @M == other.N
+      			c = Matriz_Densa.new(@N, other.M)
+      			i = 0
+      			while(i < @N)
+        			j = 0
+        			while(j < other.M)
+	      				k = 0
+					#if (get(i,j).is_a? Fraccion)	      				
+					#	c.set(i, j, Fraccion.null
+					#else
+					#	c.set(i, j, 0)
+					#end
+	      				while(k < @M)
+	        				c.set(i, j, get(i, k) * other.get(k,j) + c.get(i,j))
+	        				k += 1
+	      				end # while k
+          				j += 1
+        			end # while j
+	    			i += 1
+      			end # while i
+    		end # while else
+  		if (c.null_percent > 0.6)
+			c = MathExpansion::Matriz_Dispersa.copy(c)
+		end
+    		c
+    end # *(other)
+
+    def max
+	m = get(0,0)
+        i = 0
+	while(i < @N)
+		j = 0
+		while(j < @M)
+			if (get(i,j) > m)
+				m = get(i,j)
+			end
+			j += 1
+		end
+		i += 1
+	end
+	m
+    end
+
+    def min
+	m = get(0,0)
+        i = 0
+	while(i < @N)
+		j = 0
+		while(j < @M)
+			if (get(i,j) < m)
+				m = get(i,j)
+			end
+			j += 1
+		end
+		i += 1
+	end
+	m
+    end # Method min
+  end # Class
+end # Module

data/lib/math_expansion/matriz_dispersa.rb

@@ -0,0 +1,250 @@
+require "./lib/math_expansion/matriz.rb"
+require "./lib/math_expansion/matriz_densa.rb"
+
+module MathExpansion
+  class Matriz_Dispersa < Matriz  
+    def reset
+      @contenido = Array.new(@N) # Array con @N filas y ninguna columna (vacio)
+      i = 0
+      while(i < @N)
+        @contenido[i] = {}
+        i += 1
+      end
+    end
+    
+    def initialize(n, m)
+      super
+      reset
+    end
+    
+    # Metodo factoria
+    def self.copy(matriz)
+      raise ArgumentError, 'Tipo invalido' unless matriz.is_a? MathExpansion::Matriz_Densa
+      obj = new(matriz.N, matriz.M)
+    
+      i = 0
+      while(i < matriz.N)
+        j = 0
+        while(j < matriz.M)
+          value = matriz.get(i,j)
+          raise RuntimeError , "No se ha definido \"null\" para la clase #{value.class}" unless value.class.respond_to? :null
+          
+          if( value != value.class.null)
+            obj.contenido[i][j] = value
+          end #endif
+          j += 1
+        end #endwhile j
+        i += 1
+      end #endwhile i
+      obj
+    end #endmethod copy
+    
+    def null_percent
+      total = @N*@M
+      no_nulos = 0
+      
+      i = 0
+      while(i < @N)
+        no_nulos += @contenido[i].size # Nunca habra elementos nulos en alguna fila
+        i += 1
+      end
+      
+      nulos = total - no_nulos
+      nulos.to_f/total.to_f
+    end #endmethod null_percent
+    
+    def get(i, j)
+      if( !(i.is_a? Fixnum) or i < 0 or i >=@N or !(j.is_a? Fixnum) or j < 0 or j >= @M)
+        return nil
+      end
+        
+      if(@contenido[i][j] != nil) # Elemento no nulo (esta en el hash)
+        return @contenido[i][j]
+      else # Elemento nulo (no esta en el Hash)
+        return 0
+      end
+    end #endmethod get
+    
+    def set(i, j, value)
+      if(!(value.class.respond_to? :null))
+        puts "Se debe definir el metodo \"null\" que devuelva un elemento nulo para la clase #{value.class}"
+        return nil
+      end
+      
+      if( !(i.is_a? Fixnum) or i < 0 or i >=@N or !(j.is_a? Fixnum) or j < 0 or j >= @M)
+        return nil
+      end
+      
+      if(value == nil or value == value.class.null)
+        @contenido[i].delete(j) # Borrar elemento (valor nulo)
+      else
+        @contenido[i][j] = value
+      end
+      
+      if(null_percent < 0.6) # Si se ha sobrepasado el número de elementos nulos, borramos el último elemento modificado
+          @contenido[i].delete(j)
+          puts "Borrado el elemento #{i},#{j} por sobrepasar el numero de elementos no nulos (Porcentaje actual: #{null_percent}"
+      end
+      
+    end #endmethod set
+    
+    def to_s
+      # Ejemplo: "Fila 0: \nFila 1: 0=>1 1=>3 \nFila 2: \n"
+      # 0 0
+      # 1 3
+      # 0 0
+      i = 0
+      output = ""
+      while(i < @N)
+        output += "Fila #{i}: "
+        @contenido[i].sort.each{|k, v| output += "#{k.to_s}=>#{v.to_s} "}
+        output += "\n"
+        i += 1
+      end
+      output
+    end
+
+    def +(other)
+      		raise ArgumentError , 'Tipo invalido' unless other.is_a? Matriz
+      		raise ArgumentError , 'Matriz no compatible' unless @N == other.N and @M == other.M
+  
+      		c = Matriz_Densa.new(@N, @M)
+      		i = 0
+		while(i < @N)
+        		j = 0
+        		while(j < @M)
+	      			c.set(i, j, get(i,j) + other.get(i,j))
+    	  			j += 1
+    			end 
+    			i += 1
+      		end
+		if(c.null_percent > 0.6)
+			c = Matriz_Dispersa.copy(c)
+		end
+      		c
+    end
+	
+    def -(other)
+      		raise ArgumentError , 'Tipo invalido' unless other.is_a? Matriz
+      		raise ArgumentError , 'Matriz no compatible' unless @N == other.N and @M == other.M
+  
+      		c = Matriz_Densa.new(@N, @M)
+      		i = 0
+      		while(i < @N)
+        		j = 0
+        		while(j < @M)
+	      			c.set(i, j, get(i,j) - other.get(i,j))
+    	  			j += 1
+    			end
+    			i += 1
+      		end
+		if(c.null_percent > 0.6)
+			c = Matriz_Dispersa.copy(c)
+		end		
+      		c
+    end
+	
+    def *(other)
+    		raise ArgumentError , 'Parametro invalido' unless other.is_a? Numeric or other.is_a? Matriz
+
+    		if(other.is_a? Numeric) # Matriz * numero
+      			c = Matriz_Densa.new(@N, @M)
+	  		i = 0
+      			while(i < @N)
+        			j = 0
+        			while(j < @M)
+         				c.set(i, j, get(i,j)*other)
+         				j += 1
+	    			end # while j
+	    			i += 1
+      			end # while i
+    		else # Matriz * Matriz
+      			raise ArgumentError , 'Matriz no compatible (A.N == B.M)' unless @M == other.N
+      			c = Matriz_Densa.new(@N, other.M)
+      			i = 0
+      			while(i < @N)
+        			j = 0
+        			while(j < other.M)
+                    k = 0
+					#if (get(i,j).is_a? Fraccion)	      				
+					#	c.set(i, j, Fraccion.null
+					#else
+					#	c.set(i, j, 0)
+					#end
+                        while(k < @M)
+	        				c.set(i, j, get(i, k) * other.get(k,j) + c.get(i,j))
+	        				k += 1
+	      				end # while k
+          				j += 1
+        			end # while j
+	    			i += 1
+      			end # while i
+    		end # while else
+		if(c.null_percent > 0.6)
+			c = Matriz_Dispersa.copy(c)
+		end  
+    		c
+    end # *(other)
+	
+	def max
+	  if(null_percent == 1.0)
+	    return nil # o return 0
+	  end
+		
+	  # Valor máximo: si todos los elementos son menores que el elemento nulo
+	  # Se devolverá el mayor elemento no nulo.
+	  max = nil
+	  
+	  # Asignar al primer valor no-nulo de la matriz
+	  i = 0
+	  while(max == nil)
+	    if(@contenido[i].size != 0)
+		  max = @contenido[i].values[0]
+		end
+		i += 1
+	  end
+	  
+	  # Iterar por todos los elementos no nulos para encontrar el maximo
+	  i = 0
+	  while(i < @contenido.size)
+	    if(@contenido[i].values.max != nil and @contenido[i].values.max > max)
+		  max = @contenido[i].values.max
+		end
+	    i += 1
+	  end
+	  
+	  max
+	end
+	
+	def min
+	  if(null_percent == 1.0)
+	    return nil # o return 0
+	  end
+		
+	  # Valor máximo: si todos los elementos son menores que el elemento nulo
+	  # Se devolverá el mayor elemento no nulo.
+	  min = nil
+	  
+	  # Asignar al primer valor no-nulo de la matriz
+	  i = 0
+	  while(min == nil)
+	    if(@contenido[i].size != 0)
+		  min = @contenido[i].values[0]
+		end
+		i += 1
+	  end
+	  
+	  # Iterar por todos los elementos no nulos para encontrar el maximo
+	  i = 0
+	  while(i < @contenido.size)
+	    if(@contenido[i].values.min != nil and @contenido[i].values.min < min)
+		  min = @contenido[i].values.min
+		end
+	    i += 1
+	  end
+	  
+	  min
+    end
+    
+  end #endclass
+end #end module