DO CONCURRENT块内嵌BLOCK使用impure函数的技术问询

1. 该用法被允许的原因：属于gfortran的检查疏漏（Bug） #

根据Fortran标准（含Fortran 2018及更早版本），DO CONCURRENT块的核心约束是块内所有执行代码必须满足无副作用要求，而impure函数的定义本身就意味着它可能包含副作用（比如修改全局变量、执行I/O操作等），因此直接在DO CONCURRENT块内调用impure函数属于违反标准的行为。

你提到的在DO CONCURRENT内部嵌套BLOCK后可调用impure函数的情况，本质是gfortran（v12.2.0）的编译器检查机制存在疏漏。BLOCK构造是DO CONCURRENT块的一部分，其内部代码仍需遵守DO CONCURRENT的约束，标准并未允许通过嵌套BLOCK绕过这一限制。这种允许的行为不符合Fortran标准要求，属于编译器Bug。

2. 规避副作用的方法及gfortran的并行化行为 #

即便编译器允许这种写法，使用impure函数仍会带来不可预测风险，副作用可能导致DO CONCURRENT并行执行时出现数据竞争或结果不一致。可行的规避方案如下：

• 重构函数为pure类型：如果impure_func实际无副作用（比如仅因调用了某些非pure函数但未产生副作用），可将其修改为pure函数，从根源上符合DO CONCURRENT的要求。
• 将函数调用移出DO CONCURRENT块：若函数必须保留impure特性，先在串行代码中调用该函数生成结果数组，再将结果传入DO CONCURRENT块处理。
• 强制串行执行：若必须在DO CONCURRENT内调用，需确认函数无副作用，同时可通过编译器选项（如-fno-parallel-do-concurrent）强制DO CONCURRENT串行执行，避免并行带来的副作用问题。

关于gfortran的并行化：当编译器检测到DO CONCURRENT块内存在impure函数调用时，通常会自动禁用该块的并行优化，退化为串行执行。但这一行为依赖编译器的检测能力，并非绝对可靠，不能依赖它规避副作用风险。

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