Microsoft MVP성태의 닷넷 이야기
Math: 50. C# - MathNet.Numerics의 Matrix(행렬) 연산 [링크 복사], [링크+제목 복사]
조회: 2565
글쓴 사람
홈페이지
첨부 파일

C# - MathNet.Numerics의 Matrix(행렬) 연산

ML.NET의 경우 기본 라이브러리에는 행렬 타입이 없습니다. 그나마 기존에 Infer.NET으로 알려진 Microsoft.ML.Probabilistic에 Matrix 타입이 제공되므로 다음과 같이 추가하면 됩니다.

Install-Package Microsoft.ML.Probabilistic

<PackageReference Include="Microsoft.ML.Probabilistic" Version="0.3.1810.501" />

연산자 재정의가 되어 있으므로 행렬 곱도 자연스럽게 연산할 수 있고,

double[] data1 = Range.Create(1.0, 10.0, 1.0).ToArray();
Matrix mat1 = new Matrix(3, 3, data1);

double[] data2 = Range.Create(1.0, 10.0, 1.0).ToArray();
Matrix mat2 = new Matrix(3, 3, data2);

Matrix mat3 = mat1 * mat2;

Console.WriteLine(mat1);
/*
30  36  42
66  81  96
102 126 150
*/

요소 간 곱셈도 별도의 SetToElementwiseProduct 메서드를 이용하면 됩니다.

mat3.SetToElementwiseProduct(mat1, mat2);
Console.WriteLine(mat3);
/*
1  4  9
16 25 36
49 64 81
*/

하지만, 기능이 많지는 않습니다. 예를 들어 역행렬을 구하는 것조차 아직 미구현이니,

// 예외 발생: System.NotImplementedException: 'The method or operation is not implemented.'
Matrix invMat = mat3.Inverse();

의사 역행렬을 구하는 기능은 껍데기 메서드조차 없는 상태입니다.




이런 빈약한 행렬 처리를 보완하려면 MathNet.Numerics 라이브러리를 함께 연동하는 것이 좋습니다.

Install-Package MathNet.Numerics

<PackageReference Include="MathNet.Numerics" Version="4.7.0" />

어차피 2차원 배열로 행렬이 이뤄지기 때문에 ML.NET의 Matrix와 MathNet의 Matrix를 다음과 같이 서로 변환할 수 있습니다.

double[,] matData = ...;
Matrix mlMatrix = new Matrix(matData);

// ML.NET의 Matrix를 MathNet Matrix로 변환
Matrix<double> mathMatrix = CreateMatrix.DenseOfArray(mlMatrix.ToArray());

// MathNet Matrix를 ML.NET Matrix로 변환
mlMatrix = new Matrix(mathMatrix.ToArray());

따라서 이제부터는 MathNet의 Matrix를 이용해 의사 역행렬, 행렬 곱, 요소 간 행렬 곱 등의 연산을 할 수 있습니다.

Matrix<double> mat1 = CreateMatrix.DenseOfColumnMajor(3, 3, data1).Transpose();
Matrix<double> mat2 = CreateMatrix.DenseOfColumnMajor(3, 3, data2).Transpose();

Matrix<double> mat3 = mat1.Multiply(mat2); // 행렬 곱
               mat3 = mat1 * mat2; // 연산자도 지원
Matrix<double> mat4 = mat1.PointwiseMultiply(mat2); // 요소 간 곱
Matrix<double> invMat = mat4.Inverse(); // 역행렬

Matrix<double> mat5 = CreateMatrix.DenseOfArray(data.ToArray2D());
Matrix<double> pseudoInvMat = mat5.PseudoInverse(); // 의사 역행렬

(첨부 파일은 이 글의 예제 코드를 포함합니다.)




[이 글에 대해서 여러분들과 의견을 공유하고 싶습니다. 틀리거나 미흡한 부분 또는 의문 사항이 있으시면 언제든 댓글 남겨주십시오.]

[연관 글]





[최초 등록일: ]
[최종 수정일: 5/24/2019 ]

Creative Commons License
이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.
by SeongTae Jeong, mailto:techsharer@outlook.com

비밀번호

댓글 쓴 사람
 




