C语言封装Eigen功能的自动存储期与动态分配问题咨询

1. 现有代码的核心问题 #

你当前编写的create_matrix实现属于未定义行为：

MatrixXf(r,c)是函数内的临时对象，属于自动存储期，生命周期仅存续于当前表达式执行结束，在create_matrix函数返回前就会被析构，你返回的是指向已销毁内存的野指针，C侧拿到后访问一定会触发异常。

2. 关于引用计数的误区 #

Eigen的普通矩阵类型（比如MatrixXf）本身没有内置引用计数机制，C也不会为普通栈对象自动添加引用计数，这个机制不管你是否跨C/C调用都不存在，不要做相关假设。

3. 动态分配的正确实现 #

你必须使用new做堆动态分配才能延长对象生命周期到函数返回后，正确写法如下：

封装接口补充 #

首先在头文件补充销毁接口，避免内存泄漏：

/* eigenwrapper.h */
#ifdef __cplusplus
extern "C" {
#endif

void* create_matrix(int r, int c);
void destroy_matrix(void* mat);
void* multiply_matrix(void* mat_a, void* mat_b);
// 其他接口声明

#ifdef __cplusplus
}
#endif

C++侧实现 #

/* eigenwrapper.cxx */
#include <Eigen/Dense>
#include "eigenwrapper.h"

extern "C" {

void* create_matrix(int r, int c) {
    // 堆分配矩阵，生命周期由手动控制
    return new Eigen::MatrixXf(r, c);
}

void destroy_matrix(void* mat) {
    // 安全转换后释放
    delete static_cast<Eigen::MatrixXf*>(mat);
}

void* multiply_matrix(void* mat_a, void* mat_b) {
    auto* a = static_cast<Eigen::MatrixXf*>(mat_a);
    auto* b = static_cast<Eigen::MatrixXf*>(mat_b);
    // 计算结果是临时对象，拷贝到堆上返回给C侧
    return new Eigen::MatrixXf((*a) * (*b));
}

}

C侧使用示例 #

/* main.c */
#include "eigenwrapper.h"

int main(void) {
    void* m1 = create_matrix(3,3);
    void* m2 = create_matrix(3,3);
    void* m3 = multiply_matrix(m1, m2);
    
    /* 执行其他矩阵操作 */
    
    // 用完必须手动释放所有矩阵，避免内存泄漏
    destroy_matrix(m1);
    destroy_matrix(m2);
    destroy_matrix(m3);
    return 0;
}

4. 注意事项 #

• C侧持有的void*只能传入你封装的C接口使用，禁止直接解引用访问内存
• 可以在运算类接口中添加维度校验等逻辑，返回错误码给C侧处理异常场景

内容的提问来源于stack exchange，提问作者Gabriella Chaos