用R语言调用Fortran子程序

2024-04-22 10:38

本文主要是介绍用R语言调用Fortran子程序,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

首先写这篇博客的原因是因为关于R调用Fortran之前摸索了两天,然后没弄通就放弃了,后来因为课题需求不得不用R调用大牛写的Fortran程序(毕竟自己写不出来),所以只能硬着头皮去调通程序,然后各种搜索资料,苦苦琢磨了两三天,最后终于在用is.loaded看是否导入成功的时候,看到了我期待已久的TRUE了,此时我的心率可以说比我跑完5公里之后的心率还要快,激动的我晚上直接奖励了自己一个大鸡腿,简直是高兴得不得了(恕我要求低,碰到这么点成就就差点喜极而泣)。

回到正题,首先需要说明,你必须跟着我做好每一步,一步都不能漏,等自己搞通了之后再考虑通过其路径,我给的方法是最直接最笨的方法(很适合没有一点计算机基础的同学,特别是看到dos窗口都觉得很神奇的同学)

1.首先毋庸置疑的是你需要安装R(当然加装Rstudio肯定更好用)和Fortran(我是在六维空间下载的,装的是Fortran6.5简版)

2.然后不要问我为什么,直接去官网下载Rtools(你可以不下,很抱歉,我只会这个方法,前提是你不用我写的方法),网址是https://cran.r-project.org/bin/windows/Rtools/,然后选择一个版本下载(我下载的版本是Rtools33.exe)。安装的时候尽量安装在C盘吧,当然我装的是D盘。重要的事情来了,你装的时候会有提示要不要修改环境变量path,这个请一定要勾选,然后就会自动把Rtools路径添加到path变量了。不过忘记勾选的朋友也不要惊慌,跟着我来修改一遍:找到计算机,然后右键进入属性,然后高级系统设置,单击高级里面的环境变量,点开之后下面的系统变量中找到path变量,单击编辑选项,在后面加上 ;C:\Rtools\bin;C:\Rtools\gcc-4.6.3\bin;  注意:我帮你把分号都写好了,你如果装在D盘就把C换成D就行了。

