KmlDssSymAnalyze?I
分析数据结构,对稀疏矩阵进行行/列重排和符号分解。
接口定义
C Interface:
int KmlDssSymAnalyzeSI(KmlSolverTask **pHandle);
int KmlDssSymAnalyzeDI(KmlSolverTask **pHandle);
int KmlDssSymAnalyzeCI(KmlSolverTask **pHandle);
int KmlDssSymAnalyzeZI(KmlSolverTask **pHandle);
参数
参数名 |
类型 |
描述 |
输入/输出 |
---|---|---|---|
pHandle |
KmlSolverTask ** |
求解器句柄,传入之前步骤的变量。 |
输入 |
返回值
返回值 |
类型 |
描述 |
---|---|---|
KMLSS_NO_ERROR |
int |
正常执行。 |
KMLSS_NULL_ARGUMENT |
int |
pHandle为空指针。 |
KMLSS_BAD_HANDLE |
int |
pHandle格式存在问题。 |
KMLSS_NONZERO_INDEXING |
int |
行索引指针不是从0开始,即ia[0]不等于0。 |
KMLSS_MISSING_DIAGONAL_ELEMENT |
int |
稀疏矩阵对角线元素全部为空。 |
KMLSS_INTERNAL_ERROR |
int |
稀疏矩阵对角线存在空元素/算法内部执行失败。 |
KMLSS_REORDERING_PROBLEM |
int |
排序算法失败。 |
依赖
#include "kml_dss.h"或#include "kml_solver.h"
示例
C Interface:
int n = 8; int ia[9] = {0, 2, 4, 6, 7, 8, 9, 10, 11}; int ja[11] = {0, 7, 1, 6, 2, 5, 3, 4, 5, 6, 7}; double a[11] = {1.0, 2.0, -2.0, 3.0, 3.0, 4.0, -4.0, 5.0, -6.0, 7.0, 8.0}; KmlSolverTask *handle; int ierr = KmlDssSymInitDI(&handle, n, a, ja, ia); if (ierr != 0) { printf("\nERROR in KmlDssSymInitDI: %d", ierr); } ierr = KmlDssSymAnalyzeDI(&handle); if (ierr != 0) { printf("\nERROR in KmlDssSymAnalyzeDI: %d", ierr); }
父主题: LDLT