• umFPACK使用调用(一)


    测试平台: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;
    }

    最终显示正确结果

  • 相关阅读:
    负数幅角的选取
    记一次py交易
    区间估计
    平方和
    正态总体 下常用结论
    每日一背
    乘积的期望
    java调用javascript
    Java Agent入门
    JavaPoet入门
  • 原文地址:https://www.cnblogs.com/kmliang/p/2934418.html
Copyright © 2020-2023  润新知