KmlDssGenAnalyze?I
分析数据结构,对稀疏矩阵进行重排和符号分解。
接口定义
C Interface:
int KmlDssGenAnalyzeSI(KmlSolverTask **pHandle);
int KmlDssGenAnalyzeDI(KmlSolverTask **pHandle);
int KmlDssGenAnalyzeCI(KmlSolverTask **pHandle);
int KmlDssGenAnalyzeZI(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, 10, 12, 14}; int ja[14] = {0, 7, 1, 6, 2, 5, 3, 4, 2, 5, 1, 6, 0, 7}; double a[14] = {1.0, 2.0, -2.0, 3.0, 3.0, 4.0, -4.0, 5.0, 4.0, -6.0, 3.0, 7.0, 2.0, 8.0}; KmlSolverTask *handle; int ierr = KmlDssGenInitDI(&handle, n, a, ja, ia); if (ierr != 0) { printf("\nERROR in KmlDssGenInitDI: %d", ierr); } ierr = KmlDssGenAnalyzeDI(&handle); if (ierr != 0) { printf("\nERROR in KmlDssGenAnalyzeDI: %d", ierr); }
父主题: LU