句芒算法（morn）详情：矩阵运算

程序员文章站 2022-08-04 14:26:56

句芒算法（morn）介绍：矩阵运算（一）因为这两天有人质疑句芒算法的运行速度，想想，我也从来没有测试过，所以被质疑也很正常，所以打算以后测一测，然后顺便介绍一下句芒算法。因为这是第一次写关于句芒...

句芒算法（morn）介绍：矩阵运算（一）

因为这两天有人质疑句芒算法的运行速度，想想，我也从来没有测试过，所以被质疑也很正常，所以打算以后测一测，然后顺便介绍一下句芒算法。

因为这是第一次写关于句芒算法的博客，先说一句，句芒算法是一个轻量级的，用纯编写的算法库。

这篇先说一下算法库里最简单的几个矩阵运算。

测试程序

先看一下测试程序，这段程序再简单不过了，就是开了一个n×n的方阵，然后矩阵里的元素赋（0，1）之间的随机值。测试了三种基本的矩阵运算，①首先，算一下这个矩阵的行列式；②如果行列式不为0的话，算一下这个矩阵的逆；③把原来的矩阵乘以逆矩阵，验证一下是不是等于单位矩阵。

#include 
#include 

#include "morn_math.h"

void main()
{
    int i,j;
    mornmatrix *mat;
    mornmatrix *inv;
    mornmatrix *mul;
    int size = 8;

    float det;

    filetime beg,end;

    printf("use morn\n");

    mat = mornmatrixcreate(size,size,null);
    inv = mornmatrixcreate(size,size,null);
    mul = mornmatrixcreate(size,size,null);

    for(j=0;jdata[j][i] = ((float)rand())/((float)rand_max);

    printf("\nmat is:\n");
    printmat(mat);

    getsystemtimeasfiletime(&beg);
    det = mornmatrixdetvalue(mat);
    getsystemtimeasfiletime(&end);
    printf("\ndet is %f\n",det);
    printf("mat determinant time use is %d\n",(end.dwlowdatetime-beg.dwlowdatetime));

    if(det==0)
        return;

    getsystemtimeasfiletime(&beg);
    mornmatrixinverse(mat,inv);
    getsystemtimeasfiletime(&end);
    printf("\ninv is:\n");
    printmat(inv);
    printf("mat invere time use is %d\n",(end.dwlowdatetime-beg.dwlowdatetime));

    getsystemtimeasfiletime(&beg);
    mornmatrixmul(mat,inv,mul);
    getsystemtimeasfiletime(&end);
    printf("\nmul is:\n");
    printmat(mul);
    printf("mat mult time use is %d\n",(end.dwlowdatetime-beg.dwlowdatetime));

    mornmatrixrelease(mat);
    mornmatrixrelease(inv);
    mornmatrixrelease(mul);
}