[论文解读] smt: a Matlab structured matrices toolbox

Michela Redivo‐Zaglia, G. Rodríguez|arXiv (Cornell University)|Mar 30, 2009

Matrix Theory and Algorithms参考文献 13被引用 4

一句话总结

该论文介绍了 smt，一个 MATLAB 工具箱，通过面向对象编程为循环矩阵和托普利茨矩阵提供优化存储和快速算术运算，实现高效计算的同时保持透明性和可扩展性。它支持结构化矩阵运算、测试矩阵、循环预条件子，以及两种用于托普利茨线性系统的快速求解器，填补了公开可用、生产就绪的结构化矩阵软件中的关键空白。

ABSTRACT

We introduce the smt toolbox for Matlab. It implements optimized storage and fast arithmetics for circulant and Toeplitz matrices, and is intended to be transparent to the user and easily extensible. It also provides a set of test matrices, computation of circulant preconditioners, and two fast algorithms for Toeplitz linear systems.

研究动机与目标

• 解决公开可用、针对循环矩阵和托普利茨矩阵等结构化矩阵的优化软件缺乏的问题。
• 通过为结构化矩阵类型提供一等支持，扩展 MATLAB 原生对稀疏矩阵的支持能力。
• 提供一个透明、用户友好且可扩展的框架，模仿 MATLAB 对稀疏矩阵的处理方式。
• 使研究人员和工程师能够使用快速结构化矩阵算法，而无需从头重新实现。
• 支持图像去模糊、信号处理和大规模线性系统求解等性能关键型应用。

提出的方法

• 该工具箱利用 MATLAB 的面向对象编程系统，定义了新的数据类型：smcirc 用于循环矩阵，smtoep 用于托普利茨矩阵。
• 重载标准算术运算符和函数，以确保与现有 MATLAB 工作流的无缝集成。
• 实现优化的存储格式，以减少内存使用并加速矩阵运算。
• 工具箱包含一个图库函数（smtgallery），用于以结构化形式生成标准测试矩阵。
• 通过 MEX 文件接口集成两种用于托普利茨系统的快速求解器——toms729（基于 Fortran）和 tlls（基于 C），以实现高性能。
• 配置函数（smtconfig）允许用户自定义显示和行为设置，同时使用警告机制以保持一致性。

实验结果

研究问题

• RQ1 如何使结构化矩阵计算在 MATLAB 中的用户透明度和效率达到与稀疏矩阵操作相当的水平？
• RQ2 在 MATLAB 环境中，实现循环矩阵和托普利茨矩阵的优化存储与算术运算的最有效方法是什么？
• RQ3 能否设计一个模块化、可扩展的工具箱，以支持循环矩阵和托普利茨矩阵之外的未来结构化矩阵类别？
• RQ4 如何将现有的快速结构化矩阵算法暴露在一个高层级、生产就绪的软件框架中？
• RQ5 为确保此类工具箱的正确性、可移植性和可维护性，需要哪些关键的软件工程实践？

主要发现

• smt 工具箱成功利用 MATLAB 的面向对象系统，实现了对循环矩阵和托普利茨矩阵的优化存储与快速算术运算。
• 该工具箱支持一系列完整操作，包括矩阵创建、算术运算和线性系统求解，性能可与手工优化代码相媲美。
• 通过集成两个基于 MEX 的求解器（toms729 和 tlls），实现了大规模托普利茨系统的高性能求解。
• 工具箱通过 smtgallery 提供了一套完整的测试矩阵，其设计模仿 MATLAB 的 gallery 函数，支持实数和复数矩阵，并配有文档说明。
• 该软件与 MATLAB 6.5 及以上版本向后兼容，包含针对 isfloat 兼容性的补丁，并为多种平台提供预编译二进制文件。
• 该工具箱为开源并公开可用，网址为 http://bugs.unica.it/smt/，未来计划扩展支持更多结构化矩阵类别。