전산 과학

컴퓨터를 사용하여 과학적 문제를 해결하는 과학자를위한 Q & A

2
Quantum Monte Carlo에 대한 혼란
내 질문은 이 참조 서에 설명 된 것처럼 QMC 방법에서 관찰 가능 항목을 추출하는 것에 관한 것 입니다. Path Integral Monte Carlo와 같은 다양한 QMC 방법의 공식 파생을 이해합니다. 그러나 하루가 끝날 무렵에도 이러한 기술을 효과적으로 사용하는 방법에 대해 여전히 혼란스러워합니다. Quantum MC 방법의 도출에 대한 기본 아이디어는 트로터 근사법을 …
3
암시 적 FEM과 명시 적 FEM의 차이점은 무엇입니까?
명시 적 FEM과 암시 적 FEM의 차이점은 정확히 무엇입니까? here 에 따르면 암시 적 또는 명시 적 시간 통합을 사용하는지 여부 만 다릅니다. 내가 읽은 한 권의 책에서 알 수 있듯이 암시 적 FEM은 질량이 노드에 집중되지 않는 곳입니다. 명시 적 및 암시 적 FEM의 정확한 정의는 무엇입니까?
1
Parareal, PITA 및 PFASST의 차이점은 무엇입니까?
Parareal, PITA 및 PFASST 알고리즘은 시간에 따른 시간 문제의 솔루션을 병렬화하기위한 도메인 간 기술입니다. 이러한 방법의 기본 원칙은 무엇입니까? 그들 사이의 주요 차이점은 무엇입니까? 하나는 다른 것을 기반으로 말할 수 있습니까? 어떻게? 그들의 응용 프로그램은 어떻습니까? "어떤 것이 더 낫습니까?"라는 질문에 대한 답이 없다는 것을 알고 있지만, 응용 분야와 검증 …
2
포트란의 옥트리 코드
나는 과학 컴퓨팅을 처음 사용합니다. Octree의 Fortran (바람직하게는 f90) 구현을 찾고 있습니다. 내 문제는 모든 상자에 N 개의 입자 (또는 라플라스 방정식과 같은 방정식을 해결하기 위해 적분 방정식 방법으로 연결할 수있는 밀도 값을 알고있는 소스)가 없을 때까지 도메인을 나누는 Octree이 필요합니다. 내가 찾을 수 있었던 것은 이 C ++ 구현 …
4
부분 Singular Value Decompositions (SVD)의 메모리 효율적인 구현
모형 축소를 위해 행렬 의 20 개 가장 큰 특이 값과 연관된 왼쪽 특이 벡터를 계산하려고합니다 . 여기서 N ≈ 10 6 및 k ≈ 10 3 입니다. 불행히도 내 행렬 A 는 아무런 구조없이 밀도가 높습니다.∈ R엔, kㅏ∈아르 자형엔,케이A \in \mathbb R^{N,k}엔≈ 106엔≈106N\approx 10^6k ≈ 10삼케이≈10삼k\approx 10^3ㅏㅏA 이 크기의 …
2
불연속 Galerkin / Poisson / Fenics
불연속 Galerkin 방법 (DG)과 다음과 같은 이산화를 사용하여 2D 푸 아송 방정식을 풀려고합니다 (PNG 파일이 있지만 업로드 할 수 없습니다, 죄송합니다). 방정식 : ∇ ⋅ ( κ ∇ T) + f= 0∇⋅(κ∇티)+에프=0\nabla \cdot( \kappa \nabla T) + f = 0 새로운 방정식 : 큐= κ ∇ T∇ ⋅ q= − …
1
암시 적 시간 이산화를 이용한 cuda 및 수치 적 방법
IMPLICIT 형식의 유한 체적 방법으로 (시간 이산화) 부분 편미분 방정식 (PDE)을 해결하는 일부 코드를 포팅하려고합니다. 결과적으로 ADI / TDMA 방식으로 처리되는 x, y, z 방향의 3 각 방정식 시스템이 있습니다. CUDA를 사용한 PDE의 암시 적 솔루션과 관련하여 아무것도 찾을 수 없습니다. CUDA에서 ADI / TDMA 체계를 구현할 수 있습니까? 어딘가에 …
3
비 압축 Navier-Stokes 용으로 제조 된 솔루션 — 발산이없는 속도 필드를 찾는 방법은 무엇입니까?
제조 솔루션 (MMS) 방법에서는 정확한 솔루션을 가정하고 방정식에서 솔루션을 대체하고 해당 소스 항을 계산합니다. 그런 다음 솔루션은 코드 확인에 사용됩니다. 비 압축 Navier-Stokes 방정식의 경우 MMS는 연속성 방정식에서 (0이 아닌) 소스 항을 쉽게 유도합니다. 그러나 모든 코드가 연속성 방정식에서 소스 용어를 허용하는 것은 아니므로 이러한 코드의 경우 분기가없는 속도 필드를 …
2
행렬 인수를 사용하여 선형 시스템 해결
우리는 표준 선형 시스템을 해결하기 위해 많은 계산 방법에 익숙합니다. 그러나 더 일반적인 (유한 치수) 선형 시스템을 풀기위한 "표준"계산 방법이 있는지 궁금합니다.Ax=b.Ax=b. Ax=b. 여기서 A 는 m 1 × n 1 행렬, B 는 m 2 × n 2 행렬, L 은 m 1 × n 1 행렬을 m 2 …
2
Crank-Nicolson 이산화에 의해 열 방정식의 최대 / 최소 원리가 유지됩니까?
Crank-Nicolson 유한 차분 법을 사용하여 1D 열 방정식을 풀고 있습니다. 열 방정식의 최대 / 최소 원칙 (즉, 최대 / 최소가 초기 조건에서 또는 경계에서 발생 함)이 이산화 된 솔루션을 유지하는지 궁금합니다. 이것은 아마도 크랭크-니콜슨이 안정적이고 수렴적인 계획이라는 사실에 의해 암시 될 것입니다. 그러나 Crank-Nicolson 스텐실에서 생성 된 행렬을 사용하여 선형 …
5
비 압축 흐름의 정의
모든 사람들이 압축 불가능한 흐름이 실제로 존재하지 않는다는 것을 알고 있기 때문에 지배 방정식을 단순화하기 위해 가정을 도입했습니다. 우리는이 가정을 간단하게 적용 할 수 없습니다. 일반적으로 마하 수 (비압축성 흐름의 경우 M <0.3), 밀도 변화 (제로 밀도 변화) 및 속도의 발산 (비압축성 흐름의 경우 0과 같음)은 흐름을 비압축성 흐름으로 정의하는 …
1
과학 컴퓨팅의 맥락에서 C ++ 11 이동 시맨틱의 영향은 무엇입니까?
C ++ 11에는 예를 들어 C ++ 03이 복사 구성 또는 복사 할당을 수행해야하는 상황에서 코드 성능을 향상시킬 수있는 이동 의미론이 도입되었습니다. 이 기사에서는 C + 11로 컴파일 할 때 다음 코드에서 5 배 속도가 향상된다고보고합니다. vector<vector<int> > V; for(int k = 0; k < 100000; ++k) { vector<int> x(1000); …
4
Matlab : 프로그래밍 방식으로 코드 실행을 안전하게 중지하는 방법이 있습니까 (예 : FORTRAN의 정지)? [닫은]
닫은. 이 질문은 주제에 맞지 않습니다 . 현재 답변을받지 않습니다. 이 질문을 개선하고 싶습니까? Computational Science Stack Exchange에 대한 주제가 되도록 질문을 업데이트하십시오 . 휴일 2 년 전 . 제목에서 알 수 있듯이 특정 위치에서 코드를 중지하고 안전하게 중지하고 싶습니다. 예를 들어 FORTRAN에서 중지 명령이있는 것처럼 명령을 찾을 수 없습니다.
10 matlab 
1
MATLAB에서 linsolve와 mldivide의 차이점은 무엇입니까?
matlab에서 linsolve 와 mldivide 는 모두 결정, 과소 결정 및 불확실한 경우 모두 선형 방정식 시스템을 푸는 데 사용됩니다. 그들의 문서를 읽었을 때, 그들 사이에 어떤 차이점이 있는지 궁금했습니다. 그들은 세 가지 경우에 거의 동일한 행렬 분해 및 삼각 알고리즘 알고리즘을 사용하고 있습니까? A에 opts 속성이있는 경우 linsolve는 mldivide보다 빠릅니다. …
1
유체 시뮬레이션의 경계 조건
비디오 게임 유체 시뮬레이션에 설명 된대로 소용돌이 입자 / "vorton"을 사용하여 2D 유체 시뮬레이션 작업을하고 있습니다. 제가 생각하기에 "이산 소용돌이 법"과 같은 것입니다. 기본적으로 정의 된 소용돌이가있는 입자 모음으로 유체를 나타내고 Biot-Savart의 법칙을 사용하여 모든 소용돌이의 유도 속도를 합산하여 한 지점에서 유체 속도를 계산합니다 (예 : , 여기서p2-p1은 샘플 지점과 …
당사 사이트를 사용함과 동시에 당사의 쿠키 정책과 개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.