1  2  3  4  5  6  7  [8]  9  10  11  12  13  14  15  ...
NoWriterDateCnt.TitleFile(s)
12080정성태12/16/2019847.NET Framework: 870. C# - 프로세스의 모든 핸들을 열람
12079정성태12/13/2019622오류 유형: 584. 원격 데스크탑(rdp) 환경에서 다중 또는 고용량 파일 복사 시 "Unspecified error" 오류 발생
12078정성태12/13/2019994Linux: 26. .NET Core 응용 프로그램을 위한 메모리 덤프 방법 [1]
12077정성태12/13/2019709Linux: 25. 자주 실행할 명령어 또는 초기 환경을 "~/.bashrc" 파일에 등록
12076정성태12/17/2019813디버깅 기술: 142. Linux - lldb 환경에서 sos 확장 명령어를 이용한 닷넷 프로세스 디버깅 - 배포 방법에 따른 차이
12075정성태12/18/2019896디버깅 기술: 141. Linux - lldb 환경에서 sos 확장 명령어를 이용한 닷넷 프로세스 디버깅
12074정성태12/11/2019585디버깅 기술: 140. windbg/Visual Studio - 값이 변경된 경우를 위한 정지점(BP) 설정(Data Breakpoint)
12073정성태12/10/20191033Linux: 24. Linux/C# - 실행 파일이 아닌 스크립트 형식의 명령어를 Process.Start로 실행하는 방법
12072정성태12/9/2019415오류 유형: 583. iisreset 수행 시 "No such interface supported" 오류
12071정성태12/9/2019656오류 유형: 582. 리눅스 디스크 공간 부족 및 safemode 부팅 방법
12070정성태12/9/20191170오류 유형: 581. resize2fs: Bad magic number in super-block while trying to open /dev/.../root
12069정성태12/19/2019822디버깅 기술: 139. windbg - x64 덤프 분석 시 메서드의 인자 또는 로컬 변수의 값을 확인하는 방법
12068정성태11/28/20191326디버깅 기술: 138. windbg와 Win32 API로 알아보는 Windows Heap 정보 분석 [2]파일 다운로드2
12067정성태11/27/2019820디버깅 기술: 137. 실제 사례를 통해 Debug Diagnostics 도구가 생성한 닷넷 웹 응용 프로그램의 성능 장애 보고서 설명 [1]파일 다운로드1
12066정성태11/27/2019581디버깅 기술: 136. windbg - C# PInvoke 호출 시 마샬링을 담당하는 함수 분석 - OracleCommand.ExecuteReader에서 OpsSql.Prepare2 PInvoke 호출 분석
12065정성태11/25/2019731디버깅 기술: 135. windbg - C# PInvoke 호출 시 마샬링을 담당하는 함수 분석파일 다운로드1
12064정성태11/25/20191042오류 유형: 580. HTTP Error 500.0/500.33 - ANCM In-Process Handler Load Failure
12063정성태11/21/2019764디버깅 기술: 134. windbg - RtlReportCriticalFailure로부터 parameters 정보 찾는 방법
12062정성태12/20/2019628디버깅 기술: 133. windbg - CoTaskMemFree/FreeCoTaskMem에서 발생한 덤프 분석 사례 - 두 번째 이야기
12061정성태11/20/2019713Windows: 165. CoTaskMemAlloc/CoTaskMemFree과 윈도우 Heap의 관계
12060정성태11/21/2019794디버깅 기술: 132. windbg/Visual Studio - HeapFree x64의 동작 분석
12059정성태11/20/2019730디버깅 기술: 131. windbg/Visual Studio - HeapFree x86의 동작 분석
12058정성태11/19/2019702디버깅 기술: 130. windbg - CoTaskMemFree/FreeCoTaskMem에서 발생한 덤프 분석 사례
12057정성태11/18/2019500오류 유형: 579. Visual Studio - Memory 창에서 유효한 주소 영역임에도 "Unable to evaluate the expression." 오류 출력
12056정성태11/18/2019921개발 환경 구성: 464. "Microsoft Visual Studio Installer Projects" 프로젝트로 EXE 서명 및 MSI 파일 서명 방법파일 다운로드1
12055정성태11/17/2019499개발 환경 구성: 463. Visual Studio의 Ctrl + Alt + M, 1 (Memory 1) 등의 단축키가 동작하지 않는 경우
1  2  3  4  5  6  7  [8]  9  10  11  12  13  14  15  ...