我试图从C调用FORTRAN函数

我的问题是:

>如果fortRoutine是我的fortran子例程的名称,那么我从C调用这个
作为fortRoutine_.如果fortRoutine只包含一个字符数组参数,那么我可以这样传递:

fortRoutine_("I am in fortran");

>在调用FORTRAN子程序时,何时应该使用pass by value和通过引用?

因为我是C的新手,我对此没有任何线索.如果可能,请建议一些很好的教程链接.

解决方法

现在要做的就是在Fortran端使用Fortran ISO C Binding.这是Fortran 2003语言标准的一部分,并且可以在许多编译器中使用;它不是具体的gcc.这个网站已经在很多答案中描述过了.作为语言标准的一部分,它是编译器和平台无关的.并且您不需要知道编译器的内部传递约定. ISO C Binding在声明Fortran子程序或函数时使用时,会导致Fortran编译器使用C调用约定,以便可以从C直接调用该过程.您不需要添加隐藏的参数或命名mangle Fortran子例程名称,即没有下划线.链接器使用的名称来自“绑定”选项.

字符串是一个困难的情况,因为技术上在C它们是字符数组,你必须匹配在Fortran.您还必须处理字符串的不同定义:C为空终止,Fortran固定长度并填充空白.该示例显示了如何工作.数字更容易数组的唯一问题是C是行主和Fortran列主要,以便多维数组被转置.

int main ( void ) {

   char test [10] = "abcd";

   myfortsub (test);

   return 0;

}

subroutine myfortsub ( input_string ) bind ( C,name="myfortsub" )

   use iso_c_binding,only: C_CHAR,c_null_char
   implicit none

   character (kind=c_char,len=1),dimension (10),intent (in) :: input_string
   character (len=10) :: regular_string
   integer :: i

   regular_string = " "
   loop_string: do i=1,10
      if ( input_string (i) == c_null_char ) then
         exit loop_string
      else
         regular_string (i:i) = input_string (i)
      end if
   end do loop_string

   write (*,*) ">",trim (regular_string),"<",len_trim (regular_string)

   return

end subroutine myfortsub

您将C编译为对象文件,并使用gfortran编译fortran并链接两者:

gcc-mp-4.6   \
         -c  \
         test_fortsub.c

gfortran-mp-4.6   \
     test_fortsub.o  \
     myfortsub.f90  \
     -o test_fortsub.exe

输出为:

>abcd<           4

从C调用FORTRAN子例程的更多相关文章

  1. 为什么我的Fortran英特尔MKL非线性求解器初始化失败?

    我正在尝试在代码中使用“英特尔MKL”非线性求解器,但它无法初始化。此外,如果我尝试给eps_stop赋值为零,它会工作并打印出错误消息。当我给它赋值1e-6时,它返回了这个错误:forrtl:severe:已分配可分配数组我在调用strnlps_init之前和之后添加了测试打印输出消息,前一条打印输出,后一条没有,所以我认为问题出在这一行。然而,这些变量中没有一个是可分配的,所以我很困惑为什么会得到这个消息。知道可能发生了什么吗?

  2. Fortran函数是否开始从输入文件的特定行读取?

    有没有办法在某一行开始读取输入文件?我目前正在从文件开头读取每一行,只是为了跳过它,但这可能需要相当长的时间,具体取决于输入文件的大小。表示大致相同):

  3. 制作Fortran DLL并从C#调用它

    我在Fortran可执行文件中有一个函数,我需要把它变成一个dll文件,所以我可以从C#程序中调用它的函数P是一个浮点数,函数返回一个浮点数在这里,我对fortran一无所知,也没有用C#调用dll,所以请再解释一下.我正在使用CompaqVisualFortran和C#2008.感谢您的时间.(如果你喜欢,你可以看到完整的代码Here[这是一个计算水和蒸汽属性的程序])以下是使用单精度浮点数的示

  4. 由于在Ubuntu下的NumPy Fortran mixup,Theano失败了

    apt-getbuild-dep…

  5. 在Ubuntu上编译Fortran netCDF程序

    好的,这里有新问题.我正在尝试编译simple_xy_wr.f90–一个netCDF示例程序–在Ubuntu上使用gfortran,我必须做一些非常愚蠢的事情;我没有太多编译Fortran的经验.首先,我已经安装了libnetcdf-dev软件包,其中包括类似的文件所以,我试图用(各种命令)编译代码我得到以下输出我认为我包括了正确的库.例如.看来__netcdf_MOD_nf90_strerror

  6. 如何编写正则表达式来匹配字符串文字,其中转义是引号字符的两倍?

    我正在编写一个使用ply的解析器,需要识别FORTRAN字符串文字.这些引号用单引号引用,转义字符是双引号.即“我不明白你的意思”是一个有效的转义FORTRAN字符串.Ply以正则表达式输入.到目前为止,我的尝试不起作用,我不明白为什么.t_STRING_LIteraL=r“'[^(”)]*’”有任何想法吗?

  7. ubuntu 14.04 (64位)下 安装 intel fortran (ifort) 免费版

    licensetype=2&productid=2486,四个勾都打上,accept.最后,填写个人信息,需要注意的是邮箱地址要填写单位的。submit之后,在新打开的页面会提供下载链接及序列号。2)安装过程可选择在线安装(较快),即只下载.sh文件,然后sudobash*.sh运行。按提示,一步一步即可完成安装。

  8. 从C调用FORTRAN子例程

    我试图从C调用FORTRAN函数我的问题是:>如果fortRoutine是我的fortran子例程的名称,那么我从C调用这个作为fortRoutine_.如果fortRoutine只包含一个字符数组参数,那么我可以这样传递:>在调用FORTRAN子程序时,何时应该使用passbyvalue和通过引用?

  9. c – 为什么使用静态“常量”而不是实际值?

    我知道我正在查看的代码已从FORTRAN转换为C/C++.FORTRAN中是否有任何理由?FORTRAN对所有子例程参数使用pass-by-reference.但是,与其他具有pass-by-reference的语言不同,它仍然允许您将“rvalues”作为参数传递.在幕后,FORTRAN编译器转换代码如下:像代码一样这正是您的C代码正在做的事情:创建变量以将指针(或引用)传递给子例程FORTRAN样式.当然,这不是惯用的C或C.通常情况下,你会按值传递一倍.

