Eu tenho uma pergunta simples que é realmente difícil para o Google (além do artigo canônico O que todo cientista da computação deve saber sobre aritmética de ponto flutuante ). Quando funções como log1pou expm1devem ser usadas em vez de loge exp? Quando eles não devem ser usados? Como as...
Alguém pode me recomendar um bom sistema de gerenciamento de fluxo de trabalho (WMS), preferencialmente em Python? Até agora, tenho usado o GNU Make, mas apresenta uma camada de complexidade que eu quero evitar. Um bom WMS deve ter os seguintes recursos: integre-se facilmente com ferramentas de...
Os engenheiros geralmente insistem em usar métodos conservadores localmente, como volume finito, diferença finita conservadora ou métodos descontínuos de Galerkin para resolver PDEs. O que pode dar errado ao usar um método que não é localmente conservador? Ok, então a conservação local é...
Recentemente, encontrei uma formulação do meta-fenômeno : " dois é fácil, três é difícil " (formulado dessa maneira por Federico Poloni), que pode ser descrito da seguinte forma: Quando um determinado problema é formulado para duas entidades, é relativamente fácil de resolver; no entanto, um...
O Journal of Computational Physics tem sido uma saída importante para a ciência da computação no passado, e eu já o publiquei antes. Para o benefício daqueles (como eu) que assinaram o boicote à Elsevier , que periódicos não pertencentes à Elsevier seriam lugares apropriados para publicar artigos...
Se eu tenho uma matriz quadrada invertível e tomo seu determinante, e acho que , isso implica que a matriz está mal condicionada?det(A)≈0det(A)≈0\det(A) \approx 0 O inverso também é verdadeiro? Uma matriz mal condicionada tem um determinante quase zero? Aqui está algo que eu tentei no Octave: a...
Antecedentes: Eu acho que posso querer portar algum código que calcule produtos vetoriais exponenciais da matriz usando um método de subespaço Krylov do MATLAB ao Python. (Especificamente, a função expmvp de Jitse Niesen , que usa um algoritmo descrito neste documento .) No entanto, eu sei que, a...
O problema Atualmente, estou trabalhando em uma simulação de Elementos Finitos de Navier Stokes e gostaria de investigar os efeitos de vários parâmetros. Alguns parâmetros são especificados em um arquivo de entrada ou através de opções de linha de comando; outros parâmetros são fornecidos como...
Como alguém que é bacharel em física, fiquei um pouco escandalizado quando comecei a trabalhar com simulações moleculares. Foi um pouco chocante descobrir que mesmo as simulações mais detalhadas e computacionalmente caras não podem reproduzir quantitativamente todo o comportamento da água a partir...
A versão 3.0 do padrão MPI excluiu formalmente a interface C ++ (ela foi descontinuada anteriormente). Embora as implementações ainda possam suportá-lo, os recursos novos no MPI-3 não possuem uma interface C ++ definida no padrão MPI. Consulte
Eu conheço e usei o f2py2e para quebrar algum código antigo do Fortran 77, mas meu entendimento é que ele não funciona com o código Fortran 95 mais recente. Eu pesquisei o que eu deveria usar e me deparei com o fwrap e o G3 f2py, nenhum dos quais parece dar qualquer explicação sobre seu estado...
Eu quero modificar uma matriz de transição quadrada densa no local, alterando a ordem de várias de suas linhas e colunas, usando a biblioteca numpy do python. Matematicamente, isso corresponde à pré-multiplicação da matriz pela matriz de permutação P e à pós-multiplicação por P ^ -1 = P ^ T, mas...
Esta pergunta foi migrada do Stack Overflow porque pode ser respondida no Computational Science Stack Exchange. Migrou há 7 anos . Existem problemas / algoritmos famosos na computação científica que não podem ser acelerados pela paralelização? Parece-me ao ler livros...
EDIT: Estou testando se algum autovalor tem uma magnitude de um ou mais. Preciso encontrar o maior autovalor absoluto de uma grande matriz esparsa e não simétrica. Eu tenho usado a eigen()função de R , que usa o QR algo do EISPACK ou LAPACK para encontrar todos os autovalores e, em seguida,...
Nas estatísticas e suas diversas aplicações, geralmente calculamos a matriz de covariância , que é positiva definida (nos casos considerados) e simétrica, para vários usos. Às vezes, precisamos do inverso dessa matriz para vários cálculos (formas quadráticas com esse inverso como a (única) matriz...
Os desenvolvedores de software têm a opção de escolher uma licença apropriada de acordo com o (s) objetivo (s) do trabalho. Alguém pode dar algumas recomendações / experiências sobre qual licença escolher para o software? Quais são os prós / contras de "doar" todo o trabalho codificado como...
Uma das principais questões com as quais temos de lidar em simulações moleculares é o cálculo de forças dependentes da distância. Se pudermos restringir as funções de força e distância para ter potências iguais da distância de separação , podemos calcular o quadrado da distância r 2 = r ⋅ r e não...
Não estou muito familiarizado com os esquemas comuns de discretização para PDEs. Eu sei que Crank-Nicolson é um esquema popular para discretizar a equação de difusão. Também é uma boa escolha para o termo de advecção? Eu sou interessante em resolver a equação Reação-Difusão-Advecção...
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á...
Existem ferramentas projetadas especificamente para compactar dados científicos de ponto flutuante? Se uma função é suave, obviamente há muita correlação entre os números que representam essa função; portanto, os dados devem ser compactados bem. Compactar / compactar dados binários de ponto...