[论文解读] A Basic Introduction on Math-Link in Mathematica

Sonam Maiti|ArXiv.org|Mar 1, 2006

Polynomial and algebraic computation参考文献 2被引用 3

一句话总结

本文提供了一份关于在 Mathematica 中使用 MathLink 集成外部程序（特别是 XL Fortran-90 代码）以加速复杂数值计算的基础指南。演示了如何使用 SetDirectory、Export、ReadList 和 Eigenvalues 等命令在 Mathematica 笔记本中编译、链接并执行外部 Fortran 程序，实现高效的数据交换与高性能计算工作流。

ABSTRACT

Starting from the basic ideas of mathematica, we give a detailed description about the way of linking of external programs with mathematica through proper mathlink commands. This article may be quite helpful for the beginners to start with and write programs in mathematica. In the first part, we illustrate how to use a mathemtica notebook and write a complete program in the notebook. Following with this, we also mention elaborately about the utility of the local and global variables those are very essential for writing a program in mathematica. All the commands needed for doing different mathematical operations can be found with some proper examples in the mathematica book written by Stephen Wolfram \cite{wolfram}. In the rest of this article, we concentrate our study on the most significant issue which is the process of linking of {\em external programs} with mathematica, so-called the mathlink operation. By using proper mathlink commands one can run very tedious jobs efficiently and the operations become extremely fast.

研究动机与目标

• 为初学者提供 Mathematica 笔记本使用与程序开发的全面介绍。
• 解释在 Mathematica 编程中局部变量与全局变量的作用及其区别。
• 演示如何使用 MathLink 将外部 Fortran 程序（特别是 XL Fortran-90）集成到 Mathematica 中，以提升计算性能。
• 说明在 Mathematica 与外部程序之间进行文件 I/O、目录管理与矩阵操作的实用技术。
• 使研究人员能够将计算密集型任务卸载到编译后的外部代码，同时在 Mathematica 内保持无缝的数据流。

提出的方法

• 使用 Mathematica 笔记本编写和执行程序，输入在单元格中处理，结果在按下 Shift+Enter 后在输出单元格中显示。
• 使用 Block 结构定义函数以封装局部变量并模块化代码，例如 sample[time_] 函数。
• 使用 XL Fortran 编译器（如 xlf 或 xlf95）将外部 Fortran-90 程序编译为可执行文件（例如 'mat'），并使用优化标志。
• 使用 SetDirectory[] 定位到包含外部可执行文件的目录，确保文件操作的路径正确解析。
• 使用 Export["mat3.dat", output] 将 Mathematica 中的矩阵数据写入文件，供外部程序使用。
• 使用 ReadList["!mat < mat3.dat", Number, RecordLists -> True] 通过 MathLink 调用外部程序，并将结果读回 Mathematica 以进行进一步处理。

实验结果

研究问题

• RQ1 如何有效将外部 Fortran 程序链接到 Mathematica 以实现高性能数值计算？
• RQ2 通过 MathLink 使用 Mathematica 集成 XL Fortran-90 源文件时，正确的编译与优化流程是什么？
• RQ3 如何使用基于文件的 I/O 和 MathLink 命令在 Mathematica 与外部程序之间交换数据？
• RQ4 局部变量与全局变量在构建与外部代码集成的稳健 Mathematica 程序中起什么作用？
• RQ5 从外部 Fortran 程序检索数据后，如何执行矩阵运算与特征值计算？

主要发现

• MathLink 允许在 Mathematica 笔记本中无缝执行外部 Fortran 程序，显著提升计算密集型任务的性能。
• 使用 Export 与 ReadList 配合 '!program < file' 语法，可可靠地调用外部程序并从编译后的程序中检索结果。
• 通过 SetDirectory[] 和 ResetDirectory[] 进行正确的目录管理，可确保程序执行期间路径解析正确。
• 在从外部程序读取数据后，可在 Mathematica 中计算矩阵运算（如 Eigenvalues），证明了端到端的数据流可行性。
• 只要使用正确的编译与链接命令，XL Fortran-90、F77 和 F95 程序通过 MathLink 集成是可行且高效的。
• 该方法支持涉及动态矩阵生成、文件 I/O 和符号计算的复杂工作流，可实现可扩展的数值研究。