A resposta impulso é utilizada como ferramenta padrão no estudo direto de sistemas concentrados, discretos e distribuídos de ordem arbitrária. Esta abordagem leva ao desenvolvimento de uma plataforma unificada para a obtenção de respostas dinâmicas. Em particular, as respostas forçadas dos sistemas são decompostas na soma de uma resposta permanente e de uma resposta livre induzida pelos valores iniciais da resposta permanente. A teoria desenvolve-se de maneira geral e direta para sistemas de nésima ordem, introduzindo-se a base dinâmica gerada pela resposta impulso na forma padrão e normalizada, sem utilizar-se a formulação de estado, através da qual reduz-se um sistema de ordem superior para um sistema de primeira ordem. Considerou-se sistemas de primeira ordem a fim de acompanhar-se os muitos resultados apresentados na literatura através da formulação de espaço de estado. Os métodos para o cálculo da resposta impulso foram classificados em espectrais, não espectrais e numéricos. A ênfase é dada aos métodos não espectrais, pois a resposta impulso admite uma fórmula fechada que requer o uso de três equações características do tipo algébrica, diferencial e em diferenças. Realizou-se simulações numéricas onde foram apresentados modelos vibratórios clássicos e não clássicos. Os sistemas considerados foram sistemas do tipo concentrado, discreto e distribuído. Os resultados da decomposição da resposta dinâmica de sistemas concentrados diante de cargas harmônicas e não harmônicas foram apresentados em detalhe. A decomposição para o caso discreto foi desenvolvida utilizando-se os esquemas de integração numérica de Adams-Basforth, Strömer e Numerov. Para sistemas distribuídos, foi considerado o modelo de Euler-Bernoulli com força axial, sujeito a entradas oscilatórias com amplitude triangular, pulso e harmônica. As soluções permanentes foram calculadas com o uso da função de Green espacial. A resposta impulso foi aproximada com o uso do método espectral. / The impulse response is employed as a standard tool for a direct study of concentrated, discrete and distributed systems of arbitrary order. This approach leads to the development o f a unified platform for obtaining dynamical responses. In particular, forced responses are decomposed into the sum of a permanent response and a free response induced by the initial values of the permanent solution. The theory is developed in a general manner for n-th order systems; being introduced the standard dynamical basis generated by the impulse response and the normalized one, without employing the state formulation, through which a higher-order system is reduced to a first-order system. In order to follow the many results found in the literature through the state space formulation, first-order systems were considered. The methods for computing the impulse response were classified into spectral, non spectral and numeric. Emphasis was given to non spectral methods, because the impulse response has a closed-form formula that requires the use of three characteristic equations of algebraic, differential and difference type. Numerical simulations were performed with classical and non classical vibrating models. The systems considered were concentrated, discrete and distributed. The decomposition results of the forced response of concentrated systems subject to harmonic and non harmonic loads were worked out in detail. The decomposition for the discrete case was developed by using the numerical integration schemes of Adams-Basforth, Strõmer and Numerov. For distributed systems was considered the Euler-Bernoulli model with an axial force subject to oscillating inputs with triangular, pulse and harmonic amplitude. The permanent solutions were computed with the spatial Green function. The impulse response was approximated with the spectral method.

