Legendre多项式程序在LMDE5下(4n+1)阶触发free()内存错误

问题现象 #

• 程序通过动态分配两个一维数组存储Legendre多项式计算结果，完成计算后释放内存。
• 仅当输入阶数为**5、9、13…即4*n+1（n>0）**时崩溃，抛出free(): invalid pointer错误并触发SIGABRT信号。
• 环境差异：Windows 11下运行完全正常，但在LMDE5（Linux Mint Debian Edition 5）环境中，将计算结果存入数组时触发问题。

崩溃时的运行输出（翻译后） #

输入Legendre多项式的阶数n。
4
迭代计算与递归计算结果相等：T
输入Legendre多项式的阶数n。
6
迭代计算与递归计算结果相等：T
输入Legendre多项式的阶数n。
5
free(): invalid pointer

程序收到信号SIGABRT：进程中止信号。

错误回溯：
#0 0x7fa540867bd0 in ???
#1 0x7fa540866e25 in ???
#2 0x7fa5406acd5f in ???
#3 0x7fa5406acce1 in ???
#4 0x7fa540696536 in ???
#5 0x7fa5406ee3a7 in ???
#6 0x7fa5406f5699 in ???
#7 0x7fa5406f6853 in ???
#8 0x55d4b8d74044 in __legendre_MOD_rleg
at src/legendre.f95:186
#9 0x55d4b8d73da2 in __legendre_MOD_rleg
at src/legendre.f95:186
#10 0x55d4b8d73da2 in __legendre_MOD_rleg
at src/legendre.f95:186
#11 0x55d4b8d74ac4 in leg
at src/legendre.f95:234
#12 0x55d4b8d74e91 in main
at src/legendre.f95:221
[1] 21146 中止 (核心已转储) ./myapp

可能的原因 #

1. 数组越界写入：Linux下glibc内存分配器对边界的检查远严于Windows平台，当阶数为4n+1时，数组大小计算可能存在误差，导致写入操作超出分配范围，破坏了malloc维护的内存管理结构，最终触发free时的错误。
2. 递归逻辑内存操作错误：回溯显示错误发生在递归函数rleg中，递归过程中可能存在数组重复释放、非法指针传递或索引计算错误，在特定递归深度（对应4n+1阶）下暴露问题。
3. Fortran索引误用：Fortran默认采用1-based索引，若代码中错误使用0-based索引，或数组大小计算时少算1项，在特定阶数下刚好触发越界。

解决步骤 #

1. 校验数组大小计算：Legendre多项式n阶包含n+1项（从x⁰到xⁿ），检查代码中动态数组的分配逻辑，确认是否在4n+1阶时错误设置数组大小为n而非n+1。
2. 添加越界检查：在数组写入的关键位置（如rleg函数的186行）添加索引合法性校验，例如：

   if (idx < 1 .or. idx > size(target_arr)) then
       print *, "Array out of bounds at index: ", idx
       stop 1
   end if
   

3. 用Valgrind定位内存问题：在LMDE5下运行valgrind ./myapp，它会精准定位内存越界的具体位置、写入地址与预期边界，直接指出错误根源。
4. 统一编译选项：检查Windows与Linux下的Fortran编译选项是否一致，在Linux下添加-Wall -Wextra -fcheck=all等严格检查选项，提前发现编译阶段的警告。
5. 梳理递归逻辑：检查rleg递归函数中数组的传递、赋值与释放逻辑，确认递归过程中没有出现指针非法复用、内存重复释放的情况。

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