[发明专利]一种光路追踪计算方法及系统

说明书

本发明提供一种光路追踪计算方法及系统，包括：采集被测激光的参数值；基于所述采集的参数值，利用高精度计算方法求解预先构建的光线参数方程，计算激光传播时间；基于所述激光传播时间计算所述激光的轨迹点。本发明提供的技术方案根据浮点稳定性计算理论，采用高精度计算方法求解光线参数方程，提高了光路追踪的准确性。

技术领域

本发明涉及科学研究计算中的高性能计算领域，具体涉及一种光路追踪计算方法。

背景技术

激光在等离子体中的能量沉积问题是激光间接驱动惯性约束聚变(ICF，InertialConfinement Fusion)中的重要物理过程，其模拟精度会影响ICF全过程数值模拟结果的正确性。在ICF整体数值模拟中，目前激光主要采用几何光路建模，在模拟中需要通过求解光线与网格的相交关系来确定光线的运动轨迹。

光路追踪算法有广泛的应用，例如无限通信系统、晶体光学、电磁辐射传播及碰撞检测、激光的传播等。在大规模并行机上实现快速高精度光路追踪计算一直是一件有挑战性的任务。计算光路与一系列(3-D或2-D)几何体的交点是光路追踪算法的核心组成部分。光的传播轨迹是由光线与几何体的交点组成的。相交点处及附近区域的物质特性决定了光线的状态和传播方向。如何提高算法的精度和效率往往因应用问题而异。在激光约束聚变问题中，激光作为聚变的源项出现在物理建模之中。具体地，通常采用几何光线近似来模拟激光在等离子体中吸收与传播。激光束被表示成多条激光光线(直线或曲线)。每条光线被赋予一定的激光能量。随着激光的传播，其携带的能量沉积到沿途的等离子体中。在每个计算网格的电子数密度线性分布的假设下，求交算法的核心就是求出光线参数方程与每个网格边线段的交点。在算法的设计上，一方面需要考虑到电子数密度为常数等可能的退化情形，另一方面需要考虑到计算机数字的有限表示位数，理论上相同数值上不同的算法可能会得到不同的数值解。由于浮点计算的局限性，很多常用的算法难以计算出理论上的精确解。而且在计算交点的算法中当判据或阈值达到计算机浮点精度的量级时，数值结果还可能会因计算机不同而不同，数值结果不同可能会引起光路传播方向的改变，或者传播路径长度的改变。

IEEE二进制浮点数算术标准(IEEE 754)是20世纪80年代以来最广泛使用的浮点数运算标准，为许多CPU与浮点运算器所采用，这个标准定义了表示浮点数的格式与反常值，一些特殊数值，例如：无穷(Inf)与非数值(NaN)，以及这些数值的“浮点数运算符”；它也指明了四种数值舍入规则和五种例外状况，包括例外发生的时机与处理方式。

因为IEEE 754标准指明了浮点数的有效存储位数，导致浮点数在计算机中运算结果的不精确性来自于浮点数的有限位数无法精确表达自然界中的一个自然数，进一步导致了一些算法在理论上正确无误，但在计算机上则可能计算结果与真实数值不同，而计算错误。

发明内容

为解决上述技术的不足，本发明提供一种光路追踪计算方法，根据浮点稳定性计算理论，采用扩展精度浮点数计算方法求解光线参数方程，解决了减法相消带来的计算错误，提高了光线参数方程求解的准确性，能精准地实现光路的追踪。

一种光路追踪计算方法，其改进之处在于，所述方法包括：

采集被测激光的参数值；

基于所述采集的参数值，利用高精度计算方法求解预先构建的光线参数方程，计算激光传播时间；

基于所述激光传播时间计算所述激光的轨迹点；

所述光线参数方程包括一元四次方程。

优选地，所述基于所述采集的参数值，利用高精度计算方法求解预先构建的光线参数方程，计算激光传播时间包括：

将一元四次光线参数方程简化为一元三次方程；

基于扩展精度浮点数采用三角函数法求解所述一元三次方程得到所述一元三次方程的根；