Na literatura de métodos numéricos, muitos trabalhos de pesquisa consistem em uma descrição de uma nova variação algorítmica, seguida de alguns problemas de teste comparando o novo método com um ou dois métodos existentes. Isso dificulta a determinação Como o novo método funcionaria em outros...
Estou interessado nos casos em que o gradiente conjugado funciona muito melhor que o método GMRES. Em geral, o CG é a escolha preferível em muitos casos de SPD (simétrico-positivo-definido) porque requer menos armazenamento e o limite teórico da taxa de convergência para o CG é o dobro do GMRES....
Estou à procura de um software de geração de malha que é gratuito e de código aberto, fornece uma interface de script sã para especificação de domínio, trabalha para geometrias complexas, pode gerar malhas 2D e 3D, Que opções eu
Eu sempre calculo as coisas manualmente, mas agora meus companheiros estão ficando desagradáveis e fazendo muitos exercícios repetitivos envolvendo apenas conectar coisas como a expressão acima. Estou particularmente interessado em software de código aberto, como Python ou R, para simplificar...
Suponha que eu tenho alguma função e quero encontrar tal que . Eu poderia usar o método Newton-Raphson. Mas isso requer que eu conheça a função derivada . Uma expressão analítica para pode estar indisponível. Por exemplo, pode ser definido por um código complicado de computador que consulta um...
Eu preciso escrever alguns dados de uma computação, que serão lidos posteriormente pelo Paraview (arquivo .vtu ou vtk). No que diz respeito ao tamanho do arquivo, devo escolher o formato ASCII ou o formato
Recentemente, comparei diferentes solucionadores não lineares do scipy e fiquei particularmente impressionado com o exemplo de Newton-Krylov no Scipy Cookbook, no qual eles resolvem uma equação de equação diferencial de segunda ordem com o termo de reação não linear em cerca de 20 linhas de...
Muitos métodos numéricos para EDPs hiperbólicas são baseados no uso de solucionadores de Riemann. Esses solucionadores são essenciais para capturar com precisão as ondas de choque. Há uma variedade desses solucionadores disponíveis para os sistemas mais bem estudados (por exemplo, solucionadores...
Os métodos numéricos para resolver EDPs (ou EDOs) se enquadram em duas grandes categorias: métodos explícitos e implícitos. Métodos implícitos permitem timestais estáveis maiores, mas exigem mais trabalho por etapa. Para PDEs hiperbólicas, o senso comum é que os métodos implícitos geralmente não...
Eu gostaria de ser independente do software comercial para o meu trabalho científico. Acho que uma dependência de pacotes comerciais como o Matlab e suas caixas de ferramentas é insatisfatória, porque não sei se terei acesso ao Matlab no futuro e porque não gosto do idioma. Portanto, estou...
Dada uma evolução PDE vocêt= A u + B uvocêt=UMAvocê+Bvocêu_t = Au + Bu onde são operadores diferenciais (possivelmente não lineares) que não se deslocam, uma abordagem numérica comum é alternar entre resolverA , BUMA,BA,B vocêt= A uvocêt=UMAvocêu_t = Au e vocêt= B u .vocêt=Bvocê.u_t = Bu. A...
Estou em um campo que não faz necessariamente muito trabalho com HPC e, quando o encontra, geralmente é o resultado de pesquisadores de outros campos explorando novas aplicações para seus métodos e afins. Principalmente o que isso significa é que ele nunca é realmente introduzido no curso dos...
Considere com quase singular, o que significa que há um autovalor de que é muito pequeno. O critério de parada usual de um método iterativo é baseado no residual e diz respeito às iterações podem parar quando com n o número da iteração. Mas no caso que estamos considerando, pode haver um grande...
Estou executando uma tarefa de classificação de texto com R e obtenho uma matriz de termo de documento com tamanho 22490 por 120.000 (apenas 4 milhões de entradas diferentes de zero, menos de 1%). Agora, quero reduzir a dimensionalidade utilizando o PCA (Principal Component Analysis). Infelizmente,...
O que mais detesto no MPI é lidar com tipos de dados (isto é, mapas de dados / máscaras) porque eles não se encaixam muito bem com o C ++ orientado a objetos. boost::mpino entanto, apenas suporta o MPI 1.1 no site: O boost :: mpi é uma interface amigável em C ++ para a Interface de transmissão...
Comecei a usar python como a linguagem de programação para realizar todas as minhas tarefas no CFD. Eu tenho uma experiência muito pequena em programação. Sou formado em engenharia mecânica e estou cursando o ensino superior em engenharia aeroespacial. Às vezes, o aspecto computacional do CFD se...
Estou usando o Valgrind + Callgrind para criar um perfil para o solucionador que escrevi. Conforme o manual do usuário do Valgrind, compilei meu código com as opções de depuração para o compilador: "Sem as informações de depuração, a melhor ferramenta que o Valgrind poderá fazer é adivinhar a...
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
Provavelmente é uma questão de nível de estudante, mas não posso fazer isso exatamente para mim. Por que é mais preciso usar grades não uniformes nos métodos numéricos? Estou pensando no contexto de algum método de diferenças finitas para o PDE da forma . E suponha que eu esteja interessado em uma...
Sou um grande defensor do desenvolvimento orientado a testes em computação científica. Sua utilidade na prática é simplesmente impressionante e realmente alivia os problemas clássicos que os desenvolvedores de código conhecem. No entanto, existem dificuldades inerentes ao teste de códigos...