3.前面两步一定要全部做好了,才能进行接下来我要介绍的步骤,否则会出现命令无效或者命令不存在等等让你焦虑的错误,不要问我怎么改或者怎么办,返回去按照我说的再做一遍吧,再次声明:我只会这样做(怪我太笨)。接下来这个步骤是给那些可能都不知道Fortran怎么用的同学。先打开Fortran,点击file按钮然后new,然后选择Fortran Dynamic Link Library,给个项目名称multiply(这里采用的例子是来自网址https://en.wikiversity.org/wiki/R_(programming_language)/Tutorials/Connecting_Fortran_and_R),然后点击ok键,会弹出一个窗口,选项有两个(第一个说是创建空白库,第二个是会自动写入空函数,具体英文我就不写了,直接选择第二个就行了),我们选择第二个(免得又要自己去创建函数,省去这个步骤),然后finish就行了。创建完之后呢我们可以在multiply files中的source files下看到multiply.f90这个函数,双击进入函数,将给你的网址例子中的代码复制进去并编译执行(这个按钮就不用我教了吧,用过C的人都知道在哪里),执行完之后一定确保不报错才行哈。

4.前面三步很关键,没完成或者报错的同学就不要继续了,返回去好好改改,看看哪步忘记了。完成第三步之后我们就找到multiply这个文件夹(在你安装的Fortran6.5文件夹下),然后复制其中的文件multiply.f90,把它放在C:\Rtools\gcc-4.6.3\bin目录下。好了,到这里,你该做的准备工作就结束了,这是非常重要的,不要漏掉每个细节,不然没法继续下一步了。

5.此时就可以利用Rtools的gfortran工具了,具体操作是这样的,先打开DOS窗口(方法是:运行(快捷键是Windows键+R)->cmd),这个应该不用我教的,相信大家都会。然后在DOS里面执行下列代码:
c:  #设置路径到c盘下,如果装在D盘上就写d:
cd Rtools\gcc-4.6.3\bin  #相当于将目录更改到C:\Rtools\gcc-4.6.3\bin中
gfortran -shared -o multiply.dll multiply.f90  #执行完之后呢不能报错,否则就不会在bin目录下生成multiply.dll文件

6.完成第五步之后几乎该做的工作都做完了,剩下的就是R的事情了,R中的代码如下:
dyn.load("C:/Rtools/gcc-4.6.3/bin/multiply.dll") #加载动态库
is.loaded("multiply") #看是否导进来了,得到TRUE才说明加载成功
a=5
b=2
.Fortran("multiply",as.integer(a),as.integer(b),c=integer(1)) #调用mutiply函数

到这里你就漂亮的完成一次R调用Fortran了,我用了最简单最笨的的办法来进行的,应该都能看得懂。你应该庆祝下自己完成了如此漂亮的突破,但是这貌似只是刚起步,因为还有一件很棘手的事情需要你处理,那就是修改Fortran代码,改成R能调用的形式(给的这个例子未免过于简单了,因为就做了一个乘法),这就需要你稍微了解Fortran的语法和数据结构了。首先声明一下,R只能调用Fortran中的子程序(即subroutine,相当于R中的function),也不要问我为什么,别人是这么设计的我也无法回答(毕竟不是我写的,否则我就不会这么费事了)。这里我觉得我不用说太多,因为不会Fortran的我说了也是白说,会Fortran的看了我接下来的提醒就会知道怎么去修改程序了。

首先对于那些函数互相调用,出现一堆common的程序我相信很多人头都大了,这个你不用想太多,只要保证你想调用的主程序里面必须也有同样common,而且common中的变量确实存在(不存在的话就需要通过参数输入了)。然后呢就是变量的维数要做到R和Fortran中的保持一致,毕竟Fortran所有变量都需要指定类型(向量就要设置维数),所以你要把R中输入到Fortran中的参数事先维数统一。其他的就没啥问题了,前提是你的Fortran程序是正确的。

最后呢在这里要感谢在网上分享心得的热心网友,虽然我都不知道自己在哪看到的具体方法,但是希望看到我写的这篇博客的朋友在学习中有什么好的收获麻烦花点时间分享出来,互相学习才能更快的进步。



这篇关于用R语言调用Fortran子程序的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/925562

相关文章

Go语言连接MySQL数据库执行基本的增删改查

《Go语言连接MySQL数据库执行基本的增删改查》在后端开发中,MySQL是最常用的关系型数据库之一,本文主要为大家详细介绍了如何使用Go连接MySQL数据库并执行基本的增删改查吧... 目录Go语言连接mysql数据库准备工作安装 MySQL 驱动代码实现运行结果注意事项Go语言执行基本的增删改查准备工作

Java调用Python脚本实现HelloWorld的示例详解

《Java调用Python脚本实现HelloWorld的示例详解》作为程序员,我们经常会遇到需要在Java项目中调用Python脚本的场景,下面我们来看看如何从基础到进阶,一步步实现Java与Pyth... 目录一、环境准备二、基础调用:使用 Runtime.exec()2.1 实现步骤2.2 代码解析三、

Go语言使用Gin处理路由参数和查询参数

《Go语言使用Gin处理路由参数和查询参数》在WebAPI开发中,处理路由参数(PathParameter)和查询参数(QueryParameter)是非常常见的需求,下面我们就来看看Go语言... 目录一、路由参数 vs 查询参数二、Gin 获取路由参数和查询参数三、示例代码四、运行与测试1. 测试编程路

Python如何调用另一个类的方法和属性

《Python如何调用另一个类的方法和属性》在Python面向对象编程中,类与类之间的交互是非常常见的场景,本文将详细介绍在Python中一个类如何调用另一个类的方法和属性,大家可以根据需要进行选择... 目录一、前言二、基本调用方式通过实例化调用通过类继承调用三、高级调用方式通过组合方式调用通过类方法/静

C#控制台程序同步调用WebApi实现方式

《C#控制台程序同步调用WebApi实现方式》控制台程序作为Job时,需同步调用WebApi以确保获取返回结果后执行后续操作,否则会引发TaskCanceledException异常,同步处理可避免异... 目录同步调用WebApi方法Cls001类里面的写法总结控制台程序一般当作Job使用,有时候需要控制

Go语言使用net/http构建一个RESTful API的示例代码

《Go语言使用net/http构建一个RESTfulAPI的示例代码》Go的标准库net/http提供了构建Web服务所需的强大功能,虽然众多第三方框架(如Gin、Echo)已经封装了很多功能,但... 目录引言一、什么是 RESTful API?二、实战目标:用户信息管理 API三、代码实现1. 用户数据

Python用Flask封装API及调用详解

《Python用Flask封装API及调用详解》本文介绍Flask的优势(轻量、灵活、易扩展),对比GET/POST表单/JSON请求方式,涵盖错误处理、开发建议及生产环境部署注意事项... 目录一、Flask的优势一、基础设置二、GET请求方式服务端代码客户端调用三、POST表单方式服务端代码客户端调用四

Go语言网络故障诊断与调试技巧

《Go语言网络故障诊断与调试技巧》在分布式系统和微服务架构的浪潮中,网络编程成为系统性能和可靠性的核心支柱,从高并发的API服务到实时通信应用,网络的稳定性直接影响用户体验,本文面向熟悉Go基本语法和... 目录1. 引言2. Go 语言网络编程的优势与特色2.1 简洁高效的标准库2.2 强大的并发模型2.

Go语言使用sync.Mutex实现资源加锁

《Go语言使用sync.Mutex实现资源加锁》数据共享是一把双刃剑,Go语言为我们提供了sync.Mutex,一种最基础也是最常用的加锁方式,用于保证在任意时刻只有一个goroutine能访问共享... 目录一、什么是 Mutex二、为什么需要加锁三、实战案例:并发安全的计数器1. 未加锁示例(存在竞态)

Python跨文件实例化、跨文件调用及导入库示例代码

《Python跨文件实例化、跨文件调用及导入库示例代码》在Python开发过程中,经常会遇到需要在一个工程中调用另一个工程的Python文件的情况,:本文主要介绍Python跨文件实例化、跨文件调... 目录1. 核心对比表格(完整汇总)1.1 自定义模块跨文件调用汇总表1.2 第三方库使用汇总表1.3 导