OpenBLAS用户手册：从编译到应用开发全指南-CSDN博客

本文链接：https://blog.csdn.net/gitblog_00730/article/details/148465528

OpenBLAS用户手册：从编译到应用开发全指南

OpenBLAS 项目地址: https://gitcode.com/gh_mirrors/ope/OpenBLAS

OpenBLAS是一个高性能的BLAS（基本线性代数子程序）库实现，广泛应用于科学计算、机器学习等领域。本文将全面介绍OpenBLAS的使用方法，包括编译安装、链接配置以及实际应用开发。

1. OpenBLAS编译指南

1.1 常规编译方法

OpenBLAS支持通过Make工具进行编译安装：

make  # 可使用-j参数并行编译，如-j4表示使用4个进程
make install

编译时OpenBLAS会自动检测CPU架构并进行优化。如需手动指定目标CPU架构，可使用TARGET参数：

make TARGET=NEHALEM  # 为Intel Nehalem CPU优化

完整的CPU架构支持列表可在项目文件TargetList.txt中查看。

1.2 交叉编译方法

交叉编译需要设置三个关键参数：

CC和FC：指定交叉工具链
HOSTCC：指定主机编译器
TARGET：明确指定目标CPU架构

示例（为ARM Cortex-A9编译）：

make CC=arm-linux-gnueabihf-gcc FC=arm-linux-gnueabihf-gfortran HOSTCC=gcc TARGET=CORTEXA9

1.3 调试版本编译

添加DEBUG=1参数可编译调试版本：

make DEBUG=1

1.4 自定义安装目录

使用PREFIX参数指定安装目录：

make install PREFIX=/custom/path

注意：安装时需保持与编译时相同的参数设置。

2. 链接OpenBLAS库

2.1 动态链接

动态链接需要指定库路径和头文件路径：

gcc -o test test.c -I/path/to/include -L/path/to/lib -lopenblas

可能需要的附加库：

多线程支持：-lpthread
LAPACK函数支持：-lgfortran

2.2 静态链接

静态链接更为简单：

gcc -o test test.c /path/to/libopenblas.a

3. 开发实战示例

3.1 CBLAS接口调用

以下示例展示如何使用CBLAS接口进行矩阵乘法运算：

#include <cblas.h>
#include <stdio.h>

int main() {
    double A[6] = {1.0,2.0,1.0,-3.0,4.0,-1.0};
    double B[6] = {1.0,2.0,1.0,-3.0,4.0,-1.0};
    double C[9] = {.5,.5,.5,.5,.5,.5,.5,.5,.5};
    
    cblas_dgemm(CblasColMajor, CblasNoTrans, CblasTrans,
                3,3,2,1,A,3,B,3,2,C,3);

    for(int i=0; i<9; i++)
        printf("%lf ", C[i]);
    printf("\n");
    return 0;
}

编译命令：

gcc -o test test.c -I/path/to/include -L/path/to/lib -lopenblas -lpthread -lgfortran

3.2 Fortran BLAS接口调用

从C代码调用Fortran风格的BLAS接口：

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

// 声明Fortran接口函数
extern void dgemm_(char*, char*, int*, int*, int*, 
                  double*, double*, int*, 
                  double*, int*, double*, 
                  double*, int*);

int main(int argc, char* argv[]) {
    if(argc < 4) {
        printf("需要输入m,n,k三个参数\n");
        return 1;
    }
    
    int m = atoi(argv[1]);
    int n = atoi(argv[2]);
    int k = atoi(argv[3]);
    
    // 矩阵初始化与计算...
    
    // 调用dgemm
    char trans = 'N';
    double alpha = 1.0, beta = 0.0;
    dgemm_(&trans, &trans, &m, &n, &k, 
           &alpha, A, &m, B, &k, &beta, C, &m);
    
    // 结果处理...
    return 0;
}