전산 과학

현재 일부 특정 알고리즘에 의해 생성 된 매우 큰 대칭 (그러나 양의 한정이 아닌) 시스템을 해결하기 위해 노력하고 있습니다. 이 행렬에는 병렬 해결에 사용할 수있는 훌륭한 블록 희소성이 있습니다. 그러나 직접 접근 방식 (예 : 다중 정면) 또는 반복적 접근 방식 (사전 조건 GMRES 또는 MINRES)을 사용 해야하는지 결정할 수 …

10 parallel-computing  linear-solver  sparse-matrix  linear-system  gmres 

차등 형식 이론에 대해 오늘 읽은 후, 2 차 유한 체적 법 (FVM)을 생각 나게하는 것이 얼마나 큰지 감명 받았습니다. 나는이 방법을 사소한 것으로 생각하거나 더 깊은 연결이 있는지 알아 내려고 고심하고 있습니다. 음, 차등 형태는 표면을 통한 유체의 흐름과 같이 2 차 FVM에 깊이 뿌리 박은 일부 개념을 일반화하는 …

10 fluid-dynamics  finite-volume 

비 부동 소수점 컴퓨팅 라이브러리 / 패키지는 거의 없습니다. 왜 거기에 부동 소수점 표현의 다양한 부정확성을 감안할 때, 문제가 발생하지 적어도 몇몇 이 높은 정확도의 가치를 고정 소수점 작업의 복잡성을 수 있습니다 필드. 고정 소수점 고유 값 솔버를 사용하는 데 큰 어려움이 있습니까? 얼마나 느리고 빠르며 부정확하고 정확합니까? 관련 : …

10 floating-point  numerics 

PETSc 패키지를 처음 사용합니다. 매트릭스 시장 형식으로 ~ 4000x4000 매트릭스 A가 있으며 PETSc가 여러 프로세서를 사용 하여이 문제를 해결하도록하고 싶습니다. 단일 프로세서에서 시스템을 해결하는 방법을 알고 있지만 매트릭스와 벡터를 다른 프로세서로 분산시키는 방법을 모르겠습니다. 이를 수행하기위한 간단한 지침이 있습니까?

10 petsc  matrix 

문제는 최대 f( x ) A x = b에 종속 최대에프(엑스) 에 따라 ㅏ엑스=비\max f(\mathbf{x}) \text{ subject to } \mathbf{Ax} = \mathbf{b} 여기서 에프( x ) = ∑엔나는 = 11 + x4나는( ∑엔나는 = 1엑스2나는)2−−−−−−−−−−√에프(엑스)=∑나는=1엔1+엑스나는4(∑나는=1엔엑스나는2)2f(\mathbf{x}) = \sum_{i=1}^N\sqrt{1+\frac{x_i^4}{(\sum_{i=1}^{N}x_i^2)^2}} , x =[ x1, x2, . . . , x엔]티∈ R엔× 1엑스=[엑스1,엑스2,...,엑스엔]티∈아르 자형엔×1\mathbf{x} …

10 optimization 

다음과 같은 최적화 문제가 있습니다. minJ,Bs.t.∑ij|Jij|MJ+BY=XminJ,B∑ij|Jij|s.t.MJ+BY=X \begin{array}{rl} \min_{J,B} & \sum_{ij} |J_{ij}|\\ \textrm{s.t.} & MJ + BY =X \end{array} 여기서 내 변수는 행렬 JJJ 와 BBB 이지만 전체 문제는 여전히 선형 프로그램입니다. 나머지 변수는 고정되어 있습니다. 이 프로그램을 자주 사용하는 선형 프로그래밍 도구에 입력하려고하면 문제가 발생합니다. 즉, 이것을 "표준"선형 프로그램 형식으로 …

10 optimization  convex-optimization  linear-programming 

실제 물리적 시스템을 나타내는 수치 모델이 있고 혼돈을 나타내는 수치 모델 (예 : 유체 역학 모델, 기후 모델)이 모델이 제대로 수행되고 있음을 어떻게 알 수 있습니까? 초기 조건의 작은 변화만으로도 개별 시뮬레이션의 출력이 크게 변하기 때문에 두 세트의 모델 출력을 직접 비교할 수 없습니다. 관측치의 초기 조건을 충분히 자세하게 알 …

