关于SolverChallenge

SolverChallenge竞赛是解法器快速算法及应用研讨会（SOLVER会议）的一项特色活动，由SOLVER会议组委会主办，每年举行一次。竞赛针对实际应用中抽取出来的线性代数问题，要求参赛队伍在一定约束条件下通过算法设计与性能调优等手段对这些系统进行数值求解，并基于不同标准从多个维度对所获结果进行评价。第三届竞赛（SolverChallenge23）将与第六届解法器快速算法及应用研讨会（SOLVER23，青海西宁）同期启动，欢迎解法器领域的同行及研究生报名参赛。

竞赛背景

线性代数方程组 Ax=b 是科学与工程计算领域面临的共性问题，求解该类方程的解法器已成为很多应用软件的性能瓶颈。当前，国产高性能计算机的峰值性能已达每秒十亿亿次量级，很快将进入百亿亿次（E级）计算时代。这些国际领先的超级计算机为我国科学与工程计算应用迈进超大规模计算时代、实现更高精细的数值模拟成为可能。然而，超大规模计算给解法器算法设计与性能优化带来巨大挑战。一方面，超大规模计算必然是求解越来越复杂的实际应用问题，更加精细的物理建模蕴含着越来越复杂的应用特征（如：多时空尺度、多物理耦合、多介质、复杂流动、强非线性强间断、长时间积分等），从而导致现有解法器算法的计算复杂度急剧增长。另一方面，随着计算机硬件性能的提升，超级计算机呈现“多级嵌套海量并行、异构众核加速”的复杂体系结构特征，导致解法器算法的实现效率急剧下降。这种“应用与机器”的双重复杂性导致解法器越来越成为实际数值模拟应用的性能瓶颈，如何设计能同时匹配实际应用特征和机器体系结构特征的快速算法与性能优化技术是当前亟待解决的关键问题。

竞赛宗旨

SolverChallenge竞赛瞄准我国重要行业应用中的解法器瓶颈问题，赛题源于实际应用的需求，并根据每年的需求和优化情况持续更新。一方面，希望通过检验参赛队伍的算法与软件对这些问题的适应性，从中筛选出合适的算法及解法器策略，为应用提供新的思路与解决途径，提升我国重要行业应用中自主软件的解法器计算能力。另一方面，希望通过实际问题的牵引，吸引更多科研人员开展应用驱动的解法器算法创新与性能优化研究，同时培养具有交叉学科背景的优秀青年研究人员和研究生。

竞赛方案与规则

竞赛方案

竞赛组委会提供10套来源于不同应用领域的线性系统（每套线性系统包括矩阵A和右端向量b、存储格式的描述、求解约束条件等），参赛选手针对这10套系统进行求解，不限制求解方法（即迭代法、直接法或混合方法均可），也不限制解法器软件包的来源。竞赛网站将提供程序中输入、输出与计时部分，并开放接口供参赛选手提交自己的执行代码（若是闭源的，则提交链接文件）。最终将使用不同的右端向量b测试，其求解效果（速度、精度、浮点运算次数等）将作为本竞赛奖项评判标准。

竞赛规则

硬件环境：竞赛指定硬件作为程序性能的测试平台，参赛选手有权使用竞赛测试平台，但使用时间将被限制，在使用时间内参赛选手可以配置环境与测试程序。在竞赛提供的平台测出的成绩会作为参赛选手的最终成绩。

指定硬件包括：

1. 商用CPU（Intel Xeon Platinum 9242, 2路共96核，频率2.3 GHz, 内存384 GB DDR4）;
2. 商用GPU（NVIDIA V100，5120 CUDA核，频率1.2 GHz，显存32 GB HBM2）;
3. 国产CPU（Huawei Kunpeng 920-6426，2路共128核，频率2.6 GHz，内存512GB DDR4）;

软件环境：准许使用任何软件包参与竞赛，但是在评分时会考虑所使用软件包的自主性和创新性。

参赛队伍：竞赛面向本科生与研究生开放，学生采用自愿组合方式组成参赛队伍。每支队伍人数不超过3人，且应指定1-2名指导教师。

文档内容规范

算法设计实现涉及到的技术栈：包含采用的硬件、系统环境、软件；
算法详细内容：包含采用的算法、参数；
总体思路：包含问题的简要分析、求解策略；
结果分析（注明使用内存的情况）；
遇到的问题。

评分依据

