以球谐函数为例,说明如何使用python调用fortran的接口。
fortran 模块
建立fortran模块sh_value, 其中主要包括连带勒让得多项式plm(l,m,x)函数,以及球谐函数ylm(l,m,theta,phi)。
使用f2py生成动态库
通过f2py生成动态库sh.so,然后直接在pyhton中使用, 命令如下
勒让得多项式
建立python脚本 plot_pl.py 画出勒让得多项式, python需要先安装 matplotlib 和 numpy 库。
运行该程序,得到勒让得多项式图像
连带勒让得多项式
plot_plm.py 画出连带勒让得多项式, python需要先安装 matplotlib 和 numpy 库。
运行该程序,得到连带勒让得多项式图像
球谐函数
plot_ylm.py 画出球谐函数, python需要先安装 mayavi 和 numpy 库
运行该程序,得到球谐函数图像
领取专属 10元无门槛券
私享最新 技术干货