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 superior em um...
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 e seria, talvez, 80% e assim...
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...
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...
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...
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 1000 com 100 valores. Ao esquadrar essas...
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
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...
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...
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...
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,...
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 5 segundos de espera...
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...
Existe uma implementação C aberta para a solução de equações quárticas: ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0 Estou pensando em uma implementação da solução da Ferrari. Na Wikipedia, li que a solução é estável em termos computacionais apenas para algumas das combinações possíveis...
Considere o seguinte problema que o termo forçado pode depender de (consulte o Edit 1 abaixo para a formulação) e e seus primeiros derivados. Esta é uma equação de onda dimensional 1 + 1. Temos dados iniciais prescritos em .Wuv=FWuv=F W_{uv} = F u,vu,vu,vWWW{u+v=0}{u+v=0}\{u+v = 0\} Estou...
Estou tentando otimizar um aplicativo MPI com um padrão de comunicação altamente assíncrono. Cada classificação possui uma lista de itens a serem calculados e envia mensagens conforme necessário, se as entradas ou saídas residirem em uma classificação diferente. Além disso, cada classificação é...
Nos últimos anos, surgiram várias bibliotecas / projetos de software que oferecem uma forma ou outra de paralelismo de memória compartilhada de uso geral orientado a dados. A idéia principal é que, em vez de escrever um código explicitamente encadeado, os programadores implementam seus algoritmos...
Qual é a maneira mais avançada de implementar funções especiais de precisão dupla? Preciso da seguinte integral: para e , que pode ser escrita em termos da função gama incompleta inferior. Aqui está minha implementação de Fortran e C: m=0,1,2,. . . t>0Fm( t ) = ∫1 10 0você2 me- t u2du = γ( m +...
Eu tenho uma função numérica f(x, y)retornando um número de ponto flutuante duplo que implementa alguma fórmula e quero verificar se ele está correto em relação às expressões analíticas para todas as combinações de parâmetros xe no yqual estou interessado. Qual é a maneira correta de comparar os...
Estou interessado em maximizar uma função , onde .θ ∈ R pf( θ )f(θ)f(\mathbf \theta)θ ∈ Rpθ∈Rp\theta \in \mathbb R^p O problema é que não conheço a forma analítica da função ou de suas derivadas. A única coisa que posso fazer é avaliar a função ponto a ponto, um valor e obtendo uma estimativa...