«floating-point» 태그된 질문

볼록하지 않은 전역 최적화 문제를 해결해야 할 몇 가지 문제가 있습니다. 현재 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 

Google에 어려운 질문이 있습니다 ( 모든 컴퓨터 과학자가 부동 소수점 산술 용지에 대해 알아야 할 표준 외에도 ). 과 대신에 log1p또는 같은 기능을 사용해야하는 경우는 언제 입니까? 언제 사용해서는 안됩니까? 이러한 기능의 다른 구현은 사용법 측면에서 어떻게 다릅니 까?expm1logexp

30 floating-point 

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

25 software  floating-point  accuracy 

Markov Chain Monte Carlo의 가장 일반적인 문제는 지수 지수가 큰 컴퓨팅 확률과 관련이 있습니다. ea1+ea2+...ea1+ea2+... e^{a_1} + e^{a_2} + ... 여기서의 구성 요소 매우 큰 매우 작은 범위 일 수있다. 내 접근 방식은 가장 큰 지수 항 하여 다음과 같습니다.aaaK:=maxi(ai)K:=maxi(ai)K := \max_{i}(a_{i}) a′=K+log(ea1−K+ea2−K+...)a′=K+log(ea1−K+ea2−K+...)a' =K + log\left( e^{a_1 - K} + …

24 monte-carlo  floating-point  statistics 

양수 유한 시퀀스가 ​​증가 합니다. 다음 두 알고리즘 중 숫자의 합을 계산하는 데 더 적합한 알고리즘은 무엇입니까?지1, z2, . . . . . 지엔z1,z2,.....znz_{1} ,z_{2},.....z_{n} s=0; for \ i=1:n s=s + z_{i} ; end 또는: s=0; for \ i=1:n s=s + z_{n-i+1} ; end 제 생각에는 오류가 점점 작아지기 때문에 …

22 algorithms  floating-point 

여러 부동 소수점 값 (단일 또는 이중)을 입력으로 사용하고 계산을 수행하며 출력 부동 소수점 값 (단일 또는 이중)을 생성하는 함수가 있다고 가정합니다. 저는 주로 MSVC 2008과 함께 일하고 있지만 MinGW / GCC 와도 함께 일할 계획입니다. C ++로 프로그래밍 중입니다. 결과적으로 얼마나 많은 오류가 발생하는지 프로그래밍 방식으로 측정하는 일반적인 방법은 …

20 floating-point  c++ 

상대적 오류 가 낮은 배정 밀도 부동 소수점에서 다음 함수를 구현하려고합니다 . logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log⁡(exp⁡(x)+exp⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) 이는 통계적 응용 프로그램에서 광범위하게 사용되어 로그 공간에 표시되는 확률 또는 확률 밀도를 추가합니다. 물론 또는 는 쉽게 오버플로 또는 언더 플로가 될 수 있습니다. 로그 공간은 처음에 언더 플로를 피하기 위해 사용되므로 …

18 floating-point  stability  numerics 

두 가지 다른 BLAS 구현이 주어지면 정확히 동일한 부동 소수점 계산을 수행하고 동일한 결과를 반환 할 것으로 기대할 수 있습니까? 또는 예를 들어 스칼라 곱을 하고 다른 하나는 계산하여 IEEE 부동 소수점에서 다른 결과를 얻을 수 산수?( ( x1와이1+x2와이2) + x삼와이삼)+ x4와이4((엑스1와이1+엑스2와이2)+엑스삼와이삼)+엑스4와이4 ((x_1y_1 + x_2y_2) + x_3y_3) + x_4y_4 ( …

17 floating-point  blas 

24 개의 Intel Xeon CPU를 포함하는 노드로 구성된 Ubuntu Linux 클러스터에서 MD (Molecular Dynamics) 코드 GROMACS 를 실행했습니다 . 내 관심의 대상은 부동 소수점 산술 정밀도에 다소 민감하다는 것이 밝혀 졌으므로 배정 밀도의 높은 계산 비용에도 불구하고 단 정밀도가 아닌 배정 밀도로 GROMACS를 실행해야했습니다. 클러스터에서 GROMACS를 배정도로 컴파일했습니다. CPU에 비해 …

17 hpc  floating-point  molecular-dynamics  gpu 

하자 x, y두 개의 부동 소수점 숫자합니다. 평균을 계산하는 올바른 방법은 무엇입니까? 순진 방법은 (x+y)/2오버 플로우가 발생할 때 수 x와 y너무 큽니다. 내 생각 엔 0.5 * x + 0.5 * y아마도 더 나은,하지만 (아마도 비효율적 인)이 곱셈을 포함, 그것은 충분히 좋은 경우에 나는 확실하지 않다. 더 좋은 방법이 있습니까? …

15 floating-point 

컴퓨터가 할 수있는 이론적 인 수의 FLOP (Floating Point Operations)를 결정하고 싶습니다. 누군가 나를 도와주세요. (컴퓨터와 다른 슈퍼 컴퓨터의 차이점을 이해하기 위해 일부 슈퍼 컴퓨터와 비교하고 싶습니다)

15 floating-point 

함수를 숫자로 구현하는 데 문제가 있습니다. 큰 입력 값에서 결과는 매우 많은 수의 매우 작은 수라는 사실로 인해 어려움을 겪습니다. 치명적인 취소가 올바른 용어인지 확실하지 않으므로 수정하십시오. 무언가 잘못되었다는 증거 : 6의 더 큰 입력에 대해 진동과 0.0의 어설 션을 어떻게 피할 수 있습니까? 내 기능은 다음과 같습니다. import numpy …

13 python  numerical-analysis  floating-point 

Allen과 Tildesley 의 액체 컴퓨터 시뮬레이션 책을 읽었습니다 . 71 페이지부터 저자는 분자 역학 (MD) 시뮬레이션에서 뉴턴의 운동 방정식을 통합하는 데 사용되는 다양한 알고리즘에 대해 설명합니다. 78 페이지에서 시작하여 저자는 MD의 표준 통합 알고리즘 인 Verlet 알고리즘에 대해 설명합니다. 그들은 말한다 : 운동 방정식을 통합하는 데 가장 널리 사용되는 방법은 …

13 floating-point  computer-arithmetic 

간단한 함수에서 FLOP의 수를 평가할 때 종종 기본 산술 연산자를 계산하는 식으로 내려갈 수 있습니다. 그러나 짝수 나누기가 포함 된 수학적 진술의 경우이를 수행 할 수 없으며 덧셈과 곱셈 만있는 함수의 FLOP 카운트와 비교할 수있을 것으로 기대합니다. 작업이 라이브러리에서 구현되면 상황이 더욱 악화됩니다. 따라서 특수 기능의 성능에 대한 합리적인 개념을 …

13 performance  complexity  floating-point 

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

11 floating-point  accuracy