Neste trabalho estudou-se a modelagem matemática de processo de polimerização em suspensão de estireno, envolvendo a equação de balanço populacional para a previsão da distribuição de tamanhos de partículas. Neste problema os termos referentes aos processos de quebra e coalescência de partículas são importantes. Foram estudados métodos de solução da equação de balanço populacional e a solução desta equação, para o problema de polimerização em suspensão, foi obtida pelo método discretização do domínio em classes, usando a técnica dos pivôs fixos. Foram investigados os efeitos no número de elementos de discretização (número de classes) e o modo de modelar a distribuição inicial de tamanhos de gotas de monômero, quando esta não é medida experimentalmente. Neste caso, adotou-se uma distribuição normal com média elevada e variância pequena, arbitrárias, e simulou-se o processo de agitação com ocorrência de quebra e coalescência, até que se atingisse uma distribuição que não variasse com o tempo. Esta distribuição era usada então como condição inicial para a simulação do processo. Os resultados da simulação foram validados por comparação com dados experimentais referentes à evolução da distribuição de tamanhos das partículas, obtidos da literatura. O primeiro conjunto de dados foi referente ao processo de dispersão de gotas de estireno em água (sem ocorrência de polimerização). O segundo conjunto de dados foi referente ao processo de polimerização em suspensão de estireno em reator de batelada de escala de laboratório, nas quais os efeitos da velocidade de agitação, quantidade de agente estabilizante e fração de fase dispersa foram estudados. Em ambos os casos, o modelo representou corretamente o comportamento dos dados experimentais. / The mathematical modeling of suspension polymerization of styrene accounting for the population balance equation to predict the dynamic behavior of particle size distribution was studied. In this process, the population balance is determined by the rates of breakage and coalescence. Different methods for solving population balance (BP) equations were studied and, for the case of suspension polymerization, the numerical solution was obtained using the discretized BP by the method of classes and the fixed pivot technique. The effect of the number of discretization elements (number of classes) was investigated. In order to model the initial particle size distribution (PSD), which was not measured experimentally, an initial normal distribution with high average and low variance was assumed and the agitation process with breakage and coalescence was simulated until a constant PSD is reached. This PSD is then used as the initial condition for the process simulation. Simulation results were validated by comparison with experimental data taken from the literature. The first data set involves the dispersion process of styrene droplets in water, without polymerization. The second data set refers to the process of suspension polymerization of styrene carried out in a laboratory-scale batch reactor, in which the effects of agitation speed, amount of stabilizer (surface-active agent) and disperse phase holdup were studied. In both cases, the model was able to represent correctly the trends of the experimental data.

As equações polinomiais são estudadas desde a antiguidade e atualmente são utilizadas, por exemplo, para modelar problemas do cotidiano nas mais variadas áreas do conhecimento. As técnicas de solução de equações polinomiais nem sempre são triviais, principalmente quando envolvem equações de alta ordem e raízes complexas. O ensino desse tema no Ensino Básico é limitado a equações de segundo ou terceiro grau e coeficientes inteiros, o que restringe a aplicação em problemas mais realistas. Assim, o objetivo deste trabalho é trazer uma contribuição aos estudantes, aos professores do Ensino Básico e aos demais interessados, apresentando um material que aborde técnicas de resolução para equação polinomial de diversas naturezas. Iniciamos por uma revisão dos números complexos e dos polinômios, suas operações e propriedades. Embasamos o trabalho com teoremas e permeamos de exemplos com um crescente grau de dificuldade. Dividimos as técnicas de resolução em analíticas e numéricas. Entre as primeiras, tratamos das relações de Girard, das fórmulas resolventes e de alguns casos particulares de equações. Entre as técnicas numéricas, estudamos o método de Newton, o método das secantes e o método de Newton-Bairstow, este último para encontrar raízes complexas. / Polynomial equations have been studied since antiquity and are currently used, for example, to model everyday problems in the most varied areas of knowledge. The solution techniques of polynomial equations are not always trivial, especially when they involve high order equations and complex roots. The teaching of this subject in Basic Education is limited to second or third degree equations and integer coefficients, which restricts the application to more realistic problems. Thus, the objective of this work is to bring a contribution to students, teachers of Basic Education and other interested parties, presenting a material that treats of resolution techniques for polynomial equation of different natures. We begin with a review of complex numbers and polynomials, their operations and properties. We support the work with theorems and permeate examples with an increasing degree of difficulty. We divide the techniques of resolution into analytical and numerical. Among the first, we deal with Girards relations, the resolvent formulas, and some particular cases of equations. Among numerical techniques, we studied the Newton method, the secant method, and the Newton-Bairstow method, the last one to find complex roots.

