Perguntas com a marcação «petsc»

Eu tenho vários problemas desafiadores de otimização global não convexa para resolver. Atualmente, uso o Optimization Toolbox do MATLAB (especificamente, fmincon()com o algoritmo = 'sqp'), o que é bastante eficaz . No entanto, a maior parte do meu código está em Python, e eu adoraria fazer a...

O que pode dar errado ao usar métodos Krylov pré-condicionados da KSP ( pacote de solucionadores lineares do PETSc ) para resolver um sistema linear esparso, como os obtidos pela discretização e linearização de equações diferenciais parciais? Que etapas posso tomar para determinar o que está...

Até onde eu sei, as duas grandes estruturas genéricas de software de ciências computacionais do Departamento de Energia dos EUA são PETSc e Trilinos . Eles parecem semelhantes à primeira vista, além das diferenças de linguagem (C versus C ++). Quais são as principais diferenças entre as duas...

Estou usando o SNES do pacote de resolução não-linear do PETSc para resolver um sistema de equações não-lineares obtido discretizando uma equação diferencial parcial. Como posso determinar por que o solucionador não está convergindo e o que posso fazer para resolver com êxito minhas...

Estou planejando usar multigrid para calcular alguns autovalores e vetores, e notei que o PETSc tem suporte de alto nível para multigrid. A documentação do PETSc diz que essa parte do PETSc não deve ser usada, pois será substituída em breve. Quais outras bibliotecas têm suporte de alto nível para...

Gostaria de saber se existe uma maneira rápida de calcular a distância euclidiana de dois vetores na oitava. Parece que não há função especial para isso, então devo usar apenas a fórmula com

Eu tenho um PETSc Mate gostaria de estimar seu número de

Onde quer que eu tenha visto, os tutoriais / documentos do PETSc etc. dizem que é útil para álgebra linear e geralmente especifica que sistemas esparsos serão beneficiados. E matrizes densas? Estou preocupado sobre como resolver para densa .AA x = bUMAx=bAx=bUMAUMAA Eu escrevi meu próprio código...

Eu estava brincando com o PETSc e percebi que quando executo meu programa com mais de um processo via MPI, ele parece correr ainda mais devagar ! Como posso verificar para ver o que está

Normalmente, escrevo código serial e, quando o faço, escrevo testes de unidade com alguma estrutura de teste no estilo xUnit (MATLAB xUnit, PyUnit / nose ou estrutura de teste C ++ do Google). Com base em uma pesquisa superficial do Google, não tenho visto muito sobre como os profissionais...

Pergunta: Quais métodos estão disponíveis para calcular com precisão e eficiência a estrutura de esparsidade de uma matriz de elementos finitos? Info: Estou trabalhando em um solucionador de equações de pressão de Poisson, usando o método de Galerkin com base em Lagrange quadrático, escrito em C,...

Meu sistema é um problema simétrico de FE com multiplicadores de lagrange (por exemplo, fluxo incompressível de Stokes): ( ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} onde é o caso típico (eu até me certifiquei de que as equações sejam numeradas para que os multiplicadores de...

Preciso resolver o mesmo sistema linear esparso (300x300 a 1000x1000) com muitos lados do lado direito (300 a 1000). Além deste primeiro problema, eu também gostaria de resolver sistemas diferentes, mas com os mesmos elementos diferentes de zero (apenas valores diferentes), ou seja, muitos sistemas...

Aqui está um uso típico de ponteiros de função em C. Eu gostaria de fazer algo semelhante no Fortran. Tenho algumas idéias, mas gostaria de saber se existe alguma maneira canônica de fazer isso. Os ponteiros de função e os contextos transmitidos pelo usuário são armazenados e depois chamados...

Quanto mais lento o petsc4py vs c / c ++ / fortran? Sei que isso dependerá significativamente do código que está sendo executado, mas e algo simples como um produto de vetor de

A compilação do PETSc com uma biblioteca externa BLAS / LAPACK afeta significativamente o desempenho em matrizes esparsas ou apenas usa essas bibliotecas para matemática densa de

Eu sou muito novo na programação GPGPU, então, desculpe-me se a pergunta não for particularmente apropriada. Pelo que entendi, a programação da GPU é um trabalho muito complicado de engenharia, quando comparado à programação usual da CPU. É preciso ter muito cuidado com questões de divergência,...

Eu sou novo no pacote PETSc. Eu tenho uma matriz A de 4000 x 4000 no formato de mercado de matriz e quero que o PETSc resolva isso usando vários processadores. Sei como resolver o sistema em um único processador, mas não sei como distribuir a matriz e os vetores entre diferentes processadores....

Basicamente, o FEM parece ser um problema praticamente "resolvido". Existem inúmeras estruturas poderosas, como Trilinos, PETSc, FEniCS, Libmesh ou MOOSE. Uma coisa que eles têm em comum: são extremamente "pesados". Primeiro, a instalação normalmente é super dolorosa. Segundo, a interface / API...

Estou mergulhando no mundo fascinante da análise de elementos finitos e gostaria de resolver um grande problema termomecânico (apenas térmico mecânico, sem retorno).→→\rightarrow Para o problema mecânico, eu já entendi da resposta de Geoff , que precisarei usar o solucionador iterativo devido ao...