t50_matriz_dispersa_densa 0.0.1

data/lib/matriz_densa.rb

@@ -0,0 +1,128 @@
+require "matriz.rb"
+
+class Matriz_densa < Matriz
+
+	# Operacion de la suma en la matriz densa. Realiza la operacion entre dos matrices densa y una densa con una dispersa.
+        def +(other)
+           
+			  
+			  case other
+			  when Matriz_densa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @f == other.f && @c == other.c
+		           r =@m
+		     
+		          #for i in 0...@f do
+						0.upto(@f) do |i|
+		                   #for j in 0...@c do
+								 0.upto(@c) do |j|
+		                         r[i][j]=@m[i][j]+other.m[i][j]
+		                   end
+		          end
+		          return r
+			  when Matriz_Dispersa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @f == other.matrix.size && @c == other.matrix[1].vector.size
+					other + self
+				
+			  else
+					raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+			  end
+        end
+
+	# Operacion de la resta en la matriz densa. Realiza la operacion entre dos matrices densa y una densa con una dispersa.
+        def -(other)
+          
+           case other
+			  when Matriz_densa
+					  raise ArgumentError, "Las matrices no son iguales en longitud." unless @f == other.f && @c == other.c
+		           rs =@m
+		     
+		          #for i in 0...@f do
+						@f.times do |i|
+		                   #for j in 0...@c do
+								 @c.times do |j|
+		                         rs[i][j]=@m[i][j]+other.m[i][j]
+		                   end
+		          end
+		          return rs
+			  when Matriz_Dispersa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @f == other.matrix.size && @c == other.matrix[1].vector.size
+					other - self
+				
+			  else
+					raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+			  end
+        end
+
+	# Operación de la multiplicación en la matriz densa. Realiza la operacion entre dos matrices densa y una densa con una dispersa.
+        def *(other)				
+		  		case other
+				when Matriz_densa
+					raise ArgumentError, "Las columnas de una matriz no coinciden con las filas de la otra." unless @c == other.f
+                z = Array.new
+                for i in 0...@f do
+                   z[i] = Array.new
+                   for j in 0...@c do
+                        z[i][j] = 0
+                   end
+                end
+                #for i in 0...@f do
+						0.upto(@f) do |i| 
+                        #for j in 0...@c do
+								0.upto(@c) do |j|
+                                    #for k in 0...@f do
+												0.upto(@f) do |k|
+                                        z[i][j] += @m[i][k] * other.m[k][j]
+                                    end
+                        end
+                    end        
+                return z
+				when Matriz_Dispersa
+					raise ArgumentError, "Las columnas de una matriz no coinciden con las filas de la otra." unless @c == other.matrix.size
+					other * self
+				else
+					raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+				end
+        end
+
+        
+        def to_s
+                for i in 0...@f do
+                        for j in 0...@c do
+                                print @m[i][j]        
+                        end
+                        puts
+                end
+                
+        end
+        
+	# Devuelve el maximo de una matriz densa
+        def max
+                
+                maximo = @m[0][0].to_f
+                for i in 0...@f do
+                        for j in 0...@c do
+                                if @m[i][j].to_f > maximo.to_f
+                                        maximo = @m[i][j]
+                                end
+                        end
+                end
+                maximo
+
+        end
+
+	# Devuelve el minimo de una matriz densa
+        def min
+
+                minimo = @m[0][0].to_f
+                for i in 0...@f do
+                        for j in 0...@c do
+                                if @m[i][j].to_f < minimo.to_f
+                                        minimo = @m[i][j]
+                                end
+                        end
+                end
+                minimo
+
+        end
+
+end

data/lib/matriz_densa.rb~

