用遗迹提升矩阵

boost matrix with relic

本文关键字：遗迹更新时间：2023-10-16

我正在尝试使用fb_t的boost矩阵，这是一个表示有限域元素的遗迹对象。以下是如何从文档中定义fb_t：

typedef uint64_t dig_t
typedef align dig_t fb_t[FB_DIGS+PADDING(FB_BYTES)/(FB_DIGIT/8)]

这是我的代码：

#include <iostream>
#include <boost/numeric/ublas/matrix.hpp>
extern "C" {
#include <relic.h>
}
#include "relic_test.h"
using namespace std;
typedef boost::numeric::ublas::matrix<fb_t> matrix;
int main(void) {
    core_init();
    fb_param_set_any();
    fb_param_print();
    matrix mat(2,2);
   core_clean();
   return 0;
}

我得到以下错误：

Compiling: main.cpp
In file included from /usr/include/boost/numeric/ublas/vector.hpp:19:0,
                 from /usr/include/boost/numeric/ublas/matrix.hpp:16,
                 from /home/Foo/main.cpp:2:
/usr/include/boost/numeric/ublas/storage.hpp: In instantiation of ‘boost::numeric::ublas::unbounded_array<T, ALLOC>::unbounded_array(boost::numeric::ublas::unbounded_array<T, ALLOC>::size_type, const ALLOC&) [with T = unsigned int [4]; ALLOC = std::allocator<unsigned int [4]>; boost::numeric::ublas::unbounded_array<T, ALLOC>::size_type = unsigned int]’:
/usr/include/boost/numeric/ublas/matrix.hpp:131:92:   required from ‘boost::numeric::ublas::matrix<T, L, A>::matrix(boost::numeric::ublas::matrix<T, L, A>::size_type, boost::numeric::ublas::matrix<T, L, A>::size_type) [with T = unsigned int [4]; L = boost::numeric::ublas::basic_row_major<>; A = boost::numeric::ublas::unbounded_array<unsigned int [4], std::allocator<unsigned int [4]> >; boost::numeric::ublas::matrix<T, L, A>::size_type = unsigned int]’
/home/Foo/main.cpp:21:19:   required from here
/usr/include/boost/numeric/ublas/storage.hpp:71:23: error: functional cast to array type ‘boost::numeric::ublas::unbounded_array<unsigned int [4], std::allocator<unsigned int [4]> >::value_type {aka unsigned int [4]}’
/usr/include/boost/numeric/ublas/storage.hpp: In instantiation of ‘static void boost::numeric::ublas::unbounded_array<T, ALLOC>::iterator_destroy(T*&) [with T = unsigned int [4]; ALLOC = std::allocator<unsigned int [4]>; boost::numeric::ublas::unbounded_array<T, ALLOC>::iterator = unsigned int (*)[4]]’:
/usr/include/boost/numeric/ublas/storage.hpp:106:25:   required from ‘boost::numeric::ublas::unbounded_array<T, ALLOC>::~unbounded_array() [with T = unsigned int [4]; ALLOC = std::allocator<unsigned int [4]>]’
/usr/include/boost/numeric/ublas/matrix.hpp:90:11:   required from here
/usr/include/boost/numeric/ublas/storage.hpp:290:13: error: request for member ‘~boost::numeric::ublas::unbounded_array<unsigned int [4], std::allocator<unsigned int [4]> >::value_type’ in ‘* i’, which is of non-class type ‘unsigned int [4]’

我不太确定错误信息是关于什么的。知道吗？

我不认为boost:：ublas支持数组作为值类型。错误消息可能是因为数组不可复制。为了解决这个问题，可以将fb_t封装在一个类中：

struct fb_t_wrapper {
    fb_t value;
};

或使用std::array

using fb_t_array = std::array<dig_t, FB_DIGS+PADDING(FB_BYTES)/(FB_DIGIT/8)>;

用遗迹提升矩阵