전산 과학

Matlab pde 도구 상자를 사용하여 2D에서 특정 타원 방정식을 풀고 있습니다. 주어진 선을 따라, 즉 솔루션을 나타내는 3D 메쉬에서 평면 슬라이스를 자르기 위해 플롯해야하지만 솔루션은 좋습니다. 도구 상자 기능을 똑똑하게 포함하는 방법을 알 수 없습니다 (즉, 삼각 메쉬에 낮은 수준의 보간이 포함되지 않음). 도움을 주셔서 감사합니다.

9 pde  matlab 

인구 통계 모델을 데이터에 더 잘 맞출 수 있도록 일부 인구 통계 모델링 소프트웨어의 최적화 프로세스를 개선하기 위해 노력하고 있습니다. 최적화 시간을 줄이고 싶습니다. 목적 함수를 평가하는 데 걸리는 시간은 입력 값에 따라 크게 다릅니다. 목적 함수를 평가하는 시간과 입력 간의 관계가 알려져 있습니다. 평가할 포인트를 선택할 때 목적 함수의 …

9 optimization 

조건 번호의 정의에서 그것을 계산하기 위해 행렬 반전이 필요한 것 같습니다. 일반 정사각형 행렬 (또는 대칭 양수 한정의 경우)이 행렬 분해를 이용하여 조건 번호를 계산할 수 있는지 궁금합니다. 더 빠른 방법.

9 linear-algebra  matrix  condition-number 

연구 작업을 위해 C ++ 계산 코드와 Octave / Matlab 간의 통신을 용이하게하는 작은 C ++ 라이브러리를 개발했습니다 (후자는 후 처리 목적으로 사용되는 경우). Sourceforge 또는 Github와 같은 무료 호스팅 사이트에서 GPL로 그러한 라이브러리를 릴리스하고 싶습니다. 그러나 솔직히 말해서 (소규모) 도서관에 더 적합한 사이트를 선택할 수있는 경험이 부족합니다. 한편으로, github는 …

9 matlab  libraries  c++ 

균일 한 그리드에서 체비 쇼프 보간 다항식의 빠른 (대략적인) 평가를 구현하는 바람직한 방법이 있습니까 (체비 쇼프 노드에서 함수 값이 제공됨)? 내 문제는 보간 다항식의 정도가 증가하면 보간이 느려진다는 것입니다. 다음과 같은 아이디어가 떠 올랐습니다. NFFT (Non-Uniform FFT) 기술을 적용하십시오. FFT를 사용하면 Chebyshev 노드에서 파생물을 계산할 수 있습니다. 처음에는 더 미세한 …

9 interpolation  fourier-analysis  fftw 

나는 해결하고 싶다 케이u = bKu=bK u = b 어디 케이KK강성 매트릭스입니다. 그러나 일부 제약 조건이 누락되어 시스템에 일부 강체 운동이 여전히 존재할 수 있습니다 (고유 값 0으로 인해). 선형 시스템을 해결하기 위해 CG를 사용하고 있기 때문에 때로는 CG가 반 양성 문제에 수렴하지 않지만 때로는 수렴 할 수 있으므로 허용되지 …

9 finite-element  iterative-method  conjugate-gradient 

FFT를 사용하여 푸리에 유사 스펙트럼 방법 (Orzag & Patterson, PRL, 1972)으로 등방성 난류의 직접 수치 시뮬레이션을 수행합니다. 난기류 커뮤니티에서 널리 사용되는 방법의 배경은 다음 코스를 참조하십시오. http://www.math.ualberta.ca/~bowman/m655/lab3d.pdf 로 구성되는 소위 거래 처리에 규칙 사용 여기서 는 시간, 는 파수, 는 최대 파수 및 는 속도의 스펙트럼 진폭입니다.2 / 32/32/3유^( k …

9 fluid-dynamics  simulation 

함수의 최소값을 찾기 위해 Broyden-Fletcher-Goldfarb-Shanno 메소드를 구현하려고합니다. 두 가지 초기 추측이 필요합니다x−1x−1x_{-1}& 및 초기 Hessian Matrix 근사값 . 에 대해 찾은 유일한 요구 사항 은 Hessian이 대칭 양수이면 이어야한다는 것 입니다. Wikipedia를 살펴보면 일반적인 초기 근사값은 (정체 행렬)입니다. 이것은 항상 좋은 초기 입니까? 내가 아닌 다른 선택 할 수 있습니다 …