@@ -0,0 +1,148 @@
+require "matriz.rb"
+
+class Matriz_densa < Matriz
+
+	# Operacion de la suma en la matriz densa. Realiza la operacion entre dos matrices densa y una densa con una dispersa.
+        def +(other)
+           
+			  
+			  case other
+			  when Matriz_densa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @f == other.f && @c == other.c
+		           r =@m
+		     
+		          #for i in 0...@f do
+						0.upto(@f) do |i|
+		                   #for j in 0...@c do
+								 0.upto(@c) do |j|
+		                         r[i][j]=@m[i][j]+other.m[i][j]
+		                   end
+		          end
+		          return r
+			  when Matriz_Dispersa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @f == other.matrix.size && @c == other.matrix[1].vector.size
+					other + self
+				
+			  else
+					raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+			  end
+        end
+
+	# Operacion de la resta en la matriz densa. Realiza la operacion entre dos matrices densa y una densa con una dispersa.
+        def -(other)
+          
+           case other
+			  when Matriz_densa
+					  raise ArgumentError, "Las matrices no son iguales en longitud." unless @f == other.f && @c == other.c
+		           rs =@m
+		     
+		          #for i in 0...@f do
+						@f.times do |i|
+		                   #for j in 0...@c do
+								 @c.times do |j|
+		                         rs[i][j]=@m[i][j]+other.m[i][j]
+		                   end
+		          end
+		          return rs
+			  when Matriz_Dispersa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @f == other.matrix.size && @c == other.matrix[1].vector.size
+					other - self
+				
+			  else
+					raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+			  end
+        end
+
+	# Operación de la multiplicación en la matriz densa. Realiza la operacion entre dos matrices densa y una densa con una dispersa.
+        def *(other)				
+		  		case other
+				when Matriz_densa
+					raise ArgumentError, "Las columnas de una matriz no coinciden con las filas de la otra." unless @c == other.f
+                z = Array.new
+                for i in 0...@f do
+                   z[i] = Array.new
+                   for j in 0...@c do
+                        z[i][j] = 0
+                   end
+                end
+                #for i in 0...@f do
+						0.upto(@f) do |i| 
+                        #for j in 0...@c do
+								0.upto(@c) do |j|
+                                    #for k in 0...@f do
+												0.upto(@f) do |k|
+                                        z[i][j] += @m[i][k] * other.m[k][j]
+                                    end
+                        end
+                    end        
+                return z
+				when Matriz_Dispersa
+					raise ArgumentError, "Las columnas de una matriz no coinciden con las filas de la otra." unless @c == other.matrix.size
+					other * self
+				else
+					raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+				end
+        end
+
+        
+        def to_s
+                for i in 0...@f do
+                        for j in 0...@c do
+                                print @m[i][j]        
+                        end
+                        puts
+                end
+                
+        end
+
+        def /(other)
+        
+                w = Array.new
+                for i in 0...@f do
+                   w[i] = Array.new
+                   for j in 0...@c do
+                        w[i][j] = 0
+                   end
+                end
+
+                for i in 0...@f do
+                        for j in 0...@c do
+                                    for k in 0...@f - 1 do
+                                        w[i][j] += @m[i][k] * other.m[i][j]
+                                    end
+                        end
+                    end
+                return w       
+        end
+        
+	# Devuelve el máximo de una matriz densa
+        def max
+                
+                maximo = @m[0][0].to_f
+                for i in 0...@f do
+                        for j in 0...@c do
+                                if @m[i][j].to_f > maximo.to_f
+                                        maximo = @m[i][j]
+                                end
+                        end
+                end
+                maximo
+
+        end
+
+	# Devuelve el minimo de una matriz densa
+        def min
+
+                minimo = @m[0][0].to_f
+                for i in 0...@f do
+                        for j in 0...@c do
+                                if @m[i][j].to_f < minimo.to_f
+                                        minimo = @m[i][j]
+                                end
+                        end
+                end
+                minimo
+
+        end
+
+end

data/lib/matriz_dispersa.rb

