'Index type not supported yet' error when doing QR factorization using Eigen and SuiteSparseQR

-

Environment

Windows 10 21H1 / MSVC C++ 14 / eigen 3.3.7 / suitesparse  5.3.0

Error Specification

Assertion failed: false && "Index type not supported yet", file C:\vcpkg\installed\x64-windows\include\eigen3\Eigen\src/CholmodSupport/CholmodSupport.h, line 93

Solution

Just comment eigen_assert on line 93 of CholmodSupport.h.

Test Code

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
#include <Eigen/Sparse>
#include <Eigen/SPQRSupport>
#include <Eigen/SparseCholesky>
#include <Eigen/OrderingMethods>
#include <vector>
#include <iostream>

typedef Eigen::SparseMatrix<double> SparseColMat;

int main(int argc, char** argv)
{

    SparseColMat m(4, 3);
    m.coeffRef(0, 0) = -1; m.coeffRef(0, 1) = -1; m.coeffRef(0, 2) = 1;
    m.coeffRef(1, 0) = 1; m.coeffRef(1, 1) = 3; m.coeffRef(1, 2) = 3;
    m.coeffRef(2, 0) = -1; m.coeffRef(2, 1) = -1; m.coeffRef(2, 2) = 5;
    m.coeffRef(3, 0) = 1; m.coeffRef(3, 1) = 3; m.coeffRef(3, 2) = 7;

    std::cout << "Here is the matrix m:\n" << m << std::endl;

    Eigen::SPQR<SparseColMat> solver1(m);
    std::cout << "QR solver1 :\n"  << solver1.matrixR() << std::endl;

    m.makeCompressed();
    Eigen::SparseQR<SparseColMat, Eigen::COLAMDOrdering<int>> solver2(m);
    std::cout << "QR solver2 :\n" << solver2.matrixR() << std::endl;

    return 0;
}

Result



댓글

댓글 쓰기

이 블로그의 인기 게시물

Linux에서 특정한 디렉토리가 차지하는 용량을 효율적이고, 빠르게 계산하는 법(Fast, efficient way to calculate directory size recursively on linux)

-
이미지
[주의] 아직 업데이트 중인 글입니다. Standard unix command 중에서 '특정 directory와 하위 파일들이 파일시스템에서 차지하는 용량의 총합'을 알아낼 수 있는 가장 대표적인 방법으로 'du -s' ~ 'du --max-depth=0'가 있다. du command는 glibc의 fts를 이용하여 file system hierarchy를 (default)preorder로 recursive하게 traversing하여 모든 파일/디렉토리에 대하여 stat system call을 호출한다. 더 자세한 구현 내용은 Unix 'coreutils' repository의 'du.c'를 읽어보면 되겠다. (사실상 위의 두줄이 전부다.) glibc의 fts module implementation을 읽어보지 않아서, 실제로 directory tree traversing이 recursive하게 이루어지는지를 알 수는 없다. (설마?) 만약, 그렇다면 pre-order traversing 을iterative하게 바꾸면 속도가 좀 빨라지긴 할 것이다 .일단 이건 논외로 하고(...) 다른 대책으로는 다음과 같은 해결방법들이 있을 수 있다. 해결방법(?) (1) 모든 subdirectory들을 recursive하게 방문하면서 계산한다. (2) file과 관련된 system call이 일어날 때마다 어딘가에 subdirectory summarized size를 caching 한다. (장난하냐?) File system '전체'의 free space에 대한 정보는 file system header에 operation이 일어날 때마다 항상 업데이트되기 때문에 df command로 (2)가 가능하지만, single directory에서는 이게 안된다(꼼수가 있긴 한데, 나중에 실험해보고 업데이트 해야겠다.) 젠장! 그런데 내가(혹은 독자가) 이 문제를 왜 고민하고 있는가? 바...
자세한 내용 보기

Proof of well-known 'Intersection Of Three Planes' formula.

-
Introduction & Motivation Recently, I found out well-known formula 'Intersection between Three Planes'  (Equation (22.31) of [1]) has no proof on literaly anywhere. Even (almost) the book where it is propsed[2]. So, I decided to figure out how this formula is derived.. (of course, experimentally it is wokrs very well.) Notation Let there be three planes on euclidean space $P_1, P_2, P_3$ which intersect in a point. If we let there normals as $n_1=(a_1, b_1, c_1), n_2=(a_2, b_2, c_2), n_3=(a_3, b_3, c_3)$, they must be linearly independent. Also, let the planes $P_1, P_2, P_3$ as the system of equation like below :  $$\begin{cases} a_1x+b_1y+c_1z+d_1=0 \\ a_2x+b_2y+c_2z+d_2=0 \\ a_3x+b_3y+c_3z +d_3=0 \end{cases}$$ Proof A. Solve problem with Linear Algebra. To find a intersection point by linear algebra, write system of equations as matrix form : $$\begin{bmatrix}a_1 & b_1 & c_1 \\ a_2 & b_2 & c_2 \\ a_3 & b_3 & c_3 \\ \end{bmatrix}\begin{bmatrix} x...
자세한 내용 보기

선형대수와 군

-
이번학기 현대대수 과목에서 부교재로 주어진 책이다. 처음에는 주교재만 열심히 보려고 했는데 선형대수를 다 까먹어서 진도를 1도 못나가게 되는 상황에 처해버린 것이다. 그래서 이 책을 Friedberg, David의 Linear Algebra 책과 같이 보려고 한다. 다른 책들은 딱딱해서 읽다가 조금씩 지치는데 이 책은 교수님이 유머를 섞어서 써놓으셨기 때문에 읽는데 재미있다. 그리고, 가끔씩(...) 던지는 연습문제들도 괜찮은 것 같다.  단, 어려운 개념을 갑가지 돌직구로 던진다던지, 2015년 개정판에서까지 국한혼용(...한자 6급한테)을 유지한다던지, 추상대수에서나 나오는 개념이나 notation을 써서 설명을 한다던지(책 이름이 선형대수와 군이니까...... 이 점은 이해 해야겠다.), 이북을 팔지 않아서 팔아프게 들고다녀야한다던지 하는 단점이 있다. 책이 하드커버인데다가 조금 두께가 있어서 스캔업체에다가 맡길까 생각중인데, 종이책을 날릴까 무서워서 고민중이다.
자세한 내용 보기