전산 과학

I는 C ++로 MATLAB에서 기존 코드를 포팅하고 해결하는 선형 시스템이 오전 (오히려 더 일반적인 형태보다 X = B를 )xA=bxA=bxA=bAx=bAx=bAx=b 행렬 는 밀도가 높고 일반적인 형태이지만 1000x1000보다 크지 않습니다. MATLAB에서 솔루션은 함수 또는 슬래시 표기법으로 찾을 수 있습니다.AAAmrdivide(b,A)x = b/A; BLAS 및 LAPACK 루틴을 사용하여 C ++ 코드에서이 문제를 어떻게 해결해야합니까? …

12 linear-solver  c++  lapack 

MATLAB을 사용하여 모든 시간 단계에서 를 풀어야하는 문제를 해결하고 있습니다 . 여기서 b 는 시간에 따라 바뀝니다. 지금은 MATLAB을 사용 하여이 작업을 수행하고 있습니다 .Ax=bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}bb\mathbf{b}mldivide x = A\b 필요한만큼 사전 계산을 수행 할 수있는 유연성이 있으므로보다 빠르고 정확한 방법이 있는지 궁금합니다 mldivide. 여기서 일반적으로 수행되는 작업은 무엇입니까? 모두 감사합니다!

12 linear-algebra  matlab  linear-solver  matrix  linear-system 

오른쪽 (300 ~ 1000)이 많은 동일한 희소 선형 시스템 (300x300 ~ 1000x1000)을 해결해야합니다. 이 첫 번째 문제 외에도 다른 시스템을 해결하고 싶지만 0이 아닌 동일한 요소 (단지 다른 값)로 일정한 희소성 패턴을 가진 많은 희소 시스템입니다. 내 행렬은 무기한입니다. 분해 및 초기화 성능은 중요하지 않지만 해결 단계의 성능은 중요합니다. 현재 …

12 linear-algebra  petsc  linear-solver  sparse-matrix 

그래서 V 사이클이 어떻게 실행되는지 이해합니다 (또는 적어도 믿습니다). Matlab에서 1-D, 재귀 버전의 V 사이클을 작성했습니다. 그러나 FMG에 대한 코드를 실행했을 때 솔루션이 수렴되지 않았습니다. 내 문제는 실제 FMG 부분을 이해하는 데 있다고 생각합니다. 내가 현재 알고있는 것은 이것입니다 : FMG 보간 직전, 내 솔루션 완화 한uuu 보간 모두 오류 …

12 multigrid 

비동기 셀룰러 오토마타로 설명 할 수있는 계산 모델 모음이 있습니다. 이 모델은 Ising 모델과 비슷하지만 약간 더 복잡합니다. 마치 그러한 모델이 CPU가 아닌 GPU에서 실행되는 것이 도움이 될 것 같습니다. 불행히도 그러한 모델을 병렬화하는 것은 매우 간단하지 않으며 어떻게 진행 해야하는지 명확하지 않습니다. 나는 그 주제에 관한 문헌이 있다는 것을 …

12 parallel-computing  monte-carlo  gpu 

변수 및 x 2 (다른 모든 상수)에 대한 다음 방정식 시스템을 풀려고합니다 .피, x1P,x1P,x_1엑스2x2x_2 A(1−P)2−k1x1=0AP2−k2x2=0(1−P)(r1+x1)4L1−P(r1+x2)4L2=0A(1−P)2−k1x1=0AP2−k2x2=0(1−P)(r1+x1)4L1−P(r1+x2)4L2=0\frac{A(1-P)}{2}-k_1x_1=0 \\ \frac{AP}{2}-k_2x_2=0 \\ \frac{(1-P)(r_1+x_1)^4}{L_1}-\frac{P(r_1+x_2)^4}{L_2}=0 나는 x 1 과 x 2에 대한 방정식 1과 2를 각각 풀고 방정식 3으로 대체 함으로써이 방정식 시스템을 단일 변수 의 단일 방정식으로 바꿀 수 있음을 알 수 있습니다. matlab의 명령을 …

12 optimization  convergence  newton-method 

비정형 셀 기반 유한 체적 CFD에서 셀 탐색을위한 효율적인 데이터 구조에 대한 조언에 관심이 있습니다. ( dolfyn CFD 코드에서) 발생한 한 가지 예 는 다음과 같습니다 (관련 세그먼트 표시) 따라서 각 셀의면 수가 저장되는 배열 NFace가 있습니다. 그런 다음 CFace 배열은 로컬-셀 얼굴 번호를 전역 얼굴 번호로 매핑합니다.\begin{listing}do ip=1,Ncel ... …

