如何解决PARI/GP中计算双二次数域基本单位范数时的‘[not a vector]’错误

问题根源 #

你的代码触发[not a vector]错误，核心原因是硬编码索引提取单位的方式不可靠：你固定用k1[8][5][1]获取基本单位，但当素数p、q足够大时，对应的二次域判别式增大，bnfinit返回的结构中单位相关字段的索引可能发生变化，直接按固定位置取值会访问到非向量类型的数据，导致报错。

修复代码 #

改用bnfinit返回结构的units字段提取基本单位（该字段专门存储单位群的生成元，排除 torsion 根单位），同时初始化空列表避免未定义变量问题：

lis4 = [];  // 提前初始化空列表，避免concat时出错
for (q=5, 200,
  for(p=q+1, 200,
    if (isprime(p) && isprime(q),
      k1 = bnfinit(y^2 - 2*p, 1);
      k2 = bnfinit(y^2 - q, 1);
      k3 = bnfinit(y^2 - 2*p*q, 1);
      // 从units字段提取第一个基本单位（非根单位）
      ep1 = k1.units[1];
      ep2 = k2.units[1];
      ep3 = k3.units[1];
      normep1 = nfeltnorm(k1, ep1);
      normep2 = nfeltnorm(k2, ep2);
      normep3 = nfeltnorm(k3, ep3);
      li = [[q,p], [normep1, normep2, normep3]];
      lis4 = concat(lis4, [li]);
    )
  )
);

关键说明 #

• bnfinit的units字段是专门用于存储单位群生成元的稳定字段，不受判别式大小影响，比硬编码索引更可靠。
• 必须初始化lis4 = []，否则第一次执行concat时会因变量未定义抛出错误。
• 若计算大判别式域时速度较慢，可尝试将bnfinit的flag参数改为2，启用更多计算优化，不影响结果正确性。

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