Este trabalho apresenta a aplicação de técnicas de processamento paralelo na resolução de equações algébricas oriundas do Método de Elementos Finitos aplicado ao Eletromagnetismo, nos regimes estático e harmônico. As técnicas de programação paralelas utilizadas foram OpenMP, CUDA e GPUDirect, sendo esta última para as plataformas do tipo Multi-GPU. Os métodos iterativos abordados incluem aqueles do subespaço Krylov: Gradientes Conjugados, Gradientes Biconjugados, Conjugado Residual, Gradientes Biconjugados Estabilizados, Gradientes Conjugados para equações normais (CGNE e CGNR) e Gradientes Conjugados ao Quadrado. Todas as implementações fizeram uso das bibliotecas CUSP, CUSPARSE e CUBLAS. Para problemas estáticos, os seguintes pré-condicionadores foram adotados, todos eles com implementações paralelizadas e executadas na GPU: Decomposições Incompletas LU e de Cholesky, Multigrid Algébrico, Diagonal e Inversa Aproximada. Para os problemas harmônicos, apenas os dois primeiros pré-condicionadores foram utilizados, porém na sua versão sequencial, com execução na CPU, resultando em uma implementação híbrida CPU-GPU. As ferramentas computacionais desenvolvidas foram testadas na simulação de problemas de aterramento elétrico. No caso do regime harmônico, em que o fenômeno é regido pela Equação de Onda completa com perdas e não homogênea, a formulação adotada foi aquela em dois potenciais, A-V aresta-nodal. Em todas as situações, os aplicativos desenvolvidos para GPU apresentaram speedups apreciáveis, demonstrando a potencialidade dessa tecnologia para a simulação de problemas de larga escala na Engenharia Elétrica, com excelente relação custo-benefício. / This work presents the use of parallel processing techniques in Graphics Processing Units (GPU) for the solution of algebraic equations arising from the Finite Element modeling of electromagnetic phenomena, both in steadystate and time-harmonic regime. The techniques used were parallel programming OpenMP, CUDA and GPUDirect, the latter for those platforms of type Multi-GPU. The iterative methods discussed include those of the Krylov subspace: Conjugate Gradients, Bi-conjugate Gradients, Conjugate Residual, Bi-conjugate Gradients Stabilized, Conjugate Gradients for Normal Equations (CGNE and CGNR) and Conjugate Gradients Squared. All implementations have made use of CUSP, CUSPARSE and CUBLAS libraries. For the static problems, the following pre-conditioners were adopted, all with parallelized implementations and executed on the GPU: Incomplete decompositions, both LU and Cholesky, Algebraic Multigrid, Diagonal and Approximate Inverse. For the time-harmonic varying problems, only the first two pre-conditioners were used, but in their sequential version and running in the CPU, which yielded a hybrid CPU-GPU implementation. The developed computational tools were tested in the simulation of electrical grounding systems. In the case of the harmonic regime, in which the phenomenon is governed by the driven, lossy wave equation, the formulation adopted was that in two potential, the ungauged edge A-V formulation. In all cases, the developed GPU-based tools showed considerable speedups, showing that this is a promising technology for the simulation of large-scale Electrical Engineering problems, with excellent cost-benefit.

Neste trabalho analisa-se o comportamento transiente de um modelo unidimensional de ablação por sublimação na sua superfície. Este estudo trata de uma proteção térmica com Teflon, que tem propriedades que são funções da temperatura e está sujeita a um fluxo de calor constante. Emprega-se a consideração de que o material sublimado é substituído por um material fictício e a equação da condução para a entalpia, com um "calor específico efetivo" que considera o calor latente. O cálculo do "calor específico efetivo" é feito a partir da variação da entalpia com a temperatura através de um esquema implícito. Sendo a condutividade térmica função da temperatura, com um alto valor para a fase fictícia, garante-se que não ocorra gradiente de temperatura neste meio para que o fluxo de calor aplicado no contorno original seja conduzido até a localização da fronteira de fases. A influência do método de cálculo da capacidade térmica efetiva é avaliada comparando-se os perfis de temperatura. Estes resultados são decorrentes das simulações numéricas para refinamento da malha, curva da condutividade térmica e intervalo de troca de fase (2DT) nas três metodologias que são: propriedades como função da temperatura do nó com uso do Método de Diferenças Finitas, propriedades como função da temperatura variando linearmente no elemento com uso do Método de Elementos Finitos de Galerkin e propriedades como função da temperatura local no elemento através do Método de Elementos Finitos de Galerkin com integração numérica pela regra de Simpson. Um procedimento implícito no tempo é adotado de maneira a resolver-se um sistema de equações algébricas através do TDMA (algoritmo de Thomas) a cada passo de tempo. Através da análise desenvolvida neste trabalho identifica-se a forma de implementação que mais aproxima-se dos resultados disponíveis na literatura.