10 testing  numerics 

CPU, 메모리, 네트워크 및 스토리지를 통한 부동 소수점 연산 및 데이터 전송의 현재 및 미래 에너지 비용을 모델링하기 위해 가장 중요한 문헌 및 슬라이드 참조를 찾고 있습니다. 이 질문을 커뮤니티 위키로 표시했으며 각 답변을 다음 형식으로 제한하는 것이 좋습니다. 제목, 저자, 위치 / 컨퍼런스 / 저널, 연도, DOI / URL, …

10 performance  architecture  exascale 

함수 를 최대화하는 데 관심이 있는데 여기서 θ ∈ R p 입니다.f(θ)f(θ)f(\mathbf \theta)θ∈Rpθ∈Rp\theta \in \mathbb R^p 문제는 함수 또는 파생 함수의 분석 형식을 모른다는 것입니다. 내가 할 수있는 유일한 방법은 값에 연결하여, 포인트 - 현명한 기능을 평가하는 것입니다 잡음 추정 얻을 F를 ( θ * ) 그 시점에서. 원하는 경우 …

10 optimization  monte-carlo  simulation 

배정도 특수 기능을 구현하는 최첨단 방법은 무엇입니까? 다음 적분이 필요합니다. 대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 

최근 몇 가지 형태의 범용 데이터 중심 공유 메모리 병렬 처리를 제공하는 여러 라이브러리 / 소프트웨어 프로젝트가 등장했습니다. 주요 아이디어는 프로그래머가 명시 적으로 스레드 된 코드를 작성하는 대신 알고리즘을 상호 종속적 인 작업으로 구현 한 다음 공유 메모리 시스템의 범용 미들웨어에 의해 동적으로 스케줄됩니다. 이러한 라이브러리의 예는 다음과 같습니다. QUARK …

10 parallel-computing  software  hpc  libraries  multicore 

나는 수치 기능이 f(x, y)구현 몇 가지 공식과 나는 그것이 모든 매개 변수 조합에 대한 분석 표현에 대한 올바른지 확인하려는 이중 부동 소수점 숫자 반환 x과 y나는에 관심이 있습니다. 무엇을 계산하고 비교하는 적절한 방법입니다 분석 부동 소수점 숫자? 두 숫자가 a및 라고 가정 해 보겠습니다 b. 지금까지 절대 abs(a-b) < …

10 floating-point 

쿼트 방정식의 해에 대한 개방형 C 구현이 있습니까? ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0 페라리 솔루션의 구현을 생각하고 있습니다. Wikipedia에서 솔루션이 계수의 가능한 부호 조합 중 일부에 대해서만 계산적으로 안정적이라는 것을 읽었습니다. 그러나 운이 좋을 수도 있습니다 ... 컴퓨터 대수 시스템을 사용하여 분석적으로 해결하고 C로 내 보내서 실용적인 솔루션을 얻었습니다. 그러나 테스트를 거친 구현이 있으면 …

10 polynomials  nonlinear-equations  roots 

비동기식 통신 패턴으로 MPI 응용 프로그램을 최적화하려고합니다. 각 순위에는 계산할 항목 목록이 있으며 입력 또는 출력이 다른 순위에있는 경우 필요에 따라 메시지를 보냅니다. 또한 각 순위는 스레드되어 있습니다 (현재는 하나의 통신 스레드와 5 명의 작업자가 있음). 코드의 다른 성능 결정 부분에 대한 타이머를 사용하여 코드를 계측하여 각 스레드에 대한 (시작, …

10 mpi  visualization 

공학을 전공하고 컴퓨터 과학 프로젝트를 진행하고 있습니다. 나는 당신의 자신의 알고리즘을 구현하는 장점과 단점에 대한 논의가 있었다 것으로 나타났습니다 여기 . 초보 대학원 생일 때 자체 알고리즘 구현과 라이브러리 사용에 대한 비용-이익 분석이 다른가?

10 software  education