Question

我正在使用Ruby C API。我需要在函数内创建一个结构，但我想我正在做一些分配错误。这是我的代码

#include <stdio.h>
#include <ruby.h>
#include <ruby/thread.h>

typedef struct {
    double *matrix;
    int nrows;
    int ncols;
}Matrix;


void createMatrix(VALUE matrix, Matrix *mat) {
    printf("In\n");
    mat->nrows = RARRAY_LEN(matrix);
    VALUE firstElement = rb_ary_entry(matrix, 0);
    mat->ncols = RARRAY_LEN(firstElement);
    printf("Matrix shape: (%d,%d)\n", mat->nrows, mat->ncols);
    int i,j;

    double *tempMat = (double *)malloc(mat->nrows * mat->ncols * sizeof(double));

    printf("Allocated\n");
    VALUE row;
    for (i=0; i<mat->nrows; i++)
    {
        row = rb_ary_entry(matrix, i);
        for (j=0; j<mat->ncols; j++)
        {
            tempMat[i * mat->ncols + j] = NUM2DBL(rb_ary_entry(row, j));
//          printf("Matrix A Element(%d,%d)=%f\n", i, j, matA[i * colsA + j]);
        }
    }
    mat->matrix = tempMat;
    for (i=0; i<mat->nrows; i++)
        {
            for (j=0; j<mat->ncols; j++)
            {
                printf("Matrix temp Element(%d,%d)=%f\n", i, j, mat->matrix[i * mat->ncols + j]);
            }
        }

    printf("Assigned\n");
    return;
}


VALUE matmat_mul(VALUE self, VALUE matrixA, VALUE matrixB)
{
    int i,j;

    Matrix *matA;
    createMatrix(matrixA, matA);

    Matrix *matB;
    createMatrix(matrixB, matB);

    return Qnil;
}

void Init_la_ruby_ext()
{
    VALUE rg = rb_define_module("RG");
    VALUE linalg = rb_define_module_under(rg, "LinearAlgebra");
    VALUE operation = rb_define_class_under(linalg, "Operation", rb_cObject);
    rb_define_method(operation, "matmat_mul", matmat_mul, 2);
}

extconf.rb文件是

require 'mkmf'
extension_name = 'la_ruby_ext'
create_makefile(extension_name)

您可以使用

运行测试

require './la_ruby_ext'
rows = 3
cols = 3
mat = Array.new(rows){Array.new(cols)}
mat[0] = [0.0, 1.0, 2.0]
mat[1] = [3.0, 4.0, 5.0]
mat[2] = [6.0, 7.0, 8.0]
operation = RG::LinearAlgebra::Operation.new 
matC = operation.matmat_mul(mat, mat.transpose)

正如您在运行测试时所看到的，代码在第二次调用函数createMatrix时会出现分段错误。这是我的理解：

我首先创建一个指向struct Matrix
然后我将此指针作为参数传递给createMatrix
我分配了一个指向2D数组的指针
创建2D数组
我将指针指向matrix-＆gt; nmatrix

有什么想法吗？这是正确的方法吗？

Answer 1

鉴于定义void createMatrix(VALUE xxx, Matrix *mat)以及在函数内使用mat->nrows = …等，createMatrix()中对matmul_mul()的调用是错误的。

你有：

VALUE matmat_mul(VALUE self, VALUE matrixA, VALUE matrixB)
{
    int i,j;

    Matrix *matA;
    createMatrix(matrixA, matA);

    Matrix *matB;
    createMatrix(matrixB, matB);

你需要：

VALUE matmat_mul(VALUE self, VALUE matrixA, VALUE matrixB)
{
    int i,j;

    Matrix matA;
    createMatrix(matrixA, &matA);

    Matrix matB;
    createMatrix(matrixB, &matB);

在当前代码中，您将未初始化的指针传递给createMatrix()，然后分配给＆＃39; random＆＃39;指针指向的内存。这很少会带来快乐。在修订后的代码中，您有一对Matrix值，并将指针传递给createMatrix()函数的指针，然后为您填写详细信息。

C - 使用struct指针作为参数在函数内分配结构

1 个答案: