测试平台:VS2008
新建工程-添加源文件和头文件
main函数
#include<iostream> #include"callumfpack.h" using namespace std; //int print(int **A) int main() { //int A[2][2]={0}; //printf("%d ",A[1][1]); /*int **A=new int *[2]; for(int i=0;i<2;i++) A[i]=new int [2]; for(int i=0;i<2;i++) for(int j=0;j<2;j++) A[i][j]=0;*/ //int A[5][5]={ // {2, 3, 0, 0, 0}, // {3, 0, 4, 0, 6}, // {0,-1,-3, 2, 0}, // {0, 0, 1, 0, 0}, // {0, 4, 2, 0, 1} //}; int n=5; int **A=new int *[n]; for(int j=0;j<n;j++) A[j]= new int[n]; for(int i=0;i<n;i++) for(int j=0;j<n;j++) scanf("%d",&A[i][j]); double b [ ] = {8., 45., -3., 3., 19.} ; double *x=new double[n]; print( A ,b, x); return 0; }
callumfpack头文件
//Data:2013-2-24 //修改了Ai Ax的类型 利用最大维数n*n来保存,可以调用正确结果 不过不知系统随机分配的值 函数没有用 //Data:2013-2-26 //调用UMFPACK包来实现求解方程组 //UMFPACK采用CSC(列压缩存储) matlab中的接口为A/b #include <stdio.h> #include <math.h> #include "umfpack.h" int print(int **A,double *b,double *x) { printf("--"); //printf("%d \n",A[1][1]); int n=5; double *null =(double *)NULL ; void *Symbolic, *Numeric ; int i,j; /* //定义矩阵A int **A=new int *[n]; for(j=0;j<n;j++) A[j]= new int[n]; for(i=0;i<n;i++) for(j=0;j<n;j++) scanf("%d",&A[i][j]); */ int *Ap=new int [n+1]; Ap[0]=0; int *Ai=new int[n*n]; for(i=0;i<n*n;i++) Ai[i]=1; double *Ax=new double[n*n]; for(i=0;i<n*n;i++) Ax[i]=1; double epsilon=0.00001; int NZnum=0;//矩阵非零元的个数 int k=0,l=0,m=0; for(j=0;j<n;j++) { for(i=0;i<n;i++) { if(abs(A[i][j])>epsilon) { Ai[l++]=i; Ax[m++]=A[i][j]; NZnum++; } } Ap[k+1]=NZnum; k++; } Ap[n]=NZnum; printf("--Ai---\n"); for (int i=0; i<n*n; i++) { printf("%d ", Ai[i]); } printf("\n"); umfpack_di_symbolic (n, n, Ap, Ai, Ax, &Symbolic, null, null); umfpack_di_numeric (Ap, Ai, Ax, Symbolic, &Numeric, null, null) ; umfpack_di_free_symbolic (&Symbolic); umfpack_di_solve (UMFPACK_A, Ap, Ai, Ax, x, b, Numeric, null, null); umfpack_di_free_numeric (&Numeric) ; for(i=0;i<n;i++) printf("x[%d]=%g\n", i, x[i]) ; return 0; }
最终显示正确结果