gsl을 사용하여 정방 행렬의 역함수를 얻는 예제 코드가 필요합니다. LU 분해

Question

누군가가 gsl 함수 gsl_linalg_LU_decomp()를 호출하고 행렬의 역행렬을 얻는 것과 관련하여 친절하게도 예제 C++ 코드를 보여줄 수 있습니까? 대단히 감사합니다!

Answer 1

실제 역변환은 필요하지 않지만 Ax = b 유형의 문제를 해결해야합니다. 그렇다면 here이라는 좋은 예가 있습니다. 데이터에 STL 컨테이너를 사용하는 경우 (예 : 그런 다음 나는 LU 분해를 확보하고 함수 gsl_linalg_LU_invert를 호출하는 링크의 예를 따르의 실제 역을해야하는 경우 std::vector는, 당신은

std::vector<double> vec(length,val); 
gsl_needs_ptr_to_double(&vec[0]); 

처럼 첫 번째 데이터 항목에 대한 포인터를 전달해야합니다. gsl 라이브러리는 GNU 프로젝트이며 일반적으로 온라인에서 문서화가 잘되어 있으므로 약간의 시간을 들여 읽어 보는 것이 좋습니다.