求解速度（到解墙上时间）：评分标准同时考虑求解的线性系统个数与性能。每个系统满分10分，能正确解出且性能最优者得满分，其他选手的成绩依次顺延，未能解出者不得分。例如：求解某一系统时，排名第一的参赛小组用时120ms、排名第二的参赛小组用时240ms、排名第三的参赛小组未解出，则第一参赛小组得10分，第二参赛小组得5分，第三参赛小组得0分。所有10组系统得分累加得到总成绩。

限制：确保求解出的解向量x满足相应的精度要求，每道题限制使用单节点计算资源。

奖项设定

竞赛专委会将根据参赛作品的创新性、完成度、文档质量、难度以及答辩表现等情况进行评奖。

基础奖项：按照求解速度与浮点运算次数得分，分别设特等奖（不超1名），一等奖（不超过15%）、二等奖（不超过30%）、三等奖若干；
特别奖项：创新奖（在算法或性能优化方法方面有独特创新）、单系统奖（对于某一个线性系统有独特贡献）等；
企业特别奖：赞助厂商可根据自己的需求评选企业特别奖。

竞赛流程

事项	时间
开启报名	5月29日
公布赛题	5月29日
报名截止	6月23日
提交程序	7月3日
关闭提交程序	~~7月21日~~ 7月24日
初评、公布入围名单	7月28日
在SOLVER23会议上进行答辩、颁奖	8月7日-8月8日

竞赛截止报名时间

2023 年 6 月 23 日

竞赛报名方式

竞赛报名请点击此处填写表格，或扫描下方二维码

第三届解法器算法及优化竞赛（SolverChallenge23）赛题

百度网盘：下载链接提取码: uryj

赛题说明

端到端时间可分为如下3部分:
1. 分析
2. 准备
3. 求解
矩阵和向量的读入、格式转换、分发至多进程、矩阵坐标排序不列入统计时间。

直接法和迭代法对于上述3步的定义略有不同。

直接法：
1. 分析：利用矩阵结构信息（即不依赖矩阵数值），做解法器预处理，包括图分析、图优化、预开辟空间、节点优化编号等操作，参考intel pardiso的phase=11（注意，pardiso phase=11可以传入矩阵的值做一定的处理，本次比赛不允许）。
2. 准备：利用矩阵数值信息做数值分解，参考intel pardiso phase=22。
3. 求解：输入右端向量b，得到解向量x，参考intel pardiso phase=33。

迭代法：
1. 分析：利用矩阵结构信息（即不依赖矩阵数值），做解法器预处理，包括预条件子构建、多重网格构建等。
2. 准备：利用矩阵数值信息做数值准备，包括预条件子数值准备，多重网格数值准备等。
3. 求解：输入右端向量b，得到解向量x。

注：若将直接法作为预条件子用在迭代法中，则直接法部分参考“直接法”的说明计入迭代法三项相应时间中。

不同题目对于时间的要求不同，请选手注意。
此外，本次比赛没有内存使用限制，但是需要选手在最终报告中注明内存使用情况。

赛题名称	solverchallenge23_01
阶数非零元数	2,097,152
非零元数	14,581,760
应用领域	激光聚变
来源与背景描述	PDE：辐射流体力学三维单温方程；网格类型：结构网格；网格规模： 128×128×128；离散格式：有限体积7点格式。
MD5	f2c3e20cccc6d7842c6a6b1bbf53b872
约束	\|\|r\|\|₂ / \|\|b\|\|₂ < 1e-10
矩阵说明	0-base 矩阵（矩阵中非零元素索引从0开始）
关注时间	端到端时间

赛题名称	solverchallenge23_02
阶数非零元数	6,291,456
非零元数	52,133,888
应用领域	激光聚变
来源与背景描述	PDE：辐射流体力学三维三温方程；网格类型：结构网格；网格规模： 128×128×128；每个网格3个自由度。离散格式：有限体积7点格式。
MD5	38118979b8d7571317b10e06ea9f6487
约束	\|\|r\|\|₂ / \|\|b\|\|₂ < 1e-10
矩阵说明	0-base 矩阵（矩阵中非零元素索引从0开始）
关注时间	端到端时间

赛题名称	solverchallenge23_03
阶数非零元数	12,582,822
非零元数	678,943,908
应用领域	结构力学
来源与背景描述	PDE：三维线弹性方程；网格类型：非结构网格；网格规模：非结构网格; 离散格式：有限元。
MD5	c92ef3ea2fe76da54cb4b401e98fe617
约束	\|\|r\|\|₂ / \|\|b\|\|₂ < 1e-8
矩阵说明	1-base 矩阵（矩阵中非零元素索引从1开始）
关注时间	端到端时间