12 fluid-dynamics  finite-volume  data-structures 

드리프트-확산 모델에서 시작하여 교육 학적 목적으로 기본 반도체 모델을 시뮬레이션하려고합니다. 상용 반도체 시뮬레이터를 사용하고 싶지는 않지만 다른 (일반적이거나 최근이거나 모호한) 모델을 배우고 있지만 상용 PDE 솔버를 사용하고 싶습니다. 그러나 간단한 1D 사례의 경우에도 드리프트-확산 모델은 여러 개의 결합 된 비선형 PDE로 구성됩니다. 전류 밀도 방정식 J p = q p …

12 pde 

나는 0에 가까운 특정 값에서 NaN을 나타내는 까다로운 적분을 처리하고 있으며 현재 ISNAN 문을 사용하여 NaN을 처리하는 까다로운 적분을 다루고 있습니다. FORTRAN의 NMS 라이브러리 (q1da 루틴-q1dax는 다르지 않음)와 C의 GSL 라이브러리 (QAGS 루틴 사용)를 사용하여이를 시도했습니다. 나는 integrand에서 NaN과 INF를 처리하도록 특별히 설계된 CQUAD (C에 대한 GSL 라이브러리의 일부)를 살펴 …

12 quadrature 

수치 안정성을 희생하지 않고 큰 희소 행렬에서 행렬 연산을 수행하는 라이브러리를 찾고 있습니다. 행렬은 1000+ x 1000+이고 행렬 값은 0과 1000 사이입니다. 인덱스 미적분 알고리즘을 수행 하므로 행렬의 행 벡터를 연속적으로 생성합니다. 각 행을 개발할 때 선형 독립성을 테스트해야합니다. 원하는 수의 선형 독립 벡터로 행렬을 채우면 행렬을 줄 바꿈 형태로 …

12 linear-algebra  algorithms  matrix 

로프 객체를 시뮬레이션하려고합니다. 내가 이해하는 제제는 스프링으로 연결된 입자 배열입니다. 이 스프링은 k 값이 매우 커서 선이 변형되지만 거의 늘어나지 않습니다. 로프는 진자 (닫힌 형태가 아님)의 일반화이므로 닫힌 함수에서는 이것을 시간의 함수로 해결하는 것이 불가능하다는 결론을 내 렸습니다. 그러면 대략적인 솔루션에 정착합니다. 잘 확장되는 알고리즘이 필요합니다. 예에서는 입자를 이동시키기 위해 …

12 simulation 

3 노드의 4 노드 요소에 다항식을 통합하고 싶습니다. FEA에 관한 몇 권의 책은 임의의 평평한 4-noned 요소에 대해 통합이 수행되는 경우를 다룹니다. 이 경우 일반적인 절차는 Jacobi 행렬을 찾아 적분 기준을 더 간단한 적분 한계 [-1; 1]가 있고 가우스 레전드 르 직교 기법이 쉽게 사용되는 정규화 된 것으로 변경하는 결정 …

12 finite-element  quadrature 

Nocedal & Wright의 Book Numerical Optimization (2006)에 따르면 부정확 한 라인 검색에 대한 Wolfe의 조건은 하강 방향 ,피pp 충분한 감소 : 곡률 조건 : ∇ f ( x + α p ) T p ≥ c 2 ∇ f ( x ) T p 대 0 < C 1 < …

12 optimization 

수치 적으로 평가하려는 일부 변수 x i 의 실제 값 함수 가 있다고 가정합니다 . 일반적으로 f에 대한 공식 은 제품, 합리성, 초월 함수 등을 포함 할 수 있으며 수치 안정성을 분석적으로 조사하는 데 오랜 시간이 걸릴 것입니다. 또는 실제로 실제로는 시간이 오래 걸릴 것입니다. 보장 된 안정성과 동등한 것이 …

12 stability 

나는 궁금해 : d u 를 해결하는 가장 좋은 알고리즘은 무엇입니까 dudt=Aududt=Au\begin{equation} \frac{du}{dt} = Au \end{equation} 여기서AAA는 실수n×nn×nn\times n행렬입니다. A는 명시 적으로 시간 의존적이지 않으며, 일반적으로 드물지만 반드시 밴딩되는 것은 아닙니다. 고유 값에는 양수가 아닌 실제 부분이 있습니다. A는 대각선 화가 가능하지만 전체 대각선 화가 계산적으로 효율적 이기에는 너무 클 …

12 linear-algebra  ode