전산 과학

C ++로 과학 코드를 작성한다고 가정 해 봅시다. 최근 동료와의 토론에서 표현 템플릿은 정말 나쁜 것일 수 있으며 특정 버전의 gcc에서만 소프트웨어를 컴파일 할 수 있다고 주장했습니다. 아마도이 문제는 이 몰락 패러디의 자막에서 암시 된 것처럼 몇 가지 과학적 코드에 영향을 미쳤습니다 . (이것들은 내가 아는 유일한 예이므로 링크입니다.) 그러나 …

24 performance  c++  automatic-differentiation  interval-arithmetic 

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 

강력한 형식의 PDE에서 FEM 형식으로 갈 때 항상 변형 형식을 지정하여이 작업을 수행해야합니다. 이를 위해 일부 (Sobolev) 공간의 요소에 강력한 형식을 곱하고 해당 지역에 통합하십시오. 나는 이것을 받아 들일 수있다. 내가 이해하지 못하는 것은 왜 Green의 공식을 사용해야합니까 (하나 또는 여러 번)입니다. 나는 주로 포아송 방정식을 사용하여 작업했습니다. 예를 들어 …

24 pde  finite-element  elliptic-pde 

제가 알 수있는 한, 미국의 에너지 에너지 전산 과학 소프트웨어 프레임 워크는 PETSc 와 Trilinos 입니다. 언어 (C와 C ++)의 차이를 넘어 언뜻 보면 비슷해 보입니다. 두 프레임 워크의 주요 차이점은 무엇이며 다른 요소를 선택하는 데 영향을 미치는 요소는 무엇입니까? (기관의 편견과 기존 인프라는 무시하십시오.)

24 petsc  software  trilinos 

실험 설명 : 라그랑주 보간에서 정확한 방정식은 NNN 점 (다항식 차 N−1N−1N - 1 ) 에서 샘플링 되고 101 점에서 보간됩니다. 여기서 NNN 은 2에서 64까지 다양하다. 매번 L1L1L_1 , L2L2L_2 및 L∞L∞L_\infty 에러 플롯이 준비된다. 함수가 등 간격 포인트에서 샘플링 될 때, 에러는 초기에 떨어지고 ( NNN 이 약 …

24 finite-difference  interpolation  error-estimation  polynomials 

내가 이해하는 것처럼 선형 방정식 시스템을 푸는 반복적 인 방법에는 두 가지 주요 범주가 있습니다. 고정 방법 (Jacobi, Gauss-Seidel, SOR, Multigrid) Krylov Subspace 방법 (콘쥬 게이트 그라디언트, GMRES 등) 나는 대부분의 고정 방법이 푸리에 오류 모드를 반복적으로 완화 (평활)함으로써 작동한다는 것을 이해합니다. 내가 알기로, Conjugate Gradient 방법 (Krylov subspace 방법)은 …

24 linear-algebra  convergence  krylov-method  conjugate-gradient 

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

23 fluid-dynamics  accuracy 

(xi,yi)(xi,yi)(x_i,y_i)y(x)y(x)y(x)xxxf(x)≡y(x)−(ax+b)f(x)≡y(x)−(ax+b)f(x) \equiv y(x) - (ax + b)x→∞x→∞x \to \inftyf′(x)f′(x)f'(x)f′′(x)f″(x)f''(x)f (x) 의 기능적 형태 f(x)f(x)f(x)가 기본 기능으로 설명 될 수있는 경우조차도 모르겠습니다 . 내 목표는 점근 적 경사 a 의 가능한 최상의 추정치를 얻는 것입니다 aaa. 명백한 조잡한 방법은 마지막 몇 가지 데이터 포인트를 선택하고 선형 회귀를 수행하는 것입니다. 물론 f(x)f(x)f(x) …

23 algorithms 

대한 아마존 EC2 컴퓨팅 클러스터 비용 $ 800 $ (1000)는 삼년의 과정을 통해 물리적 CPU 코어 당 (듀티 사이클에 따라 다름). 우리의 마지막 하드웨어 획득 과정에서, 랩에서 Amazon의 클러스터와 비슷한 48 개의 코어 하드웨어를 코어 당 약 300 달러에 집어 들었습니다. 여기에 뭔가 빠졌습니까? 분자 역학 시뮬레이션과 같은 높은 CPU …

23 cloud-computing  hpc 

1 차원과 더 높은 차원의 고진동 적분을 임의의 정밀도로 근사화하는 데있어 최첨단 기술은 무엇입니까?

23 quadrature  precision  numerical-analysis 

컴퓨터 과학 연구에서 재현성은 점점 더 중요 해지고 있습니다. (예를 들어, 과학 에서 Roger Peng의이 기사를 참조하십시오 . 다른 기사와 웹 사이트도 알고 있습니다.) 그러나 저널 기사 (또는 온라인)에 얼마나 많은 정보를 포함시켜야하는지 명확하지 않습니다. 나의 전산 연구는 재현 가능합니다 (지적 재산권 계약과 같은 다른 장애물이 없다고 가정). 거기에 지침이 …

23 software  publications  documentation  reproducibility 

이 질문은 Computational Science Stack Exchange에서 답변을 얻을 수 있기 때문에 Cross Validated에서 마이그레이션 되었습니다. 7 년 전에 이주했습니다 . 나는 스스로 재활 약 / 분석을 수행 할 수 있고 온라인 등에서 찾은 일부 MATLAB 스크립트를 연구 / 수정할 수 있도록 MATLAB 프로그래밍을 배우고 싶습니다. 그러나 문제는 MATLAB을 감당할 수 …

23 software  matlab  languages  octave 

다음 방정식에서 유한 차분 이산 법은 무엇입니까? ∂ρ∂티+ ∇ ⋅ ( ρ u ) = 0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0 ? 우리는 1D 사례를 취할 수 있습니다. ∂ρ∂티+ d디엑스( ρ u ) = 0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} + \frac{d}{dx}\left(\rho u\right)=0 어떤 이유로 든 내가 찾을 수있는 모든 계획은 …

22 finite-difference  fluid-dynamics 

현재로서는이 질문이 Q & A 형식에 적합하지 않습니다. 답변, 사실, 참고 자료 또는 전문 지식을 통해 답변이 뒷받침 될 것으로 예상되지만이 질문은 토론, 논쟁, 여론 조사 또는 광범위한 토론을 요구할 것입니다. 이 질문을 개선하고 다시 열 수 있다고 생각 되면 도움말 센터 를 방문하여 안내를 받으십시오 . 휴일 칠년 전에 …

22 languages 

양수 유한 시퀀스가 ​​증가 합니다. 다음 두 알고리즘 중 숫자의 합을 계산하는 데 더 적합한 알고리즘은 무엇입니까?지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