c++多线程不适合直接并行求解潮流/暂态方程本身，因节点变量强耦合、无法独立计算；但适合任务级（如区域划分、时间断面）和数据级（如雅可比矩阵分块组装、稀疏lu分解）并行。
直接说结论：C++多线程不适合直接并行求解大规模电力系统潮流/暂态方程本身，但适合在「任务级」和「数据级」做合理切分——比如按区域电网划分、按时间断面并行、或对雅可比矩阵块进行多线程组装与稀疏LU分解的调度。
为什么不能直接用 std::thread 对每个节点方程开一个线程
电力系统潮流方程（如 NR 法）本质是非线性耦合代数方程组，节点电压 V_i 和 θ_i 相互依赖，无法独立计算；暂态仿真中 dδ/dt、dω/dt 也强耦合于全网功率平衡。强行让每个线程算一个节点，会立刻遇到：race condition、inconsistent state、迭代不收敛。
- 典型错误现象：NR 迭代残差不下降，甚至发散；
std::vector::operator[]在多线程写同一块内存时触发heap-use-after-free（ASan 可捕获） - 真正可并行的是「重复结构化操作」：比如对 5000 个发电机模型同时调用
update_state()，只要它们不共享状态变量（如全局system_frequency），就安全 - 注意
std::thread创建开销约 10–100 μs，若单个计算粒度 std::async +std::launch::deferred或线程池
雅可比矩阵组装与稀疏求解的线程安全切分
大规模系统（>10k 节点）的雅可比矩阵是稀疏块状结构，非零元集中在导纳矩阵相关行。组装过程可并行，但必须避免写冲突：
- 按行分片：每个线程负责一组节点对应的雅可比行（例如线程 0 算节点 0–999，线程 1 算 1000–1999），前提是每行只由一个线程写 —— 需预分配好
std::vector<:vector>></:vector>的行指针，且禁止跨行更新 - 禁用
std::map或std::unordered_map动态插入非零元；改用std::vector<:tuple int double>></:tuple>先收集三元组，最后单线程压缩存储（CSR 格式） - 求解器层面：不要自己手写 LU，优先用线程安全的
SuiteSparse/UMFPACK（需编译时加-DNGPL=ON）或Intel MKL PARDISO，它们内部已做 OpenMP 分区；启用前确认链接了mkl_intel_thread而非mkl_sequential
时间域并行（Parareal）在暂态仿真中的落地要点
对长时间尺度暂态（如 10 秒、50Hz 采样 → 500k 步），传统串行积分（trapezoidal 或 IMEX-RK）瓶颈明显。Parareal 把时间轴切成若干段，用粗粒度解器（如隐式欧拉）快速生成初始猜测，再用细粒度解器（如 4 阶 RK）多线程修正 —— 这才是 C++ 多线程真正能发挥优势的场景：
- 关键约束：各时间段的初值必须来自上一段的终值，因此需显式同步点；推荐用
std::barrier（C++20）或std::condition_variable+ 计数器 - 粗解器输出必须是完整状态向量（
std::vector<double></double>），不能只存 delta；否则细解器无法启动 - 实测发现：当时间分段数 > 物理核心数 × 2 时，加速比开始下降；建议分段数设为
std::thread::hardware_concurrency() * 1.5 - 错误配置示例：
pardiso_handle被多个线程共用未加锁 → 触发PARDISO error -10（internal error）
最易被忽略的一点：电力系统模型里大量使用 std::shared_ptr 管理拓扑对象（如 Bus、Line），而 std::shared_ptr 的引用计数操作不是无锁的 —— 若 10 个线程同时析构同一对象，会卡在 __atomic_fetch_add_8。解决方法是提前明确所有权，用 std::unique_ptr + std::move 传递，或用索引（int bus_id）替代指针访问。
C++免费学习笔记（深入）：立即使用
在学习笔记中，你将探索 C++ 的入门与实战技巧！