Made available in DSpace on 2016-04-27T16:57:17Z (GMT). No. of bitstreams: 1 Antonia Zulmira da Silva.pdf: 1733925 bytes, checksum: 3460611fce63b2ca5913c2f811ada3ac (MD5) Previous issue date: 2012-05-14 / Secretaria da Educação do Estado de São Paulo / The purpose of this investigation was to find evidence of indicators of algebraic thinking development for the topic 'First-degree algebraic equations' from the mathematics Teacher's Manual adopted by public middle schools in São Paulo state, Brazil, and thus provide a written contribution to this teaching material. The investigation sought to answer the following research questions: Do the activities proposed in the topic 'First-degree algebraic equations' from the mathematics Teacher s Manual for the third quarter of the eighth grade enable teachers to foster the development of algebraic thinking among students? If so, which indicators are most evident? The definition used for indicators of algebraic thinking development drew on Fiorentini, Miorim, and Miguel (1993) and Fiorentini, Fernandes and Cristóvão (2005) with regard to aspects of algebraic thinking and on Ursini et al. (2005) concerning use of variables. Concurrently, the so-called multimeanings of equations, as defined by Ribeiro and Machado (2009), were taken into account. Desk research, as defined by Lüdke and André (1986), was the method selected for the study. Of the twelve indicators of algebraic thinking development investigated, nine were detected in the activities examined. The results obtained showed that these activities enable teachers to foster the development of algebraic thinking among students. A final, stand-alone section summarizes the theoretical framework adopted and includes a chart of the algebraic thinking indicators investigated, in addition to a synthetic view of the analyses providing evidence of these indicators in the activities. This summarized section is also available in CD-ROM format / O presente estudo teve por objetivo evidenciar indicadores de desenvolvimento do pensamento algébrico no tópico Equações algébricas de primeiro grau do Caderno do professor de Matemática adotado na docência do Ensino Fundamental da rede pública do Estado de São Paulo, com a finalidade de escrever um produto que contribuísse com esse material. O objetivo se desdobrou nas seguintes questões de pesquisa: As atividades presentes no tópico Equações algébricas de primeiro grau do Caderno do professor de Matemática do terceiro bimestre do oitavo ano do Ensino Fundamental possibilitam que o professor conduza os alunos ao desenvolvimento do pensamento algébrico? Em caso afirmativo, que indicadores são priorizados? Para definir os indicadores de desenvolvimento do pensamento algébrico, tomamos como referências sobre o pensamento algébrico Fiorentini, Miorim e Miguel (1993) e Fiorentini, Fernandes e Cristóvão (2005) e, a respeito do uso das variáveis, Ursini et al. (2005). Ao mesmo tempo, investigamos os multissignificados das equações, segundo Ribeiro e Machado (2009). Para a condução da pesquisa, utilizamos o método de análise documental, conforme Lüdke e André (1986). Dentre os doze indicadores de desenvolvimento do pensamento algébrico considerados, nove foram evidenciados nas atividades analisadas. Os resultados permitiram concluir que as atividades analisadas possibilitam que o professor conduza os alunos a desenvolver o pensamento algébrico. O produto deste trabalho contém referências aos elementos teóricos do trabalho, um quadro com os indicadores do pensamento algébrico utilizados nas análises e a síntese das análises das atividades, evidenciando os indicadores do pensamento algébrico. Esse produto está anexado a esta dissertação e também encontra-se disponível em CD-ROM

Submitted by JOCIMAR MONTANHA null (jocimarmontanha@hotmail.com) on 2017-02-17T00:51:02Z No. of bitstreams: 1 Dissertação_Jocimar_Corrigida.pdf: 4981816 bytes, checksum: 3df5599a820e86618dced98b217c53ec (MD5) / Approved for entry into archive by LUIZA DE MENEZES ROMANETTO (luizamenezes@reitoria.unesp.br) on 2017-02-21T19:45:29Z (GMT) No. of bitstreams: 1 montanha_j_me_prud.pdf: 4981816 bytes, checksum: 3df5599a820e86618dced98b217c53ec (MD5) / Made available in DSpace on 2017-02-21T19:45:29Z (GMT). No. of bitstreams: 1 montanha_j_me_prud.pdf: 4981816 bytes, checksum: 3df5599a820e86618dced98b217c53ec (MD5) Previous issue date: 2017-02-03 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Este trabalho tem por objetivo principal apresentar uma sugestão de como introduzir e contextualizar os conceitos de números complexos, utilizando como motivações áudios, vídeos e software, além de outras atividades complementares sugeridas. Os áudios tratam dos números complexos através de uma história livremente inspirada no livro O Médico e o Monstro, do escritor escocês Robert Louis Stevenson. Os vídeos mostram uma maneira divertida e curiosa de olhar para os números complexos contando um pouco sobre sua história. O software tem a finalidade de estudar as transformações geométricas no plano (translação, rotação, dilatação e contração), utilizando os conceitos e operações de números complexos, propriedades e características geométricas. Este material faz parte da coleção M3 - Matemática Multimídia da Universidade Estadual de Campinas, e serviu como base para organizarmos o nosso trabalho. Outro software utilizado é o GeoGebra que servirá de suporte para a realização das soluções das demais atividades sugeridas. / This work has as main objective to present a suggestion of how to introduce and contextualize the concepts of complex numbers, using as motivation audios, videos and software, and other complementary activities suggested. Audios deal with complex numbers through a story loosely inspired by the book The Doctor and Monster, the Scottish writer Robert Louis Stevenson. The videos show a fun and funny way to look at the complex numbers telling a little about their history. The software aims to study the geometric transformations in the plane (translation, rotation, expansion and contraction), using the concepts and operations of complex numbers, geometric properties and characteristics. This material is part of the M3 - Multimedia Mathematics collection of the State University of Campinas, and served as a basis for organizing our work. Another software used is GeoGebra that will be used to support the solutions of the other suggested activities.

