전산 과학

재현 가능한 계산 연구는 다른 논문에서 계산 논문의 결과를 생성하는 데 필요한 코드를 다른 연구원들이 이용할 수 있도록하여이 코드를 실행하여 해당 논문의 결과를 재현 할 수 있도록하는 것을 목표로합니다. 두 논문은 내가 사용에 대한 사내 자동 차별화 패키지 (라고하고 있어요 : 내 연구의 재현을 모두하고 싶습니다,하지만 난하다가 조금 실행 해요 …

19 software  publications  documentation  reproducibility 

주어진 연산 / 알고리즘의 1, 2, 3 차원 버전에 대해 매우 유사한 코드를 작성하는 경우가 종종 있습니다. 이 모든 버전을 유지 관리하는 것은 지루할 수 있습니다. 간단한 코드 생성은 잘 작동하지만 더 나은 방법이 있어야한다고 생각되는 것 같습니다. 작업을 한 번 작성하고 더 높거나 더 작은 차원으로 일반화하는 비교적 간단한 …

19 algorithms  software  fortran 

원하는 상태 와이0와이0y_0 과 정규화 매개 변수 β∈ Rβ∈아르 자형\beta \in \mathbb R 을 고려하여 , 기능 을 최소화하기 위해 상태 와이와이y 와 제어 유유u 를 찾는 문제를 고려하십시오.Ay=u입니다. y,y0,u∈RnA∈Rn×n12∥ y− y0∥2+ β2∥ u ∥212”와이−와이0”2+β2”유”2\begin{equation} \frac{1}{2} \| y - y_0 \|^2 + \frac{\beta}{2} \| u \|^2 \end{equation}Y= u .ㅏ와이=유.\begin{equation} Ay …

19 optimization  reference-request  time-integration 

n × n G n n GAAA 는 대칭 양수 한정 (SPD) 희소 행렬입니다. 는 희소 대각선 행렬입니다. 은 크고 ( > 10000) 이 아닌 숫자 는 일반적으로 100 ~ 1000입니다.n×nn×nn \times nGGGnnnnnnGGG L D L TAAA 는 cholesky 형식에서 로 인수 분해되었습니다 .LDLTLDLTLDL^T 가 가 될 때 과 효율적 …

19 linear-algebra 

반복적 방법을 사용 대형 선형 시스템 A x = b 의 솔루션의 경우, 사전 조건을 도입하는 것이 종종 중요합니다. 예를 들어 대신 M - 1 ( A x = b )를 해결 하십시오. 여기서 M 은 시스템의 왼쪽 사전 조건에 사용됩니다. 일반적으로, 우리는 M - 1 ≈ A - 1을 …

19 linear-algebra  iterative-method  preconditioning 

PDE 솔버에 대한 몇 가지 문헌을 읽는 동안 나는 오늘 의사 시간 단계 라는 용어를 만났습니다 . 일반적인 용어 인 것 같지만 좋은 정의 나 소개 기사를 찾지 못했습니다. 따라서 의사 시간 스텝핑이란 무엇이며 일반적으로 어떻게 사용됩니까?

18 pde  time-integration 

FMM (Fast Multipole Method)의 대부분 (모두?) 구현에서 octree는 관련 도메인을 분해하는 데 사용됩니다. 이론적으로 octree는 간단한 체적 경계를 제공하며, 이는 FMM의 O (n) 런타임을 증명하는 데 유용합니다. 이 이론적 근거를 넘어서 다른 트리 또는 트리 데이터 구조보다 Octree를 사용하면 이점이 있습니까? 셀이 바로 인접한 이웃을 알기 때문에 octree를 사용하면 상호 …

18 algorithms 

상대적 오류 가 낮은 배정 밀도 부동 소수점에서 다음 함수를 구현하려고합니다 . 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 

비 구조적 쿼드 메시를 2D로 생성하는 데 가장 적합한 (확장 성 및 효율성) 알고리즘은 무엇입니까? 좋은 구조화되지 않은 쿼드 메쉬 생성기를 어디에서 찾을 수 있습니까? (오픈 소스 선호)

18 mesh-generation  computational-geometry 

웨이블릿 방법을 PDE에 적용하는 방법을 배우고 싶지만 불행히도이 주제에 대해 배울 수있는 좋은 자료는 없습니다. 웨이블릿에 대한 많은 소개는 예를 들어, 바람직하게는 적은 웨이블릿의 중첩에 의해 신호를 조립하는 보간 이론에 초점을 둔 것으로 보인다. PDE에 대한 응용은 때때로 그 주제에 대해 깊이 들어 가지 않고 언급됩니다. 나는 WFT를 보았지만 그 …

18 pde  wavelet 

f2py현대의 포트란과 함께 사용 하고 싶습니다 . 특히 다음 기본 예제를 작동 시키려고합니다. 이것은 내가 생성 할 수있는 가장 작은 유용한 예입니다. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8, intent(out) :: z(:) ! Variable Declarations ! real*8, allocatable :: y(:) integer …

18 python  fortran 

압력, 밀도, 에너지, 온도 및 농도와 같은 물리량은 항상 양수 여야하지만 숫자 방법은 때때로 솔루션 프로세스 중에 음수 값을 계산합니다. 방정식이 복잡하거나 무한한 값 (일반적으로 코드 충돌)을 계산하기 때문에 이것은 좋지 않습니다. 이러한 양이 양을 유지하기 위해 어떤 수치 방법을 사용할 수 있습니까? 다음 중 가장 효율적인 방법은 무엇입니까?

18 pde  fluid-dynamics  hyperbolic-pde 

이 게시물을 개선하고 싶습니까? 인용과 답변이 올바른 이유에 대한 설명을 포함하여이 질문에 대한 자세한 답변을 제공하십시오. 세부 사항이없는 답변은 편집하거나 삭제할 수 있습니다. AMR ( Adaptive Mesh Refinement )은 PDE의 수치 솔루션에서 광범위하게 변화하는 공간 스케일 문제를 처리하기위한 일반적인 기술입니다. 구조화 된 그리드에 AMR에 대해 어떤 범용 라이브러리가 있습니까? 라이브러리가 …

18 pde  libraries  parallel-computing  adaptive-mesh-refinement 

GPU (CUDA)에서 제한되지 않은 비선형 최적화 문제를 해결하려고합니다. 목적 함수는 부드러운 비선형 함수이며 기울기는 분석적으로 계산하기에 상대적으로 저렴하므로 수치 근사로 신경 쓸 필요가 없습니다. 대부분의 fp32 수학 연산 으로이 문제를 해결하고 싶습니다 (다양한 이유로). 비선형 최적화 방법은 반올림 오류에 대해 더 강력하지만 성능이 우수합니까? (예를 들어 켤레 그라디언트 / 쿼시 …

18 optimization  cuda 

재귀 적으로 다음과 같이 정의 된 행렬 방정식이 있다고 가정합니다. A[n] = inverse([1 - b[n]A[n+1]]) * a[n] 그런 다음 A [1]에 대한 방정식은 지루한 재 계산을 피하는 매우 효율적인 방법이있는 연속 분수와 유사합니다 (일부 예는 "수치 레시피"참조). 그러나 계수 b [n] 및 a [n]을 행렬로 만드는 유사한 방법이 있는지 궁금합니다. …

18 algorithms