@@ -0,0 +1,276 @@
+require 'Fraccion.rb'
+require 'matriz.rb'
+require "matrix"
+
+class SparseVector
+    attr_reader :vector
+
+    def initialize(h = {})
+        @vector = Hash.new(0)
+        @vector = @vector.merge!(h)
+    end
+
+    def [](i)
+        @vector[i]
+    end
+
+    def to_s
+        @vector.to_s
+    end
+end
+
+#Clase que se encarga de crear el hash
+class SparseMatrix < Matriz
+
+    attr_reader :matrix
+
+    # Inicializador propio de la clase
+    def initialize(h = {})
+        @matrix = Hash.new({})
+        for k in h.keys do
+            @matrix[k] = if h[k].is_a? SparseVector 
+                                h[k]
+                        else
+                              @matrix[k] = SparseVector.new(h[k])
+                            end
+        end
+    end
+
+    # Obtener un valor de la matriz dada una posicion
+    def [](i)
+        @matrix[i]
+    end
+
+    def col(j)
+        c = {}
+        for r in @matrix.keys do
+            c[r] = @matrix[r].vector[j] if @matrix[r].vector.keys.include? j
+        end
+        SparseVector.new c
+    end
+
+        def comprobar
+                contar = 0
+                contartotal = 0
+                result = 0
+           for r in @matrix.keys do
+                         for j in @matrix[r].vector.keys do
+                                contartotal = contartotal + 1
+                                if @matrix[r].vector[j] == 0 
+                                        contar = contar + 1
+                                end
+                        end
+      end
+                result = (contar * 100) / contartotal
+                if result < 60 
+                         "La matriz no es dispersa"
+                end
+        end
+
+        def mostrar
+        for r in @matrix.keys do
+                                  for j in @matrix[r].vector.keys do
+                                        print "#{@matrix[r].vector[j]}  "  
+                                end
+                                puts
+        end
+        end
+
+	# Operacion de la suma en la matriz dispersa. Realiza la operacion entre dos matrices dispersa y una dispersa con una densa    
+        def +(other)
+			  case other
+			  when Matriz_Dispersa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @matrix.size == other.matrix.size && @matrix[1].vector.size == other.matrix[1].vector.size
+					resul = Array.new
+					 @matrix.keys.times do |r|
+					 r = r + 1
+					  #for r in @matrix.keys do
+						resul_f = Array.new
+						 1.upto(@matrix[r].vector.keys) do |j|  
+						 #for j in @matrix[r].vector.keys do       
+					  		#print "  #{@matrix[r].vector[j] + other.matrix[r].vector[j]}  "
+							resul_f << @matrix[r].vector[j] + other.matrix[r].vector[j]
+					  	end
+					  		#puts
+						resul << resul_f
+					end
+					return resul
+			  when Matriz_densa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @matrix.size == other.f && @matrix[1].vector.size == other.c
+					resul = Array.new
+					1.upto(other.f) do |r|
+					 #for r in 1..other.f do
+						resul_f = Array.new
+						1.upto(other.c) do |j|
+						 #for j in 1..other.c do
+               		#print "  #{@matrix[r].vector[j] + other.m[r - 1][j - 1]}"
+							resul_f << @matrix[r].vector[j] + other.m[r - 1][j - 1]
+               	end
+               	#puts
+						resul << resul_f
+        			end
+					return resul
+			  else
+					raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+			  end
+        end
+	
+	# Operacion de la resta en la matriz dispersa. Realiza la operacion entre dos matrices dispersa y una dispersa con una densa
+        def -(other)
+			  case other
+			  when Matriz_Dispersa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @matrix.size == other.matrix.size && @matrix[1].vector.size == other.matrix[1].vector.size
+					resul = Array.new
+		     		#for r in @matrix.keys do
+					 1.upto(@matrix.keys) do |r|
+						resul_f = Array.new
+					  	#for j in @matrix[r].vector.keys do  
+						 1.upto(@matrix[r].vector.keys) do |j|      
+					  		#print "  #{@matrix[r].vector[j] - other.matrix[r].vector[j]}  "
+							resul_f << @matrix[r].vector[j] - other.matrix[r].vector[j]
+					  	end
+					  		#puts
+						resul << resul_f
+					end
+					return resul
+			  when Matriz_densa
+					raise ArgumentError, "Las matrices no son iguales en longitud." unless @matrix.size == other.f && @matrix[1].vector.size == other.c
+					resul = Array.new
+					#for r in 1..other.f do
+					1.upto(other.f) do |r|
+						resul_f = Array.new
+               	#for j in 1..other.c do
+						1.upto(other.c) do |j|
+               		#print "  #{@matrix[r].vector[j] - other.m[r - 1][j - 1]}"
+							resul_f << @matrix[r].vector[j] - other.m[r - 1][j - 1]
+               	end
+               	#puts
+						resul << resul_f
+        			end
+					return resul
+			  else
+					raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+			  end
+
+        end
+
+	# Operacion de la multiplicacion en la matriz dispersa. Realiza la operacion entre dos matrices dispersa y una dispersa con una densa
+        def *(other) 
+        
+				case other
+				when Matriz_Dispersa
+					 raise ArgumentError, "Las columnas de una matriz no coinciden con las filas de la otra." unless @matrix[1].vector.size == other.matrix.size
+                sumatotal = 0
+                mul = Array.new                
+
+               #for i in @matrix.keys do
+					1.upto(@matrix.keys) do |i|
+                        mul_f = Array.new
+                     #for j in @matrix[i].vector.keys do
+							 1.upto(@matrix[i].vector.keys) do |j|
+                                #for k in @matrix[i].vector.keys do
+											 1.upto(@matrix[i].vector.keys) do |k|
+                                    suma = @matrix[i].vector[k] * other.matrix[k].vector[j]
+                                    sumatotal = sumatotal + suma;
+                                end
+                                mul_f << sumatotal
+                                sumatotal = 0
+                    end
+							mul << mul_f
+                end
+=begin        
+                for r in @matrix.keys do
+                   for j in @matrix[r].vector.keys do
+                      print "  #{mul[r][j]}  "
+                   end
+                   puts
+        			end
+=end
+					return mul
+			  when Matriz_densa
+					raise ArgumentError, "Las columnas de una matriz no coinciden con las filas de la otra." unless @matrix[1].vector.size == other.f
+					sumatotal = 0
+                mul = Array.new                
+
+                #for i in @matrix.keys do
+					  1.upto(@matrix.size) do |i|
+                        mul_f = Array.new
+                     #for j in @matrix[i].vector.keys do
+							 1.upto(@matrix[i].vector.size) do |j|
+                           #for k in @matrix[i].vector.keys do
+									1.upto(@matrix[i].vector.size) do |k|
+                              suma = @matrix[i].vector[k] * other.m[k - 1][j - 1]
+                              sumatotal = sumatotal + suma;
+                           end
+                           mul_f << sumatotal
+                           sumatotal = 0
+                    end
+							mul << mul_f
+                end 
+=begin       
+                for r in @matrix.keys do
+                   for j in @matrix[r].vector.keys do
+                      print "  #{mul[r][j]}  "
+                   end
+                   puts
+        			end
+=end
+					return mul
+			  else
+				  raise TypeError.new("No se puede pasar #{other.inspect} a Matriz")
+			  end
+        end
+
+	# Obtención del maximo de una matriz dispersa
+
+        def max
+                   
+                maximo = @matrix[1].vector[1]
+           
+                #for i in @matrix.keys do
+						@matrix.size.times do |i|
+								i = i + 1
+                        #for j in @matrix[i].vector.keys do
+								#@matrix[i].vector.keys do |j|
+								1.upto(@matrix[i].vector.size) do |j|
+                           if @matrix[i].vector[j] > maximo
+                              maximo = @matrix[i].vector[j]
+                      end
+                   end
+                end
+					 return maximo
+        end
+
+		  def [](i,j)
+			
+				return @matrix[i].vector[j]
+			
+		  end
+
+	# Obtención del maximo de una matriz dispersa
+        def min
+
+           minimo = @matrix[1].vector[1]
+                
+                #for i in @matrix.keys do
+						1.upto(@matrix.size) do |i|
+                        #for j in @matrix[i].vector.keys do
+								@matrix[i].vector.size.times do |j|
+								j = j + 1
+                           if @matrix[i].vector[j] < minimo
+                              minimo = @matrix[i].vector[j]
+                      end
+                   end
+                end
+					return minimo
+
+        end
+
+        def coerce(other)
+                [other,self]
+        end
+end
+
+class Matriz_Dispersa < SparseMatrix
+
+end