Este trabalho apresenta a aplicação de técnicas de processamento paralelo na resolução de equações algébricas oriundas do Método de Elementos Finitos aplicado ao Eletromagnetismo, nos regimes estático e harmônico. As técnicas de programação paralelas utilizadas foram OpenMP, CUDA e GPUDirect, sendo esta última para as plataformas do tipo Multi-GPU. Os métodos iterativos abordados incluem aqueles do subespaço Krylov: Gradientes Conjugados, Gradientes Biconjugados, Conjugado Residual, Gradientes Biconjugados Estabilizados, Gradientes Conjugados para equações normais (CGNE e CGNR) e Gradientes Conjugados ao Quadrado. Todas as implementações fizeram uso das bibliotecas CUSP, CUSPARSE e CUBLAS. Para problemas estáticos, os seguintes pré-condicionadores foram adotados, todos eles com implementações paralelizadas e executadas na GPU: Decomposições Incompletas LU e de Cholesky, Multigrid Algébrico, Diagonal e Inversa Aproximada. Para os problemas harmônicos, apenas os dois primeiros pré-condicionadores foram utilizados, porém na sua versão sequencial, com execução na CPU, resultando em uma implementação híbrida CPU-GPU. As ferramentas computacionais desenvolvidas foram testadas na simulação de problemas de aterramento elétrico. No caso do regime harmônico, em que o fenômeno é regido pela Equação de Onda completa com perdas e não homogênea, a formulação adotada foi aquela em dois potenciais, A-V aresta-nodal. Em todas as situações, os aplicativos desenvolvidos para GPU apresentaram speedups apreciáveis, demonstrando a potencialidade dessa tecnologia para a simulação de problemas de larga escala na Engenharia Elétrica, com excelente relação custo-benefício. / This work presents the use of parallel processing techniques in Graphics Processing Units (GPU) for the solution of algebraic equations arising from the Finite Element modeling of electromagnetic phenomena, both in steadystate and time-harmonic regime. The techniques used were parallel programming OpenMP, CUDA and GPUDirect, the latter for those platforms of type Multi-GPU. The iterative methods discussed include those of the Krylov subspace: Conjugate Gradients, Bi-conjugate Gradients, Conjugate Residual, Bi-conjugate Gradients Stabilized, Conjugate Gradients for Normal Equations (CGNE and CGNR) and Conjugate Gradients Squared. All implementations have made use of CUSP, CUSPARSE and CUBLAS libraries. For the static problems, the following pre-conditioners were adopted, all with parallelized implementations and executed on the GPU: Incomplete decompositions, both LU and Cholesky, Algebraic Multigrid, Diagonal and Approximate Inverse. For the time-harmonic varying problems, only the first two pre-conditioners were used, but in their sequential version and running in the CPU, which yielded a hybrid CPU-GPU implementation. The developed computational tools were tested in the simulation of electrical grounding systems. In the case of the harmonic regime, in which the phenomenon is governed by the driven, lossy wave equation, the formulation adopted was that in two potential, the ungauged edge A-V formulation. In all cases, the developed GPU-based tools showed considerable speedups, showing that this is a promising technology for the simulation of large-scale Electrical Engineering problems, with excellent cost-benefit.