9 optimization  algorithms 

작은 (일반적으로 60x60보다 작은) 비대칭 행렬의 모든 양의 (고유 값이 양수) 양의 고유 쌍을 찾아야하는 문제가 있습니다. 고유 값이 특정 임계 값보다 작은 경우 계산을 중지 할 수 있습니다. 고유 값이 실제라는 것을 알고 있습니다. 최고의 성능을 짜기 위해 사용할 수있는 알고리즘에 대한 제안이 있으십니까? 이러한 분해를 수천 번해야하므로 속도가 …

9 performance  eigensystem  gpu 

내 프로그래밍의 대부분은 내가 사용하는 C의 일회용 연구 코드입니다. 가까운 공동 작업자 이외의 다른 코드를 배포 한 적이 없습니다. 과학 저널에 게시하는 알고리즘을 개발했습니다. 기사의 온라인 보충 자료에 소스 코드와 실행 코드를 제공하고 싶습니다. 동료가 알고리즘을 일반화하여 C ++ (ack!)로 작성해야하고 작은 조밀 한 선형 시스템을 풀어야한다고 요구했습니다. 알고리즘에 대한 …

9 linear-solver  libraries  c++  c 

몇 가지 교과서를 통해, 나는 줄을 검색하는 동안 처음에 최소값을 브라케팅하는 문제가 (적어도 학부 텍스트에서) 나중에 생각되는 경향이 있음을 알았습니다. 이러한 유형의 문제에 대해 잘 정립 된 기술 또는 모범 사례가 있습니까? 아니면 솔루션이 일반적으로 응용 프로그램에 종속됩니까? 누구든지 주제에 대한 참조를 추천 할 수 있습니까?

9 optimization  reference-request 

나는 매우 다른 자원들과 연구들과의 대화에서 수치 부분 미분 방정식에서 고정밀 계산에 대한 수요가 증가하고 있다는 인상을 받았습니다. 여기서 높은 정밀도는 표준 64 비트 배정도보다 더 높은 정밀도를 의미합니다. 이 주제의 최신 기술이 궁금합니다. 비교를 위해, 예를 들어, 멀티 코어 방법, 대규모 병렬화 또는 GPU 컴퓨팅을 목표로하는 수치 적 PDE에는 …

9 reference-request  precision 

Answer 는 혼돈 동기화에서 결합 발진기에 대한 조건부 Lyapunov 지수 (CLE)를 계산하는 소프트웨어를 제공 합니다. 그러나 따르기가 어렵고 플롯의 그래픽 출력이 없습니다 (그리고 C에는 더 복잡합니다). 누구나 결합되지 않은 시스템에 매우 적합한 LET Toolbox 를 수정하는 방법을 알고 있지만 CLE을 수용하기 위해 동기화 된 시스템으로 작업하는 방법을 이해하지 못합니다. CLE의 …

9 matlab  ode  simulation  quantum-mechanics 

주어진 o 포인트가 주어진 알고리즘이 있는지 알고 싶습니다. 각도가 이고 주어진 경우 "perimeter"에 더 가까운 엔벨로프를 계산하면 각도가 볼록-헐을 계산합니다. ".α=0α=0\alpha = 0α>0α>0\alpha > 0 그리고 점 세트의 교차하지 않는 둘레의 정의가있는 경우이 경우 가 클 때 결과 다각형이됩니다 .αα\alpha 문제의 또 다른 관점은 대한 최소 주변 솔루션 (볼록-헐)과 (정규화 …

9 algorithms  computational-geometry 

복잡한 목적 함수를 최소화하고 싶지만 볼록한지 확실하지 않습니다. 볼록하지 않다는 것을 증명하는 멋진 알고리즘이 있습니까? 물론 알고리즘은이를 증명하지 못할 수 있는데,이 경우 볼록인지 아닌지 알 수 없으며 이것이 정상입니다. 목적 함수가 볼록한지 여부를 분석적으로 결정하기 위해 많은 시간을 소비하기 전에 볼록성을 배제하려고합니다. 예를 들어 문제를 볼록한 것으로 알려진 표준 형식으로 …

9 convex-optimization