随机推荐

  1. 从C到C#的zlib(如何将byte []转换为流并将流转换为byte [])

    我的任务是使用zlib解压缩数据包(已接收),然后使用算法从数据中生成图片好消息是我在C中有代码,但任务是在C#中完成C我正在尝试使用zlib.NET,但所有演示都有该代码进行解压缩(C#)我的问题:我不想在解压缩后保存文件,因为我必须使用C代码中显示的算法.如何将byte[]数组转换为类似于C#zlib代码中的流来解压缩数据然后如何将流转换回字节数组?

  2. 为什么C标准使用不确定的变量未定义?

    垃圾价值存储在哪里,为什么目的?解决方法由于效率原因,C选择不将变量初始化为某些自动值.为了初始化这些数据,必须添加指令.以下是一个例子:产生:虽然这段代码:产生:你可以看到,一个完整的额外的指令用来移动1到x.这对于嵌入式系统来说至关重要.

  3. 如何使用命名管道从c调用WCF方法?

    更新:通过协议here,我无法弄清楚未知的信封记录.我在网上找不到任何例子.原版的:我有以下WCF服务我输出添加5行,所以我知道服务器是否处理了请求与否.我有一个.NET客户端,我曾经测试这一切,一切正常工作预期.现在我想为这个做一个非托管的C客户端.我想出了如何得到管道的名称,并写信给它.我从here下载了协议我可以写信给管道,但我看不懂.每当我尝试读取它,我得到一个ERROR_broKEN_P

  4. “这”是否保证指向C中的对象的开始?

    我想使用fwrite将一个对象写入顺序文件.班级就像当我将一个对象写入文件时.我正在游荡,我可以使用fwrite(this,sizeof(int),2,fo)写入前两个整数.问题是:这是否保证指向对象数据的开始,即使对象的最开始可能存在虚拟表.所以上面的操作是安全的.解决方法这提供了对象的地址,这不一定是第一个成员的地址.唯一的例外是所谓的标准布局类型.从C11标准:(9.2/20)Apointe

  5. c – 编译单元之间共享的全局const对象

    当我声明并初始化一个const对象时.两个cpp文件包含此标头.和当我构建解决方案时,没有链接错误,你会得到什么如果g_Const是一个非const基本类型!PrintInUnit1()和PrintInUnit2()表明在两个编译单元中有两个独立的“g_Const”具有不同的地址,为什么?

  6. 什么是C名称查找在这里? (&amp;GCC对吗?)

    为什么在第三个变体找到func,但是在实例化的时候,原始变体中不合格查找找不到func?解决方法一般规则是,任何不在模板定义上下文中的内容只能通过ADL来获取.换句话说,正常的不合格查找仅在模板定义上下文中执行.因为在定义中间语句时没有声明func,并且func不在与ns::type相关联的命名空间中,所以代码形式不正确.

  7. c – 在输出参数中使用auto

    有没有办法在这种情况下使用auto关键字:当然,不可能知道什么类型的.因此,解决方案应该是以某种方式将它们合并为一个句子.这可用吗?解决方法看起来您希望默认初始化给定函数期望作为参数的类型的对象.您无法使用auto执行此操作,但您可以编写一个特征来提取函数所需的类型,然后使用它来声明您的变量:然后你就像这样使用它:当然,只要你重载函数,这一切都会失败.

  8. 在C中说“推动一切浮动”的确定性方式

    鉴于我更喜欢将程序中的数字保留为int或任何内容,那么使用这些数字的浮点数等效的任意算术最方便的方法是什么?说,我有我想写通过将转换放在解析的运算符树叶中,无需将表达式转化为混乱是否可以使用C风格的宏?应该用新的类和重载操作符完成吗?解决方法这是一个非常复杂的表达.更好地给它一个名字:现在当您使用整数参数调用它时,由于参数的类型为double,因此使用常规的算术转换将参数转换为double用C11lambda……

  9. objective-c – 如何获取未知大小的NSArray的第一个X元素?

    在objectiveC中,我有一个NSArray,我们称之为NSArray*largeArray,我想要获得一个新的NSArray*smallArray,只有第一个x对象…

  10. c – Setprecision是混乱

    我只是想问一下setprecision,因为我有点困惑.这里是代码:其中x=以下:方程的左边是x的值.1.105=1.10应为1.111.115=1.11应为1.121.125=1.12应为1.131.135=1.14是正确的1.145=1.15也正确但如果x是:2.115=2.12是正确的2.125=2.12应为2.13所以为什么在一定的价值是正确的,但有时是错误的?请启发我谢谢解决方法没有理由期望使用浮点系统可以正确地表示您的帖子中的任何常量.因此,一旦将它们存储在一个双变量中,那么你所拥有的确切的一

返回
顶部