Intel MKL 라이브러리를 사용하여 복수의 오른쪽 (rhs) 벡터가있는 선형 방정식 시스템 (A*x = b)을 풀고 있습니다. rhs 벡터는 비동기 적으로 생성되며 별도의 루틴을 통해 생성되므로 한꺼번에 해결할 수 없습니다. 프로그램을 신속하게 처리하기 위해 각 스레드가 단일 rhs 벡터를 해결할 책임이있는 다중 스레드 프로그램이 사용됩니다. 행렬 A은
중/대형 치수 (1000-5000 방정식)의 선형 시스템을 해결하기 위해 Armadillo C++ 라이브러리를 사용하고 있습니다. I 다른 선형 시스템 AX를 해결하고 있으므로 는 = B, A는 항상 동일한 및 B의 변경 인가 , I는 LU 싶습니다 번만를 인수 분해하고 다른과 LU 인수 분해를 재사용 비. 불행히도 나는 아르마딜로에서 이런 종류의 작업을
에 선형 방정식의 시스템 설정 I I 해결하고자 다음 식 가지고 H*b0 = M(Q+1)b(Q+1)+l+M'B 는 unknowns가 b0, b(q+1) 및 B됩니다. 알려진 행렬의 크기는 다음과 같습니다. H=(42 x 42)
M(Q+1) = (42 x 21-P)
l = (42 x 1)
M' = (42 x 4)
그래서 벡터를 찾는 방법을 알고
mxn 행렬을 얻고 QR이 그것을 분해하는 프로그램을 작성하려고합니다. 아직 완료되지 않았지만 문제가 발생했습니다. 나는 그들이 Q1라는 무엇 위키 피 디아 http://en.wikipedia.org/wiki/QR_decomposition A=[12,-51,4;6,167,-68;-4,24,-41]
, Q2에서의 예 내 프로그램을 실행 해봤 ... 내가
나는 특별한 데이터를 사용 할 수있는 구현을 찾고 스파 스 매트릭스 A와 B 그러므로 내가 의 일반화 된 특이 값 분해 (GSVD)을 계산 할 스파 스 행렬에 대한 구조. 내가 구현 한 유일한 구현 (here)은 Fortran 77로 작성된 LAPACK 패키지의 일부입니다. 꽤 잘 작동하지만 불행히도 스파 스 매트릭스를 처리 할 수 없습니다.
나는 사용자가 그린 몇 가지 포인트가 있습니다. 나는이 점들 사이에서 중요한 점을 찾고 싶다. 필자가 정의한 중요한 점은 포인트의 방향이 갑자기 바뀌는 지점입니다. 예를 들어 손으로 그려지는 'Z'는 2 개의 중요한 점을 가져야합니다. 나는 인접한 포인트 사이에서 계산 각을 시도했지만 이것이 원하는 결과를 얻지 못했습니다. 그리고 기울기의 변화를 계산하는
데이터 프레임의 열에서 그룹 단위 평균을 계산하고 뺄 수있는 병렬화 가능 코드 (plyr 및 doMC)를 작성하려고합니다. plyr 구문을 올바르게 수정하는 데 어려움을 겪고 있습니다. 여기 을위한 루프 작업과 스크립트입니다 data = data.frame(x = rnorm(100),y = rnorm(100),ID = round(runif(100)*10))
대용량 대칭 행렬 (~ 800MB)을 사용하는 선형 대수학을하고 있는데, 여러 가지 분해를 시도하고 있습니다. 현재 LDLT를 구현 중이며 L^T로 행렬의 위쪽 삼각형을 덮어 써서 메모리 요구량을 반으로 줄이려합니다. Eigen의 문서는 이것을 수행하기 위해 어쨌든 나열되어 있지 않지만 내부 네임 스페이스에 숨겨진 방법이 있습니다. 방법는이 구조체에서 찾을
로프 물리 시뮬레이션을 C++ (스프링 물리학)으로 만드는 중입니다. 나는 몇 개의 "뼈"(위치, 무게 등이있는 일부 입자 [질량])로 구성된 기본 로프 이동을 성공적으로 구현했습니다. 각 파티클의 위치가 계산 될 때, 나는 스킨 된 메쉬 뼈대의 해당 뼈에 각각을 할당합니다. 내가 원하는 것은 그것의 자 입자로 회전되도록 입자의 회전을 계산하는 것이다 (원
스칼라에 대한 일반적인 수학/통계 프레임 워크에 대해 묻는 질문은 questions입니다. 저는 하나의 특정 문제, 즉 거대한 선형 시스템을 해결하는 것에 관심이 있습니다. 본질적으로 나는 scipy.sparse.linalg.spsolve과 동등한 것을 찾고 있습니다. 현재 breeze-math은 ScalaNLP Breeze입니다.이 라이브러리 컬렉션의 초