基于线程池的parfor中调用geometryFromEdges为何出现解析错误？

问题根源 #

• PDE工具箱的线程兼容性限制：MATLAB的线程池（Thread-based parpool）与进程池（Process-based parpool）环境差异显著。pde.GeometricModel这类底层类的实现未考虑线程安全场景，在共享内存的线程并行环境中，类的元数据初始化或全局状态操作会出现冲突，导致系统无法正确识别父类定义。
• 共享内存中的类加载冲突：线程池内所有线程共享主进程的工作区与类定义，而PDE工具箱部分类在初始化时会修改全局状态，多线程同时操作会干扰类的正常加载流程，触发“找不到父类”的错误。
• 官方隐含限制：部分MATLAB工具箱（包括PDE Toolbox）对线程并行的支持有限，仅推荐在进程并行环境中使用相关功能。

可行解决方案 #

• 改用进程并行（推荐）：既然进程并行能正常运行，直接删除parpool('threads')语句，使用默认的进程池并行。每个worker进程拥有独立的地址空间，不会出现共享内存导致的类加载冲突。
• 预加载PDE类：在启动线程池前，先在主进程中完整执行一次PDE模型初始化流程，提前加载所需的类定义到共享内存中，避免线程内重复加载时出现错误。示例代码如下：

  delete(gcp('nocreate'));
  % 主进程预加载PDE相关类
  model = createpde();
  R1 = [3,4,0,1,1,0,1,1,0,0]';
  sf = 'R1';
  ns = char('R1')';
  gm = R1;
  g = decsg(gm,sf,ns);
  geometryFromEdges(model, g);
  clear model R1 sf ns gm g % 清理临时变量
  
  parpool('threads')
  parfor i = 1:1
      model = createpde();
      R1 = [3,4,0,1,1,0,1,1,0,0]';
      sf = 'R1';
      ns = char('R1')';
      gm = R1;
      g = decsg(gm,sf,ns);
      geometryFromEdges(model, g);
  end
  

• 升级MATLAB版本：旧版本的PDE Toolbox对线程并行的支持不完善，升级到最新版MATLAB可能修复此类兼容性问题。

内容的提问来源于stack exchange，提问作者Brazilian Cérebro