Made available in DSpace on 2016-04-27T16:58:48Z (GMT). No. of bitstreams: 1 Adriano de Morais Martins.pdf: 735599 bytes, checksum: 20cee7a87d8b09965aa11f3b1d3d7646 (MD5) Previous issue date: 2008-11-07 / Secretaria da Educação do Estado de São Paulo / The present study makes a synthesis of brazilian s researches whose focus are in Fundamental schools, that are about algebraic equations. The analysis material constituted of nine dissertations published on the Mathematics Educations between 1998 and 2004. The selections of these essays were done according to some researches instruments and pre-defined criterions. The methodology characterizes itself as a documental study denominated qualitative metanalysis which makes a systematic revision of searches set, aiming at realization of a synthesis of these productions. Here, is presented one of these study confront, more specifically among their objectives referential, methodological theoretical, conclusions and suggestions indicated to futures researches. The results confirms that almost all the analyzed researches , in most part, don t mention the last ones. They present some worries about the equations meaning. There is a big dispersion in the choice of the theoretical indicators and/or theoretical methodological. The study indicate results that can offer subsidies and new researches / O presente estudo faz uma síntese de pesquisas brasileiras voltadas ao Ensino Fundamental ciclo II, que tratam das equações algébricas. O material de análise constitui-se de nove dissertações publicadas entre os anos de 1998 a 2004 na área de Educação Matemática. A seleção desses trabalhos se deu por meio da busca de dissertações e teses segundo alguns instrumentos de pesquisas e critérios pré - determinados. Metodologicamente esse estudo se caracteriza como um estudo documental denominado metanálise qualitativa na qual se procura fazer uma revisão sistemática de um conjunto de pesquisas, visando à realização de uma síntese dessas produções. Aqui é apresentado um confronto desses estudos, mais especificamente entre seus objetivos, referenciais teóricos, teórico-metodológicos, conclusões e sugestões deixadas para futuras pesquisas. Os resultados atestam que as pesquisas analisadas, em sua maioria, não citam as anteriores já publicadas, apresentam uma preocupação com o significado equações e que existe uma grande dispersão nas escolhas dos indicadores teóricos e/ou teórico metodológicos. Com isso o presente estudo aponta resultados que possam oferecer subsídios a novas pesquisas, tais como a carência de alguns pontos a serem pesquisados dentro desse tema de ensino

Estudamos a otimização do método SOR clássico, para a resolução de um sistema linear Ax = b, com A não-singular, a partir dos resultados de Young [55, 57] e Varga [50, 51] para matrizes de blocos p-cíclicas consistentemente ordenadas. Num primeiro nível, a otimização refere-se à escolha do parâmetro de relaxação do SOR que produz a maior velocidade de convergência, e, num segundo nível, à escolha da p-ciclicidade que apresenta o melhor desempenho com os valores ótimos do parâmetro, e damos ênfase ao caso 2-cíclico. Além disso, descrevemos a otimização do parâmetro em três generalizações: a) num relaxamento das condições sobre o espectro da matriz de Jacobi associada a A; b) no método SOR para matrizes singulares; c) num novo método SOR, que substitui a decomposição A = D - L - U, onde D, L e U são a diagonal de A, a parte triangular inferior estrita de A e a parte triangular superior estrita de A, pela A = D - P - Q, onde P pertence a uma classe de matrizes constru ída a partir das matrizes-escada. Descrevemos também a aplicação do caso singular às cadeias de Markov, comentamos a computação paralela aplicada ao SOR, e apresentamos diversas simulações relativas à otimização desse método. / We study the optimization of the classic SOR method for solving a linear system Ax = b, where A is a nonsingular p-cyclic consistently ordered block matrix, based on the discoveries of Young [55, 57] and Varga [50, 51]. In a first levei, the optimization refers to the choice of the SOR relaxation parameter, which produces the greatest convergence speed and, in a second levei, to the p-cyclicity that presents the best performance with the optimal parameter values and emphasize the 2- cyclic case. Moreover we describe three SOR generalizations concerning optimization: a) by weakening the conditions on the spectrum of Jacobi matrix associated with A; b) by considering the SOR method for singular matrices; c) by approaching a new SOR, that replaces the splitting A = D - L - U, where O, L and U are the diagonal of A, the strict lower triangular part of A and the strict upper triangular part of A. respectively, by this one A = D - P - Q, where P is a stair matrix or a matrix even more general than a stair matrix. We also describe the application of the singular case to Markov chains, discuss parallel computing applied to SOR method, and present severa! simulations regarding the optimization of that method.