Ciência computacional

Estou iniciando meus estudos de pós-graduação em engenharia e trabalhando em projetos de ciências computacionais. Percebi que houve alguma discussão sobre as vantagens e desvantagens de implementar seus próprios algoritmos aqui . A análise de custo-benefício da implementação de seus próprios algoritmos versus o uso de bibliotecas é diferente quando …

10 software  education 

Eu sou muito novo na programação GPGPU, então, desculpe-me se a pergunta não for particularmente apropriada. Pelo que entendi, a programação da GPU é um trabalho muito complicado de engenharia, quando comparado à programação usual da CPU. É preciso ter muito cuidado com questões de divergência, lado a lado, alocação …

10 parallel-computing  petsc  gpu  cuda 

Estou resolvendo o sistema de reação-difusão de Turing com o seguinte código C ++. É muito lento: para a textura de 128x128 pixels, o número aceitável de iterações é 200 - o que resulta em 2,5 segundos de atraso. Preciso de 400 iterações para obter uma imagem interessante - mas …

10 pde  stiffness 

Estou tentando diagonalizar algumas matrizes densas e mal condicionadas. Na precisão da máquina, os resultados são imprecisos (retornando valores próprios negativos, os vetores próprios não possuem as simetrias esperadas). Eu mudei para a função Eigensystem [] do Mathematica para tirar proveito da precisão arbitrária, mas os cálculos são extremamente lentos. …

10 linear-algebra  parallel-computing  eigensystem  precision  dense-matrix 

No momento, estou trabalhando em um método de solução de equações diferenciais chamado colocação de base-spline. Estou com problemas para criar um método para construir um spline de ordem arbitrária, com o relacionamento com o condição inicial e estou tendo problemas até para começar com esse problema, pois é recursivo …

10 c++  b-spline 

Como trabalho na faculdade, fiz uma simulação em 1d. A declaração do problema era resolver o problema do tubo de choque 1d envolvendo gás ideal compressível como fluido de trabalho. Para esse problema, resolvi o sistema de equações de Eulers usando o solucionador Riemann de Roe. Quero saber, para resolver …

10 fluid-dynamics 

A compilação do PETSc com uma biblioteca externa BLAS / LAPACK afeta significativamente o desempenho em matrizes esparsas ou apenas usa essas bibliotecas para matemática densa de matriz?

10 petsc  sparse-matrix  lapack  blas 

A multiplicação da matriz (Mat * Mat e Mat * Vec) é dimensionada com o número de não-zeros ou com o tamanho da matriz? Ou alguma combinação dos dois. Que tal com forma. Por exemplo, eu tenho uma matriz 100 x 100 com 100 valores ou uma matriz 1000 x …

10 linear-algebra  performance  sparse-matrix 

Eu tenho que fazer uma transformação de coordenadas entre dois sistemas de referência (eixos). Para isso, três matrizes ( ) devem ser multiplicadas devido ao uso de alguns eixos intermediários. Pensei em duas abordagens para resolver isso:3 × 33×33\times3 Método # 1 : Fazendo a multiplicação diretamente, ou seja, vf= …

10 matrix  performance  matlab 

Estou tentando combinar itens. Dado um conjunto de itens, posso classificar em uma escala de 0 a 100 de quão semelhantes eles são entre si. Por exemplo, se o item é leite e o item também é leite, a semelhança entre e seria 100%. Se for refrigerante, a similaridade entre …

10 regression  data-analysis 

Estou tentando escrever uma implementação SVM completa em Python e tenho alguns problemas ao calcular os coeficientes de Lagrange. Primeiro, deixe-me reformular o que entendo no algoritmo para garantir que estou no caminho certo. Se x1,x2,...,xnx1,x2,...,xnx_1, x_2, ..., x_n é um conjunto de dados e yi∈{−1,1}yi∈{−1,1}y_i \in \{-1, 1\} é …

10 optimization  machine-learning  support-vector-machines  quadratic-programming 

Deixe- Ser um triângulo e deixe f ser uma função suave em T .TTTfffTTT Podemos usar ponto médio quadratura , em que X M é a-ponto médio de T .∫fdx ≈ | T| ⋅f( xM)∫fdx≈|T|⋅f(xM)\int f dx \approx |T|\cdot f(x_M)xMxMx_MTTT Você pode me fornecer (uma referência para) fórmulas de ordem …

10 quadrature 

Para otimização, da Wikipedia : Na ciência da computação, a metaheurística designa um método computacional que otimiza um problema, tentando iterativamente melhorar uma solução candidata em relação a uma determinada medida de qualidade. As metaheurísticas fazem poucas ou nenhuma suposição sobre o problema que está sendo otimizado e podem pesquisar …

10 optimization  heuristics 

Eu estou querendo saber como encontrar os autovalores de alguma matriz esparsa em determinado intervalo [a, b] pelo método iterativo. Para meu entendimento pessoal, é mais óbvio usar o método do subespaço de Krylov para encontrar os valores próprios extremos e não os interiores.

10 linear-algebra 

Seja um domínio Lipschitz delimitado poligonalmente convexo em , seja .R 2 f ∈ L 2 ( Ω )ΩΩ\OmegaR2R2\mathbb R^2f∈L2(Ω)f∈L2(Ω)f \in L^2(\Omega) Então a solução do problema de Dirichlet em , em tem uma solução única em e está bem posicionada, ou seja, para constante , temos .Ω traço u …

10 pde  finite-element  error-estimation