Perguntas com a marcação «quadrature»

Eu tenho vários problemas desafiadores de otimização global não convexa para resolver. Atualmente, uso o Optimization Toolbox do MATLAB (especificamente, fmincon()com o algoritmo = 'sqp'), o que é bastante eficaz . No entanto, a maior parte do meu código está em Python, e eu adoraria fazer a otimização em Python …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

Preciso avaliar numericamente a integral abaixo: ∫∞0 0s i n c′( x r ) r E( r )----√dr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr onde , e . Aqui é a função de Bessel modificada do segundo tipo. No meu caso particular, tenho , e .x∈R+λ,κ,ν>0Kλ=0,00313κ=0,00825ν=0,33E( r ) = r4( λ κ2+ …

25 matlab  quadrature  special-functions 

Qual é o estado da arte na aproximação de integrais altamente oscilatórias em uma dimensão e dimensões superiores à precisão arbitrária?

23 quadrature  precision  numerical-analysis 

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 

Gostaria de saber se existe uma maneira rápida de calcular a distância euclidiana de dois vetores na oitava. Parece que não há função especial para isso, então devo usar apenas a fórmula com sqrt?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Ouvi de maneira anedótica que quando alguém está tentando fazer numericamente uma integral do formulário ∫∞0f(x)J0(x)dx∫0∞f(x)J0(x)dx\int_0^\infty f(x) J_0(x)\,\mathrm{d}x com suave e bem-comportado (por exemplo, não seja altamente oscilatório, não singular, etc.), ajudará a precisão a reescrevê-lo comof(x)f(x)f(x) 1π∫π0∫∞0f(x)cos(xsinθ)dxdθ1π∫0π∫0∞f(x)cos⁡(xsin⁡θ)dxdθ\frac{1}{\pi}\int_0^\pi \int_0^\infty f(x) \cos(x\sin\theta) \,\mathrm{d}x\,\mathrm{d}\theta e execute a integral interna numericamente primeiro. Não …

15 quadrature  special-functions 

Estou com uma certa frustração com a maneira como o matlab lida com a integração numérica versus o Scipy. Observo as seguintes diferenças no meu código de teste abaixo: A versão do Matlab é executada em média 24 vezes mais rápido que o meu equivalente em python! A versão do …

13 matlab  python  quadrature  numba 

Eu tenho uma função tal que seja finita e quero aproximar essa integral. ∫ R 3 f ( x , y , z ) d Vf( x , y, z)f(x,y,z)f(x,y,z) ∫R3f( x , y, z) dV∫R3f(x,y,z)dV\int_{R^3} f(x,y,z)dV Eu estou familiarizado com regras de quadratura e aproximações de integrais de monte …

13 monte-carlo  quadrature 

Apenas as duas primeiras seções desta longa questão são essenciais. Os outros são apenas para ilustração. fundo Quadraturas avançadas, como Newton-Cotes, Gauß-Legendre e Romberg, parecem ser principalmente destinadas a casos em que é possível provar a função com precisão, mas não integrar analiticamente. No entanto, para funções com estruturas mais …

12 reference-request  quadrature 

Existem várias famílias de métodos para quadratura numérica. Se eu tiver uma classe específica de integrandos, como seleciono o método ideal? Quais são as perguntas relevantes a serem feitas sobre o integrando (por exemplo, é suave? Possui singularidades?) E o problema computacional (por exemplo, tolerância a erros, orçamento computacional)? Como …

12 quadrature 

Estou lidando com uma integral complicada que exibe NaNs em certos valores próximos de zero e, no momento, estou lidando com eles de maneira bastante grosseira, usando uma instrução ISNAN que define o integrando como zero quando isso ocorre. Eu tentei isso com a biblioteca NMS em FORTRAN (a rotina …

12 quadrature 

Quero integrar uma expressão polinomial sobre um elemento de 4 nós em 3D. Vários livros sobre FEA cobrem o caso em que a integração é executada sobre um elemento plano arbitrário de 4 neros. O procedimento usual neste caso é encontrar a matriz de Jacobi e usá-lo é determinante para …

12 finite-element  quadrature 

Eu tenho uma integral incorreta (bidimensional) I=∫AW(x,y)F(x,y)dxdyI=∫AW(x,y)F(x,y)dxdyI=\int_A \frac{W(x,y)}{F(x,y)}\,\mbox{d}x\mbox{d}y onde o domínio de integração é menor que , mas ainda mais restrito por . Como e são suaves ex = [ - 1 , 1 ] y = [ - 1 , 1 ] F ( x , y ) > …

12 monte-carlo  quadrature 

Em deste curso avançado sobre aplicações da teoria de função complexa em um ponto em um exercício integral altamente oscilatório I(λ)=∫∞−∞cos(λcosx)sinxxdxI(λ)=∫−∞∞cos⁡(λcos⁡x)sin⁡xxdxI(\lambda)=\int_{-\infty}^{\infty} \cos (\lambda \cos x) \frac{\sin x}{x} d x deve ser aproximado para grandes valores de usando o método do ponto de sela no plano complexo.λλ\lambda Devido à sua natureza …

11 quadrature  integration  oscillations 

Digamos que eu tenho uma função que desejo integrar em um tetraedro . Se fosse arbitrário, a quadratura de Gauss seria uma boa solução, mas por acaso sei que é harmônico. Quanto a quadratura de Gauss pode ser acelerada usando essas informações?f:R3→Rf:R3→Rf : \mathbf{R}^3 \to \mathbf{R}T⊂R3T⊂R3T \subset \mathbf{R}^3ffffff Por exemplo, …

11 quadrature