다양한 차원의 PDE


14

필자는 PDE에 대한 근사 솔루션을 찾는 대부분의 방법이 차원 수에 따라 확장 성이 낮으며 Monte Carlo가 ~ 100 차원을 요구하는 상황에 사용된다는 것을 알고 있습니다.

~ 4-10 차원에서 PDE를 효율적으로 수치 적으로 해결하는 좋은 방법은 무엇입니까? 10 ~ 100?

몬테카를로 외에 차원의 수에 맞게 확장 할 수있는 방법이 있습니까?


1
해결하려는 문제의 종류에 대해 좀 더 자세한 정보를 제공하는 데 도움이 될 수 있습니다. 전산 과학에서 다루는 대부분의 PDE는 최대 4 차원 (시간 + 3 차원) 인 경향이 있습니다. 변수가 공간 또는 시간 변수입니까, 아니면 다른 종속 항목이 있습니까?
aeismail

1
공간 변수. 양자 역학에서 당신이 밀도 함수 이론 또는 하 트리 - 폭 상태에서 사용하는 근사값을하지 않으려면, 파동 함수는 차원 여기서 n은 전자의 수입니다. 따라서 작은 원자와 분자조차도 올바르게 처리하려면 많은 치수가 필요합니다. 3nn
Dan

1
솔루션에 대해 알고 싶은 정보에 따라 다릅니다. 전자파 함수에 대한 모든 세부 사항을 알고 싶어하는 사람은 거의 없습니다 . 따라서 계산 기술을 실제로 원하는 정보에 테일러해야합니다. n
Arnold Neumaier

1
전자식 슈뢰딩거 방정식의 Monte Carlo 솔루션에 대한 참조를 100 차원으로 인용하십시오.
Arnold Neumaier

참조가 없습니다. QCD에 사용되는 많은 차원의 시뮬레이션에 대해서만 들었습니다. 나는 4-5 차원의 Schroedinger 시뮬레이션 만하고 싶지만 몬테 카를로 이외의 것이 차원의 수와 잘 비례하는지, 그리고 100이 점근 적 스케일링을 얻기 위해 멋지고 큰 둥근 숫자처럼 보 였는지 궁금했습니다.
Dan

답변:


13

다차원에서 기본 또는 직교 (많은 경우에 MC를 대체 할 수 있음)를 제공하는보다 체계적인 방법은 희소 격자의 것입니다 . 이것은 희소 격자 의 것입니다. 차원 이 해상도 N 의 지수가되기보다는 차원 2d .Nd

이것은 일련의 1 차원 규칙 Q 1 l 을 결합한 Smolyak 구적법을 통해 수행됩니다.Ql1 됩니다.

Qnd=ln(Qi1Qi11)Qmi+1d1

이것은 공간에서 높은 혼합 차수가 제거 텐서 곱 구적 공간과 동일 합니다. 이것이 충분히 심각한 방식으로 수행되면 복잡성이 크게 향상 될 수 있습니다. 그러나,이를 수행하고 근사를 양호하게 유지하기 위해서는 용액의 규칙 성이 혼합 유도체를 충분히 소멸시켜야한다.

스파 스 그리드가 같은 것들에 대한 GRIEBEL 그룹에 의해 죽음에 구타 한 구성 공간에서 슈뢰딩거 방정식 다른 고차원적인 것들에 . 응용 프로그램에서 사용되는 기본 기능은 중첩 할 수있는 한 일반적 일 수 있습니다. 예를 들어, 평면파 또는 계층 기반이 일반적입니다.

자신을 코딩하는 것도 매우 간단합니다. 그러나 내 경험상 실제로 이러한 문제를 해결하는 것은 매우 어렵습니다. 좋은 튜토리얼 이 있습니다.

빠르게 죽는 파생 상품을 특징으로하는 특화된 Sobolev 공간에 솔루션이 존재하는 문제의 경우 희소 그리드 접근 방식으로 인해 더 큰 결과를 얻을 수 있습니다.

Acta Numerica 검토 논문, 고차원 파라 메트릭 및 확률 론적 PDE의 스파 스 텐서 이산화를 참조하십시오 .


희소 격자가 적용되지 않는 잘 알려진 예가 있습니까?
MRocklin

1
당신은 진정으로 규칙이 필요합니다. 또한 QM과 같이 불쾌한 고차원 교두가있는 경우 조심해야합니다. 나는 그것이 아니라는 것을 (심지어 증거로) 인정하기 시작하는 스파 스 그리드 파벌에 대한 몇 가지 이야기를 듣고 훨씬 더 몬테 - 카를로보다,하지만 좋은 참조를 찾을 수 없습니다.
Peter Brune

당신이 언급 한 슈 레딩 거를위한 희소 격자에있는 종이는 2 개의 전자만을 취급합니다. 이 방법으로 몇 개의 전자가 실제로 다루기 쉬운가?
Arnold Neumaier

6

일반적으로 일반 그리드가 3 차원 또는 4 차원 문제를 넘어 설 수없는 이유를 쉽게 이해할 수 있습니다. d 차원에서 좌표 방향 당 최소 N 포인트를 원하면 N ^ d를 얻습니다. 전반적으로 포인트. 1d에서 비교적 좋은 기능을 수행하더라도이를 해결하려면 적어도 N = 10 그리드 포인트가 필요하므로 전체 포인트 수는 10 ^ d가됩니다. 예를 들어 d를 넘어 서기 어려운 최대 컴퓨터에서도 = 9, 아마도 그 이상을 넘어 가지 않을 것 입니다 . 솔루션 함수에 특정 속성이있는 경우 일부 상황에서 스파 스 그리드가 도움이 될 수 있지만 일반적으로 차원 저주의 결과에 따라 생활하고 MCMC 방법을 사용해야합니다.


MCMC는 무엇을 의미합니까?
Dan


2

d=4,...,100d=100,101,... . 따라서 치수가 4에서 100으로 증가하더라도 몬테 카를로는 어쨌든 더 빠를 것입니다.


2
O(N)107

Ck,α
당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.