赛题名称	solverchallenge23_04
阶数非零元数	4,951,702
非零元数	612,096,612
应用领域	生物医疗
来源与背景描述	PDE：三维流体-相场-刚体耦合方程；网格类型：非结构四面体网格；网格规模：618,962网格点；离散格式：P1-P1压力稳定化，带刚体约束条件。
MD5	08416e32dfb1bc5713ebc61da27fbbd5
约束	\|\|r\|\|₂ / \|\|b\|\|₂ < 1e-10
矩阵说明	1-base 矩阵（矩阵中非零元素索引从1开始）
关注时间	端到端时间

赛题名称	solverchallenge23_05
阶数非零元数	3,514,495
非零元数	634,840,119
应用领域	生物医疗
来源与背景描述	PDE：三维流体-刚体耦合方程；网格类型：非结构四面体网格；网格规模：82,615网格点；离散格式：Taylor-Hood有限元，带刚体约束条件。
MD5	f4a6e95eadeb36dc383639c66890d723
约束	\|\|r\|\|₂ / \|\|b\|\|₂ < 1e-10
矩阵说明	1-base 矩阵（矩阵中非零元素索引从1开始）
关注时间	端到端时间

赛题名称	solverchallenge23_06
阶数非零元数	92,275,200
非零元数	1,742,860,800
应用领域	天气预报
来源与背景描述	PDE：全球天气预报Helmholtz方程; 网格类型：拟结构网格（东西周期边界，南北极网格跨极点中心对称）; 网格规模：1440×720×89；（同时提供相同问题的两个不同分辨率矩阵供参考： 1.网格规模：360x180x89；行数：5767200；非零元素个数：108928800。 2. 网格规模：720x360x89；行数：23068800；非零元素个数：435715200）离散格式：有限差分19点格式。
MD5	eb317190392353e8867094888b57c8d4
约束	\|\|r\|\|₂ < 1e-6
矩阵说明	1-base 矩阵（矩阵中非零元素索引从1开始）
关注时间	求解

赛题名称	solverchallenge23_07
阶数非零元数	1,765,074
非零元数	28,533,594
应用领域	电子学系统
来源与背景描述	PDE：滤波器电磁场仿真时谐麦克斯韦方程；离散格式：一阶Nedelec有限元，每条边一个自由度；特点：离散矩阵为复型系统。
MD5	a4902ea5ac570d78e8c240630585eac4
约束	\|\|r\|\|₂ / \|\|b\|\|₂ < 1e-8
矩阵说明	1-base 矩阵（矩阵中非零元素索引从1开始）
关注时间	端到端时间

赛题名称	solverchallenge23_08
阶数非零元数	5,809,680
非零元数	194,348,610
应用领域	集成电路
来源与背景描述	电路谐波平衡分析; 电路工艺：T65；特点：离散矩阵为复型系统。
MD5	dca1af651886114ea20568bddf0f8ac7
约束	\|\|r\|\|_max < 1e-8
矩阵说明	1-base 矩阵（矩阵中非零元素索引从1开始）
提示	请认真分析矩阵元素结构和数值分布，多为块矩阵；可以考虑迭代法，预处理子一般选为将每个块矩阵变成对角矩阵（或三对角、五对角等），或变成更小的对角块矩阵；也可以考虑直接法，注意填入元个数。该题块大小15x15，块类型：满阵/对角阵。
关注时间	准备+求解

赛题名称	solverchallenge23_09
阶数非零元数	7,785,500
非零元数	117,784,837
应用领域	集成电路
来源与背景描述	集成电路瞬态仿真方程；电路类型：DAC接口；电路工艺：7nm。
MD5	f7e9979e3aa665c1686e3cc474612d44
约束	\|\|r\|\|₂ / \|\|b\|\|₂ < 1e-8
矩阵说明	1-base 矩阵（矩阵中非零元素索引从1开始）
提示	一般用直接法
关注时间	准备+求解

赛题名称	solverchallenge23_10
阶数非零元数	418,718
非零元数	921,708,586
应用领域	集成电路
来源与背景描述	射频电路谐波平衡仿真方程；电路类型：VCO；电路工艺：28nm；特点：离散矩阵为复型系统。
MD5	c9a8230de6e60fd830d071bb73a7cbda
约束	\|\|r\|\|₂ / \|\|b\|\|₂ < 1e-8
矩阵说明	1-base 矩阵（矩阵中非零元素索引从1开始）
提示	请认真分析矩阵元素结构和数值分布，多为块矩阵；可以考虑迭代法，预处理子一般选为将每个块矩阵变成对角矩阵（或三对角、五对角等），或变成更小的对角块矩阵；也可以考虑直接法，注意填入元个数。
关注时间	准备+求解

第三届线性解法器算法与性能优化竞赛