Question

来源：

//$ DO DECLARATION
#define VEC_MAX_SIZE 600000
mpz_t Vec[ VEC_MAX_SIZE ];

//$ DO INITIALIZATION
mpz_array_init( Vec[ 0 ], VEC_MAX_SIZE, 2000 ); // Stack overflow error - 3rd param

问题：

问题是我不知道如何正确初始化mpz_t数组，在这种情况下是Vec。错误是堆栈溢出。我认为它是由第3个参数引起的。这是mpz_array_init的原型/定义：

#define mpz_array_init __gmpz_array_init
__GMP_DECLSPEC void mpz_array_init __GMP_PROTO ((mpz_ptr, mp_size_t, mp_size_t));

错误：

Unhandled exception at 0x008e8377 in TestMPIR.exe: 0xC00000FD: Stack overflow.

注意：

我在Windows 7 64位上使用mpir 3.1.1。

这个来自mpir的例子很好用：

mpz_t Vec[ 20000 ];
mpz_array_init( Vec[ 0 ], 20000, 512 );

解答：

感谢@Zeta他/她的精彩回答。如果我错了，请纠正我。我只想在这里发布，以便其他新手可以从中学习。

#define MAX_VEC_SIZE 600000

malloc方法（我猜应该是C风格）：

mpz_t *Vec;

if( ( Vec = ( mpz_t *)::malloc( sizeof( mpz_t ) * MAX_VEC_SIZE ) ) == 0 ) {
    ::perror( "Out of memory" );
    return false;
}

//> Initialize Vec loop ( ... )

::free( Vec );

C ++风格？

mpz_t *VecEx = new mpz_t[ MAX_VEC_SIZE ];

//> Initialize VecEx loop ( ... )

delete [ ] VecEx;

这个答案基于@ Zeta的答案（接受的答案）。

Answer 1

对于如此大的向量（“堆栈溢出”），堆栈上的空间不足。你必须动态分配它：

mpz_t * Vec = malloc(VEC_MAX_SIZE * sizeof(mpz_t));

Vec大小约为6MiB。堆栈对于这样的对象来说不够大。请注意，GMP还提供C++ interface。

C ++ mpir：如何初始化mpz_t数组？

1 个答案: