Li essa resposta e percebi que tenho usado a diferença entre iterações sucessivas para definir um critério de parada para um método iterativo de encontrar autovalores / vetores. Quais são os bons critérios de parada para métodos iterativos que convergem para autovalores e...
Eu tenho um conjunto de arquivos que consiste em pontos selecionados aleatoriamente em um conjunto de dados, cada arquivo pertencendo a uma classe específica. Cada linha desses arquivos contém as coordenadas no espaço n do ponto. Gostaria de comparar as distribuições no espaço n de cada um desses...
Quais métodos numéricos estão disponíveis para encontrar o ponto fixo de um operador que atua nas funções ? Eu estou procurando a função para a qual .UMAUMAAf: [ a , b ] → [ a , b ]f:[uma,b]→[uma,b]f : [a,b] \rightarrow [a,b]fffA f= fUMAf=fAf = f Detalhes essenciais: Minha função é na verdade uma...
Lembro-me fracamente das minhas primeiras palestras "numéricas" de que os solucionadores lineares iterativos para geralmente exigem que quando seja decomposto comoAx=bAx=bAx=bAAA A=D+MA=D+MA=D + M onde D é uma matriz diagonal e tem diagonal zero, os elementos de devem ser dominantes sobre as...
Estou pensando em usar redes neurais artificiais (RNA) para prever as taxas de reação no meu fluido, em vez de resolver o sistema completo de ODEs rígidas. Algumas pessoas do meu laboratório já fizeram algum trabalho para não começar do zero, mas estou tendo problemas com meus aplicativos. Acho que...
O cálculo dos modos próprios de uma membrana semicircular se reduz ao seguinte problema de valor próprio ∇2u=k2u,∇2u=k2u,\nabla^2u=k^2u\;, onde a região de interesse é um semicírculo definido por e φ ∈ [ 0 , π ] .r ∈ [ 0 , 1 ]r∈[0,1]r\in[0,1]& Phi; ∈ [ 0 , π]φ∈[0,π]\varphi\in[0,\pi] É...
Estou analisando o desempenho dos solucionadores lineares do PETSc. Como eu entendo, speedup = Tempo sequencialTempo Paralelo.acelerar=Tempo SeqüencialTempo Paralelo.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Eu sei que a execução do código paralelo em um processador...
Quais são os "prós e contras" desses dois métodos de refinamento de malha? Ambos parecem ser os métodos predominantes. Naturalmente, posso imaginar que o refinamento global vermelho é relativamente fácil de implementar e atraente para a pesquisa em métodos multigrid, enquanto o NVB é usado na...
Estou tentando resolver um problema de mínimos quadrados não lineares de um lado com restrições lineares, ou seja, o problema: minx∑mi=1ri(x) s.t Ax≤bminx∑i=1mri(x) s.t Ax≤b\min_{\mathbf{x}} \quad \sum^m_{i=1} \mathbf{r}_i(\mathbf{x}) \qquad \text{ s.t } \quad A\mathbf{x} \leq...
Temos um código de refinamento de malha adaptativa (AMR) para resolver a equação de onda elástica com interfaces de falha por atrito (baseado no Chombo para aqueles que estão interessados). Uma das coisas que percebemos é que nossos resultados estão sendo fortemente afetados pela presença do limite...
Gostaria de testar a influência dos parâmetros de solvatação nos modelos implícitos de solvatação e saber quais códigos estão disponíveis gratuitamente como programas independentes para dobrar proteínas de pequenas proteínas e que usam abordagens de minimização de energia em vez de dinâmica, por...
Eu criei um spline para ajustar meus dados em python usando: spline=scipy.interpolate.UnivariateSpline(energy, fpp, k=4) A equação que quero usar envolve um somatório entre n = 2 en = infinito, onde n é a ordem do diferencial em um ponto Eo. No entanto, usando; UnivariateSpline.__call__(spline,...
Estou brincando com implementações simples de algoritmos para encontrar a cobertura mínima de vértices / correspondência de cardinalidade máxima em gráficos bipartidos. No entanto, pareço ter problemas para pesquisar alguns conjuntos de dados de teste que eu poderia usar para verificar, se...
Normalmente, quando eu proponho um gerador seqüencial de números aleatórios em C, uso a chamada srand(time(NULL)) então use rand() mod N para obter um número inteiro aleatório entre 0 e N-1. No entanto, quando faço isso em paralelo, as chamadas ao tempo (NULL) são tão próximas uma da outra...
Estou procurando uma biblioteca C ++ gratuita e conveniente para gráficos que incluem um solucionador para o problema do máximo subgrafo comum (MCS). Estou ciente da Boost Graph Library e do LEMON , mas nenhuma delas inclui um solucionador de
Eu sei que devo usar uma tolerância para comparar números de ponto flutuante. Mas para comparar vetores, posso pensar em 3 soluções possíveis correspondentes a diferentes métricas de distância: Compare os componentes de cada vetor individualmente: os vetores são iguais se todos os 3 estiverem...
Qual é uma maneira comum e boa de expressar dicionários (= mapas) no pseudocódigo? Ou seja, estruturas de dados que basicamente permitem armazenar valores para chaves, iterar sobre todos os pares chave / valor, testar a inclusão de uma determinada chave etc. Eu tenho algo parecido com o seguinte...
Suponha que você resolva um sistema linear por um método iterativo, por exemplo, gradientes conjugados ou iteração de Richardson. Então você tenta resolver um sistema linear que é ligeiramente perturbado na matriz e no lado direito, digamos, .˜ A ˜ u = ˜ fA u = fAu=fAu = fUMA~você~= f~A~u~=f~\tilde...
Quais são alguns bons recursos (livros, artigos, sites) sobre interseção de polígonos e algoritmos de
Escrevo códigos em Fortran e C para vários algoritmos matriciais. No entanto, quando eu perfil meus códigos usando o VTune, geralmente encontro alguma terminologia que não consigo apreciar completamente. Existe um bom recurso para aprender a criar perfis do zero para um nível bastante...