전산 과학

내가 사용해야하는 프로젝트가 차 필드 형태의 구체적 숫자를a + b− 3−−−√a+b−3a + b \sqrt{-3} 와 a , b ∈ Qa,b∈Qa,b \in \mathbb{Q}. 예를 들어 다음은 아이젠 슈타인 정수 의 소수입니다 . 세이지를 사용하고 싶지 않습니다. 통합 할 데이터 형식을 직접 작성하고 싶습니다 numpy. PARI는 유용하지만 Python과 호환되지 않습니다. 이 객체에 …

9 linear-algebra  precision 

많은 적응 형 FEM 라이브러리는보다 고급 메쉬 데이터 구조를 사용하여 노드, 모서리, 삼각형, 사면체 등 추가 / 제거를 처리합니다. 예를 들어, p4est 라이브러리는 적응 형 메시 개선을 위해 옥트리 데이터 구조를 사용합니다. 스태틱 메시에서 계산에 사용되는 옥트리는 종종 없습니다. 적응 형 FEM에 대한 선형 대수 측의 변화는 무엇입니까? 내가 생각할 …

9 linear-algebra  finite-element  sparse-matrix 

나는 쉬운 병렬화가 DG 방법의 장점 중 하나라는 것을 항상 들어 왔지만, 왜 그런 이유가 연속 Galerkin에도 적용되지 않는지는 알 수 없습니다.

9 finite-element  parallel-computing  discontinuous-galerkin 

FEM과 XFEM의 주요 차이점은 무엇입니까? 우리는 언제 FEM의 XFEM을 사용해서는 안됩니까? 다시 말해, 새로운 문제를 만났을 때 어떤 문제를 사용하는지 어떻게 알 수 있습니까?

9 finite-element  numerical-analysis  adaptive-mesh-refinement  numerics 

큰 입방 고유 값 문제가 있습니다. (ㅏ0+ λㅏ1+λ2ㅏ2+λ삼ㅏ삼) x =0(ㅏ0+λㅏ1+λ2ㅏ2+λ삼ㅏ삼)엑스=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. 선형 고유 값 문제로 변환 하여이 문제를 해결할 수는 있지만 시스템은 .삼2삼23^2 ⎡⎣⎢−ㅏ0000나는000나는⎤⎦⎥⎡⎣⎢엑스와이지⎤⎦⎥= λ⎡⎣⎢ㅏ1나는0ㅏ20나는ㅏ삼00⎤⎦⎥⎡⎣⎢엑스와이지⎤⎦⎥,[−ㅏ0000나는000나는][엑스와이지]=λ[ㅏ1ㅏ2ㅏ삼나는000나는0][엑스와이지],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ 0 & \mathbf{I} & 0 \\ 0 & 0 & \mathbf{I} \end{bmatrix} …

9 c++  eigensystem  eigenvalues 

double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } 위의 내용은 trapezia를 사용하여 func()한도 사이 …

9 c++  performance 

Gauss-Seidel 대신 SOR을 수행 할 가치 가 있는지에 대한 간단한 경험 법칙이 있습니까? (및 Realxation 매개 변수를 추정하는 가능한 방법ωω\omega) 나는 매트릭스 를 보거나 매트릭스가 나타내는 특정 문제에 대한 지식 만으로 의미 합니까? 나는이 질문들에 대한 답을 읽고 있었다 : 연속적인 이완법 (SOR) 방법을 최적화하기위한 휴리스틱이 있습니까? 하지만 너무 정교합니다. …

9 linear-solver  iterative-method  heuristics 

이 질문은 comp-sci에서 약간 주제가 아닐 수 있습니다. 필요한 경우 어디에 맞는지 제안하십시오. 문제는 모든 시뮬레이션 실행을 효율적으로 관리하는 방법에 관한 것입니다. 예를 들어, 시뮬레이션에는 특정 제안 된 값 범위에서 정의되어야하는 2 개의 매개 변수를 수정해야합니다. 예를 들어 시뮬레이션 결과와 실험 데이터를 비교하여 두 매개 변수 쌍으로 생성 된 더 …

