«c++» 태그된 질문

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 

SIAM Journal of Scientific Computing의 "표현 템플릿 표현 : 현재 방법론의 성능 분석"은 "블레이즈"선형 대수 라이브러리를 참조합니다. 이전에 들어 본 적이 없으며 온라인 참조를 찾을 수 없습니다. 명백한 구글 검색은 위의 논문을 돌려주고있다. 이 라이브러리는 무엇이며 어디에서 더 배울 수 있습니까?

12 linear-algebra  reference-request  c++ 

이 QuasiMonteCarlo방법을 사용하여 3 또는 4 차원으로 일부 적분을 수행하는 Mathematica 프로그램이 있습니다. 문제는 HPC 클러스터에서 사용할 수있는 최대 작업 시간으로 이러한 계산 중 일부를 완료 할 수없는 시점까지 실행하는 데 시간이 오래 걸린다는 것입니다. 그래서 C ++로 프로그램을 다시 작성하는 것을 고려하고 있습니다. 나는 GSL 문서를 보았고 유사 시퀀스 …

12 libraries  monte-carlo  c++  mathematica 

암시 적 숫자 체계를 사용하여 물리적 문제를 해결하고 있습니다. 이것은 3 각형 행렬로 선형 방정식을 푸는 데 도움이됩니다. 이 알고리즘 을 Wikipedia에서 코딩 했습니다 . 이 유형의 방정식을 최적화 된 방식으로 해결할 수있는 효율적인 라이브러리가 있는지 궁금합니다. 중요한 점은 시스템 매개 변수가 변경 될 때만 행렬 자체가 변경되므로 좋은 성능 …

12 linear-algebra  libraries  c++ 

오일러 방법을 Runge-Kutta 4 차 순서로 대체하여 일정한 중력 크기가 아닌 자유 낙하 운동을 결정하려면 어떻게해야합니까 (예 : 지상 10,000km에서 자유 낙하)? 지금까지 나는 오일러 방법으로 간단한 통합을 작성했습니다. while() { v += getMagnitude(x) * dt; x += v * dt; time += dt; } x 변수는 현재 위치를, v는 …

11 numerical-analysis  c++  runge-kutta  integral-equations 

12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ 나는 현재 armadillo 라이브러리로 이것을하고 있지만 너무 느립니다. 문제는 1 조의 행렬에 대해이 작업을 수행해야하며 두 결정 요인을 계산하는 것이 내 프로그램의 병목 현상이라는 것이 밝혀졌습니다. 따라서 두 가지 질문이 있습니다 크기를 알면 결정자를 더 빨리 계산하는 데 사용할 수있는 트릭이 있습니까? …

11 linear-algebra  matrix  c++  graph-theory  c 

나는 1967 년 Bulirsch & Stoer (Numerische Mathematik, 9, 271-278)가 발표 한 ALGOL 프로그램의 C ++ 적응 인 수치 적분 (quadrature)을위한 작은 서브 루틴을 가지고 있습니다. 더 현대적인 (적응) 알고리즘으로 업그레이드하고 그러한 (무료) C ++ 라이브러리가 있는지 궁금합니다. 나는 GSL (C)로 보였지만 끔찍한 API가 제공됩니다 (숫자가 좋을 수도 있지만). 다른 …

10 c++  quadrature 

Windows 기계를 사용하여 C ++로 작성된 천체 역학 / 우주 역학 응용 프로그램에서 Runge-Kutta 8 차 방법 (89)을 사용하고 싶습니다. 따라서 누구나 문서화되고 무료로 사용할 수있는 좋은 라이브러리 / 구현을 알고 있는지 궁금합니다. 컴파일 문제가 예상되지 않는 한 C로 작성해도 괜찮습니다. 지금 까지이 라이브러리 (mymathlib)를 찾았습니다 . 코드는 괜찮아 보이지만 …

10 c++  reference-request  libraries  c  runge-kutta 

C ++ 11에는 예를 들어 C ++ 03이 복사 구성 또는 복사 할당을 수행해야하는 상황에서 코드 성능을 향상시킬 수있는 이동 의미론이 도입되었습니다. 이 기사에서는 C + 11로 컴파일 할 때 다음 코드에서 5 배 속도가 향상된다고보고합니다. vector<vector<int> > V; for(int k = 0; k < 100000; ++k) { vector<int> x(1000); …

10 finite-element  performance  c++ 

현재 나는 기초 스플라인 배열이라는 미분 방정식 해결 방법을 연구하고 있습니다. 내가 겪고있는 것은 B k + 1 i ( x ) = x − x i 의 관계로 임의의 순서 스플라인을 만드는 방법을 만드는 것입니다 초기 조건 B 1 i (x)={ 1 인 x k + i + 1 …

10 c++  b-spline 

기본적으로 FEM은 거의 "해결 된"문제인 것 같습니다. Trilinos, PETSc, FEniCS, Libmesh 또는 MOOSE와 같은 수많은 강력한 프레임 워크가 존재합니다. 그들이 공통적으로 가지고있는 것 : 그들은 매우 "무거운"것입니다. 첫째, 설치는 일반적으로 매우 고통 스럽습니다. 둘째, 인터페이스 / API가 두껍고 무겁습니다. 전체 아이디어를 해당 라이브러리의 사고로 변환해야합니다. 또한 특별한 요구 사항이나 기존 …

9 finite-element  python  c++  petsc  fenics 

큰 입방 고유 값 문제가 있습니다. (ㅏ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 

그래서 저는 파이썬 학습을 귀찮게해야하는지에 대해 토론했습니다. 교수님과의 대화에서 Matlab은 학계 에 관한 한 응용 수학 / 계산 과학에 사용되는 공통 언어 인 것 같습니다 . 반면에 산업 내 교수 (특히. 업계에서 일한 사람들) 말했다 학습 C ++는 가장 안전한 길이다. 나는 파이썬을 귀찮게 해야하는지 또는 당분간 내가 아는 것 …

9 matlab  python  c++ 

현재 matlab "fmincon"함수로 구현 된 비선형 제약 최소화 문제를 해결하려고합니다. 내 기대는 최소화 (fun1, x0, uB, lB, fun2)입니다. 여기서 x0은 초기 상태이고 fun1은 최소화 해야하는 함수, uB는 상한, lB는 하한 및 fun2는 비선형 평등의 벡터를 제공하는 함수입니다 http://www.mathworks.com/help/optim/ug/fmincon.html에 설명 된 / inequalitiesnonlcon 기능으로. 이러한 벡터는 반복을 통해 변경됩니다 (x_n, n …

9 c++  nonlinear-programming