Quero publicar o pacote de software que escrevi para meu trabalho de pós-graduação, juntamente com um artigo descrevendo o pacote. Até agora, estive considerando licenças liberais como o BSD. No entanto, agora estou ciente de que a publicação do seu código pode envolver um acordo com os termos de...
Fundo: Atualmente, estou executando uma grande quantidade de experimentos de variação de parâmetros. Eles estão sendo executados no Python 2.6+, usando numpy. Esses experimentos levarão cerca de 2 semanas para serem executados. Aproximadamente, estou variando três parâmetros (variáveis...
Eu codifiquei um problema do Stokes Flow usando elementos finitos e estou no processo de verificar se ele funciona. Só não tenho certeza de qual taxa de convergência devo esperar enquanto refino a malha globalmente. Conheço problemas escalares usando funções de base linear que eu esperaria...
Estou entrando no meu terceiro ano do meu programa de doutorado e até agora meu código de pesquisa (problemas inversos numéricos / imagens médicas / processamento de imagens / etc) consiste principalmente de scripts e funções desorganizados do MATLAB, com alguns C ++ e CUDA MEX'd rotinas lançadas...
Durante minha peregrinação no Mathematica.se , notei gradualmente que um certo tipo de problema de solução de equações diferenciais está nos "incomodando" o tempo todo, ou seja, o problema do valor limite (BVP) das equações diferenciais ordinárias não lineares (ODEs). O método de disparo , usado...
Ao escrever código MPI 3.0, que deve passar argce argvà MPI_Initchamada ou não, e por quê? EDIT: Como a resposta de Bill Barth levantou algumas perguntas, quero fazer algumas observações: A aprovação argc/ argvnão é necessária desde o MPI 1.1. A pergunta é especificamente sobre por que deveria /...
dada uma matriz simétrica , e uma matriz arbitrária X ∈ R n × n , e um vetor v ∈ R n × 1 , é possível calcular a seguinte expressão no tempo O ( n 2 ) ?Y∈Rn×nY∈Rn×nY \in \mathbb{R}^{n \times n}X∈Rn×nX∈Rn×nX \in \mathbb{R}^{n \times n}v∈Rn×1v∈Rn×1v \in \mathbb{R}^{n \times...
Estou procurando um método ou algoritmo para produzir um valor que descreva como duas malhas diferentes são geometricamente, mas que têm topologias diferentes. Um exemplo seria alguns dados CAD que tiveram uma malha FE aplicada e a peça depois que ela foi fabricada e, em seguida, digitalizada em...
Agora que o Formulize / Eureqa começou a cobrar US $ 2500 por ano por usá-lo e ter aleijado a versão de teste, alguém sabe de alguma substituição que possa fazer coisas semelhantes, como encontrar uma equação com base em pontos de dados? Era gratuito, mas agora está fora da faixa de preço de...
Ultimamente, tenho lido muitos artigos sobre o ADMM e também tentei resolver vários problemas usando-o, nos quais foi muito eficaz. Ao contrário de outros métodos de otimização, não consigo ter uma boa intuição de como e por que esse método é tão eficaz (é claro, já vi análises de convergência em...
Em um artigo recente da SIAM News , há um longo artigo descrevendo uma organização sistemática dos elementos finitos, apropriadamente apelidada de Tabela Periódica dos Elementos Finitos . É realmente muito fascinante ver como a classificação pode ser realizada através do cálculo externo de...
Para um protetor de tela bobo que estou tentando desenvolver, gostaria de gerar aleatoriamente uma matriz 2D livre de divergência de vetores 2D e usá-la para gerar um gráfico de convolução integral de linha. Ouvi que uma maneira de fazer isso é gerar ruído aleatório e depois projetar o componente...
Eu tenho um vetor que pode ser decomposto no espaço eigens do operador esparso hermitiano :MVVVMMM V=∑ivim^iV=∑ivim^iV = \sum_i v_i \hat{m}_i Existe uma maneira de encontrar o (o próprio vetor próprio) que corresponde ao maior (em magnitude)?vim^im^i\hat{m}_iviviv_i Eu quero essencialmente os...
Fiquei me perguntando se existe algum teorema que me permita colocar um limite superior no erro introduzido pela omissão de pequenos elementos da matriz antes da diagonalização. Vamos supor que temos uma matriz grande, cujos elementos da matriz variam de a . Se eu definisse todos os elementos da...
Tentei ir às fontes primárias para entender como usar as tabelas de Butcher para simplificar a álgebra que preciso fazer ao usar a série Taylor para encontrar a ordem de precisão de um esquema, por exemplo. No entanto, talvez por falta de antecedentes relevantes, achei particularmente difícil...
É um truque bastante conhecido evitar divisão no cálculo de raízes quadradas e aplicar o método de Newton para encontrar 1 / x--√1/x1/\sqrt{x} , e provavelmente melhor conhecido, usando o método de Newton para encontrarrecíprocos sem divisão. Ao resgatar um encadeamento StackOverflow Semeando a...
Eu tenho uma pergunta bastante específica em relação ao número da condição. Eu corro simulações FEM que possuem várias escalas de comprimento, o que resulta em uma enorme disparidade entre as entradas maiores e as menores da minha matriz. O número da condição pode chegar a 10 ^ 15 em algumas...
Existe alguma biblioteca / caixa de ferramentas que tenha implementação de SVD incremental no MATLAB. Eu implementei este documento , é rápido, mas não funciona bem. Eu tentei isso, mas neste erro também se propaga rapidamente (dentro de atualização de 5 a 10 pontos, o erro é...
As equações de movimento para um sólido elástico são dadas por ∇⋅σ+f=ρu¨σ=Cεε=12(∇u+[∇u]T)∇⋅σ+f=ρu¨σ=Cεε=12(∇u+[∇u]T)\begin{align} &\nabla \cdot \boldsymbol{\sigma} + \mathbf{f} = \rho \ddot{\mathbf{u}}\\ &\boldsymbol{\sigma} = \mathbb{C}\boldsymbol{\varepsilon}\\ &\boldsymbol{\varepsilon} =...
Suponha que tenhamos um problema de valor inicial no formato onde x 0 ∈ R n é conhecido exatamente (ou seja, com precisão ilimitada) e podemos avaliar eficientemente f : R n → R n com qualquer precisão. Ou seja, temos uma caixa preta que, dado um vetor x ∈ R n e um número inteiro M , retorna uma...