Source: http://docplayer.es/1575150-Revista-investigacion-operacional-vol-35-no-3-201-214-2014.html
Timestamp: 2017-04-24 11:45:43+00:00

Document:
REVISTA INVESTIGACIÓN OPERACIONAL VOL. 35, NO. 3, , PDF
REVISTA INVESTIGACIÓN OPERACIONAL VOL. 35, NO. 3, , 2014
Download "REVISTA INVESTIGACIÓN OPERACIONAL VOL. 35, NO. 3, 201-214, 2014"
Javier Pinto Torregrosa
1 REVISTA INVESTIGACIÓN OPERACIONAL VOL. 35, NO. 3, , 2014 W_ILU_ GMRES: UN DISEÑO EN PARALELO Yisleidy Linares Zaila, Marta L. Baguer Díaz Romañach, Ángela León Mecías y Universidad de La Habana, Facultad de Matemática y Computación, San Lázaro y L, Vedado 10400, Habana 4, Cuba ABSTRACT This paper proposes a parallel design of W_ILU_GMRES strategy, which is an iterative technique for solving systems of linear equations with dense and large matrices. First basic concepts related to parallelism are explained. Then we propose a parallel design, which uses multi-core technology and therefore has shared memory architecture. As a language for developing W_ILU_GMRES was chosen FSharp. Some experiments were performed in order to measure the quality of the parallel implementation. KEY WORDS: parallel computing, preconditioning, parallel wavelet compression, Krylov subspaces methods. MSC: 68W10 RESUMEN En este trabajo se propone un diseño en paralelo de la estrategia W_ILU_GMRES, la cual es una técnica iterativa para resolver sistemas de ecuaciones lineales con matriz densa y de grandes dimensiones. Primeramente se explican conceptos básicos relacionados con el paralelismo. Después se presenta el diseño en paralelo que utiliza tecnología multi-core y por tanto está implementado sobre una arquitectura de memoria compartida. Como lenguaje para el desarrollo de W_ILU_GMRES se escogió FSharp. Se realizaron algunos experimentos para medir la calidad de la implementación en paralelo. 1. INTRODUCCIÓN La solución de sistemas de ecuaciones lineales de grandes dimensiones con matrices densas por métodos directos trae aparejado un alto costo computacional, de ahí que se empleen como vía alternativa técnicas iterativas, especialmente las basadas en los subespacios de Krylov (Y. Saad, 2000), que aplicadas de conjunto con precondicionadores permiten acelerar la convergencia. En un trabajo reciente (L. Villarín, A. León, M. L. Baguer, Y. Linares, 2012), se propuso la construcción de un precondicionador para GMRES, al cual denominamos W_ILU_GMRES. La construcción del precondicionador está basada en la factorización LU incompleta de la aproximación dispersa no singular de la matriz densa. El diseño e implementación del algoritmo fue presentado en forma secuencial. Como resultado de la experimentación numérica realizada en (Y. Saad, 2000), se obtuvo que cuando la cantidad de elementos no nulos en la aproximación dispersa crece, la calidad del precondicionador también crece, sin embargo, su costo de obtención crece en el mismo orden. Además, la estrategia de selección automática del umbral que se propone se basa en un algoritmo de ordenación, con costo, donde es la dimensión de la matriz, que se impone como cota inferior de la propuesta. Como también se explica en (L. Villarín, A. León, M. L. Baguer, Y. Linares, 2012), para grandes dimensiones se requieren aproximaciones dispersas con una cantidad de elementos no nulos mayor que, entonces la obtención de la factorización LU Incompleta (ILU) tiene un costo mayor que, dominando el costo total de la estrategia. Por lo tanto, al aumentar la cantidad de elementos no nulos de la matriz, aumenta el costo de la obtención de la factorización ILU. Para matrices dispersas, obtener el precondicionador ILU tendría un costo, con lo que no se podría mejorar la solución vía métodos directos. La interrogante sería, cómo encontrar una alternativa para superar esta limitante? En los últimos años el número de aplicaciones que involucran el uso del paralelismo ha sido incrementado en gran medida. Existe una amplia variedad de algoritmos cuyo objetivo es paralelizar tanto precondicionadores como técnicas iterativas que resuelven sistemas de ecuaciones lineales de grandes dimensiones. Sin embargo, la mayoría de estos, véase (R. D. D. Cunha, T. Hopkins, 1994), (J. Erhel, 1995), (M. Benzi, W. Joubert, G. Mateescu, 1999) están desarrollados para plataformas y arquitecturas específicas tales como clústeres y sistemas con memoria distribuida respectivamente. En la revisión bibliográfica realizada no se encontró trabajo alguno acerca del diseño en paralelo de estos algoritmos implementados sobre las nuevas tecnologías, como la tecnología multi-core. 2012 En el presente trabajo se propone un diseño en paralelo de W_ILU_GMRES, estructurado en tres secciones además de la introducción. En la segunda sección se introducen algunos conceptos esenciales relacionados con el desarrollo de programas en paralelo, así como las principales características del hardware utilizado y del lenguaje FSharp, este último seleccionado para la implementación de la estrategia. En la tercera sección se presenta el diseño en paralelo así como la estructura de datos utilizada para la representación de las matrices dispersas, que constituye el núcleo fundamental del trabajo. El análisis de los costos de cada una de las etapas y los resultados de los experimentos realizados aparecen en la última sección. 2. COMPUTACIÓN EN PARALELO. ALGUNAS DEFINICIONES 2.1. Hardware En este trabajo se utilizó la tecnología multi-core. Esta consiste en computadoras cuyo procesador contiene más de un núcleo en el mismo paquete físico. El objetivo de este diseño es lograr que varias tareas sean ejecutadas simultáneamente, y por tanto mejorar el desempeño de todo el sistema. Esta tecnología, según la Taxonomía de Flynn 1, encontrada en (B. Barney, 2009) se sitúa en el modelo MIMD (Multiple Instruction Multiple Data) donde diferentes procesadores ejecutan diferentes instrucciones sobre distintos segmentos de datos (Figura 1) Por otra parte, atendiendo al acceso a la memoria RAM, la tecnología multi-core presenta una arquitectura de memoria compartida (Shared Memory)(B. Barney, 2009) Software Figura 1Modelo MIMD El primer paso para diseñar un algoritmo en paralelo es identificar aquellas partes que pueden ejecutarse simultáneamente. Es preciso enfocarse en los puntos críticos, es decir, aquellas secciones que requieren un alto consumo de CPU e identificar las partes que no se pueden paralelizar. Posteriormente se pasa al diseño de la solución, el cual está dividido en fases entre las que se encuentran: el particionado, la comunicación entre los diferentes procesadores y la fase de sincronización. La primera es la etapa donde se define la forma en la que se fragmenta el problema para poder diseñar una solución del mismo en paralelo. En la segunda se determina el modo en el que los procesadores trasmitirán la información. En la última fase se analiza la dependencia de las distintas tareas garantizando que no se viole el orden de precedencia. Para una explicación más detallada véase(i. Foster, 1995). Otro aspecto importante a tener en cuenta es la eficiencia del programa en paralelo. El factor determinante en el costo de un algoritmo en paralelo no es el número de procesadores sino el porciento del problema que puede escribirse en paralelo. Así lo demuestra la ley de Amdahl 2, formulada en Michael J. Flynn (nacido el 20 de mayo de 1934 en Nueva York) es un profesor emérito de la Universidad Stanford estadounidense, con estudios en ingeniería electrónica y ciencias de la computación. En 1995 recibió el premio Harry H. Goode Memorial Award por sus contribuciones al área del procesamiento de la información. 2 Eugene MyronAmdahl (nacido el 16 de noviembre de 1922, Flandreau, Dakota del Sur) norteamericano de origen noruego, arquitecto computacional y una de las personalidades más importantes y excéntricas en la historia de la informática y la computación. Fundó cuatro compañías tecnológicas en diferentes ámbitos, la mayoría de las cuales, se fundaron con el objetivo de 2023 por Gene Amdahl, la cual establece que la porción del problema que no se puede paralelizar constituye un límite en la velocidad de procesamiento, véase(i. Foster, 1995). Una forma práctica de medir el aumento de esta velocidad es mediante la fórmula: donde es el tiempo de ejecución del programa secuencial y es el tiempo de ejecución en paralelo con procesadores. y son funciones que dependen de la dimensión de la matriz. Esta fórmula introduce un concepto importante en la programación en paralelo: la escalabilidad. Un programa en paralelo se dice escalable cuando al aumentar el número de procesadores aumenta la velocidad de ejecución Microsoft y el paralelismo La industria de los microprocesadores se ha inclinado hacia la producción de procesadores multi-core. Debido a este cambio, Visual Studio 2010 y.net 4.0, incluyen la biblioteca Task Parallel Library (TPL) y un nuevo lenguaje funcional y orientado a objetos llamado FSharp (F#) que hace la programación paralela y asíncrona más sencilla en un ambiente amigable. Más detalles acerca de cada una de estas bibliotecas se encuentran en(microsoft Developers, 2010) FSharp F# combina las cualidades de la programación funcional, tales como código claro y conciso, con las bibliotecas, herramientas, interoperabilidad y modelo de objetos de la plataforma.net. Es el único lenguaje con esta característica, es decir, el único capaz de sacar lo mejor de ambos paradigmas (el funcional y el orientado a objetos). Como ventajas significativas de este lenguaje se tienen que a través del mismo se puede extender la estrategia a una aplicación MPI 3 lo cual es de gran importancia teniendo en cuenta que se lograría incrementar mucho más la dimensión de los sistemas de ecuaciones lineales. También en FSharp las estructuras de datos son inmutables, lo cual permite un mayor aprovechamiento del paralelismo. Otro factor importante es que su uso no solo se restringe a Windows, este también puede utilizarse en Linux, Apple Mac OS X and OpenBSD. Para una argumentación más profunda véanse (J. Harrop, 2008),(R. Pickering, 2007), (D. Syme, 2007) Factorización LU Incompleta (ILU) Este proceso se divide en dos pasos fundamentales, el primero es la obtención de los factores y y el segundo es la aplicación de esta factorización, es decir, la sustitución hacia adelante y retrógrada. Existen varias formas de realizar esta factorización LU Incompleta, siguiendo diferentes estrategias de creación de relleno, ver (Y. Saad, 2000); aquí se utiliza aquella que no incorpora nuevas entradas no nulas durante su obtención, conocida como ILU (0). Este algoritmo está basado en la eliminación gaussiana y su única diferencia con la LU completa es que solo utiliza los elementos distintos de cero de la matriz a factorizar. Las optimizaciones para el indizado de matrices sparse hacen que el costo de este algoritmo dependa de la cantidad de elementos no nulos que tenga la matriz. Un aporte del trabajo es la modificación que se realiza a este algoritmo en cuanto a la forma en que se actualizan los elementos. La implementación usual realiza el proceso por filas, o sea, en el momento en el que se pasa a procesar la fila, la fila ya está completamente lista. En la Figura 2, se representa este mecanismo. Las componentes con tonalidad más oscura son aquellas que ya tienen su resultado final. Las más claras, son las que serán modificadas utilizando el multiplicador. Nótese como no está listo para utilizarse como multiplicador hasta que no fue actualizado en el paso anterior por. El inconveniente de esta alternativa es la sincronización, pues para utilizar a, antes este tuvo que ser actualizado por, y lo mismo ocurre con. La modificación que se propone es realizar el proceso por columnas, verfigura 3. competir con otras compañías donde él mismo había trabajado o incluso creado anteriormente, y que ya había abandonado. Un ejemplo de esto es la famosa empresa informática IBM que sigue en auge actualmente, y que abandonó para hacerle la competencia. 3 Del inglés Message Passing Interface, plataforma que permite la ejecución de algoritmos en paralelo en una arquitectura de memoria distribuida. 203 4 Figura 2 Proceso de factorización ILU (0) usual Con este cambio se resuelve el problema de la sincronización, pues no depende de en ningún momento. Esta nueva implementación es equivalente a la anterior, la diferencia es que no son las filas las que se irán terminando, sino las columnas. Cuando se pase a la columna la ya estará lista. De esta forma se pueden procesar las filas paralelamente Figura 3 Factorización ILU (0) modificada En la Figura 4 se puede apreciar en qué consiste el proceso en paralelo. La matriz se divide en bloques filas y se modifican con su multiplicador correspondiente. Nótese cómo la actualización de las filas utilizando los multiplicadores que están en la columna es independiente, por tanto no se tiene el problema de sincronización presente en el diseño usual. Otra observación importante es que no se pasa a la columna hasta que todos los procesadores terminaron su tarea en la columna, es decir, hasta que todas las filas fueron actualizadas con su multiplicador correspondiente. En este punto es donde se tiene en cuenta la comunicación entre los procesadores, cada uno comunica sus modificaciones y se redistribuyen las filas.. Figura 4 Factorización ILU (0) en paralelo 204 5 Sustitución hacia adelante y retrógrada En la aplicación del precondicionador el análisis es muy similar, ver (Figura 5). Figura 5 Sustitución hacia adelante y retrógrada usual La modificación que fue introducida aquí es, por tanto, muy parecida a la que se realizó para la obtención del precondicionador. La actualización de las se hizo por las columnas, es decir, en el paso se calcula la componente y se modifican, para la calculada, las componentes, con para la sustitución hacia adelante y para la retrógrada, ver Figura 6. De esta forma en el paso se actualizan todas las componentes que dependen de la calculada simultáneamente. Figura 6 Sustitución hacia adelante y retrógrada modificada En la Figura 7se muestran dos pasos en paralelo de la sustitución hacia adelante y retrógrada respectivamente. En ambos casos el paralelismo se aplica en la actualización de las. No se procesa la columna hasta que las actualizaciones de las correspondientes a la columna no hayan concluido evitando cualquier problema de sincronización. Figura 7 Sustitución hacia adelante y retrógrada en paralelo 205 6 3. W_ILU_GMRES EN PARALELO El diseño en paralelo se desarrolló sobre la tecnología multi-core, el tipo de particionado utilizado fue Descomposición del Dominio (Domain Decomposition), el cual está centrado en los datos. Cada procesador trabaja, en paralelo, con una porción de los mismos. Como se explica en (L. Villarín, A. León, M. L. Baguer, Y. Linares, 2012), dado un sistema de ecuaciones lineales, aplicar W_ILU_GMRES significa: a) Representar en una base Wavelet. b) Escoger un umbral de truncamiento respecto al cual se obtiene una aproximación dispersa. c) Obtener el precondicionador ILU (0). d) Resolver el sistema disperso mediante GMRES precondicionado. En lo que sigue se explicará el diseño en paralelo de cada una de estas etapas Transformada Wavelet Bidimensional en paralelo El diseño consiste en dividir la matriz en una cierta cantidad de bloques, dependiente del número de procesadores disponibles. De esta manera cada procesador se encargará de aplicar la transformada a su bloque correspondiente. Teniendo en cuenta que se utilizó la transformada estándar, el algoritmo primero divide a la matriz en bloques por filas, cada procesador transforma la fila que le corresponde hasta el nivel deseado y luego pasa a la. Cuando un procesador termina con una fila computa otra que no haya sido procesada, de esta forma un procesador se detiene solo cuando no queden filas. Una vez que se terminan de calcular todas las filas, se realiza el mismo procedimiento por las columnas de la matriz resultante, ver (Figura 8 y Figura 9).. Figura 8 Wavelet por filas, en paralelo Cuando se espera a que todos los procesadores terminen con sus bloques filas, para comenzar con las columnas, se está teniendo en cuenta la sincronización y comunicación entre estos, pues aquí es donde se manifiesta una dependencia de las tareas y de los datos. Figura 9 Wavelet por columnas, a la matriz resultante, en paralelo 206 7 3.2. Selección del umbral de truncamiento La matriz se convierte en un vector de componentes que se ordena para obtener el menor que constituirá el umbral a partir del cual se obtendrá la aproximación dispersa de la matriz densa original. Una descripción detallada de la estrategia, se puede ver en (L. Villarín, A. León, M. L. Baguer, Y. Linares, 2012). El proceso de conversión a vector es muy fácil de paralelizar (Figura 10). Figura 10 Proceso de conversión de matriz a vector, en paralelo Cada procesador se encarga de sus filas correspondientes y una vez que todos terminaron, se tiene el vector con todos los datos. Se realizó una implementación en paralelo del algoritmo QuickSort, véase (T. H. Cormen, C. E. Rivest, R. L. Leiserson, 2009), pues es, por sus características, uno de los que más aprovecha el paralelismo debido fundamentalmente a su estrategia de divide y vencerás. La idea de este algoritmo es dividir el vector (o porción del vector) a ordenar en dos partes (no necesariamente de igual tamaño) donde cualquier elemento de la primera es menor que cualquier otro de la segunda. Esta división se realiza usando un valor del mismo vector a ordenar (elemento pivote). Una vez realizada esta operación, se le aplica el algoritmo recursivamente a cada una de las dos partes obtenidas. El diseño en paralelo que se propone aprovecha la implementación secuencial, pues el paralelismo se aplica a cada uno de los llamados recursivos. Un detalle interesante es que con un vector de tamaño relativamente pequeño ( ) es más eficiente, en cuanto a tiempo de ejecución, utilizar el algoritmo secuencial, pues resulta más costoso crear dos nuevos hilos de ejecución, que ordenar el segmento de datos de forma secuencial GMRES en paralelo GMRES es un algoritmo iterativo de naturaleza secuencial, véase (Y. Saad, 2000), luego aquí solo se pueden paralelizar las operaciones que se realizan en una iteración. Las operaciones de mayor peso son el proceso de Arnoldi-Householder y la resolución del sistema. Arnoldi-Householder se basa en operaciones x y x Ambas constituyen ejemplos clásicos de métodos en paralelo. En la multiplicación x y x secuencial el cálculo de cada una de las componentes es independiente del resto, por tanto la implementación en paralelo de estas operaciones se realiza directamente. Tanto la matriz como el vector columna se dividen en bloques filas, (Figura 11 y Figura 12) y cada procesador se encarga de calcular su parte totalmente independiente de los otros. Esta idea es bastante similar a la encontrada en (G. Golub, C. F. Van Loan, 1996), solo que Golub tuvo en cuenta un conjunto de factores que hoy, con la evolución de la tecnología, son transparentes al programador. Nótese que en el caso de x para calcular la componente se requiere de la fila y de todas las componentes del vector por el que se multiplica y en el caso de x para calcular la fila se necesita de la componente del vector columna. Esto pudiera suponer un problema de acceso a memoria, pero como estos valores no son modificados en ningún momento, los procesadores pueden acceder simultáneamente sin provocar ningún error en el resultado. 2078 Para la resolución del sistema, con matriz triangular superior, basta hacer la sustitución retrógrada. Figura 11 Multiplicación matriz x vector en paralelo Figura 12 Multiplicación vector x vector en paralelo 3.4. Estructura Sparse En el momento en el que se obtiene una aproximación dispersa de la matriz densa, se hace necesario utilizar alguna estructura que almacene los elementos distintos de cero y permita implementar las operaciones básicas aprovechando la característica dispersa. Existen distintas estructuras de datos que resuelven este problema, tales como CRS (Compressed Row Storage) y JDS (Jagged Diagonal Storage), (S. Hossian, 2006)y(G. Haase, M. Liebmann, G. Plank, 2007). Todos estos patrones de almacenamiento han sido diseñados para ajustarse a las características de problemas específicos. Cualquiera de ellos pudo haber sido utilizado, pero cuando se tiene en cuenta que es preciso realizar eficientemente: el recorrido tanto por filas, como por columnas, la multiplicación x y x tanto secuencial como en paralelo y las operaciones relacionadas con el precondicionador ILU (0), resultó más conveniente diseñar una estructura propia afín con esta estrategia. El diseño que se propone se basa en listas doblemente enlazadas, constituidas por nodos con los campos: 2089 , donde se almacenará el valor del elemento en la matriz que este representa., indica la fila correspondiente a., corresponde a la columna de., referencian a los elementos correspondientes a la izquierda, derecha, arriba y abajo respectivamente. Se tienen dos arreglos unidimensionales ( : izquierda-derecha y : arriba-abajo) de tamaño (dimensión de la matriz). tiene en la posición una referencia al primero y último nodo correspondiente a la columna, mientras que tiene en la posición una referencia al primero y último nodo correspondiente a la fila (Figura 13). Por otra parte se implementaron enumeradores que garantizan el recorrido tanto por las filas como por las columnas. Las operaciones de x y x se realizan de manera muy similar al procedimiento a seguir con las matrices densas. Figura 13 Estructura de datos para la matriz dispersa Para optimizar las operaciones en ILU (0) se implementó otro enumerador que itera a partir del primer elemento que constituya un multiplicador en una columna, pues en el procedimiento de ILU (0) se tendrían que recorrer elementos innecesarios (Y. Saad, 2000), véase la Figura 13. Como solución a este problema se implementaron dos nuevos enumeradores uno que itera por todos los multiplicadores y otro que itera por la fila a partir de cierta columna. De esta manera el diseño cumple con todos los requisitos mencionados anteriormente. 4. EXPERIMENTACIÓN NUMÉRICA Y ANÁLISIS DE RESULTADOS Las experimentaciones están orientadas a medir la escalabilidad del diseño en paralelo, ya que es así como se evalúa la eficiencia del programa. Se realizaron cuatro experimentos fundamentales. En el primero y el segundo se analizan los algoritmos que intervienen en la obtención de la aproximación dispersa y del precondicionador ILU (0) respectivamente. El tercer experimento está dedicado al estudio del comportamiento de GMRES en paralelo y por último se realiza un análisis de la estrategia en su conjunto. Se utilizó la matriz de Lehmer, que es simétrica y definida positiva, y haciendo uso de su definición implícita se puede variar la dimensión tanto como se desee Comportamiento de la obtención de la aproximación dispersa Este primer experimento se diseñó con el objetivo de analizar el comportamiento de la velocidad de ejecución en la obtención de la aproximación dispersa de la matriz original (aplicación de la transformada Wavelet bidimensional estándar y selección del umbral de truncamiento), para uno, dos y cuatro 20910 procesadores haciendo variar la dimensión de la matriz y el nivel de resolución de la transformada Wavelet. La Figura 14 y la Figura 15 muestran el buen aprovechamiento del paralelismo. Nótese cómo el speedup se acerca bastante al número de procesadores, lo que indica la escalabilidad del diseño. Figura 14 Speedup variando la dimensión de la matriz Figura 15 Speedup variando el nivel de resolución La conclusión fundamental es que se puede aumentar la dimensión de la matriz y/o el nivel de resolución y el speedup se mantendrá constante. Esto induce a pensar que incrementar la dimensión no afecta la velocidad de ejecución. Por otra parte, según la teoría llegará un momento en que aumentar la cantidad de procesadores no supondrá un aumento de la velocidad. Asumiendo que la porción del programa que realmente se ejecuta en paralelo es (P 0.98) y aplicando las fórmulas correspondientes que aparecen en (I. Foster, 1995) a este ejemplo, se obtiene la Tabla 1 donde se aprecia la dependencia de la velocidad de ejecución con respecto al número de procesadores. Tabla 1Análisis del número de procesadores contra parte en paralelo N P = Este análisis no tiene en cuenta los costos adicionales, como la comunicación entre procesadores y por tanto su interés es esencialmente teórico. Tanto es así que el resultado práctico obtenido para dos procesadores es de aproximadamente 1.82, véase la Figura 15. Lo que sí se puede afirmar y comprobar en la práctica (disponiendo de la tecnología adecuada) es que un aumento de la cantidad de procesadores no siempre redunda en una disminución significativa del tiempo de ejecución. 21011 4.2. Comportamiento de ILU (0) En esta segunda parte de la experimentación se hicieron dos pruebas, una haciendo variar la dimensión de la matriz y otra cambiando la cantidad de elementos no nulos. Figura 16 Speedup variando la dimensión en ILU (0) Figura 17 Speedup variando el umbral en ILU (0) Esta etapa constituye el nodo crítico de la estrategia, tanto para el diseño secuencial como para el paralelo. Ambos gráficos así lo demuestran. En la implementación en paralelo el problema radica en que cuando la matriz es muy dispersa, contiene muy pocos elementos por fila, además de que, como se vio, en cada fila no se procesan todas las componentes. Esto trae como consecuencia que se queden procesadores ociosos y exista una diferencia notable entre el momento en que termina el primero y el último. Nótese cómo en lafigura 16 al aumentar la dimensión de la matriz y mantener fijo el umbral esta se fue haciendo más dispersa, por lo que se afectó la escalabilidad, y en la Figura 17 al aumentar la cantidad de elementos no nulos se aumenta la escalabilidad. Este aspecto conlleva a una conclusión interesante. En el caso de la implementación secuencial, se tiene que a mayor cantidad de elementos nulos, se obtiene un menor costo de ILU (0). Como ILU (0) es quien domina el costo total de la estrategia, esto obliga a que se necesite una matriz bien dispersa. Por otra parte no se puede utilizar una aproximación con muy pocos elementos distintos de cero porque la probabilidad de que la matriz sea singular es muy alta, entonces hay que buscar un umbral que equilibre el costo de ILU (0) y la condición de no singularidad. Con la implementación en paralelo es cierto que el costo se mantiene. Sin embargo, resultó que mientras menos dispersa es la matriz el diseño es más escalable, es decir, aprovecha mejor los procesadores disponibles. Esto brinda la posibilidad de aumentar un poco la cantidad de elementos no nulos y por ende encontrar una mejor aproximación de la matriz Comportamiento de GMRES Figura 18 GMRES variando la dimensión de la matriz 21112 Con este experimento se pretende mostrar la eficiencia de los algoritmos que se encargan de efectuar las operaciones de multiplicación x y x en paralelo. Para ello se varió la dimensión de la matriz y se analizó si un aumento de la dimensión influía en la escalabilidad del algoritmo. En la Figura 18 se observa que el speedup se mantiene constante. Este resultado es el esperado, pues las operaciones paralelizadas son los algoritmos clásicos que sacan gran provecho del paralelismo, y prácticamente todas las operaciones que se hacen en una iteración se basan en estos Comportamiento de la estrategia general Para analizar la estrategia en su conjunto se corrieron todas las partes integradas variando la dimensión de la matriz y el umbral de truncamiento. Figura 19 Variación de la dimensión en W_ILU_GMRES Figura 20 Variación del umbral en W_ILU_GMRES En las figuras se observa que la obtención del precondicionador ILU (0) no afecta considerablemente el speedup de la estrategia general, pues aumentar la dimensión de la matriz o el umbral de truncamiento no supuso una variación de la escalabilidad, lo que mantiene el número teórico de procesadores hasta el cual se obtiene una mejora considerable de la velocidad de ejecución. Figura 21 Estrategia general 21213 5. CONCLUSIONES En el presente trabajo se propuso una implementación de W_ILU_GMRES en paralelo, que supera desde el punto de vista computacional la implementación secuencial presentada en (L. Villarín, A. León, M. L. Baguer, Y. Linares, 2012). Se realizó un diseño modular y reutilizable que permite incorporar otros precondicionadores, otras bases wavelet y/o métodos iterativos. La implementación en paralelo de la transformada Wavelet, que permitió su aplicación a cada fila y columna de manera totalmente independiente, y la elección de Quicksort para la selección del umbral de truncamiento propiciaron que para dos y cuatro procesadores se alcanzara un speedup de aproximadamente 1.7 y 3.8 respectivamente. En el análisis del precondicionador ILU (0) los resultados ilustran un mejor comportamiento cuando se tiene una mayor cantidad de elementos no nulos en la aproximación dispersa. Se mostró que ILU (0) continúa siendo el nodo crítico de la estrategia, pues afecta la escalabilidad. Con dos procesadores para un 10%, 25% y 40% de elementos no nulos se obtienen speedup entre 1.3 y 1.4 y utilizando cuatro procesadores para los mismos porcientos de dispersión el speedup toma valores entre 2.2 y 2.3 aproximadamente. No obstante, se logró minimizar el problema existente en la implementación secuencial donde para disminuir el costo de ILU (0) se necesita una matriz bien dispersa, lo que implica una mayor probabilidad de que la matriz sea singular. En la evaluación de GMRES en paralelo se obtuvo un speedup de aproximadamente 1.7 y 3.7 para dos y cuatro procesadores respectivamente, lo cual es un resultado alentador. En la estrategia general W_ILU_GMRES en paralelo se obtuvo para dos y cuatro procesadores un speedup de aproximadamente 1.6 y 3.5 respectivamente, lo cual indica que la influencia de la obtención del precondicionador ILU (0) no fue trascendental. Los resultados muestran la tendencia a una disminución del tiempo de ejecución con el aumento del número de procesadores, por lo cual sería recomendable experimentar con un incremento de los mismos. Además queda abierta la incorporación de nuevas técnicas de precondicionamiento y la extensión del diseño a una aplicación MPI para contar con una herramienta multiplataforma. RECEIVED FEBRUARY, 2013 REVISED APRIL, 2014 REFERENCIAS 1 BARNEY,.(2009): Introduction to Parallel Computing [En línea]. - Enero de Marzo de https://computing.llnl.gov/tutorials/parallel_comp/. 2 BENZI, M.,W. JOUBERT and G. MATEESCU(1999): Numerical Experiments with parallel orderings for ILU preconditioners. Electronic Transactions on Numerical Analysis, 8, CORMEN, T. H., C. E. RIVEST and R. L. LEISERSON(2009): Introducion to Algorithm. MIT Press, Cambridge. 4 CUNHA, R.D.D: and T. HOPKINS(1994): A parallel implementation of the restarted GMRES iterative method for nonsymmetric systems of linear equations. Advances in Computational Mathematics. - Canterbury 2, ERHEL, J.(1995): A parallel GMRES version for General Sparse Matrices.Electronic Transactions on Numerical Analysis. 3, FOSTER, I. (1995): Designing and Building Parallel Programs: Concepts and Tools for Parallel Software Engineering. Addison-Wesley, Boston. 7 GOLUB, G. and C. F. VAN LOAN(1996): Parallel Matrix Computations,Third Edition. : Johns Hopkins University Press, Maryland. 8 HAASE, G., M. LIEBMANN and G. PLANK(2007): A Hilbert-order multiplication scheme for unstructured sparse matrix. International Journal of Parallel, Emergent and Distributed Systems, 22, HARROP, J. (2008): F# for scientist.wiley-interscience, Chichester. 10 HOSSIAN, S. (2006): On Efficient storage of Sparse Matrix. The Proceedings of ICCSE, 2005, MICROSOFT DEVELOPERS(2010): What's New in Visual F# Microsoft Developer Network. - MSDN, msdn.microsoft.com/en-us/library/vstudio/dd PICKERING, R. (2007): Foundations of F#. Apress, California. 13 SAAD, Y. (2000): Iterative Methods for Sparse Linear Systems. (2 nd edition), SIAM, Philadelphia. 21314 14 SYME, D.(2007): Expert F#. Apress, California. 15 VILLARÍN, L., A. LEÓN, M. L. BAGUER, y Y. LINARES (2012): GMRES precondicionado con Wavelets. Un algoritmo de selección del umbral para la obtención del patrón de dispersión.i Revista Investigación Operacional, 33, Mostrar más
Gradiente conjugado MSc. Miguel Vargas-Félix miguelvargas@cimat.mx http://www.cimat.mx/~miguelvargas 07/09/11 1/23 Contenido Contenido Estructura de una matriz rala a partir de conectividades Compressed Más detalles CLUSTER FING: PARALELISMO de MEMORIA DISTRIBUIDA
CLUSTER FING: PARALELISMO de MEMORIA DISTRIBUIDA SERGIO NESMACHNOW Centro de Cálculo, Instituto de Computación FACULTAD DE INGENIERÍA, UNIVERSIDAD DE LA REPÚBLICA, URUGUAY CONTENIDO Introducción: arquitecturas Más detalles Tópicos de implementación en paralelo
Apéndice C Esquema PGP. Implementación en un entorno paralelo. Para poder abordar eficazmente la simulación computacional de grandes problemas estructurales, es necesario utilizar estrategias numéricas Más detalles 4. Programación Paralela
4. Programación Paralela La necesidad que surge para resolver problemas que requieren tiempo elevado de cómputo origina lo que hoy se conoce como computación paralela. Mediante el uso concurrente de varios Más detalles Capítulo 5: Comparación con ScaLAPACK
Capítulo 5: Comparación con ScaLAPACK En este capítulo se presentan dos aspectos importantes en cuanto a la validez y utilización de los aportes de esta tesis: 1) Aplicación de los principios de paralelización Más detalles Desarrollo de un cluster computacional para la compilación de. algoritmos en paralelo en el Observatorio Astronómico.
Desarrollo de un cluster computacional para la compilación de algoritmos en paralelo en el Observatorio Astronómico. John Jairo Parra Pérez Resumen Este artículo muestra cómo funciona la supercomputación Más detalles T E C N O L O G Í A OPTIMIZACIÓN DE MATERIALES MEDIANTE PATRONES DE CORTE EFICIENTE. Aplicación. a la INDUSTRIA
OPTIMIZACIÓN DE MATERIALES MEDIANTE PATRONES DE CORTE EFICIENTE Aplicación a la INDUSTRIA de la construcción 1 El presente estudio propone el uso de un algoritmo comúnmente utilizado en la rama de investigación Más detalles SISTEMAS DE ECUACIONES LINEALES
SISTEMAS DE ECUACIONES LINEALES INTRODUCCIÓN En el presente documento se explican detalladamente dos importantes temas: 1. Descomposición LU. 2. Método de Gauss-Seidel. Se trata de dos importantes herramientas Más detalles Módulo: Modelos de programación para Big Data
Program. paralela/distribuida Módulo: Modelos de programación para Big Data (título original: Entornos de programación paralela basados en modelos/paradigmas) Fernando Pérez Costoya Introducción Big Data Más detalles Novedades de software/new softwares
REVISTA INVESTIGACIÓN OPERACIONAL VOL., 3, No. 3, 275-28, 2 Novedades de software/new softwares ALGORITMO PARA LA GENERACIÓN ALEATORIA DE MATRICES BOOLEANAS INVERSIBLES P. Freyre*, N. Díaz*, E. R. Morgado** Más detalles Algoritmo para resolver exactamente sistemas de ecuaciones lineales con coeficientes enteros
Miscelánea Matemática 43 (2006) 7 132 SMM Algoritmo para resolver exactamente sistemas de ecuaciones lineales con coeficientes enteros Daniel Gómez-García Facultad de Ingeniería Universidad Autónoma de Más detalles construcción de programas Prof. Eliana Guzmán U.
Unidad II. Metodología para la construcción de programas Prof. Eliana Guzmán U. Semestre: A-2015 Introducción Resolver un problema con una computadora conduce a la escritura de un programa y a su ejecución. Más detalles Prácticas de programación en C con MinGW Developer Studio
Prácticas de programación en C con MinGW Developer Studio MinGW Developer Studio es un entorno de desarrollo integrado (IDE) para la programación en lenguaje C gratuito y cómodo de usar. Se ejecuta en Más detalles Descomposición de dominios
Descomposición de dominios Miguel Vargas 27/10/10 1/29 Contenido Contenido Solución de ecuaciones diferenciales con descomposición de dominios Dominios sin traslape, complemento de Schur Método alternante Más detalles Trabajo Práctico Final. Transformada Rápida de Fourier con Paralelismo
Trabajo Práctico Final Transformada Rápida de Fourier con Paralelismo Sistemas Operativos II Facultad de Ciencias Exactas, Físicas y aturales UC Mauricio G. Jost 2009 Índice 1. Introducción 1 2. Marco Más detalles Facultad de Ingeniería ISSN: 0121-1129 revista.ingenieria@uptc.edu.co. Universidad Pedagógica y Tecnológica de Colombia. Colombia
Facultad de Ingeniería ISSN: 0121-1129 revista.ingenieria@uptc.edu.co Universidad Pedagógica y Tecnológica de Colombia Colombia Amézquita-Mesa, Diego Germán; Amézquita-Becerra, Germán; Galindo-Parra, Omaira Más detalles COMPARACIÓN DEL RENDIMIENTO COMPUTACIONAL ENTRE DIFERENTES METODOLOGÍAS DE PROCESAMIENTO EN PARALELO PARA FEA VÍA ANSYS 14.5
Second International Conference on Advanced Mechatronics, Design, and Manufacturing Technology - AMDM 2014 1 COMPARACIÓN DEL RENDIMIENTO COMPUTACIONAL ENTRE DIFERENTES METODOLOGÍAS DE PROCESAMIENTO EN Más detalles MineSight y la Tecnología Informática de última Generación
MineSight y la Tecnología Informática de última Generación Así como se ha modificado la tecnología en el escritorio y en las estaciones de trabajo de los últimos años, del mismo modo ha ido variando el Más detalles UNIVERSIDAD AUTÓNOMA DE YUCATÁN FACULTAD DE MATEMÁTICAS MISIÓN
UNIVERSIDAD AUTÓNOMA DE YUCATÁN FACULTAD DE MATEMÁTICAS MISIÓN Formar profesionales altamente capacitados, desarrollar investigación y realizar actividades de extensión, en Matemáticas y Computación, así Más detalles TEMA 4: CALCULO NUMERICO DE AUTOVALORES
Lino Alvarez - Aurea Martinez METODOS NUMERICOS TEMA 4: CALCULO NUMERICO DE AUTOVALORES 1 INTRODUCCION La determinación de autovalores y autovectores de una matriz cuadrada A de orden n es un problema Más detalles GPU IMPLEMENTATIONS OF SCHEDULING HEURISTICS FOR HETEROGENEOUS COMPUTING ENVIRONMENTS
GPU IMPLEMENTATIONS OF SCHEDULING HEURISTICS FOR HETEROGENEOUS COMPUTING ENVIRONMENTS MAURO CANABÉ SERGIO NESMACHNOW Centro de Cálculo, Facultad de Ingeniería Universidad de la República, Uruguay GPU IMPLEMENTATIONS Más detalles Cómo ordenar una lista de números?
Cómo ordenar una lista de números? Germán Ariel Torres Resumen. Este trabajo trata acerca de métodos y técnicas usadas para el ordenamiento eficiente de listas de números. Se analizan los métodos básicos, Más detalles Resumen sobre Grafcet RESUMEN SOBRE GRAFCET
RESUMEN SOBRE GRAFCET INDICE GENERAL 1 INTRODUCCIÓN 1 2 METODOLOGÍA GRAFCET: CONCEPTOS BÁSICOS 2 2.1 ETAPAS 2 2.1.1 ACCIONES ASOCIADAS A LAS ETAPAS 4 2.2 TRANSICIÓN 5 2.2.1 RECEPTIVIDAD ASOCIADA A LA TRANSICIÓN Más detalles ALN - Formatos dispersos
ALN - Formatos dispersos In. Co. Facultad de Ingeniería Universidad de la República Temario Matrices dispersas Motivación e historia Formatos estáticos Formatos dinámicos Otros formatos Versión 1.0 2 Motivación Más detalles Paralelización de la factorización LDL T usando el lenguaje de programación paralela ZPL
REVISTA INGENIERÍA UC. Vol. 15, N o 2, 72-80, 2008 Paralelización de la factorización LDL T usando el lenguaje de programación paralela ZPL Edwin Vargas, Enrique Flores, Demetrio Rey Lago Instituto de Más detalles Programación Distribuida
PROGRAMA DE ESTUDIO Programa Educativo: Área de Formación : Licenciatura en Ciencias Computacionales Integral profesional Programa elaborado por: Programación Distribuida Horas teóricas: 1 Horas prácticas: Más detalles Implementación de la transformada wavelet discreta para imágenes en un FPGA
Implementación de la transformada wavelet discreta para imágenes en un FPGA Madeleine León 1, Carlos A. Murgas 1, Lorena Vargas 2, Leiner Barba 2, Cesar Torres 2 1 Estudiantes de pregrado de la Universidad Más detalles Ejemplos de conversión de reales a enteros
Ejemplos de conversión de reales a enteros Con el siguiente programa se pueden apreciar las diferencias entre las cuatro funciones para convertir de reales a enteros: program convertir_real_a_entero print Más detalles CAPITULO 1. Introducción a los Conceptos Generales de Bases de Datos Distribuidas
CAPITULO 1 Introducción a los Conceptos Generales de 1.1 Preliminares Las empresas necesitan almacenar información. La información puede ser de todo tipo. Cada elemento informativo es lo que se conoce Más detalles Paralelización de métodos iterativos. Parallelization of iterative methods
Paralelización de métodos iterativos Germán Larrazábal Departamento de Computación Facultad de Ciencias y Tecnología (FACYT), Universidad de Carabobo Avenida Montes de Oca, No. 120-267, Edificio FACYT, Más detalles Resolución de problemas en paralelo
Resolución de problemas en paralelo Algoritmos Paralelos Tema 1. Introducción a la computación paralela (segunda parte) Vicente Cerverón Universitat de València Resolución de problemas en paralelo Descomposición Más detalles Programación híbrida en arquitecturas cluster de multicore. Escalabilidad y comparación con memoria compartida y pasaje de mensajes.
Programación híbrida en arquitecturas cluster de multicore. Escalabilidad y comparación con memoria compartida y pasaje de mensajes. Fabiana Leibovich, Armando De Giusti, Marcelo Naiouf, Laura De Giusti, Más detalles MATRICES PRODUCTO DE MATRICES POTENCIAS NATURALES DE MATRICES CUADRADAS
Tema 1.- MATRICES MATRICES PRODUCTO DE MATRICES POTENCIAS NATURALES DE MATRICES CUADRADAS Fundamentos Matemáticos de la Ingeniería 1 Un poco de historia Lord Cayley es uno de los fundadores de la teoría Más detalles Factorizaciones de Cholesky, matrices definidas. semidefinidas positivas.
Factorizaciones de Cholesky, matrices definidas y semidefinidas positivas Héctor Manuel Mora Escobar Universidad Central, Bogotá hectormora@yahoo.com Junio de 2011 1 Introducción Este documento presenta, Más detalles Las matrices tienen un número cada vez mas creciente de aplicaciones en la solución de problemas en Ciencia y Tecnología.
Aplicaciones de las Matrices a la Solución de Problemas de Redes Eléctricas Resumen Se muestra como obtener, sistemas de ecuaciones lineales que permitan calcular intensidades de corrientes en los ramales Más detalles SEDO: SOFTWARE EDUCATIVO DE MATEMÁTICA NUMÉRICA. Lic. Maikel León Espinosa. mle@uclv.edu.cu
EDU082 Resumen SEDO: SOFTWARE EDUCATIVO DE MATEMÁTICA NUMÉRICA Lic. Maikel León Espinosa mle@uclv.edu.cu Departamento de Ciencia de la Computación Universidad Central Marta Abreu de Las Villas Carretera Más detalles Computación paralela y distribuida
Computación paralela y distribuida MSC-0302 Nombre de la asignatura: Computación paralela y distribuida Línea de trabajo: Tecnologías web Tiempo de dedicación del estudiante a las actividades de: DOC TIS Más detalles PFC- Aplicaciones Web para trabajo colaborativo:
PFC- Aplicaciones Web para trabajo colaborativo: Aplicación para Control de una Integración de S.I. 2º Ciclo Ingeniería Informática Curso 2011-2012 Consultor : Fatos Xhafa Autor : Miguel Angel Pineda Cruz Más detalles DISEÑO DE UN ALGORITMO PARALELO PARA CODIFICACIÓN DE VIDEO MPEG4, SOBRE UN CLUSTER DE COMPUTADORAS PERSONALES.
DISEÑO DE UN ALGORITMO PARALELO PARA CODIFICACIÓN DE VIDEO MPEG4, SOBRE UN CLUSTER DE COMPUTADORAS PERSONALES. Autores: Abelardo Rodríguez L.*, Manuel Pérez M. **, Alberto González T. **, José Hernández Más detalles ALGEBRA LINEAL. Héctor Jairo Martínez R. Ana María Sanabria R.
ALGEBRA LINEAL Héctor Jairo Martínez R. Ana María Sanabria R. SEGUNDO SEMESTRE 8 Índice general. SISTEMAS DE ECUACIONES LINEALES.. Introducción................................................ Conceptos Más detalles LECCIÓN 8: CIRCUITOS Y ALGORITMOS DE MULTIPLICACIÓN DE ENTEROS
ESTRUCTURA DE COMPUTADORES Pag. 8.1 LECCIÓN 8: CIRCUITOS Y ALGORITMOS DE MULTIPLICACIÓN DE ENTEROS 1. Circuitos de multiplicación La operación de multiplicar es mas compleja que la suma y por tanto se Más detalles Sistemas de Numeración
UNIDAD Sistemas de Numeración Introducción a la unidad Para la mayoría de nosotros el sistema numérico base 0 aparentemente es algo natural, sin embargo si se establecen reglas de construcción basadas Más detalles High Performance Computing and Architectures Group
HPCA Group 1 High Performance Computing and Architectures Group http://www.hpca.uji.es Universidad Jaime I de Castellón ANACAP, noviembre de 2008 HPCA Group 2 Generalidades Creado en 1991, al mismo tiempo Más detalles CAPÍTULO 1 Instrumentación Virtual
CAPÍTULO 1 Instrumentación Virtual 1.1 Qué es Instrumentación Virtual? En las últimas décadas se han incrementado de manera considerable las aplicaciones que corren a través de redes debido al surgimiento Más detalles PROGRAMACIÓN BÁSICA DE LA COMPUTADORA. 1 Introducción. Tabla 1: Instrucciones MIPS
PROGRAMACIÓN BÁSICA DE LA COMPUTADORA 1 Introducción Un sistema de computadora total incluye tanto circuitería (hardware) como programación (software). El hardware consta de los componentes físicos y todo Más detalles SERVICIO RECOMENDADOR DE RECURSOS DIGITALES PARA EL SISTEMA DSPACE
DESARROLLO DE UN SERVICIO RECOMENDADOR DE RECURSOS DIGITALES PARA EL SISTEMA DSPACE 1. INTRODUCCIÓN Desde hace dos décadas, las bibliotecas han efectuado con éxito una serie de transformaciones en su filosofía, Más detalles Problemas Resueltos de Desigualdades y Programación Lineal
Universidad de Sonora División de Ciencias Exactas y Naturales Departamento de Matemáticas. Problemas Resueltos de Desigualdades y Programación Lineal Para el curso de Cálculo Diferencial de Químico Biólogo Más detalles GUÍA DOCENTE. Curso 2014-2015 1. DESCRIPCIÓN DE LA ASIGNATURA. Ingeniería Informática en Sistemas de Información Doble Grado: Módulo: Módulo 6
1. DESCRIPCIÓN DE LA ASIGNATURA Grado: Ingeniería Informática en Sistemas de Información Doble Grado: Asignatura: Ingeniería del Sotware II Módulo: Módulo 6 Departamento: Deporte e Informática Año académico: Más detalles Métodos Iterativos para Resolver Sistemas Lineales
Métodos Iterativos para Resolver Sistemas Lineales Departamento de Matemáticas, CCIR/ITESM 17 de julio de 2009 Índice 3.1. Introducción............................................... 1 3.2. Objetivos................................................ Más detalles Universidad Politécnica de Madrid
Universidad Politécnica de Madrid Escuela Técnica Superior de Ingenieros de Minas Departamento de Matemática Aplicada y Métodos Informáticos x r 1 x 1 x 2 = x * d 1 * d 1 f(x)=cte. r Resolución de sistemas Más detalles Informática y Programación Escuela de Ingenierías Industriales y Civiles Grado en Ingeniería en Ingeniería Química Curso 2010/2011
Módulo 1. Fundamentos de Computadores Informática y Programación Escuela de Ingenierías Industriales y Civiles Grado en Ingeniería en Ingeniería Química Curso 2010/2011 1 CONTENIDO Tema 1. Introducción Más detalles 1. INVERSA DE UNA MATRIZ REGULAR
. INVERSA DE UNA MATRIZ REGULAR Calcular la inversa de una matriz regular es un trabajo bastante tedioso. A través de ejemplos se expondrán diferentes técnicas para calcular la matriz inversa de una matriz Más detalles Computación Científica en Paralelo
Computación Científica en Paralelo Luis Miguel de la Cruz luiggix@gmail.com www.dci.dgsca.unam.mx/lmcs Unidad de Investigación en Cómputo Aplicado DGSCA-UNAM. Posgrado en Ciencia e Ingeniería de la Computación Más detalles Tema 3. Espacios vectoriales
Tema 3. Espacios vectoriales Estructura del tema. Definición y propiedades. Ejemplos. Dependencia e independencia lineal. Conceptos de base y dimensión. Coordenadas Subespacios vectoriales. 0.1. Definición Más detalles Fundamentos Matemáticos II Electrónicos Curso 2006-07 2
Tema 2.- MATRICES!ESPACIO VECTORIAL!PRODUCTO DE MATRICES!POTENCIAS NATURALES DE MATRICES CUADRADAS Fundamentos Matemáticos II Electrónicos Curso 2006-07 1 Un poco de historia Lord Cayley es uno de los Más detalles ETSIINGENIO 2009 DIBUJO DE GRAFOS MEDIANTE ALGORITMOS GENÉTICOS
ETSIINGENIO 2009 DIBUJO DE GRAFOS MEDIANTE ALGORITMOS GENÉTICOS EtsiIngenio Inteligencia Artificial 1 Raposo López Alejandro Sánchez Palacios Manuel Resumen dibujo de grafos mediante algoritmos genéticos Más detalles Computación Científica en Paralelo
Computación Científica en Paralelo Métodos de Descomposición de Dominio Luis Miguel de la Cruz luiggix@gmail.com www.dci.dgsca.unam.mx/lmcs Unidad de Investigación en Cómputo Aplicado DGSCA-UNAM. Posgrado Más detalles Las cuatro operaciones. En la. Escuela Básica. por. Francisco Rivero Mendoza
Las cuatro operaciones En la Escuela Básica por Francisco Rivero Mendoza 1 Conociendo los números Antes de pasar a estudiar los correspondientes algoritmos de la suma y la resta, es preciso desarrollar Más detalles PROGRAMACION VECTORIAL RAFAEL VALDÉS VALDAZO UO196581 ÁNGEL MARÍA VILABOA PÉREZ UO197092 BLOQUE PARALELAS 4º INFORMÁTICA UNIVERSIDAD DE OVIEDO
PROGRAMACION VECTORIAL RAFAEL VALDÉS VALDAZO UO196581 ÁNGEL MARÍA VILABOA PÉREZ UO197092 BLOQUE PARALELAS 4º INFORMÁTICA UNIVERSIDAD DE OVIEDO INTRODUCCIÓN Ligada al concepto de vector Decodifica instrucciones Más detalles BALANCEO DE CARGA ACADÉMICA EN EL DISEÑO DE UN CURRÍCULUM BASADO EN COMPETENCIAS
BALANCEO DE CARGA ACADÉMICA EN EL DISEÑO DE UN CURRÍCULUM BASADO EN COMPETENCIAS Martín G. Solar Monsalves Facultad de Ciencias de la Ingeniería Universidad Austral de Chile, Valdivia, Chile. Casilla 567, Más detalles Matrices. Definiciones básicas de matrices. www.math.com.mx. José de Jesús Angel Angel. jjaa@math.com.mx
Matrices Definiciones básicas de matrices wwwmathcommx José de Jesús Angel Angel jjaa@mathcommx MathCon c 2007-2008 Contenido 1 Matrices 2 11 Matrices cuadradas 3 12 Matriz transpuesta 4 13 Matriz identidad Más detalles UNIVERSIDAD ALBERT EINSTEIN FACULTAD DE INGENIERIA
UNIVERSIDAD ALBERT EINSTEIN FACULTAD DE INGENIERIA Estudio de las herramientas TOAD y DBArtisan para la administración e integración de bases de datos relacionales. PREVIA OPCION AL TÍTULO DE: INGENIERO Más detalles Uso de un motor de restricciones bajo dispositivos Android
Uso de un motor de restricciones bajo dispositivos Android Gonzalo Hernández 1, Camilo Villota Ibarra 2, James Muñoz Coronel 3, Harold Muñoz Muñoz 4 Universidad de Nariño, Facultad de Ingeniería, Departamento Más detalles Cómputo en paralelo con OpenMP 1
Cómputo en paralelo con OpenMP 1 Miguel Vargas-Félix miguelvargas@cimat.mx http://www.cimat.mx/~miguelvargas CIMAT, September 2, 2015 1/34 Velocidad de los procesadores de escritorio Recientemente, la Más detalles DESARROLLO DE COMPONENTES PARA LA INTEGRACIÓN DEL PORTAL CORPORATIVO DEL CITI CON LA BPMS BIZAGI
DESARROLLO DE COMPONENTES PARA LA INTEGRACIÓN DEL PORTAL CORPORATIVO DEL CITI CON LA BPMS BIZAGI Informe de Práctica Profesional de 4to Año, Ingeniería Informática Autor: Manuel Alejandro Aguilar Díaz Más detalles CAPÍTULO 5 ANÁLISIS DE CONVERGENCIA DEL MÉTODO BINOMIAL AL MODELO DE BLACK & SCHOLES
CAPÍTULO 5 ANÁLISIS DE CONVERGENCIA DEL MÉTODO BINOMIAL AL MODELO DE BLACK & SCHOLES Para la valuación de opciones hay dos modelos ampliamente reconocidos como son el modelo binomial y el modelo de Black Más detalles 270068 - PAP - Programación y Arquitecturas Paralelas
Unidad responsable: 270 - FIB - Facultad de Informática de Barcelona Unidad que imparte: 701 - AC - Departamento de Arquitectura de Computadores Curso: Titulación: 2014 GRADO EN INGENIERÍA INFORMÁTICA Más detalles Ingeniería de Software
Ingeniería de Software MSDN Ingeniería de Software...1 Ingeniería del Software_/_ Ingeniería y Programación...1 Análisis de Requerimientos...2 Especificación...3 Diseño...4 Desarrollo en Equipo...5 Mantenimiento...6 Más detalles Espacio afín. Transformaciones afines y movimientos
Capítulo Espacio afín. Transformaciones afines y movimientos. Espacio afín y espacio afín métrico Definición. El espacio afín (tridimensional) está constituido por los siguientes elementos. El espacio Más detalles Unidad 7 Aplicación de máximos y mínimos
Unidad 7 Aplicación de máimos y mínimos Objetivos Al terminar la unidad, el alumno: Interpretará el concepto de ingreso y costos marginal. Aplicará la función de ingresos en problemas de maimización. Aplicará Más detalles RAID. Los detalles de las características segunda y tercera, cambian según los distintos niveles RAID. RAID 0 no soporta la tercera característica.
RAID Como se dijo anteriormente, el ritmo de mejora de prestaciones en memoria secundaria ha sido considerablemente menor que en procesadores y en memoria principal. Esta desigualdad ha hecho, quizás, Más detalles 1. Información General
Programación de Sistemas Distribuidos E.T.S. Ingeniería Informática Dpto. Estadística, I.O. y Computación 1. Información General Curso: 2004/2005 Código de la asignatura: 151059113 Plan: 2003 Cuatrimestre: Más detalles Ejemplos de optimización para Kepler. 1. Balanceo dinámico de la carga. Contenidos de la charla [18 diapositivas]
Ejemplos de optimización para Kepler Manuel Ujaldón Nvidia CUDA Fellow Dpto. Arquitectura de Computadores Universidad de Málaga Contenidos de la charla [18 diapositivas] 1. Balanceo dinámico de la carga. Más detalles Energía Renovable 9REN WONDERWARE GARANTIZA EL DESEMPEÑO Y MEJORA EL RENDIMIENTO DE LAS PLANTAS DE GENERACIÓN DE ENERGÍAS RENOVABLES DE 9REN
ENERGÍA RENOVABLE Energía Renovable 9REN WONDERWARE GARANTIZA EL DESEMPEÑO Y MEJORA EL RENDIMIENTO DE LAS PLANTAS DE GENERACIÓN DE ENERGÍAS RENOVABLES DE 9REN Madrid, España. Las energías renovables son Más detalles GUÍA DOCENTE. Curso 2015-2016. Ingeniería Informática en Sistemas de Información Doble Grado: Mantenimiento de Bases de Datos Módulo: Módulo 6
1. DESCRIPCIÓN DE LA ASIGNATURA Grado: Ingeniería Informática en Sistemas de Información Doble Grado: Asignatura: Mantenimiento de Bases de Datos Módulo: Módulo 6 Departamento: Deporte e Informática Año Más detalles Métodos Computacionales 0302390
Métodos Computacionales 0302390 Fecha de actualización Tue, 28 Jul 2015 06:28:33-050 Usuario que realiza la actualización Facultad Autorización Vicedecano Última versión del curso 1 Número de Acta del Más detalles Una red es un conjunto de computadoras interconectadas entre sí con el. propósito de compartir archivos y periféricos Completando esta definición
REDES RED Una red es un conjunto de computadoras interconectadas entre sí con el propósito de compartir archivos y periféricos Completando esta definición podemos añadir que una red es un sistema de comunicaciones Más detalles Una fórmula para la pendiente
LECCIÓN CONDENSADA 5.1 Una fórmula para la pendiente En esta lección aprenderás cómo calcular la pendiente de una recta dados dos puntos de la recta determinarás si un punto se encuentra en la misma recta Más detalles DISEÑO DE COMPONENTES DE SOFTWARE *
DISEÑO DE COMPONENTES DE SOFTWARE * NOTAS DEL CURSO Ingeniería de Software I DRA. MARIA DEL PILAR GÓMEZ GIL INAOEP * Resumen del capítulo 10 de libro de [Pressman 2010] V:18-11-2008 (c) P. Gomez-Gil, INAOE. Más detalles Introducción a la Programación en MATLAB
Introducción a la Programación en MATLAB La programación en MATLAB se realiza básicamente sobre archivos M, o M-Files. Se los denomina de esta forma debido a su extensión.m. Estos archivos son simple archivos Más detalles Herramientas Informáticas I. Software: Clasificación y Funcionalidad Facultad de Ciencias Económicas y Jurídicas Universidad Nacional de La Pampa
Herramientas Informáticas I Software: Clasificación y Funcionalidad Facultad de Ciencias Económicas y Jurídicas Universidad Nacional de La Pampa 2013 Introducción La clasificación del Software permitirá Más detalles TALLER 04 SOFTWARE GESTION DE PROYECTOS Y SOFTWARE GESTIO CONTROL DE VERSIONES ROBIN TAPIERO HARRY CUESTA MARÍA SUÁREZ
TALLER 04 SOFTWARE GESTION DE PROYECTOS Y SOFTWARE GESTIO CONTROL DE VERSIONES ROBIN TAPIERO HARRY CUESTA MARÍA SUÁREZ Universidad Distrital Francisco José de Caldas Facultad Tecnológica Sistematización Más detalles Una Aplicación del Cálculo Matricial a un Problema de Ingeniería
Divulgaciones Matemáticas Vol. 9 No. 221, pp. 197 25 Una Aplicación del Cálculo Matricial a un Problema de Ingeniería An Application of Matrix Calculus to an Engineering Problem P. R. Almeida Benítez Dpto. Más detalles PROCESAMIENTO DIGITAL DE IMÁGENES MEDIANTE EL USO DE UN FPGA Y LENGUAJE VHDL
PROCESAMIENTO DIGITAL DE IMÁGENES MEDIANTE EL USO DE UN FPGA Y LENGUAJE VHDL N. E. Chávez Rodríguez*, A. M. Vázquez Vargas** *Departamento de Computación **Departamento de Procesamiento Digital de Señales Más detalles Programando el robot Lego Mindstorms NXT 2.0 con Java LeJOS
Programando el robot Lego Mindstorms NXT 2.0 con Java LeJOS Lina Ernestina Arias Hernández a, Valeria Rubí Hernández Cisneros a y Luis Orlando Aguirre Fuentes a a Instituto Tecnologico de la Laguna., elearias@hotmail.com, Más detalles Proyecto de Normalización Automática de Base de Datos
Proyecto de Normalización Automática de Base de Datos Lic. Beatriz Steimberg * Resumen En el primer cuatrimestre del año 2003 se encaró el proyecto de Normalización Automática de Base de Datos. El objetivo Más detalles Compiladores y Lenguajes de Programación. Maria de Guadalupe Cota Ortiz
Compiladores y Lenguajes de Programación Maria de Guadalupe Cota Ortiz Organizaciones que rigen las normas para estandarización de Lenguajes de Programación IEEE (Instituto de Ingenieros Eléctricos y Electrónicos) Más detalles Capacidad de procesamiento del compilador Python para el Sistema Operativo Windows y Linux Palabras Clave:
Capacidad de procesamiento del compilador Python para el Sistema Operativo Windows y Linux Stiven Unsihuay, Paulo Pereira, Norma León unsihuay_carlos@hotmail.com, paulopereiraa10@gmail.com, nleonl@usmp.pe Más detalles Tema 3 Resolución de Sistemas de Ecuaciones Lineales
Tema Resolución de Sistemas de Ecuaciones Lineales Índice Introducción 2 Método de Gauss 2 Resolución de sistemas triangulares 22 Triangulación por el método de Gauss 2 Variante Gauss-Jordan 24 Comentarios Más detalles Departamento de Matemática Aplicada FUNDAMENTOS DE MATEMATICAS. Ingeniería Química (Curso 2005-06) Álgebra Lineal Práctica 3
1. Matrices en Matlab Departamento de Matemática Aplicada FUNDAMENTOS DE MATEMATICAS. Ingeniería Química (Curso 2005-06) Álgebra Lineal Práctica 3 Para introducir una matriz en Matlab se procede de la Más detalles Panorama del curso Métodos Numéricos I
Panorama del curso Métodos Numéricos I Egor Maximenko ESFM del IPN 2014 Egor Maximenko (ESFM del IPN) Métodos Numéricos I 2014 1 / 35 Contenido 1 Propósito y programa del curso, software y literatura 2 Más detalles SITEMAS OPERATIVOS PROFESORA ELIZABETH ARIAS INFORMATICA SECCION F 3 SEMESTRE YORDAN ALEXANDER MARIN CUBIDES
REPUBLICA BOLIVARIANA DE VENEZUELA MINISTERIO DEL PORDER POPULAR PARA LA EDUCACION INSTITUTO UNIVERSITARIO POLITECNICO SANTIAGO MARIÑO EXTENSION MARACAIBO SITEMAS OPERATIVOS PROFESORA ELIZABETH ARIAS INFORMATICA Más detalles Programación en Capas.
Programación en Capas. Ricardo J. Vargas Del Valle Universidad de Costa Rica, Ciencias de Computación e Informática, San José, Costa Rica, 506 ricvargas@gmail.com Juan P. Maltés Granados Universidad de Más detalles SBDA-GRID: implementación de un servicio de base de datos en una arquitectura grid
SBDA-GRID: implementación de un servicio de base de datos en una arquitectura grid Nancy Ocotitla Rojas, Axel E. Moreno Cervantes, Rogelio Castro González Escuela Superior de Cómputo-IPN, Av. Juan de Dios Más detalles Espacios generados, dependencia lineal y bases
Espacios generados dependencia lineal y bases Departamento de Matemáticas CCIR/ITESM 14 de enero de 2011 Índice 14.1. Introducción............................................... 1 14.2. Espacio Generado............................................ Más detalles Modelado de relaciones existentes en un equipo de proyecto de software Modeling relationships in a software project team
Modelado de relaciones existentes en un equipo de proyecto de software Modeling relationships in a software project team Rafael Rodríguez-Puente 1, Eliana B. Ril-Valentin 2 1 Departamento de Técnicas de Más detalles Sybase IQ Servidor analítico con arquitectura basada en columnas
Sybase IQ Servidor analítico con arquitectura basada en columnas www.sybase.es Sybase IQ Descripción Tener acceso a toda la información de que dispone su organización, con el fin de analizarla no es hoy Más detalles SQL Server Reporting Services 2000 y 2005
SQL Server Reporting Services 2000 y 2005 Este artículo presenta SQL Server Reporting Services. Se presentan primero las funcionalidades de la primera versión de esta herramienta y después las mejoras Más detalles UN ENTORNO A MEDIDA PARA EL DISEÑO Y LA SIMULACIÓN DE MAQUINARIA POR COMPUTADOR
UN ENTORNO A MEDIDA PARA EL DISEÑO Y LA SIMULACIÓN DE MAQUINARIA POR COMPUTADOR Manuel González y Javier Cuadrado Departamento de Ingeniería Industrial II, Campus de Esteiro, 15403 Ferrol Universidad de Más detalles 2017 © DocPlayer.es Política de privacidad | Condiciones del servicio | Feedback

References: resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 Resolución 

Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución