Ciência computacional

Existe um método para resolver sistemas lineares da forma onde é uma matriz SPD fixa e são matrizes diagonais positivas?O(n3+n2k)O(n3+n2k)O(n^3+n^2 k)kkk(Di+A)xi=bi(Di+A)xi=bi(D_i + A) x_i = b_iAAADiDiD_i Por exemplo, se cada é escalar, basta calcular o SVD Uma . No entanto, isso se divide em D geral devido à falta de …

21 linear-algebra  algorithms  performance  complexity 

Não consigo encontrar nenhuma literatura sobre algoritmos que possa ser usada para resolver um problema de atribuição generalizada (GAP) muitos para muitos, ou seja, modelos em que não apenas mais tarefas podem ser atribuídas a um agente, mas vários agentes também podem ser usados. atribuídos a uma tarefa (APs um …

20 algorithms  data-analysis 

Estou procurando uma biblioteca de tensores C ++ que suporte código independente de dimensão. Especificamente, preciso executar operações ao longo de cada dimensão (até 3), por exemplo, calcular uma soma ponderada. As dimensões são um parâmetro de modelo (e, portanto, uma constante em tempo de compilação). Outra restrição é que …

20 c++  libraries  tensor 

Foi sugerido que este poderia ser um lugar melhor para esta pergunta do que o Mathematics Stack Exchange, onde eu a fiz antes . Suponha que se tenha uma função de caixa preta que possa ser avaliada em qualquer lugar (barato) em um intervalo especificado e não tenha ruído (exceto …

20 interpolation  signal-processing 

Para um projeto em que estou trabalhando (em PDEs hiperbólicas), gostaria de ter uma idéia mais detalhada do comportamento, analisando alguns números. No entanto, não sou um programador muito bom. Você pode recomendar alguns recursos para aprender a codificar efetivamente esquemas de diferenças finitas no Scientific Python (outros idiomas com …

20 python  finite-difference  reference-request  hyperbolic-pde 

Fiquei muito surpreso quando comecei a ler algo sobre otimização não convexa em geral e vi declarações como esta: Muitos problemas práticos de importância são não convexos e a maioria dos problemas não convexos são difíceis (se não impossíveis) de resolver exatamente em um tempo razoável. ( fonte ) ou …

20 optimization  nonconvex 

Existe uma maneira de usar algum pacote Python estabelecido (por exemplo, SciPy) para definir minha própria função de densidade de probabilidade (sem dados anteriores, apenas ), para que eu possa fazer cálculos com ele (como obter a variância da variável aleatória contínua)? É claro que eu poderia, digamos, SymPy ou …

20 python  statistics  probability 

Eu tenho um aplicativo que pode ser trivialmente paralelizado, mas seu desempenho é em grande parte limitado a E / S. O aplicativo lê uma única matriz de entrada armazenada em um arquivo com tamanho geralmente de 2 a 5 GB (mas espero que esse número aumente no futuro). Um …

20 io  parallel-computing  mpi 

Comecei a estudar o MEF de forma mais estruturada em comparação com o que costumava fazer durante meus cursos de graduação. Estou fazendo isso porque, apesar de poder usar o "FEM" em softwares comerciais (e outros não comerciais), gostaria de realmente entender as técnicas subterrâneas que suportam o método. É …

20 finite-element 

Suponha que eu tenha uma função que tome como entrada vários valores de ponto flutuante (simples ou duplo), faça alguma computação e produza valores de ponto flutuante de saída (também únicos ou duplos). Estou trabalhando principalmente com o MSVC 2008, mas também planejo trabalhar com o MinGW / GCC. Estou …

20 floating-point  c++ 

Estou com um pouco de dificuldade em tentar entender um artigo. O artigo usa o método espectral para resolver um autovalor proveniente de um sistema de ODEs acoplados. Escreverei apenas uma equação agora, porque é suficiente para chegar ao cerne da (s) minha (s) pergunta (s). A equação é V[ …

19 eigenvalues  polynomials  spectral-method 

A maioria dos métodos numéricos para quadratura trata o integrando como uma função de caixa preta. E se tivermos mais informações? Em particular, que benefício, se houver, podemos derivar de conhecer as primeiras derivadas do integrando? Que outras informações podem ser valiosas? Para derivativos em particular: estimativas de erro para …

19 quadrature  error-estimation 

Estou pensando em usar o OpenFOAM para resolver fluxos internos básicos em CFD. Qual é a melhor maneira de começar, e alguém poderia me indicar uma boa referência on-line para responder a quaisquer perguntas que eu possa ter quando mergulhar? Ouvi dizer que é bastante nicho a partir de agora …

19 fluid-dynamics  openfoam 

A pesquisa reproduzível em computação visa disponibilizar a outros pesquisadores o código necessário para gerar os resultados em um artigo computacional, para que eles possam executar esse código para reproduzir os resultados nesse artigo. Gostaria de tornar toda a minha pesquisa reproduzível, mas estou enfrentando alguns problemas: alguns trabalhos em …

19 software  publications  documentation  reproducibility 

Costumo encontrar-me escrevendo código muito semelhante para uma, duas e três versões dimensionais de uma dada operação / algoritmo. Manter todas essas versões pode se tornar tediosa. A geração simples de código funciona razoavelmente bem, mas parece que deve haver uma maneira melhor. Existe uma maneira relativamente simples de escrever …

19 algorithms  software  fortran