9 monte-carlo  data-management 

엄격하게 볼록한 제한되지 않은 최적화 문제하자 고유의 최소값을 표시하고 x_0 에 지정된 초기 근사치 X_ \ 텍스트 {옵트}. 우리는 벡터가 호출 X \ 엡실론 의 솔루션 닫기 \ mathcal {O}를 하는 경우 \ 시작 {식} \ FRAC {|| X - X _ {\ 텍스트 {옵트}} || _2} {|| x_0 - …

9 convex-optimization  complexity 

특정 클래스 수치 문제를 해결할 작은 MATLAB 패키지를 작성 중입니다. 알고리즘에는 3 단계가 있으며 사용자는 각 단계마다 5 가지를 선택할 수 있습니다. 나는 전체 문제를 사용하여 구현했다.≈ 20≈20\approx 20기능 및 3 개의 스위치 케이스 (각 알고리즘 단계마다 하나씩). 그것은 잘 작동하지만 더 많은 일 (5 개 이상의 선택과 한 단계 …

9 matlab 

Computational Science에 대한 메일 링리스트 나 Google 그룹 이 있는지 궁금합니다 . 질문과 답변 만하는 대신 질문을 토론 할 수 있습니다. 사실, 저는 PDE의 병렬 컴퓨팅과 수치 솔루션에 더 관심이 있습니다. 그러나 나는이 지역 사람들이 무엇을 어떻게하는지 모르겠습니다. 이 지역의 로드맵을 파악하기 위해 그들의 논문 만 읽을 수 있습니다. 안내 …

9 parallel-computing  reference-request 

지난 며칠 동안이 오류를 디버깅하려고했지만 아무도 진행하는 방법에 대한 조언이 있는지 궁금합니다. 셀 센터에 알 수없는 것이 정의되고 셀면의 플럭스가있는 불균일 유한 체적 메쉬에서 단계 전하 분포 (정전기 / 반도체 물리학의 일반적인 문제)에 대한 푸 아송 방정식을 풀고 있습니다. 0 = (ϕ엑스)엑스+ ρ ( x )0=(ϕx)x+ρ(x) 0 = (\phi_x)_x + …

9 discretization  finite-volume  poisson 

Fortran에서 프로그래밍하는 방법을 배우는 가장 좋은 방법이 무엇인지 궁금합니다. 나는 책을 찾고 있는데, 나에게 괜찮은 예는 "Fortran 95/2003 Explained" 와 "Fortran For Scientists and Engineers" 였습니다. 너희들이 어떻게 생각하는지 모르겠어? 내가 궁금해하는 또 다른 것은 사용해야 할 컴파일러입니다. Windows 7 컴퓨터 (64 비트)에서 작업합니다. 너희들이 나를 도울 수 있기를 바랍니다! …

9 fortran 

3 개의 연립 방정식에 유한 차분 법을 적용하고 있습니다. 두 방정식은 결합되지 않지만 세 번째 방정식은 다른 두 방정식에 모두 연결됩니다. 방정식의 순서를 에서 로 바꾸면 계수 행렬이 대칭이된다는 것을 알았습니다 .( x , y, z)(x,y,z)(x, y, z)( x , z, y)(x,z,y)(x, z, y) 이렇게하면 어떤 이점이 있습니까? 예를 들어, …

9 finite-difference  symmetry 

선형 결함 수정 (LDC)과 완전 근사법 (FAS)을 모두 사용하여 V-Cycle 멀티 그리드 솔버를 구현했습니다. 내 문제는 다음과 같습니다. LDC를 사용하면 잔차가 사이클 당 ~ 0.03 배 줄어 듭니다. FAS 구현은 선형 계수와도 수렴하지만 그 계수는 ~ 0.58에 불과합니다. 따라서 FAS는 사이클 수의 약 20 배가 필요합니다. 대부분의 코드는 공유되지만 유일한 …

9 multigrid