一种求解机组组合问题的内点半定规划gpu并行算法.pdfVIP

第 33 卷 第 7 期 电 力 自 动 化 设 备 Ｖｏｌ．33 Ｎｏ．7 ２０ 13 年 7 月 Ｅｌｅｃｔｒｉｃ Ｐｏｗｅｒ Ａｕｔｏｍａｔｉｏｎ Ｅｑｕｉｐｍｅｎｔ Jul. ２０ 13 一种求解机组组合问题的内点半定规划 GPU 并行算法 张宁宇 高 山 赵 欣 ， ， （ 东南大 学 电气工程 学院 ，江 苏 南京 210096 ） 摘 要 ： 针 对 内点法求解机 组组合 问题 的半定规 划 （SDP ）模 型 时大规模 线性 方程 组计 算 时 间太长的 问题 ，提 出 一种基 于 图形 处理 器 （GPU ）的 Krylov 子 空 间并行 算 法 。 该算 法采 用预条件 处理 的拟 最 小残 差 法 （QMR 法 ）， 并 以矩 阵分块技 术为基础 在 存储格 式下使 用 实现 并行预 处理矩 阵的计算 通 ， CSR GPU Incomplete Cholesky 。 过 对不 同规模 线性 方程 组 的计 算分析表 明 ，与传统的 直接 法相 比 ， 并行 算法具有速度和存储优 Cholesky QMR 势 ，可获得 良好 的并行加速 比。 10~100 机 6 个 系统 的仿 真结果也表 明 ，该 SDP 并行 内点法在 减 少计算 时 间的 同时可求得 近似 最优 解 。 关键词 ： 机 组组合 ； 半定规 划 ； GPU ； QMR ； 不完全 Cholesky 分解 ； 并行 算 法 ； Krylov ； 线性规 划 中图分 类号 ： TM 76 文献标识码 ： Ａ DOI ： 10.3969／ j .issn.1006-6047.20 13.07.022 0 引言 Gauss 消元 、Cholesky 分 解 、QR 分 解 等 ，将 原 方 程 组 化 为三 角 或 三对 角 等 容 易求 解 的形 式 ，然 后 通 过 回 半 定 规 划 SDP （SemiDefinite Programming ）是 线 代 或追赶等方法得 到方程组 的解 。 该方法一般 能得 性规划 的一种推广 ，它是在满 足约束 “对称矩 阵 的仿 到准确解 ，但 由于 固有 的前 推 回代 的特点 ，使得很 难 射 组 合 半 正定 ”的条 件 下使 线 性 函数 极 大 （极 小 ）化 实 现并 行 化 ，因而所 需 的 内存 与计 算 量 均很 大 。 迭 的问题 。 目前 为止 ，已有 学 者 对 SDP 在 电力 系统优 代 法 可分 为古 典 迭代 法 和 Krylov 迭 代 法 两 大类 ，其 化 调 度 中 的应 用 作 了初 步 研 究 ，文 献 ［1］使 用 SDP 中 ，古 典 迭 代 法 有 Jacobi 迭 代 法 、Gauss鄄Seidel 迭 代 算 法来 求解 中期 水 火 电力 系统 调 度 问题 ，得 到 了较 法 、SOR 等 。 目前 很 少 用 古 典 迭代 法 直 接 求 解