«accuracy» 태그된 질문

볼록하지 않은 전역 최적화 문제를 해결해야 할 몇 가지 문제가 있습니다. 현재 MATLAB의 Optimization Toolbox (특히 fmincon()algorithm = 사용 'sqp')를 사용하고 있습니다. 그러나 내 코드의 대부분은 Python이며 Python에서도 최적화를하고 싶습니다. 파이썬 바인딩과 경쟁 할 수있는 NLP 솔버가 fmincon()있습니까? 반드시 비선형 평등과 불평등 제약을 다룰 수있다 사용자가 Jacobian을 제공하도록 요구하지 않습니다. …

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 

실제 변수의 실제 기능이 주어지면 IEEE 754 산술이 장착 된 기계에서 모든 입력에 대한 함수를 계산하기 위해 숫자로 정확한 코드를 자동으로 생성 할 수있는 소프트웨어가 있습니까? 예를 들어 평가할 실제 함수가 다음과 같은 경우 소프트웨어는 계산 정확도의 손실을 피하기 위해 특정 입력 세트에 대한 치명적인 취소 및 출력 테이블 조회를 …

25 software  floating-point  accuracy 

CFD에 대한 많은 수치 적 접근은 임의로 높은 순서로 확장 될 수 있습니다 (예를 들어, 불연속 Galerkin 방법, WENO 방법, 스펙트럼 차이 등). 주어진 문제에 대해 적절한 정확도 순서를 어떻게 선택해야합니까?

23 fluid-dynamics  accuracy 

선형 시스템을 해결하기 위해 행렬을 뒤집는 것은 시스템을 직접 해결하거나 LU, Cholesky 또는 QR 분해를 사용하는 것만 큼 정확하고 효율적이지 않기 때문에 좋은 생각이 아니라는 것을 알고 있습니다. 그러나 실제 사례로는 이것을 확인할 수 없었습니다. 이 코드를 시도했습니다 (MATLAB) M = 500; A = rand(M,M); A = real(expm(1i*(A+A.'))); b = …

16 matrix  accuracy  inverse 

대부분의 수학 라이브러리에는 여러 버전의 로그 함수가 있습니다. 대부분의 경우 우리는 그것들이 완벽하다고 가정하지만 실제로 상당히 많은 숫자는 특정 숫자의 정밀도를 제공합니다. 일부 기능의 경우 수치 적으로 더 안정적인 변형이 있습니다. 예를 들어, Fortran, R, Java 및 C는 (소수의 x 값에 대해 더 높은 정밀도를 제공하는) Math.log1p컴퓨팅 log(1.0+x)및 대응을 위해 …

11 floating-point  accuracy 

유형의 방정식을 풀려고합니다. ( − ∂2∂엑스2− f( X ) ) ψ(X)=λψ(X)(−∂2∂엑스2−에프(엑스))ψ(엑스)=λψ(엑스) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) 여기서 에프( x )에프(엑스)f(x) 는 가장 작은 N 고유 값과 고유 벡터에 대해 간단한 극점을 갖습니다 . 경계 조건은 ψ ( 0 ) = 0 및 ψ ( R …

11 finite-difference  ode  accuracy 

배정도 특수 기능을 구현하는 최첨단 방법은 무엇입니까? 다음 적분이 필요합니다. 대m=0,1,2,. . . 및t>0, 이는 불완전한 감마 함수의 하에서 작성 될 수있다. 다음은 Fortran 및 C 구현입니다.에프미디엄( t ) = ∫10유2 M이자형− t u2디u = γ( m + 12, t )2 톤m + 12Fm(t)=∫01u2me−tu2du=γ(m+12,t)2tm+12 F_m(t) = \int_0^1 u^{2m} e^{-tu^2} d …

10 efficiency  accuracy  special-functions 

비 물리적 모드를 도입하지 않고 8 가우스 포인트 미만의 육면체 유한 요소에 대해 2 차 정확도를 얻을 수 있습니까? 단일 중앙 가우스 포인트는 비 물리적 전단 모드를 도입하고 8 가우스 포인트의 표준 대칭 배열은 사면체 이산화에 비해 비쌉니다. 편집 : 누군가가 방정식을 요구했습니다. 내가 관심있는 방정식은 동적 또는 준 정적 …

10 finite-element  accuracy 

통합하려고합니다 ∫10t2n+2exp(αr0t)dt∫01t2n+2exp⁡(αr0t)dt\int^1_0 t^{2n+2}\exp\left({\frac{\alpha r_0}{t}}\right)dt 이것은 간단한 변환입니다 ∫∞1x2nexp(−αr0x)dx∫1∞x2nexp⁡(−αr0x)dx\int^{\infty}_1 x^{2n}\exp(-\alpha r_0 x)dx 사용 t=1xt=1xt = \frac1{x}부적절한 적분을 수치 적으로 근사하기가 어렵 기 때문입니다. 그러나 이로 인해 새로운 정수를 거의 0으로 평가하는 문제가 발생합니다. 간격이 길이가 1이므로 적절한 구적 노드 수를 얻는 것이 매우 쉽습니다 (따라서 비슷한dtdtdt 아주 작게 만들 수 있습니다), …

9 quadrature  accuracy